anchor.diff:      case(1) ! Uniform distribution
         elens_geo_norm(j) = (elens_r2(j)+elens_r1(j))*(elens_r2(j)-elens_r1(j))
       case(2) ! Gaussian distribution
  -      elens_geo_norm(j) = exp_mb(-0.5*(elens_r1(j)/elens_sig(j))**2)&
  -                         -exp_mb(-0.5*(elens_r2(j)/elens_sig(j))**2)
  +      elens_geo_norm(j) = exp_mb(-(((elens_r1(j)/elens_sig(j))*(elens_r1(j)/elens_sig(j)))/two)) &
  +                         -exp_mb(-(((elens_r2(j)/elens_sig(j))*(elens_r2(j)/elens_sig(j)))/two))
       case(3) ! Radial profile
         elens_radial_fr1(j) = polinterp( elens_r1(j), &
               elens_radial_profile_R(0:elens_radial_profile_nPoints(elens_iRadial(j)),elens_iRadial(j)), &
  
positive.code:   use mod_common, only : beta0, napx
    use mod_common_main, only : xv1, xv2, yv1, yv2, moidpsv, rvv
    use mathlib_bouncer
    use numerical_constants, only : zero, one, two
    use mod_utils, only : polinterp
    use crcoall, only : lerr
  
  
negative.code:   use mod_utils, only : polinterp
    use mod_common, only : bez, kz, e0f, nucm0, fort3
    use mod_settings, only : st_debug
    use numerical_constants, only : c1e3, two
    use crcoall, only : lout, lerr
  
    integer j, jj, jguess
  