diff --git a/src/PHYEX/micro/modd_param_lima.f90 b/src/PHYEX/micro/modd_param_lima.f90
index 3b8b56a4a5cb1476f0875babf18fcd80d250a40a..61fd573f42c169d41cb4d217e604e399ea4efe13 100644
--- a/src/PHYEX/micro/modd_param_lima.f90
+++ b/src/PHYEX/micro/modd_param_lima.f90
@@ -710,7 +710,7 @@ ENDIF
!
IF(LLCHECK) THEN
CALL CHECK_NAM_VAL_CHAR(KLUOUT, 'CPRISTINE_ICE_LIMA', CPRISTINE_ICE_LIMA, &
- 'PLAT', 'COLU', 'BURO')
+ 'PLAT', 'COLU', 'BURO','YPLA','YCOL','YBUR','YDRO','YHCO','YHBU')
CALL CHECK_NAM_VAL_CHAR(KLUOUT, 'CHEVRIMED_ICE_LIMA', CHEVRIMED_ICE_LIMA, &
'GRAU', 'HAIL')
diff --git a/src/PHYEX/micro/modd_param_lima_cold.f90 b/src/PHYEX/micro/modd_param_lima_cold.f90
index 3134e583f585419a5121808a9ac0fd302a1b9e89..cc3e3adc1e5639286c37689f9569905947dc8501 100644
--- a/src/PHYEX/micro/modd_param_lima_cold.f90
+++ b/src/PHYEX/micro/modd_param_lima_cold.f90
@@ -53,12 +53,14 @@ REAL :: XLBEXS,XLBS,XNS ! Snow/agg. distribution parameters
!
REAL :: XAI,XBI,XC_I,XDI ,XF0I,XF2I,XC1I ! Cloud ice charact.
REAL :: XF0IS,XF1IS ! (large Di vent. coef.)
+REAL :: XDELTAI, XGAMMAI ! cloud ice area-diameter parameters
REAL :: XAS,XBS,XCS,XDS,XCCS,XCXS,XF0S,XF1S,XC1S ! Snow/agg. charact.
!
REAL :: XLBDAS_MIN, XLBDAS_MAX ! Max values allowed for the shape parameter of snow
REAL :: XFVELOS ! Wurtz - snow fall speed parameterizaed after Thompson 2008
REAL :: XTRANS_MP_GAMMAS ! Wurtz - change between lambda value for MP and gen. gamma
!
+REAL :: XREFFI ! constante for ice crystal effective radius for ecRad
!
!-------------------------------------------------------------------------------
!
@@ -146,6 +148,8 @@ REAL, POINTER :: XLBEXI => NULL(), &
XNS => NULL(), &
XAI => NULL(), &
XBI => NULL(), &
+ XGAMMAI => NULL(), &
+ XDELTAI => NULL(), &
XC_I => NULL(), &
XDI => NULL(), &
XF0I => NULL(), &
@@ -166,6 +170,7 @@ REAL, POINTER :: XLBEXI => NULL(), &
XLBDAS_MAX => NULL(), &
XFVELOS => NULL(), &
XTRANS_MP_GAMMAS => NULL(), &
+ XREFFI => NULL(), &
XFSEDRI => NULL(), &
XFSEDCI => NULL(), &
XFSEDRS => NULL(), &
@@ -286,6 +291,8 @@ IF(.NOT. ASSOCIATED(XLBEXI)) THEN
XNS => PARAM_LIMA_COLD%XNS
XAI => PARAM_LIMA_COLD%XAI
XBI => PARAM_LIMA_COLD%XBI
+ XGAMMAI => PARAM_LIMA_COLD%XGAMMAI
+ XDELTAI => PARAM_LIMA_COLD%XDELTAI
XC_I => PARAM_LIMA_COLD%XC_I
XDI => PARAM_LIMA_COLD%XDI
XF0I => PARAM_LIMA_COLD%XF0I
@@ -306,6 +313,7 @@ IF(.NOT. ASSOCIATED(XLBEXI)) THEN
XLBDAS_MAX => PARAM_LIMA_COLD%XLBDAS_MAX
XFVELOS => PARAM_LIMA_COLD%XFVELOS
XTRANS_MP_GAMMAS => PARAM_LIMA_COLD%XTRANS_MP_GAMMAS
+ XREFFI => PARAM_LIMA_COLD%XREFFI
XFSEDRI => PARAM_LIMA_COLD%XFSEDRI
XFSEDCI => PARAM_LIMA_COLD%XFSEDCI
XFSEDRS => PARAM_LIMA_COLD%XFSEDRS
diff --git a/src/PHYEX/micro/mode_ini_lima_cold_mixed.f90 b/src/PHYEX/micro/mode_ini_lima_cold_mixed.f90
index e8774847097e970a599f89335616b835b54b3c94..0f6c055af6426f00ebb763dc1a7bd4675aec260b 100644
--- a/src/PHYEX/micro/mode_ini_lima_cold_mixed.f90
+++ b/src/PHYEX/micro/mode_ini_lima_cold_mixed.f90
@@ -130,6 +130,8 @@ REAL :: ZBOUND_RDSF_RMIN ! XDCRLIM*Lbda_r : lower bound used in the tabulated
REAL :: ZRATE_R ! Geometrical growth of Lbda_r in the tabulated function
REAL :: ZKHI_LWM ! Coefficient of Lawson et al. (2015)
!
+REAL :: ZRHOIW ! ice density
+!
!-------------------------------------------------------------------------------
!
!
@@ -151,19 +153,73 @@ SELECT CASE (CPRISTINE_ICE_LIMA)
XBI = 2.5 ! Plates
XC_I = 747. ! Plates
XDI = 1.0 ! Plates
+ XGAMMAI = 0.096
+ XDELTAI = 1.83
XC1I = 1./XPI ! Plates
CASE('COLU')
XAI = 2.14E-3 ! Columns
XBI = 1.7 ! Columns
XC_I = 1.96E5 ! Columns
XDI = 1.585 ! Columns
+ XGAMMAI = 0.659
+ XDELTAI = 2.0
XC1I = 0.8 ! Columns
CASE('BURO')
XAI = 44.0 ! Bullet rosettes
XBI = 3.0 ! Bullet rosettes
XC_I = 4.E5 ! Bullet rosettes
XDI = 1.663 ! Bullet rosettes
+ XGAMMAI = 0.062
+ XDELTAI = 1.81
XC1I = 0.5 ! Bullet rosettes
+ CASE('YPLA')
+ XAI = 0.745 ! Plates_from Yang et al (2013)
+ XBI = 2.47 ! Plates_from Yang et al (2013)
+ XC_I = 63. ! Plates_from Yang et al (2013)
+ XDI = 0.68 ! Plates_from Yang et al (2013)
+ XGAMMAI = 0.096
+ XDELTAI = 1.83
+ XC1I = 1./XPI ! Plates_from Yang et al (2013)
+ CASE('YCOL')
+ XAI = 261.102 ! Columns_from Yang et al (2013)
+ XBI = 2.99 ! Columns_from Yang et al (2013)
+ XC_I = 671 ! Columns_from Yang et al (2013)
+ XDI = 0.62 ! Columns_from Yang et al (2013)
+ XGAMMAI = 0.659
+ XDELTAI = 2.0
+ XC1I = 0.8 ! Columns_from Yang et al (2013)
+ CASE('YBUR')
+ XAI = 1.268 ! Bullet rosettes_from Yang et al (2013)
+ XBI = 2.59 ! Bullet rosettes_from Yang et al (2013)
+ XC_I = 128 ! Bullet rosettes_from Yang et al (2013)
+ XDI = 0.72 ! Bullet rosettes_from Yang et al (2013)
+ XGAMMAI = 0.062
+ XDELTAI = 1.81
+ XC1I = 0.5 ! Bullet rosettes_from Yang et al (2013)
+ CASE('YDRO')
+ XAI = 1.268 ! Droxtals_from Yang et al (2013)
+ XBI = 2.59 ! Droxtals_from Yang et al (2013)
+ XC_I = 128 ! Droxtals_from Yang et al (2013)
+ XDI = 0.72 ! Droxtals_from Yang et al (2013)
+ XGAMMAI = 0.673
+ XDELTAI = 2.0
+ XC1I = 0.5 ! Droxtals_from Yang et al (2013)
+ CASE('YHCO')
+ XAI = 217.586 ! Hollow_Columns_from Yang et al (2013)
+ XBI = 2.99 ! Hollow_Columns_from Yang et al (2013)
+ XC_I = 641 ! Hollow_Columns_from Yang et al (2013)
+ XDI = 0.63 ! Hollow_Columns_from Yang et al (2013)
+ XGAMMAI = 0.659
+ XDELTAI = 2.0
+ XC1I = 0.8 ! Hollow_Columns_from Yang et al (2013)
+ CASE('YHBU')
+ XAI = 1.258 ! Hollow_Bullet rosettes_from Yang et al (2013)
+ XBI = 2.61 ! Hollow_Bullet rosettes_from Yang et al (2013)
+ XC_I = 147 ! Hollow_Bullet rosettes_from Yang et al (2013)
+ XDI = 0.73 ! Hollow_Bullet rosettes_from Yang et al (2013)
+ XGAMMAI = 0.061
+ XDELTAI = 1.81
+ XC1I = 0.5 ! Hollow_Bullet rosettes_from Yang et al (2013)
END SELECT
!
! Note that XCCI=N_i (a locally predicted value) and XCXI=0.0, implicitly
@@ -337,6 +393,10 @@ ZCONC_MAX = 1.E6 ! Maximal concentration for falling particules set to 1 per cc
!XLBDAG_MAX = ( ZCONC_MAX/XCCG )**(1./XCXG)
!XLBDAH_MAX = ( ZCONC_MAX/XCCH )**(1./XCXH)
!
+! constante for ecRad effective radius
+ZRHOIW = 0.917
+XREFFI = (3*XAI/(2*ZRHOIW*10**3*XGAMMAI)*MOMG(XALPHAI,XNUI,XBI)/MOMG(XALPHAI,XNUI,XDELTAI))*1E6
+!
!-------------------------------------------------------------------------------
!
!