anchor.diff:  C                   Data type - real array
   C                   Dimension WDOT(*) at least KK, the total number of
   C                   species.
  -c       March 22/2010--Luwi Oluwole:
  +c       March 22/2010--Luwi Oluwole: 
   c                      jsolv= -1 ===> load reduced models only; then return
  -c                      jsolv = 0 ===> don't calculate reaction rates;
  +c                      jsolv = 0 ===> don't calculate reaction rates; 
   c                                     only return reduced model species (binary) array
   c                      jsolv = 1 ===> calculate reaction rates
   C
  
positive.code:       common /saveNUR/ NUR, modelSize
  	character line*12
  C     Columns of NUR(,) contain all models including full model
  C     Columns of rangeLo and rangeHi contain valid ranges of reduced
  C     models beginning in column 2
  C     all reduced models are stored in order of increasing size
  C     <O.Oluwole--6/26/2008>
  C     Aerodyne Research, Inc.
  
negative.code: 	if(iCall.eq.2) goto 007
  C     On first call, read and save model library
  	call readModels(lin,kk,ii,rangeLo,rangeHi,nModels)
  	close(lin)	
  	iCall = 2
  C     <modelLib: signal completion of model library upload>
  	open(78,file='donemodels',status='new')
  	close(78)
  	if(jsolv .eq. -1) then
          return ! Luwi--4/23/2010--this call was for loading models only
          endif        
  C     </modelLib>
    007	continue
  C
  