anchor.diff:              call check( nf90_inq_varid(ncid, "canwat", temp_id), trim(err)//"canwat" )
               varid(31)=temp_id
           endif
  +        call check( nf90_inq_varid(ncid, "ta2m", temp_id), trim(err)//"ta2m" )
  +        varid(34)=temp_id
  +        call check( nf90_inq_varid(ncid, "hus2m", temp_id), trim(err)//"hus2m" )
  +        varid(35)=temp_id
  +        call check( nf90_inq_varid(ncid, "u10m", temp_id), trim(err)//"u10m" )
  +        varid(36)=temp_id
  +        call check( nf90_inq_varid(ncid, "v10m", temp_id), trim(err)//"v10m" )
  +        varid(37)=temp_id
           
       end subroutine setup_varids
       
  
positive.code:             call check( nf90_put_att(ncid,temp_id,"units","W m-2"))
              varid(29)=temp_id
      
              call check( nf90_def_var(ncid, "shs", NF90_REAL, dimtwo_time, temp_id), trim(err)//"shs" )
              call check( nf90_put_att(ncid,temp_id,"standard_name","sensible_heat_flux"))
              call check( nf90_put_att(ncid,temp_id,"long_name","Sensible heat flux"))
              call check( nf90_put_att(ncid,temp_id,"positive","up"))
              call check( nf90_put_att(ncid,temp_id,"units","W m-2"))
  
negative.code:         if (options%physics%landsurface==kLSM_NOAH) then
              call check( nf90_inq_varid(ncid, "rlus",  temp_id), trim(err)//"rlus" )
              varid(29)=temp_id
              call check( nf90_inq_varid(ncid, "hfss",   temp_id), trim(err)//"hfss" )
              varid(23)=temp_id
              call check( nf90_inq_varid(ncid, "hfls",  temp_id), trim(err)//"hfls" )
              varid(24)=temp_id
  