anchor.diff:  ! !         write(*,*) "snow=",snow
   !   end do
   ! end program
  -    
  
positive.code:                 enddo
              enddo
          endif
          
          ! SEDIMENTATION for snow
          if (maxval(qs)>SMALL_VALUE) then
              fall_rate = snow_fall_rate
  
negative.code:         real,intent(inout) :: qc,q
          real,intent(in) :: conversion, qcmin
          real::delta
  
          if (qc > qcmin) then
              delta = qc-(qc*conversion)
          else
              delta = 0
          endif
  
          if (delta<qc) then
              qc = qc-delta
              q = q+delta
  