anchor.diff:        CALL READMOD(IFL,N,ND,TEFF,R,NP,P,Z,ENTOT,VELO,GRADI,RSTAR,VDOP,NF,
        $             XLAMBDA(1 : NF),FWEIGHT(1 : NF),AKEY(1 : NF),
        $             ABXYZ,NATOM,MODHEAD,JOBNUM,LBLANK)
  +
         close(IFL)
  -
  -      allocate(entot(ND), enlte(ND))
  -      allocate(T(ND), RNE(ND))
  -      allocate(XJC(ND), XJCARR(ND, NF))
  -      allocate(R(ND), TAUROSS(ND))
  -      allocate(VELO(ND), GRADI(ND))
  -
  -      allocate(HTOT(ND), GTOT(ND), XTOT(ND), ETOT(ND))
  -
   
         IFL = 3; open(IFL, file = 'POPNUM', STATUS='OLD')
         !***MH readpop also reads number of depth points, temperature, rne
  
positive.code:       NP = read_param('NP')
        NF = read_param('NF')
  
        print*, 'fioss here: ', ND, NP, NF
  
        stop
  
        IFL = 3; open(IFL, file = 'MODFILE', STATUS='OLD'); rewind IFL
  
  
negative.code:       !        points in order to be able to add other lines, it is a
        !        good idea to make NFODIM quite a bit larger
  
        real*8, allocatable, dimension(:) ::       R
        real*8, allocatable, dimension(:) ::       entot, enlte
        real*8, allocatable, dimension(:) ::       T
  