anchor.diff:                                             *odzq*DT*onstep(3))
            enddo
   
  +         if (rs(kts).gt.R1*10.) &
            pptsnow = pptsnow + sed_s(kts)*DT*onstep(3)
         enddo
   
  
positive.code:             N0_c = 1.0D-18 * Nt_c*ocg1 * lamc**cce(1)
              sum1 = 0.0d0
              sumn2 = 0.0d0
              do n = 1, nbc
                 y = Dc(n)*1.0D6
                 vol = massc(n)*orho_w
                 prob = 1.0D0 - DEXP(-120.0D0*vol*5.2D-4 * Texp)
  
negative.code:          pptgraul = 0.
           pptice = 0.
           RAINNCV(i,j) = 0.
           IF ( PRESENT (snowncv) ) THEN
              SNOWNCV(i,j) = 0.
           ENDIF
           IF ( PRESENT (graupelncv) ) THEN
              GRAUPELNCV(i,j) = 0.
           ENDIF
           SR(i,j) = 0.
  
           do k = kts, kte
  