anchor.diff:                  call check( nf90_put_att(ncid,temp_id,"long_name","Snow number concentration"))
                   call check( nf90_put_att(ncid,temp_id,"units","cm-3"))
                   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"))
  +                varid(6)=temp_id
               endif
               ! need to modify dimids for staggered grids. 
               dimids(1)=xu_id
  
positive.code:                                  domain%snow + domain%snow_bucket*kPRECIP_BUCKET_SIZE, &
                                   start_two_D), trim(filename)//":snow" )
          
          if ((options%physics%microphysics==kMP_THOMPSON).or.(options%physics%microphysics==kMP_MORRISON)) then
              call check( nf90_put_var(ncid, varid(16), &
                                       domain%graupel + domain%graupel_bucket*kPRECIP_BUCKET_SIZE, &
                                       start_two_D), trim(filename)//":graupel" )
  
negative.code:         call check( nf90_put_att(ncid,temp_id,"units","kg m-2"))
          varid(15)=temp_id
          
          if     ((options%physics%microphysics==kMP_THOMPSON)        &
              .or.(options%physics%microphysics==kMP_MORRISON)        &
              .or.(options%physics%microphysics==kMP_WSM6))           then
              call check( nf90_def_var(ncid, "graupel", NF90_REAL, dimtwo_time, temp_id), trim(err)//"graupel" )
              call check( nf90_put_att(ncid,temp_id,"standard_name","graupel_amount"))
              call check( nf90_put_att(ncid,temp_id,"long_name","Combined large scale and convective graupel (accumulated)"))
  