anchor.diff:  
               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,"non_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"))
  
positive.code:         ! 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,"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"))
  
negative.code:             varid(4)=temp_id
  
              call check( nf90_def_var(ncid, "qs", NF90_REAL, dimids, temp_id), trim(err)//"qs" )
              call check( nf90_put_att(ncid,temp_id,"standard_name","mass_fraction_of_snow_in_air"))
              call check( nf90_put_att(ncid,temp_id,"long_name","Snow ice 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"))
  