anchor.diff:                  rain=rain+snowfall
               enddo
           endif
  -        do i=1,nz
  -			! convert mixing ratio back to specific humidity
  -			qv(i)=qv(i)/(qv(i)+1)
  -		enddo
  +!         do i=1,nz
  +! 			! convert mixing ratio back to specific humidity
  +! 			qv(i)=qv(i)/(qv(i)+1)
  +! 		enddo
       end subroutine mp_simple
   
   
  
positive.code:         integer::i,j
          
  !       calculate these once for every call because they are only a function of dt
          cloud2snow=exp(-1.0*snow_const*dt)
          cloud2rain=exp(-1.0*rain_const*dt)
          !$omp parallel private(i,j,t),copyin(cloud2rain,cloud2snow,snow_melt,snow_evap,rain_evap),&
          !$omp shared(p,th,pii,qv,qc,qs,qr,rain,snow,dz),&
          !$omp firstprivate(dt,nx,ny,nz)
          allocate(t(nz))
  
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.)
  