anchor.diff:              if (qc>SMALL_VALUE) then
                   if (t>freezing_threshold) then
                       ! convert cloud water to rain drops
  -                    call cloud2hydrometeor(qc,qr,cloud2rain)
  +                    call cloud2hydrometeor(qc,qr,cloud2rain,rain_cloud_init)
   !                   if (qs>SMALL_VALUE) then
   !                       ! it is above freezing, so start melting any snow if present
   !                       call phase_change(p,t,qs,100.,qr,L_melt,dt/snow_melt_time)
  
positive.code:         integer::i,j
          
  !       calculate these once for every call because they are only a function of dt
          cloud2snow=exp(-1.0*snow_const*dt)
          cloud2rain=exp(-1.0*rain_const*dt)
          !$omp parallel private(i,j,t),copyin(cloud2rain,cloud2snow,snow_melt,snow_evap,rain_evap),&
          !$omp shared(p,th,pii,qv,qc,qs,qr,rain,snow,dz),&
          !$omp firstprivate(dt,nx,ny,nz)
          allocate(t(nz))
  
negative.code: !                   endif
                  else
                      ! convert cloud water to snow flakes
                      call cloud2hydrometeor(qc,qs,cloud2snow,snow_cloud_init)
                      
                  endif
              endif
  