anchor.diff:      INTEGER                                   :: monday(13)=(/0,31,28,31,30,31,30,31,31,30,31,30,31/)
       INTEGER                                   :: mondayi(13)
       INTEGER                                   :: my1,my2,my3, tot_valid
  -    
  +
   ! CLWRF-UC June.09  (Copy from share/wrf_tsin.F)
       IF ( .NOT. READtrFILE ) THEN
          READtrFILE= .TRUE.
  -       
  +
          INQUIRE(FILE='CAMtr_volume_mixing_ratio', EXIST=exists)
  -       
  +
          IF (exists) THEN
             !iunit = get_unused_unit()
             iunit = io_newunit()
  
positive.code:        ENDDO
      ELSE  ! There is at least 1 valid year below yearIN
         IF (gas(yr_loc) < 0.) THEN
            
            ! Find the closest valid year, look for higher year first
            IF (any(gas(yr_loc:tot_years) > 0)) THEN
               DO idata=yr_loc+1, tot_years
  
negative.code:              CALL interpolate_CAMgases(yr, julian, nyrm, njulm, yr1, yr2, nyrp, njulp, cfc11r, cfc11vmr)
            endif
         endif
  
         IF (cfc11vmr < 0. .or. found_yearIN == 0) THEN
            CALL orig_val("CFC11",model,cfc11vmr)
         ELSE
            cfc11vmr=cfc11vmr*1.e-12
         ENDIF
      ENDIF
  
      IF (PRESENT(cfc12vmr)) THEN
         cfc12vmr = -9999.999
         if (found_yearIN /= 0) then
  