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