anchor.diff:          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"))
           call check( nf90_put_att(ncid,temp_id,"coordinates","lon lat"))
           varid(23)=temp_id
  
positive.code:             call check( nf90_def_var(ncid, "hfgs", NF90_REAL, dimtwo_time, temp_id), trim(err)//"hfgs" )
              call check( nf90_put_att(ncid,temp_id,"standard_name","upward_heat_flux_at_ground_level_in_soil"))
              call check( nf90_put_att(ncid,temp_id,"long_name","Ground Heat Flux"))
              call check( nf90_put_att(ncid,temp_id,"positive","up"))
              call check( nf90_put_att(ncid,temp_id,"units","W m-2"))
              call check( nf90_put_att(ncid,temp_id,"coordinates","lon lat"))
              varid(25)=temp_id
  
negative.code:             call check( nf90_def_var(ncid, "rlds", NF90_REAL, dimtwo_time, temp_id), trim(err)//"rlds" )
              call check( nf90_put_att(ncid,temp_id,"standard_name","surface_downwelling_longwave_flux_in_air"))
              call check( nf90_put_att(ncid,temp_id,"long_name","Longwave downward radiation energy flux at the surface"))
              call check( nf90_put_att(ncid,temp_id,"units","W m-2"))
              call check( nf90_put_att(ncid,temp_id,"coordinates","lon lat"))
              varid(19)=temp_id
  