anchor.diff:          varid(32)=temp_id
           call check( nf90_inq_varid(ncid, "snow",temp_id), trim(err)//"snow" )
           varid(15)=temp_id
  -        if ((options%physics%microphysics==kMP_THOMPSON).or.(options%physics%microphysics==kMP_MORRISON)) then
  +        
  +        if ((options%physics%microphysics==kMP_THOMPSON)    &
  +        .or.(options%physics%microphysics==kMP_MORRISON)    &
  +        .or.(options%physics%microphysics==kMP_WSM6))       then
               call check( nf90_inq_varid(ncid, "graupel", temp_id), trim(err)//"graupel" )
               varid(16)=temp_id
           endif
  
positive.code:                 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
              endif
              ! need to modify dimids for staggered grids. 
              dimids(1)=xu_id
  
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)"))
  