ner wall if =0
                # phi(PF)=phintewi*te(ix,0) on PF wall if =1
  iphibcwo=0    # same for outer wall
  isutcore=2    # =1, set dut/dy=0 on iy=0 (if iphibcc=0)
                # =0, toroidal angular momentum=lzcore on iy=0 (iphibcc=0)
  isnewpot=1.0
  rnewpot=1.0
  cfnus_i=1.0;cfnus_e=1.0       #include collisionality in drift effects 
  isybdrywd=1   # use diffusive flux only on y boundary
lfililut=200;lenpfac=150;lenplufac=150
