anchor.diff:  C              8-PARAMETER SRI FORM
   C
                     XP = 1.0/(1.0 + PRLOG**2)
  -                  FC = ((FPAR(4,N)*EXP(-FPAR(5,N)/T)
  +                  FC = ((FPAR(4,N)*EXP(-FPAR(5,N)/T) 
        1                   + EXP(-T/FPAR(6,N))) **XP)
        2                  * FPAR(7,N) * T**FPAR(8,N)
   C
  
positive.code: 	print*,'nsteps= ',nsteps
  	print*,'rmaxfull= ',rmaxfull
  	print*,'reading models.chem...'
  
  	open(1, file='models.chem', status='old')
  	read (1,*) ii, kk!Luwi-9/2/10-kfuel not needed/used!, kfuel
  C     closed in subroutine findModel
  
negative.code: 		  elseif(line(:6).eq.'REGION') then
  C     model size not specified; never look for model size
  		      ireadSize = 0
  		  endif	      
  	      endif
  	   endif
  	   if(line(:6).eq.'REGION') then
  