anchor.diff:      end subroutine
   
   
  -
  -    subroutine cloud2hydrometeor(qc,q,conversion)
  +    subroutine cloud2hydrometeor(qc,q,conversion,qcmin)
           implicit none
           real,intent(inout) :: qc,q
  -        real,intent(in) :: conversion
  +        real,intent(in) :: conversion, qcmin
           real::delta
           
  -        delta=qc-(qc*conversion)
  +		if (qc > qcmin) then
  +	        delta=qc-(qc*conversion)
  +		else
  +			delta=0
  +		endif
  +		
           if (delta<qc) then
               qc=qc-delta
               q=q+delta
  
positive.code: !                   endif
                  else
                      ! convert cloud water to snow flakes
                      call cloud2hydrometeor(qc,qs,cloud2snow)
                      
                  endif
              endif
  
negative.code: !   dt=20.0
  !   rain=0
  !   snow=0
  !	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.)
  