anchor.diff:                  call check( nf90_put_var(ncid, varid(8),  reshape(domain%nice,  output_shape, order=zlast), start_three_D),trim(filename)//":nice" )
                   call check( nf90_put_var(ncid, varid(38), reshape(domain%ngraupel, output_shape, order=zlast), start_three_D),trim(filename)//":ngraupel" )
                   call check( nf90_put_var(ncid, varid(39), reshape(domain%nsnow,  output_shape, order=zlast), start_three_D),trim(filename)//":nsnow" )
  +            elseif (options%physics%microphysics==kMP_WSM6) then
  +                call check( nf90_put_var(ncid, varid(6),  reshape(domain%qgrau, output_shape, order=zlast), start_three_D),trim(filename)//":qgraupel" )
               endif
               output_shape(1)=output_shape(1)+1
               call check( nf90_put_var(ncid, varid(9),  reshape(domain%u,     output_shape, order=zlast), start_three_D),    trim(filename)//":u" )
  
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:                 varid(38)=temp_id
                  call check( nf90_inq_varid(ncid, "nsnow", temp_id), trim(err)//"nsnow" )
                  varid(39)=temp_id
              elseif (options%physics%microphysics==kMP_WSM6) then
                  call check( nf90_inq_varid(ncid, "qg", temp_id), trim(err)//"qg" )
                  varid(6)=temp_id
              endif
              call check( nf90_inq_varid(ncid, "w",  temp_id), trim(err)//"w" )
              varid(11)=temp_id
  