anchor.diff:          call check( nf90_def_var(ncid, "hfls", NF90_REAL, dimtwo_time, temp_id), trim(err)//"hfls" )
           call check( nf90_put_att(ncid,temp_id,"standard_name","surface_upward_latent_heat_flux"))
           call check( nf90_put_att(ncid,temp_id,"long_name","Latent 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(24)=temp_id
  
positive.code:             call check( nf90_def_var(ncid, "rlus", NF90_REAL, dimtwo_time, temp_id), trim(err)//"rlus" )
              call check( nf90_put_att(ncid,temp_id,"standard_name","surface_upwelling_longwave_flux_in_air"))
              call check( nf90_put_att(ncid,temp_id,"long_name","Longwave upward radiative energy flux from the surface"))
              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(29)=temp_id
  
negative.code:             call check( nf90_def_var(ncid, "rsds", NF90_REAL, dimtwo_time, temp_id), trim(err)//"rsds" )
              call check( nf90_put_att(ncid,temp_id,"standard_name","surface_downwelling_shortwave_flux_in_air"))
              call check( nf90_put_att(ncid,temp_id,"long_name","Shortwave 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(18)=temp_id
  