anchor.diff:  ! 					endif
   				else
   					! convert cloud water to snow flakes
  -					call cloud2hydrometeor(qc,qs,dt/cloud2snow_time)
  +					call cloud2hydrometeor(qc,qs,cloud2snow)
   					
   				endif
   			endif
  
positive.code: 			excess=qv-qvsat
  ! 			temperature change if all vapor is converted
  			deltat=excess*vapor2temp
  ! 			Now calculate new saturated mixing ratio at the new hypothetical tempera
  			qvsat=sat_mr(t+deltat*0.5,p)
  			excess=qv-qvsat
  			t=t+(excess*vapor2temp)
  			qv=qv-excess
  			qc=qc+excess
  			
  ! 		if unsaturated anc clouds exist, evaporate clouds
  		else if (qc>0) then
  			excess=qvsat-qv
  			if (excess<qc) then
  
negative.code: 		implicit none
  		real, intent(inout) :: p,t,qv,qc,qr,qs
  		real,intent(in)::dt
  		real :: qvsat,L_evap,L_subl,L_melt
  		
  		L_melt=LH_liquid !kJ/kg
  		L_evap=LH_vapor !kJ/kg
  		L_subl=L_melt+L_evap
  		
  		!convert cloud water to and from water vapor
  		call cloud_conversion(p,t,qv,qc,qvsat,dt)
  		! if there are no species to process we will just return
  		if ((qc+qr+qs) >SMALL_VALUE) then
  			if (qc>SMALL_VALUE) then
  				if (t>freezing_threshold) then
  					! convert cloud water to rain drops
  					call cloud2hydrometeor(qc,qr,cloud2rain)
  ! 					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)
  