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: 		real,intent(in) :: conversion
  		real::delta
  		
  		delta=qc*conversion
  ! 		if (delta>0) then
  ! 			print*,"Converting to hydrometeor!", delta
  ! 		else
  ! 			print*, " no conversion???",qc,conversion
  ! 		endif
  		if (delta<qc) then
  			qc=qc-delta
  			q=q+delta
  
negative.code: 			excess=qv-qvsat
  ! 			temperature change if all vapor is converted
  			deltat=excess*vapor2temp
  ! 			Now calculate new saturated mixing ratio at the new hypothetical temperature
  			qvsat=sat_mr(t+deltat*0.5,p)
  			excess=qv-qvsat
  			t=t+(excess*vapor2temp)
  			qv=qv-excess
  			qc=qc+excess
  			
  ! 		if unsaturated and clouds exist, evaporate clouds
  		else if (qc>0) then
  			excess=qvsat-qv
  			if (excess<qc) then
  