anchor.diff:  !    integer(kind=im),         intent( in)  :: i, j
   !    type(time_type), intent( in) :: time
   !    integer(kind=im), dimension(8) :: seed
  -!    
  +!
   !    ! Local variables
   !    integer(kind=im) :: year, month, day, hour, minute, second
  -!    
  -!    
  +!
  +!
   !    call get_date(time, year, month, day, hour, minute, second)
   !    seed = (/ i, j, year, month, day, hour, minute, second /)
   !  end function constructSeed
  
positive.code: !
  !  name     type     purpose
  !  ----   : ----   : ---------------------------------------------
  !fracrefao: real    
  !fracrefbo: real    
  ! kao     : real     
  ! kbo     : real     
  ! selfrefo: real     
  ! forrefo : real     
  !-----------------------------------------------------------------
  
        integer(kind=im), parameter :: no14 = 16
  
negative.code: !     of H2O to CO2.  The first index in the array, JS, runs
  !     from 1 to 10, and corresponds to different gas column amount ratios,
  !     as expressed through the binary species parameter eta, defined as
  !     eta = gas1/(gas1 + (rat) * gas2), where rat is the
  !     ratio of the reference MLS column amount value of gas 1
  !     to that of gas2.  The second index, JT, which
  !     runs from 1 to 5, corresponds to different temperatures.  More
  !     specifically, JT = 3 means that the data are for the corresponding
  !     reference temperature TREF for this  pressure level, JT = 2 refers
  !     to the TREF-15, JT = 1 is for TREF-30, JT = 4 is for TREF+15, and
  !     JT = 5 is for TREF+30.  The third index, JP, runs from 13 to 59 and
  !     refers to the corresponding pressure level in PREF (e.g. JP = 13 is
  