anchor.diff:          ! else we are just storing surface fields, create 2D fields for a lot of the 3D variables, but only store the lowest model layer
           else
               call check( nf90_def_var(ncid, "qv", NF90_REAL, dimtwo_time, temp_id), trim(err)//"qv" )
  -            call check( nf90_put_att(ncid,temp_id,"standard_name","water_vapor_mixing_ratio"))
  +            call check( nf90_put_att(ncid,temp_id,"non_standard_name","mass_fraction_of_water_vapor_in_air"))
               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,"coordinates","lon lat"))
  
positive.code:             varid(3)=temp_id
  
              call check( nf90_def_var(ncid, "qr", NF90_REAL, dimids, temp_id), trim(err)//"qr" )
              call check( nf90_put_att(ncid,temp_id,"standard_name","mass_fraction_of_rain_with_air"))
              call check( nf90_put_att(ncid,temp_id,"long_name","Rain water content"))
              call check( nf90_put_att(ncid,temp_id,"WARNING","Could be mixing ratio, not mass fraction w/thompson scheme"))
              call check( nf90_put_att(ncid,temp_id,"units","kg kg-1"))
  
negative.code:                 varid(8)=temp_id
  
                  call check( nf90_def_var(ncid, "ngraupel", NF90_REAL, dimids, temp_id), trim(err)//"ngraupel" )
                  call check( nf90_put_att(ncid,temp_id,"non_standard_name","number_concentration_of_graupel_particles_in_air"))
                  call check( nf90_put_att(ncid,temp_id,"long_name","Graupel number concentration"))
                  call check( nf90_put_att(ncid,temp_id,"units","cm-3"))
                  call check( nf90_put_att(ncid,temp_id,"coordinates","lon lat"))
                  varid(38)=temp_id
  
                  call check( nf90_def_var(ncid, "nsnow", NF90_REAL, dimids, temp_id), trim(err)//"nsnow" )
                  call check( nf90_put_att(ncid,temp_id,"non_standard_name","number_concentration_of_snow_particles_in_air"))
                  call check( nf90_put_att(ncid,temp_id,"long_name","Snow number concentration"))
                  call check( nf90_put_att(ncid,temp_id,"units","cm-3"))
                  call check( nf90_put_att(ncid,temp_id,"coordinates","lon lat"))
  