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: 
              if (options%lt_options%blocked_flow) then
                  call check( nf90_def_var(ncid, "froude", NF90_REAL, dimtwo_time, temp_id), trim(err)//"froude" )
                  call check( nf90_put_att(ncid,temp_id,"standard_name","froude_number"))
                  call check( nf90_put_att(ncid,temp_id,"long_name","Froude Number"))
                  call check( nf90_put_att(ncid,temp_id,"description", "Tendency for topographic blocking to occur = U / (Z * N)"))
                  call check( nf90_put_att(ncid,temp_id,"units","(m s-1) (m s-1)-1"))
  
negative.code:             call check( nf90_def_var(ncid, "clt", NF90_REAL, dimtwo_time, temp_id), trim(err)//"clt" )
              call check( nf90_put_att(ncid,temp_id,"standard_name","cloud_area_fraction"))
              call check( nf90_put_att(ncid,temp_id,"long_name","Fractional cloud cover"))
              call check( nf90_put_att(ncid,temp_id,"units","1"))
              call check( nf90_put_att(ncid,temp_id,"coordinates","lon lat"))
              varid(22)=temp_id
          endif
  