anchor.diff:  
           elseif (options%physics%microphysics==kMP_WSM3) then
               if (this_image()==1) write(*,*) "    WSM3 Microphysics"
  -            call wsm3init(rhoair0,rhowater,rhosnow,cliq,cpv, .True.)
  +            call wsm3init(rhoair0,rhowater,rhosnow,cliq,cpv, allowed_to_read=.True.)
               precip_delta=.True.
           endif
   
  
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.
  