anchor.diff:        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)
  
positive.code:       INTEGER, PARAMETER, PRIVATE:: ntb_i1 = 55
        INTEGER, PARAMETER, PRIVATE:: ntb_t = 9
        INTEGER, PRIVATE:: nic1, nic2, nii2, nii3, nir2, nir3, nis2, nig2, nig3
        INTEGER, PARAMETER, PRIVATE:: ntb_arc = 7
        INTEGER, PARAMETER, PRIVATE:: ntb_arw = 9
        INTEGER, PARAMETER, PRIVATE:: ntb_art = 7
  
negative.code:       LOGICAL, OPTIONAL, INTENT(IN) :: diagflag
        INTEGER, OPTIONAL, INTENT(IN) :: do_radar_ref
        CHARACTER*256:: mp_debug
        integer :: OMP_GET_NUM_THREADS, OMP_GET_THREAD_NUM
  
  !+---+
       !$OMP PARALLEL DEFAULT(PRIVATE) FIRSTPRIVATE(ids,ide,jds,jde,kds,kde,ims,ime,jms,jme,          &
       !$OMP kms,kme,its,ite,jts,jte,kts,kte,dt_in,itimestep,diagflag,has_reqc, has_reqi, has_reqs)   &
       !$omp shared(tcg_racg,tmr_racg,tcr_gacr,tmg_gacr,tnr_racg,tnr_gacr)            &
       !$OMP SHARED(RAINNCV,RAINNC,SNOWNCV,SNOWNC,GRAUPELNCV,GRAUPELNC,SR,w,th,pii,p,dz,qv,qc,        &
       !$OMP qi,qr,qs,qg,ni,nr,nc,nwfa,nifa,nwfa2d,refl_10cm,re_cloud,re_ice,re_snow)
       ! parameter list : Nt_c,TNO,rho_g,av_s,bv_s,fv_s,av_g,bv_g,EF_si,Ef_ri
  
        i_start = its
        j_start = jts
  