anchor.diff:                  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"))
  +                call check( nf90_put_att(ncid,temp_id,"coordinate","lat lon"))
                   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"))
  +                call check( nf90_put_att(ncid,temp_id,"coordinate","lat lon"))
                   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"))
  +                call check( nf90_put_att(ncid,temp_id,"coordinate","lat lon"))
                   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,"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,"coordinate","lat lon"))
                   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,"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,"coordinate","lat lon"))
                   varid(39)=temp_id
               elseif (options%physics%microphysics==kMP_WSM6) then
                   call check( nf90_def_var(ncid, "qg", NF90_REAL, dimids, temp_id), trim(err)//"qg" )
                   call check( nf90_put_att(ncid,temp_id,"standard_name","mass_fraction_of_graupel_in_air"))
                   call check( nf90_put_att(ncid,temp_id,"long_name","Graupel 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(6)=temp_id
               endif
               ! need to modify dimids for staggered grids.
  
positive.code:         call check( nf90_put_att(ncid,temp_id,"long_name","Sensible heat flux"))
          call check( nf90_put_att(ncid,temp_id,"positive","up"))
          call check( nf90_put_att(ncid,temp_id,"units","W m-2"))
          varid(23)=temp_id
  
          call check( nf90_def_var(ncid, "hfls", NF90_REAL, dimtwo_time, temp_id), trim(err)//"hfls" )
  
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)
  