anchor.diff:          if (allocated(this_precip)) then
               deallocate(this_precip)
           endif
  +        if (allocated(this_snow)) then
  +            deallocate(this_snow)
  +        endif
       end subroutine mp_finish
   end module microphysics
  
positive.code:     integer :: update_interval
      real*8 :: last_model_time
      ! temporary variables
      real,allocatable,dimension(:,:) :: SR, last_rain, last_snow, this_precip,refl_10cm
  
  
      ! microphysics specific flag.  If it returns the current hourly precip (e.g. Morrison), then set this to false.
  
negative.code:             allocate(this_precip(ims:ime,jms:jme))
              this_precip=0
          endif
          if (.not.allocated(this_snow)) then
              allocate(this_snow(ims:ime,jms:jme))
              this_snow=0
          endif
  
          if (.not.allocated(refl_10cm)) then
              allocate(refl_10cm(ims:ime,jms:jme))
              refl_10cm=0
  