anchor.diff:                                rainncv = this_precip,                                & ! not used outside thompson (yet)
                                 sr = SR,                                              & ! not used outside thompson (yet)
                                 snow = domain%accumulated_snowfall%data_2d,           &
  +                              snowncv = this_snow,                                  &
                                 has_reqc=0, has_reqi=0, has_reqs=0,     &
                                 ids = ids, ide = ide,                   & ! domain dims
                                 jds = jds, jde = jde,                   &
  
positive.code:         if (allocated(this_precip)) then
              deallocate(this_precip)
          endif
      end subroutine mp_finish
  end module microphysics
  
negative.code:     integer :: update_interval
      real*8 :: last_model_time
      ! temporary variables
      real,allocatable,dimension(:,:) :: SR, last_rain, last_snow, this_precip, this_snow,refl_10cm
  
  
      ! microphysics specific flag.  If it returns the current hourly precip (e.g. Morrison), then set this to false.
  