anchor.diff:          call check( nf90_put_att(ncid,lev_id,"standard_name","height"))
           call check( nf90_put_att(ncid,lev_id,"long_name","height above ground level"))
           call check( nf90_put_att(ncid,lev_id,"description","model level height (AGL)"))
  +        call check( nf90_put_att(ncid,lev_id,"positive","up"))
           call check( nf90_put_att(ncid,lev_id,"units","m"))
   
           if (calendar==GREGORIAN) then
  
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, "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,"units","W m-2"))
          call check( nf90_put_att(ncid,temp_id,"coordinates","lon lat"))
          varid(24)=temp_id
  