anchor.diff:                  excess = excess*(-1) !DEBUG
               endif
           enddo
  -        if (iteration==5) then
  +        if (iteration==15) then
               qv = sat_mr(pre_t,pressure)
               temperature = pre_t
               qc = pre_qc
           endif
  -        ! if (temperature<170) then
  -        !     temperature = 170
  -        ! else if (temperature>350) then
  -        !     temperature = 350
  -        ! endif
   
           qc = max(qc,0.)
   
  
positive.code:                     if (qv(i)<qvsat) then
                          if (qr(i)>SMALL_VALUE) then
                              ! evaporate rain
                              call phase_change(pressure(i),temperature(i),qr(i),qvsat,qv(i),L_evap,1.0) !cloud2rain)
                          endif
                      endif
                  enddo
              enddo
          endif
  
          ! SEDIMENTATION for snow
  
negative.code:             a = 17.2693882
              b = 35.86
          endif
  
          e_s = 610.78 * exp(a * (temperature - 273.16) / (temperature - b)) !(Pa)
  
          ! alternate formulations
  