anchor.diff:              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"))
  +            call check( nf90_put_att(ncid,temp_id,"coordinate","lat lon"))
               varid(5)=temp_id
   
               ! these should only be output for thompson microphysics
  
positive.code:                 call check( nf90_put_att(ncid,temp_id,"long_name","Graupel 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"))
                  varid(6)=temp_id
  
                  call check( nf90_def_var(ncid, "nr", NF90_REAL, dimids, temp_id), trim(err)//"nr" )
                  call check( nf90_put_att(ncid,temp_id,"standard_name","number_concentration_of_rain_particles_in_air"))
                  call check( nf90_put_att(ncid,temp_id,"long_name","Rain number concentration"))
                  call check( nf90_put_att(ncid,temp_id,"units","cm-3"))
                  varid(7)=temp_id
  
                  call check( nf90_def_var(ncid, "ni", NF90_REAL, dimids, temp_id), trim(err)//"ni" )
                  call check( nf90_put_att(ncid,temp_id,"standard_name","number_concentration_of_ice_crystals_in_air"))
                  call check( nf90_put_att(ncid,temp_id,"long_name","Cloud ice number concentration"))
                  call check( nf90_put_att(ncid,temp_id,"units","cm-3"))
                  varid(8)=temp_id
              elseif (options%physics%microphysics==kMP_MORRISON) then
                  call check( nf90_def_var(ncid, "qg", NF90_REAL, dimids, temp_id), trim(err)//"qg" )
  
negative.code:             call check( nf90_put_att(ncid,temp_id,"standard_name","grid_northward_wind"))
              call check( nf90_put_att(ncid,temp_id,"long_name","Grid relative northward wind"))
              call check( nf90_put_att(ncid,temp_id,"units","m s-1"))
              call check( nf90_put_att(ncid,temp_id,"coordinate","lat lon"))
              varid(10)=temp_id
              dimtwo(2)=y_id
  
  