anchor.diff:  							rain(i,j),snow(i,j),&
   							dt,dz(i,:,j),nz,((i==(nx/2+20)).and.(j==2)))
   				th(i,:,j)=t/pii(i,:,j)
  -! 				print*, maxval(t)
  -! 				if (minval(qc(i,:,j))<0) then
  -! 					print*,i,j,minval(qc(i,:,j)), "cloud"
  -! 				endif
  -! 				if (minval(qv(i,:,j))<0) then
  -! 					print*,i,j,minval(qv(i,:,j)), "vapor"
  -! 				endif
  -! 				if (minval(qr(i,:,j))<0) then
  -! 					print*,i,j,minval(qr(i,:,j)), "rain"
  -! 				endif
  -! 				if (minval(qs(i,:,j))<0) then
  -! 					print*,i,j,minval(qs(i,:,j)), "snow"
  -! 				endif
   				
   			enddo
   		enddo
  
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: ! 					endif
  				else
  					! convert cloud water to snow flakes
  					call cloud2hydrometeor(qc,qs,cloud2snow)
  					
  				endif
  			endif
  