anchor.diff:  
         LOGICAL, PARAMETER, PRIVATE:: iiwarm = .false.
         LOGICAL, PRIVATE:: is_aerosol_aware = .false.
  +      !$omp threadprivate(is_aerosol_aware)
         LOGICAL, PARAMETER, PRIVATE:: dustyIce = .true.
         LOGICAL, PARAMETER, PRIVATE:: homogIce = .true.
   
  
positive.code:       LOGICAL, OPTIONAL, INTENT(IN) :: diagflag
        INTEGER, OPTIONAL, INTENT(IN) :: do_radar_ref
        CHARACTER*256:: mp_debug
  
  !+---+
  
        i_start = its
        j_start = jts
  
negative.code:       REAL, PARAMETER, PRIVATE:: mu_g = 0.0
        REAL, PARAMETER, PRIVATE:: mu_i = 0.0
        REAL, PRIVATE:: mu_c
        !$omp threadprivate(mu_c)
  
  !..Sum of two gamma distrib for snow (Field et al. 2005).
  !.. N(D) = M2**4/M3**3 * [Kap0*exp(-M2*Lam0*D/M3)
  