anchor.diff:              call check( nf90_put_att(ncid,temp_id,"standard_name","water_vapor_mixing_ratio"))
               call check( nf90_put_att(ncid,temp_id,"long_name","Water Wapor Mixing Ratio"))
               call check( nf90_put_att(ncid,temp_id,"units","kg kg-1"))
  +            call check( nf90_put_att(ncid,temp_id,"coordinate","lat lon"))
               varid(1)=temp_id
   
               call check( nf90_def_var(ncid, "qc", NF90_REAL, dimids, temp_id), trim(err)//"qc" )
               call check( nf90_put_att(ncid,temp_id,"standard_name","cloud_liquid_water_mixing_ratio"))
               call check( nf90_put_att(ncid,temp_id,"long_name","Cloud liquid water content"))
               call check( nf90_put_att(ncid,temp_id,"units","kg kg-1"))
  +            call check( nf90_put_att(ncid,temp_id,"coordinate","lat lon"))
               varid(2)=temp_id
   
               call check( nf90_def_var(ncid, "qi", NF90_REAL, dimids, temp_id), trim(err)//"qi" )
               call check( nf90_put_att(ncid,temp_id,"standard_name","cloud_ice_mixing_ratio"))
               call check( nf90_put_att(ncid,temp_id,"long_name","Cloud ice content"))
               call check( nf90_put_att(ncid,temp_id,"units","kg kg-1"))
  +            call check( nf90_put_att(ncid,temp_id,"coordinate","lat lon"))
               varid(3)=temp_id
   
               call check( nf90_def_var(ncid, "qr", NF90_REAL, dimids, temp_id), trim(err)//"qr" )
  
positive.code:             call check( nf90_put_att(ncid,temp_id,"standard_name","grid_eastward_wind"))
              call check( nf90_put_att(ncid,temp_id,"long_name","Grid relative eastward wind"))
              call check( nf90_put_att(ncid,temp_id,"units","m s-1"))
              varid(9)=temp_id
              dimtwo(1)=x_id
  
  
negative.code:         call check( nf90_put_att(ncid,temp_id,"standard_name","surface_temperature"))
          call check( nf90_put_att(ncid,temp_id,"long_name","Land surface skin temperature"))
          call check( nf90_put_att(ncid,temp_id,"units","K"))
          call check( nf90_put_att(ncid,temp_id,"coordinate","lat lon"))
          varid(28)=temp_id
  
          ! diagnosed surface fields (e.g. T2m, U10)
  