anchor.diff:      ! Uniform distribution
       elens_geo_norm(ielens(iElem)) = (elens_r2(ielens(iElem))+elens_r1(ielens(iElem)))*&
                                       (elens_r2(ielens(iElem))-elens_r1(ielens(iElem)))
  -  else if(elens_type(ielens(iElem)) == 2) then
  -    ! Gaussian distribution
  -    elens_geo_norm(ielens(iElem)) = exp_mb(-0.5*(elens_r1(ielens(iElem))/elens_sig(ielens(iElem)))**2)&
  -                                    -exp_mb(-0.5*(elens_r2(ielens(iElem))/elens_sig(ielens(iElem)))**2)
  +!   else if(elens_type(ielens(iElem)) == 2) then
  +!     ! Gaussian distribution
  +!     elens_geo_norm(ielens(iElem)) = exp_mb(-0.5*(elens_r1(ielens(iElem))/elens_sig(ielens(iElem)))**2)&
  +!                                     -exp_mb(-0.5*(elens_r2(ielens(iElem))/elens_sig(ielens(iElem)))**2)
     end if
   
     if(st_debug) then
  
positive.code:   case("UNIFORM")
      elens_type(ielens(iElem)) = 1
    case("GAUSSIAN")
      elens_type(ielens(iElem)) = 2
      if(nSplit < 8) then
        write(lout,"(a,i0)") "ELENS> ERROR Expected at least 8 input parameters for GAUSSIAN, got ",nSplit
        iErr = .true.
        return
      end if
    case("CHEBYSHEV")
      elens_type(ielens(iElem)) = 3
      if(nSplit < 8) then
        write(lout,"(a,i0)") "ELENS> ERROR Expected at least 8 input parameters for CHEBYSHEV, got ",nSplit
        iErr = .true.
        return
      end if
    case default
      write(lout,"(a)") "ELENS> ERROR Elens type '"//trim(lnSplit(2))//"' not recognized. Remember to use all UPPER CASE."
      iErr = .true.
  
negative.code:     iErr = .true.
      return
    end if
    if(elens_r1(ielens(iElem)) <= zero) then
      write(lout,"(a)") "ELENS> ERROR R1<=0!"
      iErr = .true.
      return
    end if
  