do k=1,km
                          tem1 = grav/rd
                          wc_in_g_m3 = 1000.0_r_kind * tem1 * cloud(n)%water_content(k) * (pr(k)/dp(k)) &
                             /tp(k) * (one + fv * qh(k))
                          x = log10(wc_in_g_m3)
                          if(n .eq. 3)  then
                             cloud(n)%effective_radius(k) =  &
                                  7.934_r_kind*x*x*x + 90.858_r_kind*x*x + 387.807_r_kind*x +679.939_r_kind
                                  cloud(n)%effective_radius(k)  = max(cloud(n)%effective_radius(k), 100.0_r_kind)
                          else if(n .eq. 4)  then
                             cloud(n)%effective_radius(k) =  &
               !                   5.7135_r_kind*x*x*x + 56.939_r_kind*x*x + 268.224_r_kind*x + 607.407_r_kind !Liu DDA_type0-hexl
               !                  4.597_r_kind*x*x*x + 45.813_r_kind*x*x + 215.811_r_kind*x + 488.7155_r_kind !Liu DDA_type1-hexs
               !                   4.037_r_kind*x*x*x + 40.231_r_kind*x*x + 189.518_r_kind*x + 429.713_r_kind !Liu DDA_type2-hexb
               !                   4.549_r_kind*x*x*x + 45.337_r_kind*x*x + 213.5704_r_kind*x + 483.641_r_kind !Liu DDA_type3-hexf
               !                   5.864_r_kind*x*x*x + 58.438_r_kind*x*x + 275.287_r_kind*x + 623.4_r_kind !Liu DDA_type4-hexp
                                 9.33_r_kind*x*x*x +  84.779_r_kind*x*x + 351.1345_r_kind*x + 691.391_r_kind !Liu DDA_type5-ros3
               !                   10.612_r_kind*x*x*x + 93.341_r_kind*x*x + 368.987_r_kind*x + 690.099_r_kind !Liu DDA_type6-ros4
               !                   10.243_r_kind*x*x*x + 90.093_r_kind*x*x + 356.149_r_kind*x + 666.079_r_kind !Liu DDA_type7-ros5
               !                   9.7725_r_kind*x*x*x + 86.063_r_kind*x*x + 340.832_r_kind*x + 638.7114_r_kind !Liu DDA_type8-ros6
               !                   12.4278_r_kind*x*x*x + 102.867_r_kind*x*x + 372.06_r_kind*x + 629.84_r_kind !Liu DDA_type9-sstr
               !                   12.636_r_kind*x*x*x + 108.2445_r_kind*x*x + 411.809_r_kind*x + 737.962_r_kind !Liu DDA_type10-sden
               !                  5.3503 _r_kind*x*x*x + 53.319_r_kind*x*x + 251.1721_r_kind*x + 568.792_r_kind !Liu DDA_mie
                                cloud(n)%effective_radius(k)  = max(cloud(n)%effective_radius(k), 100.0_r_kind)
                          endif
                       enddo