anchor.diff:                  if (temperature>freezing_threshold) then
                       ! convert cloud water to rain drops
                       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(pressure,temperature,qs,100.,qr,L_melt,dt/snow_melt_time)
  -! !                         write(*,*) "Snow Melt",temperature
  -!                   endif
  +                    if (qs>SMALL_VALUE) then
  +                        ! it is above freezing, so start melting any snow if present
  +                        call phase_change(pressure,temperature,qs,100.,qr,L_melt,cloud2rain)
  +                    endif
                   else
                       ! convert cloud water to snow flakes
                       call cloud2hydrometeor(qc,qs,cloud2snow,snow_cloud_init)
  
positive.code: !           qv(i) = qv(i)/(1-qv(i))
              call mp_conversions(pressure(i),temperature(i),qv(i),qc(i),qr(i),qs(i),dt)
          enddo
  
          ! SEDIMENTATION for rain
          if (maxval(qr)>SMALL_VALUE) then
              fall_rate = rain_fall_rate
  
negative.code:                     if (qv(i)<qvsat) then
                          if (qr(i)>SMALL_VALUE) then
                              ! evaporate rain
                              call phase_change(pressure(i),temperature(i),qr(i),qvsat,qv(i),L_evap,cloud2rain/(2*nint(cfl)))
                          endif
                      endif
                  enddo
              enddo
  
          endif
  
          ! SEDIMENTATION for snow
  