anchor.diff:              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_def_var(ncid, "hfss", NF90_REAL, dimtwo_time, temp_id), trim(err)//"hfss" )
  +            call check( nf90_put_att(ncid,temp_id,"standard_name","surface_upward_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"))
  
positive.code:             call check( nf90_put_att(ncid,temp_id,"units","kg m-2"))
              varid(31)=temp_id
          endif
  
      
          ! End define mode. This tells netCDF we are done defining metadata.
          call check( nf90_enddef(ncid) )
          
  
negative.code:             call check( nf90_put_var(ncid, varid(30), domain%snow_swe,     start_two_D),  trim(filename)//":snow_swe" )
              call check( nf90_put_var(ncid, varid(31), domain%canopy_water, start_two_D),  trim(filename)//":canopy_water" )
          endif
          call check( nf90_put_var(ncid, varid(34), domain%T2m, start_two_D),  trim(filename)//":t2m" )
          call check( nf90_put_var(ncid, varid(35), domain%Q2m, start_two_D),  trim(filename)//":hus2m" )
          call check( nf90_put_var(ncid, varid(36), domain%u10, start_two_D),  trim(filename)//":u10m" )
          call check( nf90_put_var(ncid, varid(37), domain%v10, start_two_D),  trim(filename)//":v10m" )
          
          last_rain=domain%rain
          ! Close the file, freeing all resources.
  