anchor.diff:  !   dt=20.0
   !   rain=0
   !   snow=0
  -!	cloud2snow=exp(-1.0*snow_const*dt)
  -!	cloud2rain=exp(-1.0*rain_const*dt)
  +!	cloud2snow=exp(-1.0*snow_formation_time_const*dt)
  +!	cloud2rain=exp(-1.0*rain_formation_time_const*dt)
   !   
   !   do i=1,10
   !       call mp_simple(p,t,qv,qc,qr,qs,rain,snow,dt,dx2,dz,nz,.True.)
  
positive.code:                 rain=rain+snowfall
              enddo
          endif
          do i=1,nz
  			! convert mixing ratio back to specific humidity
  			qv(i)=qv(i)/(qv(i)+1)
  		enddo
      end subroutine mp_simple
  
  
  
negative.code: !                   endif
                  else
                      ! convert cloud water to snow flakes
                      call cloud2hydrometeor(qc,qs,cloud2snow,snow_cloud_init)
                      
                  endif
              endif
  