From 795df703faca4ab7dea6c3e709f3388b5f9d31f4 Mon Sep 17 00:00:00 2001
From: Marie TAUFOUR <taum@nuwa.aerologie.net>
Date: Wed, 20 Dec 2023 12:11:49 +0100
Subject: [PATCH] M. Taufour 20/12/2023: Add optic propreties of ice crystals
 from database of Yang et al. 2013 in agreement with the PSD of LIMA alpha3
 and nu3 and the parameters of the crystals shapes

---
 .../data/YBUR_ice_scattering_rrtm.nc          | Bin 0 -> 16067 bytes
 .../data/YCOL_ice_scattering_rrtm.nc          | Bin 0 -> 16067 bytes
 .../data/YDRO_ice_scattering_rrtm.nc          | Bin 0 -> 16067 bytes
 .../data/YHBU_ice_scattering_rrtm.nc          | Bin 0 -> 16067 bytes
 .../data/YHCO_ice_scattering_rrtm.nc          | Bin 0 -> 16067 bytes
 .../data/YPLA_ice_scattering_rrtm.nc          | Bin 0 -> 16067 bytes
 .../ecrad-1.4.0_mnh/ifs/radiation_setup.F90   |   4 +-
 .../radiation/radiation_cloud_optics.F90      |  49 +++++-
 .../radiation/radiation_config.F90            |  35 +++--
 .../radiation/radiation_ice_optics_shapes.F90 | 142 ++++++++++++++++++
 10 files changed, 206 insertions(+), 24 deletions(-)
 create mode 100644 src/LIB/RAD/ecrad-1.4.0_mnh/data/YBUR_ice_scattering_rrtm.nc
 create mode 100644 src/LIB/RAD/ecrad-1.4.0_mnh/data/YCOL_ice_scattering_rrtm.nc
 create mode 100644 src/LIB/RAD/ecrad-1.4.0_mnh/data/YDRO_ice_scattering_rrtm.nc
 create mode 100644 src/LIB/RAD/ecrad-1.4.0_mnh/data/YHBU_ice_scattering_rrtm.nc
 create mode 100644 src/LIB/RAD/ecrad-1.4.0_mnh/data/YHCO_ice_scattering_rrtm.nc
 create mode 100644 src/LIB/RAD/ecrad-1.4.0_mnh/data/YPLA_ice_scattering_rrtm.nc
 create mode 100644 src/LIB/RAD/ecrad-1.4.0_mnh/radiation/radiation_ice_optics_shapes.F90

diff --git a/src/LIB/RAD/ecrad-1.4.0_mnh/data/YBUR_ice_scattering_rrtm.nc b/src/LIB/RAD/ecrad-1.4.0_mnh/data/YBUR_ice_scattering_rrtm.nc
new file mode 100644
index 0000000000000000000000000000000000000000..cfcaf7ff2ecada6346d89231a840851ffa998bc8
GIT binary patch
literal 16067
zcmeHt2UHZ<)^3xUCg&&!Dxd-q3@8YE_Aa6#IZ4hLnkIBNDquh{=bUpI<CsT9Fk-+o
z<}`>A;|Qi<1paD)nfuSZZ@v4!_2#bk*0NUb-c{%9I<@ONd!MR0=bP%|;iaZ*tSrmp
zak{ql&h4aoM~Opch4;iCRt5NY1k0-KlkNGfOpe=G)pg2WwX!TnMrxv|DgAFykp6Qx
z3p*5tBj2^d?Izumk^g4Q?(vcZscC8Ph2jCS9VT4s4uj4^jzN50^1%2!L0){aKv0mv
z{&(1m<vQ#+YMuH|<>r_k@y&g;Au=+$(sEfD85#9XgDx$nEA?2bhgqj{<{Sea$Do|G
zQuuZ+a>^a|Qq}4Eyw3CWOun5UC(rpS23^OhHd0eL4wu7|nlgU;kPwazTO`FsdS4#C
zBR%~?d;|RBLm0lFCr34kg@ippq2mO9{z%Vn_QN{uLo+f9&Dmd~d0uJ}<KbnVl9`*D
zSC}ctGfyfmGS3qf{aGbSEXYhu%1N~|@7x>*a-<H)vhAyLRHY8Z`==Hqr=$(^$jr@5
z;b^k)&~c0A-<)Hw$!=D8rF;JwL(h&%PM7|7F%0aChyPzObaMWmVb~Gf*VoQ9?!a)U
zjEsd8I$1UlyW+14!@nMPVru>`;!;@}v2xZOz4hx@p0>12IT?<jRC73n8XdYTLS_Hg
z5&G@oavWVYAlL-y3J6A;VVQY}MX81VJuKMt>B`40ExnZVHdSLfVpoHWZnmjEL)*_c
zB$Q*sT1cbbfEDSuBTdAv7P|B|usR1NY1uy>OOsq>w$qjX#%Dy=?RKTv5VlUbspd$N
zO^&0%agcft?CIqh?CH;V#((p`iEUEqK}up#Vqt30pRW9oLmb|JO99ryZyvDcbm>3F
z8tLWO?d#~s@_*%#Zs!UAIgh$L+$!*nl6urJWLfm-@YAP5b>%mE)D;>s|J%?|?-+~X
z{h7c5L4ImMQ6?J{(m|{zE7EQfr4~p>Ah-Uu?4ZOBCvF4nhS)jUW~8JIa&xqE>~H5-
z`r>#gYxsTpe*pn+Ccy7ow0wI&{~6nziN|@@>7=x4rQ2W2yB@b)<z3I~j=|_pd&Axx
z=XZ|yf3(w8-t{8U)z0-k9rj)S37vyj*YS7#9f9u%d`I9r0^bq%j=*;Wz9aA*f$s=>
zN8mdG-x2tZz`u{c487w_OtXZk^&HG>daKU-BzjHDE}o+~DZ5EuV<n|+m`Rof<LPr=
zF5UkUO>0$0lF(upIXxIiEw8L;y?##$us5cDS~_%Lg(|ro;L-_S8Ct*Q3tz#<7M?kM
zk&<QwjkX^6^^h5wR-59N$K9}JiZLbz7@^t45Isg3VAdFYC|2vC%s>xb<+{lHsDn`w
z9c<<5Ky9@)yoYJSs8I`(muR8CzZQaag<|5eIM`JtgT|yI`qxb4Zq0_#+Z^=#B^PbW
z@?hX7fW<ZeT-5TRo{*1~oANR2aXy9{2yu!L;$fx`&!!2{zL}L8AyV!L5%pGxuvv$g
zIcGwdUp^Qyp2K7qE#=!ZVC!K@x7$EQpXbpm<4F|qw2-PPkuF6A(be4yJ-RuVE(rV4
zI(G{?tYku^+jPlSMx9ie73pT0ER9;&!r$ufo}Vro4*TwYu!lG7RJ@^j(F42J5cV7L
zkv<-%=eXnPAy<s99Ex)ZE|B+gMzrq`c!Ukc&Ri#O=MDn5Y9O>f48Xkp0}v<akGxZk
z_^j^;j(q~QJW7Ima2lS@NyjuM6K?ymka{#5ic5147nTdZ_qlM8<Y8k|9!4e$aJWu@
z5l;Cyy)qy2GC~B!2yu6Z5MMqD5#UsS{rLq*-&BBi+8Y^P9Y1Cc-<rvWEF&Lqk0vfW
zO3Rxz($9Mr(3lI8DQ=R8ZbT+i<&Y4%kmyOg7-!1AYDev*mUOwFDNR14N1CnbwEZEE
z4l{C;H1I3``N1ZBh-WEs){eoztD|A<SAu6|5}bV}!n;ew5L6W*cYh)3_7uSPpb(?a
z<iq!Y0NsAi!&<XE1P;$d?dTkw*`1Boud=YlCJVBHOf1`<0WVGl@^ePO;`2y|e+WYK
z_)wgF6^^ydk+`!u3jK6q5FQ$f8yRt6V&ZYgDFL_I67cR|A|@mzVQEVeE>28_VtX>y
z=B41#^%R^Pkczq)shC-xij4`C%(-9Fm~C2SOm(mvb93lz%KCVO;<Gl=wU6`Zq4{LW
zUQkSj)+Ukn{9rm;=Ry5%51}4WcC=T^lGa=>A%l25lKnxQzRcrMRjVx3?ERhpGt<Ox
z&Mm^o`T}fRoQsqBSzvaiq4r!d4m2j<Ek6z~uSTPJ+9)hfi@?LgFdP^ig2_9A5cN3#
z%HjSvan=vlhmXL{2fmnD=7TsJZ@h2t#MDh5_<5EO>P!Q$|4s;6Y9g@hbrhz|i^cSv
z@pxUD2sOoIWR6OK#>7<2o}Gqi6Vl-vlL3Y9nRs<R6A5WqnAV(ywxVo2c%6+8F*%re
zCI?^o<|3wUIWu@&4D+U?7xQ|bER$Mtj~dS&r4fTRQQxBr=-Tngw0N<ImW)rP9iu|%
z{TNSjEO4f2PwZ&^Y)jhUZc3=tBP$gRYJAC~e*SV~<Mowqa<+*-_h>18zcCi3=f*(y
z)M)g)UIK%s64=N}@SwLCiJ>Aatt>`#O%d#6im=?P5cYEmu)JP~RRe|K%*@BQMgi{n
z3b1lV9=hx1p>1p~nqK5U+rR^ZC;4FeFh5lG4ZxsnLAYEMf@kBxu$mKroqmy6RyYb5
z$422@dNg>0VzA<Q49-`?!c8j<rE}t7tPl@HVLTq)h=<Og1T31JfTd3p;5~mIbFn#&
z@zb(orZmVij&=`dKwuRO$lXjkcP}KGR!OG6ipg$E3e`n~QvE0|;yJm{zIpbv#>R@)
zUNEIoKKkS}UW58aD$!4S<)}igm7gwY=Koq>ieHwG#lW#+pj|Q=ql!!L+jt2gmWwg$
zxCnP%7sISqF>DfxFl|#IUN#lrj%NXGuM$Gonh$<rKHi=az;Lht*pi1IO!Ck)KNs87
zb74D}LGwayeCRs@hfMs@ZA~Bqhk`L_Y$ztlhod<*0+*&lf;V{-=EX!oS1%gX2cogY
zD+Z>uF*xcLi*Z|Hp`sRtn#?%VSI2>L;?XTJ9xEOmU~WaFGqZMEGqWZrGAoN7(VOZT
zigMgWMXwjrmbj_(xqUQ6`li#Njp6hu$%nk;T&c{~k;0bsrhWQmWFa!3f`ytiFI<^q
zR?Ab|p*H>$jgS278|64~tQ?7(%HcPg)sk|^C6$Bx%5l@F9Cba)VX0P*lV8d(<w+UR
zuaqJGXc-P}DTC+2GIX0z275sniet)f-?t1wgUfKqq6|gaW!QI%4`#b3KGpf)YU~J{
zpFa{B83BlI3B>f+VBDD&0>cTRcrh{zhrfiucUm~|co7H{N1*yn1hzXy;>5T}OsI`S
zu<|I3@g9Zdaid^za1<7W?qcMkl9(eidNVaWxlD@lFQjt)IL%tQh1PCfM61H4P?ARp
z-Mo-WSrK8h+|`=~afVX5w*#fsSy4+&ce1&tPhY=ikk%O`>T4uVmB(B8>b*bk@BdZ~
z!%yYd)>Mx4cjd5nU5;lj$}#?FIaC_TQT4DKtL~L!<?V9pVDYQHR*uHY<<P7vhwJ%r
z2+x*d`>Argttm(U<K++?Eyuyb<@kJ{9C|DGs4?=y;%+`jUp5?{&iUczZT{GcKm^QT
z$E1zH_%byFO1_~meie#nQ5cNgg<(!)ICA!f<Edf<P6b7vY(WH~uSQ^+QY5~*MZ!uD
z$&N$YnKvWinVs7$nbXhY8TG;klsln{ZZFtO{jV;hhZ`!1V<w@kKc~>kvQX+h#*1DD
zxscos_M}Eur1Z#?hQ#Vq<^l~`nyExqC*^2F_g23D>}GyzdpX)VrJdFPDAMv4R(@sW
z0V^k2*}zH#D=Dmavtq-FCM%!X%5k5SYF5^>Qo%}aFdqqvJutM=8`j#x;TGlx{l5NK
zzA*rL4S_Ix6$G~n!8p(2;$$0&$b+FswhzNEi^H(1ISk!A!|{1$IO;BkqfI3OXP5}6
zN+QsIYXmBX?_eI=B`~>zt(dts3e12959sYLRkY;EW?G`Vh%)Y1Qng4z!;Mm@{ctGd
z?)D=2i7xc@mOU-cvZ5gh-D!TQKGp8gppz9!wBWuR@%&r)caAjkIh=9$yJYO8VK8Z7
z@Zqx;#?SDD|Gp8pBN>UB_x@OG9|*glLC7=-M)lob6pRhQ9d0N@Wud5f7K%6u!>svX
zh`b+$((d6%i4KS5f^Za_3&*72!|^e2J2TrYo-y8P$!xqM&nz1Ifb8d2QN_~Dly-3;
z_1aZQQ-dXRNiLPNR)tdOOfOm&<3f?U?TN!!(V>T?)Do#rV`gj6lN2R#Ixa^`3|sk_
z@%LiFs{yV)B5eNX2?wDMewZ^H%fkI|{Q5|Y(+@y@i$L7t24U2hAWX{)#+r}8$j=Lb
z%dHTo3<$-Ficn}=3dK9kFr@j1p}aB-`p3eslFcv6pdHM7*97K_qZM=9K!HL10~++K
zibg-)OtyxLsOniI9i1$pY|~WAJ{n5vfAS*#MJ{Cg-kxSpvLZvX?sRs6K5aRzK_Ls3
z$n~WhMW?j#_tgEpn5evP!xKwB410S(X@?h1Huxashv9Ir_CuR+Bo<Bbhh<g(T>1v$
z-GxBNL<T|WY7pQOjH8Q!(Y-MksRKfwDh`3){t&pm3&By#P*{u##j(O2OfSC##@pYD
zDX>yt#=m|*&sj{aecDVszeQwaK81dMETPtfRGRZ5lv=NPQR)g8@_uPgS438nqt=}s
zROnNmBN}vXt`fPuk)u_Gt$dS)zZa98CC+$v-xb=AfWircZ8tp-vfB$IhWOyjBwt|p
zaA-~(0kO9q#=Z4J^VpHdVc&bZCHX^8?GIJG0K7;HK*5#(Tx<-$7pp+*ObA4;#ev9K
zbB!t8K88_z;=t^>qQ=CnZKSHYv*d5Ki@w$`r!PBZ(94K%WN4j3^>?GFkES2BGoTN;
z1L^z_Yl<rFLA_TR(X;8=G{2`Rt@Y&6>PH+J^~-0z>h(ah9dgH<GyM_%-WDH&dn3fC
z7u-&p;p;;a1eO?MjK~lM2la8TzaE}H(!stX+PHpL3!Gb;C{x!&d5i|$9aD#igF2?~
zR)de58t8^9&d9}J*rs^A^h$;PlT17}$c3;(06kwJ`m8R1^~yqgPA$R_m11CFG5i!o
zs7({$+(8j+T19v<M2wj^ViYeG<5IO4OCO6di+z{W>MlVCL5K4XM6u2V@^7pW>G$=&
zOLB_0?f4r18%a(b^p2d)CG(8VumArMz3bchei!GcvjM4UjbBmxCxfKEZ%37Hmqphy
ze<#6j<u{5OhZFfv^|u{-eWdODy__~DP9Rg$^3Q$OMXKu3|8C#C|4;Q@-*5Cr4(I=<
z@4j&`{zcy{-LtqtiaFO<+8r`~^p#>#VqOZ%t?c4<bdfK+3_8h|IpQu7#<9z{Z2tDn
zW65;QZg~!Wbe@IOX(=<VEAhHI@2|)Ib7tJce$uu%oMBoWx{Lnww=(0T87A$ij%5AV
znJHNv9<X$_uANTZV~0Y+KXzUh(WayGa>RcSXF7?oUBsGyO^nSE&+Uvb2ji>vzV+q*
zRf3-M>fdNYybkxJbRK8l&`9Z({z}jr)?xA=An0}7bLRy)Z^xLv{8t|R=LmXT9(57)
zcKF`kWqqK6*>ot4v5X$VJif2P81{cpSClVMK=CfB$X!MkpHC&DvJ#pzEQKy!3#O;O
z9#qL0LMCPX=wpfn_0Ke+RU`H2#|Cv8EU!erERv&ay;lCkpPKoXKlDP9s4wQU_QAxd
zmJnI>fKpFW-1*4}FV5&AYltq6nP}r(ktX~cHQ+Q%9o7rgU~i>{OV3r&enAx<PO9R;
z301s4u8LC!RIzE3D$M4pBK=b!xGvEcE|&ni;v_`(NJYSxbj<b1f>%p6`h@2qxq_W)
zoW@SMq~{~kREW#Fgy?NqfXaylFt}2HMs;?IhAG5!K_Q$M6+-54A(Zd1`(CZ|M%<S<
z3~T1y-foPKstWVC;v?zxxk<-Fhp6=0N~%h!q>AP&s>%)^+;*Xoy*4C2$DEEh8&aQY
zEvhS2rHcECl)7DxTGHD2k^4ULd%6C~uk0}#8^S}O{~{F2tJzn>wjd}Q1>m#V2&{SQ
z3Fo`OB5P-?%I=T0!M2#7-W%5(dLeANInL^v!8oitWR<&P^cPd`Kbc})qbVXzo8sIg
zQ%rbb3oUtXcydPKK#U(Imjz(Xv``FN9*G3gXmEdz!Oy4TAYQ>v<FT0aaZSOHJE=G@
zEFG!B41@_YG2S-|+Hx#@o3fGVoP#knIao6+7t6*^W6JK0V;V+>G1+TvnARQ&%*Gv;
z>D}(lG~Is|Rm~`+IMXaTxG9F_)dUgtdDGLIZe(OIkglz@rifiVs6Nw(mOa%XhaM_a
zsiHu6X6<~5`$vA>D+k=$Gyr@T2Lx2uAohJvl-3y_#6b-npA@jhKn~6e<WRd`0h3L5
z$X8NE^-vXTI--KgA5~Bvtb#A+l@aBxj7uk#ke;Rlr#?!^QcyzuYeH$yAXNMoiF&PA
zc;zPI9HpYUPX^SIv#`5o4t5mh!Xi5lJ4^*o+9bdWoqQY-<ioToAGg?P=?ktx1m+7d
zZY4WKe_DvcFNCn;6`(BjM`m=-N@k{~C$m&hn^_t6h*S^#L>`CMQ@s5G;wem_4DB&=
zCM$!o??sSjfG>HhbS2JIM>=z*Hw_$OMjvJvQ1d2DicM9f;oIfu{gpO;xYbAgUJWDM
zyK01*i+b4qLKR6p<e_}G-918s1K+-~P)(A@lrIX{^N@=*!xbTtC}P1pMJ(K@h>zD4
z;n}8$!?rwVNAobbl7|U<*!?>^xVQ08v(pbtUyZ`{SBaQ?Ee#gRS-AEg8;!kk5&1b6
zcW36oza<X?M+&fHwg4wj3((ReK)q2uIIj87jmpQ^qI^uBo{w#-^AWKpAIdds`Q?0M
zyqL`#`<lxP4RvIy0+kt!xCW{nRYS$Mw^HaYi|P8rskFOz3|**6r+SSDdU4!`Y$9E0
z>=;LS)6kpLBF%`u+kpC?)1-lulxh7rc~Vhq=l6>I$d75LcbD1wyZi9T94N$bFzE;f
zo1V#FXQUjq?v%&-vkI87g^THNirD>05#vjE82XwA<6tEO>`+2UyAq!HDr4tNWn8|X
z40{C?6uYS4LADCU*LlNlVFG#$PDhqX7ADMLr<RLy;BK6Y;|p_fon5;4;+zNdqCBLm
z&x81U9y^VahfQ1o`gIqeq`v^YJOo%1B7j$-0CBkjm`eoMHC_Of%pj)8%YixYs|K^g
zvYBdUU8Na=4w0wN2Ffv9KxTU;)1-P4NzNowdv*wY3h<<PD$cao&yGG`wIIKA6N+xo
zB@=UXYLekmOpYw+FK*#4w|~#~w%X%9px~zarJ--!V@g`w?~a$jXm5GE(dEL<RT1S)
zim;92!D|^0_Yd>HoaZ6-IuAcQ;32S)htpqpNLE!sT`whsxGTXnS_u^;N{C;ggbC4M
zpvZI>u}e24W%<}Rr~nf(3Zd&>1eNPWxI)EPyPz0VSBf#IwHUj)i?GE>gty)zNZ9Fc
zie*1@l0>j#KYQxaMQF@m>oP@1&lVxkm0?^;80Pe){!G=4?o56fkFovwl<vN$rlor8
z$t-aOoy-zbLU26IIX#lD%os++TqjC2=|}uNmbCerDe*t*)53|G<T6^BuHKNRpuX+=
zrs9wM<Yl&qTR9L<uM9zZs571moG|FA15PdN3$-;pU}<EGK}|YfR%@W2ni>+Ws9-~i
z3XY6YK@_X!#;9PzKozLJRYvYSW#rf@BmcA#GE$VFu`C)-Y7-FHBMmA1Y^Yt%LwmLm
z!#xV&y{ibf*kz;V`$W*m6XPd231*Z^pxhuqqEiVX%1SWpXbH|Wl_0L?Xz=|;V~}t(
zF3lZ{h@GP`Gjb~9YF^4*ML6^Jbzeqp8i(mK`V`ImvXX2*RFYp$3Ca6sk?iGoQpgLV
zlU~CqYKuGZX#llu??ZvEW^}z>pH|-1Aa8bhAmF_m4Oz}kne_a?=bL!KV67J>Z}Y?n
z0{rI=fPSnc^!FM<WUPf_E$T4V)<iP9%!zI~Sm&yPWqI0|tD=Po#TxMWO%2a$RI&Ds
z3anUcgh!MS8>Ng5vdZ9A$KaR1cw9W4jH#<Lu!vofY0_rL4hta${9FLD-bMJePcfdn
zD#o^nB4oFT;Fv1L>|<iwSC&8~K!P6AB~YrCK=f7uX9IR=WoQZ3M3tbaayz5&WICfX
z&5t=eOM~%wafSA)ZzHt@OKAAZS(KDDnT*7vsbpq4Ex8d!_mjNHJk^=%%IzqwpCyI=
zZbD@%_2}3Obvk{DM+xR~)UEz^{`-U`ey*%HMg%da&US%wsx>Od8set0I$qYPB0Eq6
zrg2&@+@yuV5>4D%uMWRks>mu(LE{8vC^RZz_*f<MG*^OFEf0o;JoK>S0S$^M*sh2p
zO|e)wI0*^nX^@+gh4fx|xbh+&8yX5=Jg^AwyA@++Z84HUL{P64LDp7`3FF1=^tBj`
zTnS8EB=Am^;Q3Su9JflKcS?dm_au17;`1WlIMXz9I<u;FFmp%PO75n$biU7zB)4`0
z36vI-uijL;e@j9QFOx}Z5kf0sJ?U+wGc8rLr+q6dX_lQSUD&Kg{a&$4e>ZqE)=Q2~
z4FAf1e4~kfwAvLR7K5N~X@LuOHBeEhfY!5I_{>#?kB%xPxv5}&sS+M2@Nk_SUt-^L
z!M(zT=2k9linvJi=Ax)K7jKNXILFFbGcF!ja-m?&MO#TCoKK`fSeuPw8wCg)SAfE(
zA~=mLX5&o+jq@VJI*TEkFNVf_F{<<=7|KY{JzawIi4v?@DuLAx5_lb!pozuf=p_jr
z-el_@NZ{FYhWYu#R3=Q{fw|W6H+rgfiM$>jrm03->5an@^1U~OG&V{o;%o|)&j=;I
zZC;claUqAj_GCQ3ib5`$(kFj?s%NKQg*i&J=bRix^lRlWUD?dPG0g>Aw>iM>pb187
zQbuMN7s544D0`^_#c?WdyrzU<>~pqtks_SE6;Y$3h|Jqu+*{5?LoOG#p<LAZaxun>
zi!2W=#C$G3xN*_w!bQ@MROr}c;gUigrYsXen^T0-!eZpdh@dGa#v3*!G@goK*I$Cj
zTnTcQNuYL&&6(RQ4sRrw&?bR~S_w`Wl|b3D1X}hb_&m4-S5Sh<4o?`1q-jjc*Y1o#
zLn{?Iy&?z8D>TFUM|wPNIX$+XL5)S@NN-&ZDKtcp(uEQ1dl|?{Z6KWuv8Lw5J;?Ha
z5&g7So90=nQVBb?cAuSE+uZP(A9o@VgVIA$``HZ}tew#OX)t`o4ubPEd&JnVxpdhI
z8<H&$)zSl>XPe<-mMQ$!b;CPTV@$0w#G)btymHq^Sr0u-?yd{*Kpi}c*GBhMTCi}B
zg+@g*ZXAt+>gE(YxtsytdD(coITz<6*yrnke5||5zLWh_fS}Ssc=jrS&5uQhv@1r#
zvSOIF6hqcugpo@`7=KrU$9iI9ju2zj7<S6dtdWV?J%%|EWXr4<8!|=ht#qpT3SHCL
zPX~3@QN@Y5B&T0N_oo+-`ke&&sXU0Z+ZgKGe+V6a)Q=)hS<tFACKTkaN2|80ll&?k
z9Z-{_x);Cm$M+4v&l88DMY|u?hudRUn>AKU>;<DZGu*T@#o5)y$TK!XV~{@N6!p+_
zMF&yqwJ~d{7FHh7#FBOm)JAKd_<}lCjZjDI9W|^dP{TZZHK=mpz+0Dqt?V~W<jNf6
zc?)pLTZor;3b5jHAuiPxVIdng5nqdOkA1%0JuO0{rWlF=VwOfAM(ZvyWNwHN@<9w=
z4GA7wNpO0o1RnkpydLnJxxYP+Ig;YQ1e`NsZf*Wd1FP#O^XxvF@nsEJ`plu)VH0WV
zMj`26iKl77fi(34(Ua)G>=JQ5I&W@487d~Eds>%l`l-{%-aJZPElcaV{mz%|6~sPg
z2Vwg)Tc{h_;nO24>~S?me5@(xs4=G08{o%FddQ5|#RUx=?75<a>MBh{-`7Bki3VgO
z>bUVj4a>6G@4{A9BrjLRb`SP@@Use<lVY)BaU$v`XQJnsT*M{j!#$(`@<$4>@Nf}a
z*l}fxfe8ENiIDt3gla~Nv`J#<9u?zMqZo3U5@<U}AU8sSGl>!y7fT?RB*EKx63h|T
zGd{_2%%mGb8QWAt=9betI*&S1dApaUyjw+w17?%m!wIymB%j`A$B{v=0Mcydlj{^G
zdYWrXCkuPgs_bs0?WjxhHmFg-Rz)(mk)@YR3;*oB7!0rQK*l%*Eqd<wZP!4&>*D}j
zPg|U@>;tuKmPlXM6TiEe!&<#N8v2-^uG|>W3P$k$(Ewk@>BBZx4;Kq{@p7&X5>9Dj
zF;5#Av06yZ3x<C~I5OYF<DN|_#w%nZ=VdnLvhVfTQv}#{F&{&I7GlfR0`!e7#FLMO
zSSu-ld1Ddg4ll;I)y4SSSd5J}>@$gdZ;M(c!Vl*}C}#01b3e<dwU#mEE_}w&Mu%DJ
z&`fB)M6>(vqZ7evh?_KrYGNkR^Z+3hhQ(7+TL4uYC-So&Oq<$m$?Rh<T2|kUbXMrn
z7p^)Dk>}CZGFhsA*utO7^@Wo%ALT>cu(om#N)sI+v%DW3ds$=k!QQBhvqU%7UU;9_
z13iwJfzR&_%|;V!*wqb{bB&=o!w3~C3{hBZfT<t#QS7Xb_-T6B`CJ$C%>8iSco_C%
z#^AuSM0o8>#l?gSEX~P+Lt8eUyXL}?@^IW(fbc6U9U(g(wZG=WKUK))q7XUkoPlw1
Y0i2E&VCdHZmYP%u+w4N9E-%D?0r|V{i2wiq

literal 0
HcmV?d00001

diff --git a/src/LIB/RAD/ecrad-1.4.0_mnh/data/YCOL_ice_scattering_rrtm.nc b/src/LIB/RAD/ecrad-1.4.0_mnh/data/YCOL_ice_scattering_rrtm.nc
new file mode 100644
index 0000000000000000000000000000000000000000..b51cdc499d6352c8d79a690fce41bef49eb95706
GIT binary patch
literal 16067
zcmeI22UHZ<*6$k-Y*GVFXp#z&qaZ;*<qx|^O%5uE0m->ZO(QbU=8Twg&SA_sVa_>+
z8OQ9Ha~c)>P76JAzkBan?|yH+d)IqwS&JX1tLjwMKL4}#sj5?ZR76Oau9mr$Dxc3&
zF5g;~1A)z3JX&^y4_P%gE+QmRRcDiG`|sP-G+K5k5Bd92OO>ahsG_g0_<u)T@t?<=
z(p>O(t&~4BEEK;~T79p~e(^I)vU76MDyqCxo2zIzH&<v`<%!ZHnciuV5=mNSNlAGY
z``=u@O0Bs*Pq(G~tLW;PoBZ8;T|X5SBgJ-A6%`dhO9f@AW~8v#z`9+FaaKGLpC^*D
zTI%21qNdf{mc|y_C$=1KZ2A2FwN|{p*Fd?h<D#gl#?#>O6;-+D;G`s;A={+TM$xVi
zaZG4zk~l6lEr|<`4&~{jGD+A0>dgnl2FHYcuOHM>zgu2@g%$ftw~}N_xsWicto-6^
zNkx8%#44jwY9%R={<TY%UY?(xQIzd&)v|bc^ArZDvhxdgItqi*VzZ^0SvlSz`NjEJ
zJbmU5%}uoWZXCNNTkIvZ8~;(m(B_>yW%;ja7~kR#|G#SJ>;FH~u-UsWE|_56tYK#r
z6<dXLs>~6U{-@ON@7G&2wfZN1sip8(HRtBudifJys5qyZ3eQwg@_43t&85;qRsYu>
z`u*Z+JR{~1EI^bF!6_mvUy?4(uK4e6!NNxwkIGVUDMd4NdN=#79`kN?s=r)2TAb94
zXU1wMyk5jmj5~@zRGvXueiy4{P*QCB+jT{dYd82VkK-bul+9I!nIGGySagaM!KTL3
z<9R47NDK`NO$?3YLesum;LA>_upld4nqHAD{mYa;V~EH9Zy~@e{B8j|Mp^zj)+jE=
zmM;TR<p0bfqm~2ya~3HrTu~C9s<5be$g*{8wlku+P{ub~Q@V!A|JF4G&0|qoEEiv1
zQkq>Z&1X)b7{uB$6x}4tE?0~|T|2t6L5U3~UA^7?+&x|MvT}U7db)ddboZ=%xTls?
z{IUGMpg=em7yaE^zTcq#^6i$u<GpAxQqi@F<?pvEZ@0?r%KLTmVDwjgQ-|i`TSojp
z&!gO~e2`F{=b~eCedYgtErXbH{YUvjfgcL|P~e9GKNR?(zz+p}DDXpp9}4_X;D-V~
z6!@XQzplW8fL)w@!%y6SgPpkgnfjc)?Q<F|K2Dc5t)+6eS(JThG-+M0p_y&+sANM5
zJ@k&Ii|2xgr|wVpB%XBbh%*)cVoPR2+mZS*b6VhGLY;*|DhU#hshJDtiz7~(+hRzg
z6(UXAqSqS>>~lAV<0ey7l#9U2F~QAo#@KV$2=bmr@PB9s*U^Sp5@3iw-U%^kyAWfB
z2yrA)h@SRBOl>m2hbsnHc`_d5b9$mWJ`K~_XTte*Hts*p#p;@TxNj>!_SQn^4=KVb
z>tdM9EyjHn2^^Co;7yfa>nRCdf01ChT?x91ORzMr1e1rBVC<X{{IRhF8V?t8uZBc$
zN55HceQdtaH{A==n6;V2+ViM#Ts=LxRY#h2d9=DqGX1_zOs7@_(vAT>WW2_WLf_a^
ziIz2a9kiggUroqyk&s%?7m(8<EmEs&P3r@apt2+q)1F7bt~LyvT|%%>m*8R@gz3El
zv1U^j^zrVDz7PG8f7}l)XFB2A3t#Y^e9=qhgY6f*(XXpF#;)^1u8SA;|J)Jhy*fgB
zuP1IKr640C4a;iNaZW!IrO8>)&B;b_d=9SK=A!9gE+jMaaI#ZA5)S7hI;sFujuzli
z*Fv0LS%^x#B5W!u!sG)*FflHMFryeDtBMi!su)}QF6Dx>E4kE|j$Bi{4rk1JNYxV#
zk=nCWR4bTHZXHHaczP8r@yQ}9|L(NHSWIJ^0%?lGmo92|plRbBNq)RNX-{cGk1v@L
ze}EDBSLoB6Ejo0xXDt?I`~*jx8pQ0aM&AMb@gblJ*54{2J}kxkh6+5$Er+ak8O)1I
zp*Oh%(X0<uOeC0HSd8r_ijd@8giDJGVPaMY@uULWGAKa*$@zGp5s7@nVCIoH^ejlg
z(Fuu|P?v=1;oY$KeK(X&>5fo`9*AAn1IPT6G5A<AZuU(<{QDGCjqQo{?NZ^eIu&1i
zdSUeTUU<@}H!{}uhM#>O48OmXoA|Pto7Kmc>tChEjb8GUOtg;CfOTtVeez67vlvZt
zR}Y|_7P&OxTr%A#jv`6hU>e)ckGfy>AP;|M3hr)8nd<FG(7~LZH5yaNQUmh2u1mj8
zu0rbR3Org_g28(WvCtq7I-|4jtV22)n)<-tSt?G8QV=$)JNA_&A$wE;mS2j4MQkkk
zJ&DGKrBTQn5DCK!F-E6GAU{1E52ayPJSh~;zlbqpM-1kFi^rFYBs}iZ9o-C*adKt~
zda0)3RzfdW5ATf=GyCA!w7$q3o(A`VesJK@u}(h&bI)fWWK<^Xe6n!qau)7NvT>s+
z8wo>mP#fLIsU?+hE)zVtFEu)xkKRM#jXOjO9<8Ex%umuhM$)CsDmoO9MVEuRQ(sFl
zUE_iHCBF1ls{;i#IMU^#?MZ)98`^Whlpa+XQFfU=1#Q-$sZMopv#7-*lb@g>s=;p4
z0a)9n8uRV?<79^_^y@6cql8Lauan}~z6#tmuE2+?a<qC-2IszInD|>M7G;;h;&lmx
zLrQQ@RD$;FB<LR{LHhnsWIv36#=%H9e~m)vooLw1je!KQ$UGg3jKnz1J{^a95%Jh|
zJRW;e5}@`f0lUT}qQ)%=LC2Gjn%NB>KXk+S5#4b|*aJZ`dO%>2jJa<&ah4W^TnD~8
zx9hnkcQN%I?cTVb9xq%$?>wi{goVS&qESkR9;Z`D$0Rz$h10xlU1{F~Z}Ok-M!$4*
zAl-4+biBxtgnLEwT53o!wR#kHTALi7*J9+UpD<)o4Ia-Q0J}-mSU$QxmX4@`?I;;M
z##MrwFU8#B71*d#0eN~k$`6&n+qVpj>r3J8TZ&uzN-#931g4KAxHVLQZZ;AK4;ABm
z|4@vUN8m#1NVL|ALeHB~C?6AzMynXOEsH_CO)PSz$6}^V9QdQ+@YWz6W0u5Y6cX_M
zmjtNyPsC?~B>cKO36@>DA>%|h6!++kRyVpsHgqF5sed8Yirl$vdRknP>>e#Wx}Ti3
zte{aem5kR7r=SZ`a(|mngYA-Ni$yp+5p|{Lh2E5$=0*|?2eM7KrjHSp6un49{recw
zulahU)2K~rFUv7}j~uI)%Q0+%9Hj%~D9DtfN3<LwA2}Abk>i7g93ijkP<yiuo~P?z
z+gJymwRP}YP>0Wx>JU4k4hA*M-%9E*F|!W-DRnT5sY76LD7HI=!@Ela60Ss`M;9^t
z#)wgINDTfPG4`5A;s`~;rZN(TcSPc?aTG2Oj)LdsC`_)7#+*;lSp8EBRKLZ5Hz*d5
zRN`=ZU>r1?;?QaE2CivFKDS=qjYGI5m!5W)j9%{}9)ATjuAf4&zQf5%O-lQ?bQ&`~
zkuo-gQSY-|Xs5)R%)H&`$}jeG-^-e2n_ANRQW1@DGNhS4dgM4-o1T4;<Lz5H9=wv{
z$a6WSK9R%gp&aw?$)S0h>2!mwUzfxCiX1~P%AtQ=juU6(`1zC^D^JL==a?Lij>uti
zm}z-Xj@|p^FyAM~K&JECrJ<N^7Y>Wx!m%_m0-2K{kaIW!#m^!j(h<YSUW{@QL%)z|
zwL**q-^3VK5Q&Dnk(iyuw0RVTE|O^Md=-u2@)(SH5rbuUu^4_M7Tr#-=VCPTIFZnm
zOPk2&uEpJ<`gZ%Ms^4<z=rV;i9UVq>(G}#ovL9KRCQ?a67#UY|p{wTJ^zptc6-=-v
zy{A^R?z{y#xr!+Gfsp#W5m0NfHa%#PqwiNa7BkK>-ZBK=<m|p9M+b%%!=2HdVaQNr
zykd23G0rjeGd3}nGG;L98G{%y#+ZAdxXKHM@x*Xw-U!Ee?Fc;bh(KL*1eWHo-$NqM
ze@O&_&P5<ePmJ?yoZfI$47<*e*t97Ud7e?&u_g*0F46e7BpUT?V^BXi2J>3Sg4?i~
z6JN;UtiC#NMP9AB!%4Sj`iwnv)_EBPyHBQX?}pNryJd7SFOB%G;%RQDQ0h|MnfU%*
zw9&$qCj4qg<HS}J1Pi*k*o3a+3CW{QK##6#(O_o#b%qXWs~aPfk<7?ulraV{hA>7k
z8W@uqGZ=Fi3mA(T%NZ*fs~PJU8yK4yS{p*qV`CV`3&OEBCLC@R?4C3^9F1GUvGZa$
zY#)c?)R%CyH;X`JR0MRUM<D!j1lmi)`0bGx{qrNS;!z|9<VE4^y(pOWjfTO+Xsif}
zLG|!eocMMI*WSj7yE;jOTatK_e2+C!!qBA@uy_*fmJg-QW@YqnM_<~R98W4MLn!iP
z0M)GTNVBH6(7B#=H1K3QdU?#8Hu#y4R6|I!>;$xIj21m(dMLC=WaKjx8jN5}V9aDF
z+P##qim{flk+Frbow19thp~@wfN_X%gmIK{oH65%5LjLb1y2wLU4&tlBn(Lt!!UD2
z80Ksa!`Pi+NIw{cE5C$cly*3bdN8e+ULTy;J$)(D$We^ARbrfTi-hi)NLaW;p?Yx?
z9+*Yrf#J_w#^MZa(kn;KFk6G0WqgzB=QYxrjHUEy)Fg5%8cM>?rF3XfUpgKZPnYI~
z(9e$p=)$^=G;O*IE$?ebvUBZd@hNk13os!Ibs^bU3rKgU7WpvSvlw#5c!r|?moe5b
z6n(8|dm}^9?t=_Pn-%Rn$vDk8%Q(llz_`e`%(%k1%1{}~p_D>UR1*T<)gef@90I3%
zA<%pjg2zuoaQ$Tne*P4KM8i-lr%?1C6AHs8p}5&A3{6+UFf}zCf?Ld&vLm4OI08~O
z?!+<w5oAW<^}?0h@~Qo}VxA*+cA+|V<G~GT+P$0ZZd*e8RVPvXry=C)RZ1V8_M!E&
z<LEXogx016P=ixPvNCs}PQTmIgS2+kKi!<NE*R6NMFymGSeKGIXi+J%UEvc;80#53
z7>f2EW}ILs+O24_vb~D7UT6HmxWV|1ag%Y2ahq|6ahGwA(L*H&!&e5eb_Ah-AB^#B
zf)Q^Wj1G+BHo?$x2nMf1Fq%TxdRZ_oZwp4Xh_HVIVUQY!EyFnM(+Yv}_z*-Gh2ra!
zP}CTO!LcC>-FhA8ZfD52T{2(pyrn*;x&Ik?85}3O$!jTFG>gve8cjF$)zC2QJjz^@
zLI)h8sqSzvC4Bay`W#Ohv)7rPU9_c6b?s>K0&}W%F`+KHLK@atKmq4taI;TmG%fSQ
zuAa`&KVl1;Tq^{<w8RZNb4<S>!qnC#csSh%hpG*+VUiF#*!$hY&IYL6s*iP{`uKKR
z53*5uc<HZ)%dZ6ZxK@CI3ISSG_rQ+&6g)iN7iqjKjPc9EgPcNS?J36S5hW;iT8i<i
z<ydyP93!hMaP^N02-Bs=-YiA;&r-ztvS-D-N}QcsiJ?0y@%lz3E`O=SUQ-#G6Lffg
zlPK2mf&6>ch~oYFUrlmKb8UVb|2HN%HK%tJwR|$qYkB+sAEkGF|GuxxIchn9qG(E+
zUHMlA$t|bmU1DVo<!Anu1iv}oQ`~sGlz-TN)k7SiIM2T~r_Gla&+Ys4kA0_1s#2E!
zYTv#3PxYPC_w+^{@BgUpzUN^4lfIj=e%erl<{IXT?oj!&uUcoMOR`vQr82ienS80N
z(2{&vRHaP9*nQ@T%fG$zSTfz4t`d(wJI_{Ov?4Q38F<Rh`}_6(d1l-YH^sSmydZ<-
zQknkrZ_SKTM3|zdnuGOSOQaMuTfoxUls{TZ>(%Oa{<ia!i8jrhS5)<9;!I0otTM6Y
zpC`r^RZVE|Fc0%{zkKh@|EmdlimQK5BjPukuSn<d7k5ihT<Px<^n#kJ{09>Blx=Rg
zAn$oio5%mmqW@fip3)*^g5GNJuWQ?TjproG%eX%h!TDb?<$NZ7qFF1>Q(ni-^vBJ)
z^m;%&*@y;G{I^1yk=L6}Uyq~uqELF%(1kAF@useA-D!@VBTbmpo?Iuip`7=oG-I<7
zy*r>!TVGn^g1ZAmBW+P+XoV9yEwJl}2+OV-;m#-_jPB`U%W*x-Zli}4+XVQzPJp+)
z1PBWeV7-R`=i0Datrp|GE<*0;qT#eIJBKb}SL#BxBnmas9<a(tMUZwM6irBj-T8DJ
zy_SivP1$JIHy0i!^3dshJ`O%Ez=bV^n2=e7Rv(IRud)~x_lj{SM1s(n5|liUV7qe(
zda$Y3==mj3-TRJ`t6jO8V}6{Ll{>e`z>F)spu*i3T&KPJwz0?edGuyeJ-u=rL{5IC
zB#OwOY3j+8I4y=$WT7-WJdg&S=|n;IJ;_?)N+~H0v|?I&a<*wpANPge-lIhPZq@@E
z7InjbvUrRvk3`2D4hN_Up1t!yPjgSq+3AK$)h_US?u49kj)?r=fc+s3c(>dhrDpbc
z+F*y5d^_A3Y>TKS8^rgw!NF-lES~I-pm7l}rZ~KKl7KO6is)-YBDOq9gxjYi%wSWF
zmkN_%CGLrSCcQB8aBnP1?2BpJ(xCH+O|zM2pwu!Gk*%_D>S7ie$7f?^Kn_lCspnpH
z9m|a`&*o-!^WX|ad@g_JFQm172U+PXp~gLv$!qlp+S7LctwBD0_>fA*&ts|G)(~3y
zp);kYdsE?NH(GSSf&R#CPqGnhXs@mraZikB<Re$CxZ#DN$?iB4X@g$sA~;OZ#hg}}
z5SjB4F@ujJTTN_C*2Edo#AGc^96Z6tg{gcrl=E>bjgM8m`Iwi@M^F_X8RPiKUCxL2
z03U->JTYaW7zw$FNZ!#MpX-y6)gcuJUiQJfCF$rrJqz~xaxgb053fq|;S^VZFJB7q
zc6cG4yeq`U_#%v+U4(u&iXb#8#@dKtES44H-TY$MXYAkxO`X7b+)LnEXIXN^zkZ~u
z;pb?f!xqYZH=ha)H&BG*U_x;TMcAg%*v|<pR>H~rd{=U~;X}{wxs%ogM>=L{L!0?+
z$uQlFl1z*#-_n2{cQwK7y~dE7&;~bI6(#4s1m>>)5_lr~Tj0jCJg6^GfqAeh+U!$>
zfsY!}7pY;qPAeppw8Dkstzghj9kYwn;kHQ~HSg77@2!Dt#Tr<)NCTNC10k<U#LNvn
zL2vtDmS;K)CuTw?I~!3qa_~Vv57%3<Df->{_$j;qt=AXev|1r9g%u)HUWfzh3i0!e
zLL{^<LX3M6&LkG$qO=IzCKusJuZ>*A>soGKuK=!miUHRv>;)}IJVBGUt)sqeXVbu?
zW2j4dEqxi7Pa{9{q&I0XWd4B2Up;`P_UTBm%Uo#uYCFn|w4$jQ7S#8S2|e2(q*J%p
zC!Jlt2byHQ37jr@8|d%wKG0_JAAtwrngT<@d6+U=1;N{x7K_!`V|XiA>Z)Vqe07-H
zYoN<K4P@xH#?!LausPitA{##DlrcZq!ADaQAMd&_-%1m+F(4K41O0HbJ`=hwS#Y|L
zg|nBkA-bJ|o=0<Wwl)uknJ%RT`4H{Q$AGW-_~ciBrlJBIn#FwPPywDiDS(Ec5Nliu
zksMKo+Vn!CP&PMlfDflP+>kr`^*wo<I!{y1Z=u;|=98&o17+_SNOAQARQqcRO;L#=
zzayYaWxnJ(+MOnDcVJWV)-<)ik_OHcQHif1-E`HX<E3o+^<!(=?zt{-;;M6j%kSI`
zTsihdV2|hzfg8iW1U{+$7Wm-|567pgpx##%3%08w)m;tiW~rfGtrd=Bw}R3BRv2Td
z4tc&h;x?-z`HMOp25KO#Rs(0ZYM{Zq8%8B$;fP8B&WQ_f{UIAiYl>jlN&=C03A&h+
zB4u|en_ex$r-fyh_Nff3Ld(&5SUJXRFUOpF<+$0p0;YBqc-*xDdy^`#JGTPc2UNhg
zAc8Zy8pkDQ`*Z87MckL1*JQr)Ah}&#Ox}hKlrl$7rPK3icWh76YltR$HKK4;f4bo5
zNli)4G|k_Z>KohA)FWo}h>do!A_F=aq)WTjYm)9S&bZd#gU=m&5a8AUq4D-e{%DB_
zcE)IXQx8i%>fqXZZS2z0#=7lV@UPOsCr>SW`CStcvo$d{Q4_CpHKB8skA-aP8DGGM
zJ@FxF$Hy7<-WZ&p2S@K5_GvE@W7y{{S(gH=s4oVuRVi{VlwpWl1*9HQ40<Ppb3-LO
zHD$O_D8snjGN`CkAq`a+E3LweMOCmlRfYMlsvtM$k3pv<agX<p;BLgma3|Lpa;0@Q
z>6_6u`n~-^vTHwy4q!OV+geE-FK1BElq9P438z~IUFq&3Z|YU!MpZ!$v}>a^*{rmr
zW35d|bEhF~*rrGQ{?I1z{4n_IhQpz>D|T&k$4?WjSbvMKskJd8WJWl5#SpyNLcF##
zK>AialuQvItx*>|2VHpW)<I&44i<LMf#{t!cJ9!|iaKpH1!?1>sy6&=`yl&B77lFd
zgCj+Kkv=ICTj%AWZ+H=`w@9!}Qi|2P$}nbkIbIB{Kw*0+v=&IQgI|fTv`U;>UJ2KS
zl_+W_gIT-`k{TKE7Rqq^kPMH1-NjkIpTHgMAIcS(2)J*@AJDZn$LM<G20Gb&E>-uc
zr_Cqiv_8F%E>Gx1+DBt4y-f%~yEEPP@uI?>u4K#ElhY|HdUD-@M)wdAuLFyr9(r`G
zN{rkO!AMN?#47dn@cU?v4Jl^09%YKeDigd<GeRZPBsa+b-Lv#DeuW--=;|SHxB!j@
z0u23`J(ncw;x|2A^tz~nVGDKO)?WvyNjlicPeZjvCV~&A;i+*Z%&YUj?_LDwM)oO~
z>CtXq8EQ|KqjG8moZO{|StLa<uM+&|N>q)m#Pdcr9q^<QAB8eZ_mN>nH%6%pTBBvi
zuszM0EFZ(A@H%l`C$+hAjZbOuyOSjBxq&nVbEwgGERC?Pqfejn=}=HA?QDpl03k<f
z0|MyUyp9z7#)Y~x+0lh*R#d#gg5s=1B!3~KvpRZIdLj^QK6F6UZcDs8XNa1wLP!P+
zG5D!ItdsPx!Bv0_z1efvUL9<V)Iq~rZT9S-4gCgf7?f(`r|#O=-B}yO*4mh=p^djs
zw2*UN3yXGWp>l6J+G%E>P%8@+>IL{}B|&BZ`z&_9j7?FMW9Q`xY?Vo2{8WmrA(c2f
zlKI3QrUBC-z(599PZ<njWZ0K2!)fLtn;K-ex<H0e8)O(U=>oScV;EPp(u0d~(dLp|
zpHkz(<23NvT5>R*MJdZilVC^<eLj>+Kj)=T3NM;M)&-O53qQKo+mrThawh&ETcUz?
zv}&|DH3>~9>#YGf*a*n#d1oBbcE)p6Q+Tgs&#h7+#s?U{vQ7`~P6C()=)!#z>wirh
zJY1&@w{mT)3D(9rkv4Q*X<^z)E!<kEh3c_dm{Fz$om4IK3DpA6OABiTXCWmf6JMfp
z;k%&-br(zUd_x&(T2&yqsRF~durVaK64fk5&cA1I;w{5e*8Z4rGSsh=;pr(EhWsvr
z=&KAq!YYK?R$-cV74o?%9QA+9t=5!r0mIxmCvQV8$^I=lr=6lduC1rytFvj4)mX}8
zexmJEK=X&C()F`3^r<UH_o4%6%bbqXaNmV&p4!o;kyg|=&w}Pyim2x+Avqc9k<sUL
z?4KEf(f5O}e{~le3-U$%!4CL!qchG7v4{Vg_Rwr?1&2*-Ah=|K%VKj3<e7ncCc;Se
ze3t+dh$b4N*9Rl4Dlh^b8KR=h5W%X3Sb4lB)~d%MYJU%$TH6Q9{W4JdDI1H_^YKPf
zh)!9>`1(PDpS()Z&$$d5FU!z-LOJf6RN%Ma70`W90fC<s)8$gEYn0;Vdp6zXREc{%
zDzRnCHLiR}5jSv5XD-6fm>aX=H6?5~MR8&lJJEB<Wx-f7)0ETprv<dhtrvaEj-|*W
z91S`hK)Xa<G`+hk8N}Pufm>Eo`HaO*rijjZ8B(VtJ@Psci%lndaNXDw{{5UG-fe@h
zsCL+S+5%>FW)SO%Aj>tz=e9<Ov|zCnF2t!B23W6cfMrwk(bq#CvoGsm?l3)c?V^Vp
ze+Y14w*Xy72r%zScf|Uopw1-&ZPIeF*sl;#ttHS?FU4$m8RpL|$I+?^7@A4(aDo&a
zUP$5KUx~NUO4u*0M6WZIsC`|DXZkWMb&;VTC&RoH8BVhMl6U+W?o!`E?(z7}oZ$r_
z_ch}sjlXw-rr%vheLQDV^olVwD!i6PCgoH7p`O&oD~6i(5<R`^PXS#zl4Ph09UNsx
zkL<0eXOIQ$+Gj#{W(mplfPf;MqHyJc7qaR);CQJ621Hq7_XkVxSDHc9Q-r+%#!y>r
zh=f`p9Htpy0edIR>ZXs2*Yz;JOb<tl^&mYYzytPf7!xFbkT1Zui@NBvR2RanNwD9Q
z3V|>S=gjgk<XjOtr<WkuqYS(u<#-^k!1*px%sV2*A$NA4sjWn-U6l|%XXA-ThEx9R
zUXv_?f0+!gM#x|@TZT?+WEj*a!_KGMxPdPUxlfvYZ0=gX@l&7B)VoKh_m$N&aQO^6
zHgpv2T~$r>0XbCBDw%B7MADSPAlkRH6K#eE&0gz7{3|weczj!`*>6T#J&mbtm;p_x
z)1{yfJut{H0$H<zu<@=hUOwr7i$2cKU1o>+o7OPd*cOj1EOA8N9O>nz2(e-B%=#t>
zax=!8Vk7K6W{CIxhS1z0#KkZntbQ|qxXu7N_69iT8-m%l<Ke%)CpyhZ!=dmjL>$k>
ztYHPvSzQFjJP9@*EkWJ2QfRC$!{Tn`xb#aoA|flWb4djjzpg-EZz=i}OHnb8`O*m~
zw!LAW@I@?khVSN<35IZY-gV(jghpIjhc~pt<1`su+CUY%=FpcHV@Xmgr|FS}WIwzY
zMebwoe>Nes%eFJE$nv5gGhAs{y*-7PTT`N~C9Rz%qNr>`>OV-226>0Uxh5F%ll>9s
z$I^>@+%e6Bjkz=Iv1y(SmgrjJr`v5YrfnOjZL~n`NOR;in89$jDIW7p@w`xk%ePE$
zH`4^=uZ=Ndv@z~D8>8cSBaF%ogTak>Z26RoyB+&rnNvEvk%d{Cb1-i!drnC&fZpyx
ys51XK_H!{R(<C_aT7vT85=>;{ae-?oLi?9u>W)%4zAHtdQyENC$`CrZ4F3fOfkv?a

literal 0
HcmV?d00001

diff --git a/src/LIB/RAD/ecrad-1.4.0_mnh/data/YDRO_ice_scattering_rrtm.nc b/src/LIB/RAD/ecrad-1.4.0_mnh/data/YDRO_ice_scattering_rrtm.nc
new file mode 100644
index 0000000000000000000000000000000000000000..f6412d1704477f6d1fe5cfda9d8068ec59c7f394
GIT binary patch
literal 16067
zcmeI22UHZx*6)WPNphATh!RvpLCi9LdpCK841%B{IWvTr!N{l>j#(6QIOd2)#hftb
zqaMsT=ZFCVh<-H!=iKk!`_{YPTkqcW-dgsW-_&&XuI^pC_O70;+7tY}ef4y#bX2vp
z1Wnub#_h1M;Yg>(6MiM@7l-(JhpXytRdx7rOiiQlR8y6|FLhJ}DoPPUL*@7VX3B4Y
zU{*sR5U4l((6CniQc?dQ%zkMlNm5c%;|q&At2T&ev}+J(+!a{F=O%f^=Sp(plO&RY
zWcJ-4U!>L`FVJf&|0=xP)1rT<uh(5g#Z-A*RYgTbzfqv6R5Mj-Y}CrGQ8{~og_giV
z!K9k~Xp5RoLt82vb)V5#-^}Jm1vPcS-`k*RU$?zdR862E&{B%>VlE;=V8RY5ZKG_L
zcR-L&a6~{zaC`*MiG2jRvCJf_LbHYn!Ca8f5BWZg@{#Eoh4$=|V4s^J<-L9FlQVKs
zatkvgx%P=NseP_Q`qwFWLP17CVs?s~edFftDNrh;%C4_3&{ZlFADkjhN>26k&dA9~
z78o*rXlSDS59QcB+2)ixV#`0;(5K;~psD;<Z5Z0<5C6Z~uxpS1nGGAf`{L44s|FkP
zQc-E6G@UAQ#3ui1vf<zFH`>(xpZKMY(qq-yHS|{Xa4loyHPuuEmda8furz2Wn><wY
zf9;_^ZmuRUWe&jtq{$(8WrStqCP-5X|GQhT@M((2rc!w;WixgAH~6jr^KN#jzg$}!
z5D_V8!6cMkZ^2NGJIX+8x<XU=!&r@jlJeN!?kj`5`LwQULU?~sQ*)cbtUEiW+;p>*
z!KNlK5Oh*%5bopa6Ydkt`^5jypewtiQiJ3KX+mL&^e<KZj3I&6e+vPo;SUX1ou=~7
zu||12wtXEMqxfeUnKoAV&uP@8;aZ7btWu+fA*+o`gP#5kMN@pUy(ZUC`QN&Re#2N4
zAIygqNb*t&q#4X9l!KT9L)lI8lmg`l<mKVY1|>F}czL>YcXM}5PfqRP<?iP0;pRT#
z>HZN+@aOjbf(88e5b+Oh`SF1M%eNZ?Pw=KuNoCh6x4%E$^tf$0-t@d~7>xdsw{&i(
z-#Fs`d7Y-?O@Absu5;0)LB8pGP~#xhwEwgGX@Q>>_-TQk7Wip_pBDIOfu9!mX@Q>>
z_-TQk7Wip_pBDJnEpYH*CEvZ*Xx_Pl7e9QeDeofuKzCfvQN`gc<W{km^aN9Bck)PT
zXPrs*J)%ik;zxOwJ?Q}FM#omPrX*V{np9>)e%0Ev%S4U#9R0$nPJGS134F*U1-ruk
zz6W+3>43o{ZDDEUg#A&iF{61a%(k~f(nwqMw6=l$BP*m_Yk@-_EYYW<B|Jx3!0Wy_
zG$PC~waN_N`kNv6gDFC0o8lB_in%o=Sbsu{R}Z7H^w~go98ZM#gcNL^l?FIwVDrFC
z<n+tJHpgr{Je!U8u{p3gn}fOzx%fkoi)Y7jQLiq6LoW#eGbKo#DM9oO30_{7;QTuY
zCZz1=cb}E;XZUV>`Z`1Y?zGo*c*iMHYrlyu>|8((6_aT8*x}@rl}az`Bgk*1H&u-5
zMm5##$!c;d^6g?tL6;54w67NZI#-2;H~Yjr|NMgM_uwA)G$0y=iJ=&u7L3!8B1~B3
zk9FOBF+|TBi%kftsSh@m3NijmZ_FFs3s-G>V%6pDC|lGGXLGtj%dZQjwev*37M(HL
z!~>Ho-Lb~G6IA+i#K@CbcrY~{^EwYeU_m_A7bZZ@FbVVGlOY_Ig3$%3a0p4mBa3t-
zoK1&nVFu1<XX3XxnHXS~h5V&isIbk(r<vJERnI{`Ne+B2<zP|QTy*(zfL{<=%*T0k
z=MT6U@!M_QP?_gx8b4<fUA9<A6;+c+-*^P|ew0dP>5-Hu@S)BZx{<h(D`_5VMeUO;
z2?9e3FV><_2UO@>)F;k8@Fh30%YAO+>JfNXG#p&JVF-OW1m8CfhWFqi^!Ap+(o}{E
zPYThtvH<sH<>Ojm9=7$DpiZ2NMZz4k@0gAAO!v#inYgQ+fi=z2(M5y#LGx6ko2H;D
zG#v685jLF;!so*w*xxe@SCYdqs(%FfxJM%NRU~@Mje<tEzHmF$7kOdPcyTTon*w6+
z@IVX>b?JxQtNS6sJQk^=VqyL)7VCrJFnD7ebhaMm_x&-LKe?bQ@AS%mzmxEa60V-4
zRU<ai+KvlI$8Zv@Ix&n2j;Bz|><F4+?@cy(-AFpRJ?Xb?MdNN+kj4Z9YHqAW!~3gH
z$^4I;+q@TCz?6HO-QJ<tVJAm&Wj-|4<>1@p47m47#j-yVA)XnJB1M0cPVI+FyZb^#
zFA^uGgu&i31cTm+F{3gNgBJSZ!W3WBm3gCM1!3{SKKR@YxLV?ct>u33J}t(o3!xA%
zjlhaKQRsLo8gWDW;r6#!>`GwQUN-=*ZV$wk`glB4Pr&QX2?)BIh}6m?jGUQ_ivv^e
z)GihGE~g@DOd3?$reo#Nbd-PE#UGE9@hyJy<U4QF=UYvDLBq06(D%6CD2C=!AODGT
zLOPV@tW2hbh2gZ*pQq<RT}g4N9SJ7dlU|?&c?b+>#!OA>_*6jCj@EFiFF)r--MGu;
zbty)Ha0F_F!?A=LhJ3%Fm=ih#Y6AviRdx~jDdf;9lVRIhDS93$g!hvIXc`t^MUQ+~
zCFa3zq6F1ja*=*F2i}G`DC4tXH#7@Nwq!za#~b@<{BfJbr0#VwmVOMvx2wS{HbM|8
z48>?>EA$UT&zoVm9ubb~hr&_PJ_1LlM?m941iD2<;_8-2EH;S3-1I0I?~j7aq%Yi4
z`(oR$UHtP;`TU2O9(=vK9$y;$jK1zYM%K+X&~TS|#K%pbO}7WrN$(_D;}u4mx)6EA
zbfLU2ZE5WhJ2r-x)8Tl1;#4##H9<f}PkrE$UOeNn>+f)CZ%5$L-QhT0H4H6|4MpPq
zA^5a!F!mfQ!Y`-f=x|krvS(84sxL&MWg%?53SgL#kDNJq=zdHB?qe=Gw9Q4kK{?pC
zARAk5WZ@l)r<UoNh>`~Y%K}g?7a`@47^ZuI@avdhoaqz-rxPLg8WM`TzlWj=h2hP{
zFnF7Vqoyz%8_tH~jZ*|vM?~Pj9}%!|hy<4(iD8wI(9ns3YGf3gyYJzpae4f>Sa;qb
zT9?0I^psZYJ4(;HmD3`{T$=y;FSP1u5iQ9{q)m!Y>OO@d)5D&0OxBh{LhVS`-JHg%
z>Qi*N2K6@)&@VkdaC=jpaf(rQxPey{&_AfaoOKFhmMPGxSb?=^3aEuC5Y<zGWvvxB
zrK`Y>nqpkJU5wHb#aLWXjL{2=p*^k`fs$g#Vv4bjD@I+rVx*ZC<5~R(41F>Jzm@vo
z)>?mPsRm-m)j)hs5@F|YHolmOF+5z1?_<U2T_MJmt74?qis4}q1dC2VNcIgvT3ir3
z^MWvVY!Ffw1mX1NAaEyw5OF&QOWp_J-R7-)sUn*{V%?EXI-<iLKl+GzemP7RPp_kS
z;d99QkMVTCOHNKI3G{>yA$NHnsyp79K2CR`oZ^;bm}f?Hp?XxUr9tDreB*rDyyxyE
zJ>`Bcz0Gy4QK0Ib0`01qO<pNb^jv}PCklLipg`GO1#nA&PuCRKbXkGI3krz-P@uyZ
z1vF17aPF7_;}0v)^?(9*_bM=9mjdqF6}VBMz_^VH%>U?z&b<Th$}<p~&jezJhX|{O
zh;aKi5oTN!Vci!IHZuFDxQe04i7`7;jB#mV#10kX<76?;FBW6KMlrfoGT68@@Uj@^
z?ufDXi5Q<gRPe&d>3q9wu6&=~+Wev659ov2Au5YoOU6#K$?Wntl9Wlww=|wS1;NDJ
zi$sS#sA`2XRjjrpm2y)WxJZ{;BsZhur|LOF{de5&Nl!S*%3GZNR|VcO3hEU&$<X}9
z+Q^7w6fwp#CNfGGQbr;ph|!(VnxW75#OmE)9A&I!OlC+Jw*LM|I2?c)c_3cg4a9c`
z5gvt$;94X?)Jzc!){78!NCcA`BHVb-#vC0ndfA8((oT#!o?`6pEyf-aBim1m&LS}a
zf|*T1#kk3Q>}uF1zBVJ37j$UP7h7uaQMc~Voeu|SZ{iwCkDWzp+mE9U`Gq8zHIRZF
zf~ah&kbDe1Xo8wEjj*t#%k51mxV<hdsaL0ZgL=+8;Vt+6;A5^(_ZH{ObU(*%Wo^r3
z%wg<eTx2|F)H8IyD`3TNVsv8kU~mk7MhGK{5yu$BNMU3$av5#<`=jVe0H(zUqQ!<l
zOurq7a&-|zP9hZeitspD#GY#+3@sC3`En61Y!D%Gs|doKBE;_&A#slg*BL$cig1dt
zW*?KO6rsxj5f)F`z`Lp?^Apdv<vrVI^7jYbrIM5TDOtLj7SAoC>66CNhtvY9+B$%)
zCW<NanHL$2bf@VhPV{7%4Q;42q27m^)027XbRq03w|7G|2dl?iyZoEn1!j-Ftlf(k
zXBcl7+N?dc3^zt^hO+J93}w5M7|J&1F=ULP3<YBh;}^yx#x%wZ#!^#%*o+In+$RBW
zY!!$)zd#(w3&fFmf#|U_5C=~MV)UgzXxt3M0Om*U9t7f-CxOs^5r}5517T7fh%Rpf
zv4}B&(Sh+BJLbsn&nV}`YZ7_!8yEh!Bn^ISuRD~hzMoD%UqwNk%1G05EIFRar$pxg
z^yHC<ei`mX#jV`Q-`$BO#M#h;ktP%{r8ynSP^XA5bzGa+YVP65M_gv>n_L95!Olio
zn6h>#?LZ8r4dNNf_7^gSF-9|#?ViF=ws|&VK4TGMDPsj=HDfKKoZ-^b51$wL<HS>c
z46q77yO03n76w3bYykA81z_H+0DPDifVPVQQ1NR3HZBc->+%3}UJ-z646l^|7|XcH
z7`Q3`7OMkL%XqOS0D~^9<G<J^@MF|nc&%s6_!kRq)0O;6`k1haq+4gw&XZ&4p+!D*
zoZg>SCyFRS%ZqAucB1uFj<iM9hPHJup&e*W;STDwy{L{eeEfzxn)rw_IB|pfjoITH
zqaAC5FGJ}IS&YF9r41%BN*Qw*%C;|KtYRqJy@9cbv4ydXv6Hczv6oTF$g1>4t&uN!
z4Dv<irM^%*?u+CqUwl6Ai|k9jICj++<u`nhd&?Ie8FI!8wjXlM7y6fdp?87RW%4$R
zJ*-YSs}sh~nKL=NoO}G?py_<44bi-GrwhMrnKmz&`Gi~wj??v5zfso}3uuJqWQy@9
zrq~nd^inUHR@V8^oywkc@|7D&`nMsY*H%O&#`Nur4jp$_rxKg5+^#p(T$o)LvREkH
zJ=g=AFM42|r~|qlcfmnpM|4(qz?C6(SRHMPE#s}354C{vc1zg*YJu_P=BT=4hQsD&
zc%5yEwwF!NHPnPXSB>#K#~2UwjM0hxN8c=o!S0Fy*!(>av)t2=J|GhhPi7-&4SUzK
z$-{K-e9ZMMfaB)^ES_D6*_Ki~A1lRy2U65}%8)Wlh9Ns;nDR=78`g42**gcN%h8aa
zBlw#{vBm}RA6X;H_v?Q($tm8o;cfiinB>%u-jUt7WS-vm_WwUh@A|R6-;{IIcmZWm
zAHPucR|d)5whgBOnk1T*`5P1b7X3(Z69{7d;r!K30shMC{CjiSx(Y)1gSG$IcTGuE
zP36DZcW?hweb@F!dZR$_f7EwBaxnf$-_6@RXN=P38dl2gQ2Dd3G!qkYlUZ(MQ*K97
z@@11iWAbHoQBxAezH?XF|Jyr{CDXn4%I)-L=e1EPt;~#T3cRMy`}_U>d1hQm2j#T|
zf<8tKWmEdozcn*X8DYwvY6#Yv#z@I*(14}0HT`HTTWx4||8G05Dbc2(^RkQnOq^*<
zjBVCf@1G~eW*3z<dRQl`svAH0^8adrp7QQL(ulMg)KI4L^azNIQQqnA6ZHBti2Mf<
z^qSh-ctgSKVYbiznMVJ)1idDWniBLj1l-){a%Cp(UpkU69~{EpFLCCR7HjioT0f>)
zmPaVP^E$E<&L+X6anxdeAw98+r+wXmY4i(FMRI3K?B+smy0oP0ycwx`>rthJ2Cd@1
zb9$xkIs2DSxq4T7v^(SsmG6$wc6Y!m<CYloyanPc&Eavy7`w?3rLp>GEzm=<t}b%o
zo8#U+9sF9NgD*XFp#4o7^G|CdcZoJ857Wl4{k2iVYooHGHcDFffqWzJU`9WLtnH7K
z6$8=Cb`TW(6HqZO5l8kW;f`(!y62{%{6rd@H8XIlbtaay%EElLY?PkPhR(1Yn3?9H
z*E}{2p(jB}fdtb|OR(w8ZK|B*%&YYF;<wN2!M9F!=JnSY@V8G?Q|I856fj~P^^2KF
zDjK8dczp&f)agreHGOH9ygS8CbEW+?t!VpDOVX@2p!S2bh@YrJ=EpyBwd-GS(c1x&
z=rBwf)E8Zn`+~wF(BfSXKCr(=j=u?EKC~;6^gH6(Z5Nm|YXkdOdxTH71;tupwxR`&
z-Lk--M01R3W(LVo6PPbI#`0f{AlPGwabucc?43?<X~n~Sq6m2pgP|@7!E1RiYE^>p
z$0reHeigxng3xh%2+l->!#E@gD!pRx>T@hS<O48zeLSq^C&1l52{-p8qo-Ob7Tcxa
zMCwF7yU$F1$*Q6J77_6o9;W=VZ_jAl`a|S*YYj<X%_7Yu<0y7=A-(85kj$#Y#QO^=
zVz)cZUhYJSqc#*?YeIwVb;;PG8I4iYbBhG;xX)vsa6@ji#k47&*xuY7gG5g7Q&^&$
z*Tv<z>R70yf_*CAg)VLCg)3&&3FoHP3RfKXB0RhFi!k@u7okgDt<cWAPT1vmov?23
zS7DV?z3}4MdZAa|H{k)D@4|DxIN|m@9yV;6<D_jEy7mgk=731FGl_wVeSchvj0d+R
z0S;4>aMe5oU*c17M3#noiRl>QmVrUHGLTW2i5VX<u{}8peU4|LzGXJ_a<j4K{eFIv
z$s)cvd-a>+=)zButMltsZqd~ll{6)9H4Si>MNXmP$hftX#ypCrj*cNT^I;zv73WE5
zPHkzajU7E}ZBD$MKIy&Dpr=v+jd!i#w1+<D0#=%0gQf)nj_TsnuPX4DeiA+nt`T;8
z^ie3f{8<<=q*nN_x=t7uT`&B$=9}>R&F?}QDnL@50Iv#FQ1DO%wPIE5+Mo(IO*J%&
zRl~bgYFPAI4KLevg{fT-7W9e2NO>&2osY*9&18g?rC`U@6gcFhAaQ*PoV%rB`n^;H
zPe{Y^R_W-qHXWB*Wgu%x1}?qIzya?}%$}KvtQ(myu+4%aTlV1R?&I~}kLOPg^5Ziu
zTk+PveWjOKmuU0)oh0A5jQT{BQtNIdR8x~nw+rKFnovZyGQH?zFLw%yaH27jZ0P7|
z6B_clIrV(6PQ$(Gxr#$?IkD&ocl`7N;Y7Q)!Xm%7LhV-Xgb~eagfFIl7H&9QC+t7|
zn@}%PfKI6@kdZ1rsH&m!IyGGIQb*-sb%diC_HS;63~LP>9Ik=9%Npp@Rui={O+@a~
zM1clj=lQ;9(<UAz*OM@Raw^KcreevP6b#c#!L-N8(D6w@-pv#?7Nw%AZW<m>OhfOw
zG~61Lj-C6`QEHrlL-83{zcK??9%kUPZ6+>;G2493<xL)V@%6FR{C+Q0-nH~Lg*hIi
z6}Q%sqvafu>=;ivH)T{8J%}hVn6#vzgS$Lv)fi{88E#8shM7`imM*<>Z$_Qc>bact
zZ@F%&PdF8yTioFDYlX{H&j`;2-xa>Be<AE3d@sE4`zPW2%39%tjbDYaTHl0*rr(9T
zRRr*q2_Sncz~;UxSaM7SK3!BXZiy-eYp6jxRSlmi)iA+G9X<oqv3I39TxJI2s(vEs
z?q*>8tsI22rK<cn5?BUEu+m9_khT&`36WrC8C%ABA;D|EJd~`+gXQ-;2>RwDe|0`Q
z-sR(U7d91MSb#3;3c%egz)kH!P<KDRz&o1XzS@sJ`QC-+Gc@@LZda-EkqR0ydp_L}
zPo%nT!)dyHI_*B$mnK&G(cYImX^pTWg|BNv3$3k5bBr+sp4B0(yK3~HcO5tF{u}P;
z_($B8U>Eou?SgB=y5ZnT4-|22F}PPtJkl~lOrineD|In$dvkp4*&MAe=s+`B2W4-y
zF?h5#j#_DB@m?)7&(MOmy%uiW(?oT-CPGJQ;z+C}`qe}O!D-mDIU9L-IcV*d4HI!T
z+?M7bJy3#*eR=5jA|Jch^uOVOLKF>=;;W?$YnRJ#*+!0%adHISm*cp55oXAX(7B=r
z7LSWC%VaRZ)64iR4@&u@Uo!cu%JzJ=m4IIne}T@o+e=eht|xL|K%RXklkJpZn!P)n
zqE7Xtpk;nEI<_Z;jdr6MuUgXzg%$O-G$w;G9hzIEM!VALI91<jF7Lb#UUv+DNezds
z)txa_><Gm<TeMwngUCu7yuM)#cQY$=GPi_dxH(K8n<9I^39616W8NDhT=z7><(Y;U
zqiqPu3<LPKH-JZ#KAw-#$AwS*uxxe;1n1JQQJ#Y4A<5XaGZjwenef?{je{$4QEi!r
z6{39T2Nxi;O(Ci;7GhSq6zv{M(K1{H|BW*2`7VR2j~u3><T$!jjv@EtC{I7aKlybj
zuj3QXha|S*6Sp<vzo@*Ryv=9Ha{M;Bcy0-eb)HVH6GoFmc`l{9_b1(7#Pt0SAsNs2
zAkPYCHXUtC!-CByY>FOPm1&U3M?iZFYq)8Kfq2^02gP?hk-OLtkv6T+OVA1zf49S#
z3|n-lvBI@#OO(V|V2*(~Qq;|G!o?JJc_t{VGDbf?V-%b<LXTu44An4#*A_#ZOgF?P
zM?);1+#lzsC*kXM7C*&_c=aL$D>h_c0DEraWanb<dkKtP^D(6@i;a(LieO$L_Si`A
zY=RWKnN8+(m0>xX9+2*kLH<I9_m*;)aB^HrmgCNin|z?=LVo{vF~7~!jxSZy;*Y$1
zN)e(H^wDV}xvpA32_GiW-7h0po6~8=mA;g9+K&v&dQ#kRH!A7ThP2LF(bg1WT6J58
z3M|!W?5jFX|421wF_taW9qx?1c`cD9w7_ClOH_+2;Mv9ux57>EO~)9E%#3h5$`D=$
z3{V?r0E3tMxUoPVDt+}4ZK;p-_w=B<Sr4<v>LG8S9&Yv2gO`;aik8P?*xMvF<(Pz+
znpCXP$b$3lY-xN0TME+2$1g?&D85|)*~mf!2&5=Zk-}@I6!zbxc;Y2PH@OUl*2tiC
zS%#6{WH{v@M=meNvq5sS-u{4p=`o8hvG!$4Tz350ml}M^(1)ZaIYQBH<s|Mlk2D@n
zpyNeDY3%3}%65*Vq9h;smfoElt6V8N-+{E<TF|hkhO}CtO*by6(%BhbxL;<x=EgvX
zEzTXF2(`lY7ZxxGvBaeua|DH&f;(c2Z|wh9Nu?n!=ow<fC<9~~8sJogJ{BeF!`nt5
zRk!rea+Mx<sU8aa^zf#o9!zU>5p`7;9of{8XO|R|PESF{eHnP}m4jpBB&f2_$Gp%2
z4Cz{kEVcw4P{_s+Hf^*UGFVNPVf7gqT4~Dhx`!P1Q{~t)RSvTXIV!5;uzfAZ1<fMp
z?^N}kVm6hZr5(+8`_P#``^k{czfw(R!_SiYq^)#2U<sZ7I*m-nj;5zaa%o>&e_FX)
zOy|`>FHU(-%57)5Z{L#S`DRpBu15izH0aYn0jcz<;ZDv@Mk_%$3=7%xZ=ElenD&9`
zA3ZSTuqW>2cf$F-uJD=A7MgpVa8ld`+q*izIl~^Ut6JiDqAj9Z*x>qGE4+Hv0;9fJ
zLdg8=p~M1XkC`LzYk$0b9|6lnQJ9m`ALC^SxIQ2S!%Na}S|<zF+4AZ3tt^c~KMzY-
z3?-WvVCR(rj2~8rm+awrs6>ja*HYXMltFia3|sEXAhDHWNjO{5zEs0UE||n?87J|X
z-#YLy&2;&iEzhZH%Sn2?e-l+XFQTfDDP-q4lJ;!Rq-lp@$Zcc*<y!S7g;OVbT;xbj
zZ>{M@z6ohxX->J{)oH-OdhSf$cihoM;n246LJUit;750cnN9}?LtJooQXAYaZ3S|)
z!~Hh4s2^;N<tA1*_}miPsx44qX@PUG=J-}&hK()FAX{LHYfh$^RB3|AgG|6RH$e;6
zei$$6i~9I@6s|~t@$(EExslD%4J0`5YaX;O<iqkn0lG6Am>5bS9VbQgBPqVHX~uI!
zGSqI8!STKf^L6ErxXUp;j6H*L<fvrNqw4<8_!j09`Q_ztd~%R8KVpd{|MKz!IudY%
z-UpV`#`Jmgu;oNr@OUT{KTaWTWF&cR@u7*EyVC`e4rIO2fr8Uo(0*GZ>akgy)K%5!
z{?jkq`}?oC0ky&SwRJD7e%TElY&t?yy)8y2w1KF~9`n^(;_XpuO#RRT=5kA<_Of99
zY>u>KGo=4!iu(qp$Q)w=e{B;)FEqx|UdD)dWP~SkjPNGR2wT|ywF7a{@SB=|-iy-E
zt27H%p}C0k&x66Xe4IH_fT?ALxYSjOVcVrRWhKL{Oc{2R%Rmp=SYjf_vfgsMi<e_@
zu^jpH<gnN*hr=;BT$8WyH{Q(UcNit`ZV~PIKQ8O=abKQNZ|MoL*}9QELl#oy)XDU<
zg@R_!$e{ZNqseHoKP?aMMLVP&>Bf~dWF2o!pPn1jpnlEi+yr$h-~W|cIQ}hn(IyVx
z%p>5>mMgnH=P+KE`AcLMOqtRN?xS5X-?%NF*gL^vR%>XawnFd8cBuVm3)3<in8aA4
z#LEhqTnktvSYqZp3tW6+j<5daQ0z9tjWGdmUmA=nyJDc3ACD`olVN!<4MruISi3qK
zYH_(ZwM&9SSMyk`=EHe(0g_r4V$9M)Y%`Lgpn&<!Q7PW($zUm#;rb{UMls)UxF*AM
zox}VihY5V-GM<mBwd5BaswIWjMXIgZK?a)3X!fb;WOj5kjp&j~AODD@JzE1Q>0@s?
zd9@QQsCOh6_U5^JtO?!N-kfZ=sneF@U%8x!x7?BYk2#~?eUT{ih2J|bB<$*jvC};;
z5^i`qqaAW`oZ-2>4QgW?FwVyw8{{pq`>72C)2v|;X9YX&7LfZ|V&@<W>?}3M^($s@
zb2r0M7DpY;OmXR1AG}!*j0n>hTr3=b{<#V0RFVwSH>vP=kPfHenQ%Cl1&x>4czHbs
rLzm{l(O-gIS0qr6%mWVQVLcm<uNCGaVQ)T8)a4_$a{)#q7vR4DnMJ1=

literal 0
HcmV?d00001

diff --git a/src/LIB/RAD/ecrad-1.4.0_mnh/data/YHBU_ice_scattering_rrtm.nc b/src/LIB/RAD/ecrad-1.4.0_mnh/data/YHBU_ice_scattering_rrtm.nc
new file mode 100644
index 0000000000000000000000000000000000000000..127f166ddd6073d17afeb31ec962628965571cd1
GIT binary patch
literal 16067
zcmeHt2UJwa)^3x7<eWtW6cGd^8Q_C`L<LDo&KZf_O#|Hx3JvB_%wYtxV+KrPjw9xr
z!x#~am@tky0{>}&nfuSZZ@v4!_2#bk*7B^bu5(V+sk7^=y?52A{rxav=rDCv6IEGd
zWlq=8-gzAH>bT?7SulJ;&9dMTLr2T1ZI`wBc27>Rv#9Hlzvim492sd9EiLJPqk{CG
z!&%UgaX9i_CyJ)hGa31BmDw}pR9<>UMq+_zfNVz<#XcPsIxjhfiQH6|L@tk;n9Ab`
z(%Aow`XafG`W*Gn{LjqOc}(Ot@zp(KWb~x>Wo2Y!G&(DE<#KvbiM4xl?-b6GW2nqA
zl(4lFzV$^;wWBWwJ7u5S`Mkc_w+H0pIe*2V>srlTT2+pt$WfM7WdeLcLpi$aEh#qA
zehu{t^bHF23l2&QWqbmBIchO1B<uqU9S;Qg1p0of@7-BHEHkUXlKmxHa?=Z$p~EcG
zvhvco1z9|<WlB+@C6`zDXOS>jkd>U0o9<}Yc{sapqyovZ?Q3w<qyi-dr5C2AWw;E@
z%F9aQXtDax(M8K|!m;mU4_j{O&;N{}Z$}}gEB|*f4C&N||6eh58~8uNutU4AuAF1i
zf#F~o8EYwYvaBL@>0cLyf4%O+)bd~GrK(h8<@$7tR?9eL9ci0#G8`jm&fyqocH~_e
zD*M0I&~Go7<LI#p!5T=HLNHPb%i<;%rWgEowP4Ms%O1OO>8qsOREz7-T}@WI*{1$f
z?Et^fFpe=>L#p+LY?AsNsUdc?(3O7!tJ5h-@BQPo)X2Nda$6tFj2PL~-7YipV9TV3
zTCUXC<T#ofC#eLZeTVsu_6=fu6TeB|#x^OHAT7BtxgfppPf`A`A&&BYGXYk@ZxXQ2
zbmf2e8tLoU<7-Kj<X=go*ZIJIPNFUeH}HnXNG0lUS=PNf<Q&nFb=fz2)ukFT|6A42
z=<r2}K}?8%$4?g&X0b{kbz)X*k`9wFT_E*9p8fl>PKk9Vo-U3aj?R5E(=uE=ogJO~
zJ35y>Jygn8{C@nufWUAjIN+PMe4C*E)a_2g<Gk(^QaZHK<FEI-rrWOjUGsH^Gy1c>
zQLm2YJ3am%?R4GmT1a%YbG3Iz{jUF{PAArN{XPGVz;^_`Bk&!8?+AQH;5!1}5%`Y4
zcLcs8@Ew8g2z*E2-$&pGr;fQXDV;fE-iHwcDl-ZzA5vrVG1C9Fj?|YfrX2Gb)aYA5
zZpXNEzb&4w{TNKA4H<gX!-FjI97%P*4V`^pPCuP9pn)Y?q;*-9E^Stzq81Llt?i3P
z>V5Fsz!rg5dqUFM0|T=xaiq-xcE`=};-DGqZ<^wTk||33O%S`x7;3V{SXN*J%QuGD
zB{swxWkXC}X@ECmfDbS9LCf^v5U!7`@o@<GH5I45GcoXSHg^7)i$N#y5V($u8KZbu
za*c-{y!d!FkB?)w_?T&skDGq^*dolwvSs<Gtk1`l#(bo<=flcKfM5p!${|2ZumEo-
zE@Aq}`Z2Q?m@&P#w34grbs8(Xm(De>CX*@i$;fjG8Rm;=<;QeN^NpZ6$4Ag=MNj&$
zsXy73*-@2e588j!m@dle(t^)w<UC7>rmT{s0GqEqwy~k;RT_wx+L1_O{1B->9EtwE
zxUh|ZelQ`n-5a)tfS?6lsIBk>RSbb>_F#x>2Vw1%fykF12%e`03?=U9e#Q+ud${3v
zxhrNoa=`|dRPbkJ;_2yZJeZjaxfglR_{2r=Wgd1<=HqFPe3;eb<G8H=BdP=_kt@Ku
z)B>E?T>z0>Au@dn!7DGsmK}xI_NWjR+C_NfQG}ieMfj3Y&1`LnXH3`hWDNGpF;foR
zqnyCQRJL#<4Y;z1mLHo&J6DyFYJLtKF^#2p>VdTMnm4JMxlws_Uy9YYrdLU3<eaBZ
zVXm4qX|)OkE|e#mOYJ@lTBUH6ABXvyOE9^#7y%h#<i&_^ZL|<^VMS0&EyUDu1yEWh
zz?I|q*zl2$KnFf@vw6s?<HG4f9*%nF!Dvn{0-xsK<FFir*JPtrCmUZ}!Z9o_8k>yb
z@Wna-a>o+kWS$Jip($wePQ}IUY0!9)2EUqgIL2n6;kOJ}R*XUA=P@Wr%f#EGnHX%I
zg)nXw_`9><@HGn$24zEDn2q9nYngEa6PW30y%_sdip*4ABh6{7ryQADY6_{QY=ufX
z?pIDvRP)HZJdW=91<||~A2Ks>Cmo3ciJPqHj;T4VHZmaXYnrreuqwT=RUlyjhn7DS
z;^j2~E?Dz$=ynd09%rImCj(*gQxO`NjA?d>*xfrG&i=7znHr644UyOy6oCa7!|*0G
z6t@*d!*Y8FEcn3)@C?ELvp^WD1mJC(KN9l75xFP^3IPe=SSKSOHWh;$(~<IX2Bv=*
zgU?1;xM`A&Rc+aLUeCTeCl_m4b8)aN50bZecpA%v<pC}}8t@R8#lwN^Jfwc%;nH9}
zR-W0yH2s*wBptA0YCkD4S3R5PPS8;rVOmEA6c<x$btPF|ET`pf^QdiTJjuKdqAQb#
z1}<|aj}H#y8)QRI#+g$Q-++pGYLQR5DlJP^pf9UA^pR7Bw=c)xMf+H+)ftOOhZ0yu
z7US+zF`5pGa8Oo+<NiXpRu!S*Wg#%U5XNf@a7d*9zvK%r>~20*jm*c|1AJ(@@)5F)
zhc5$o(Av$#3f)l<D~H3VHWDrSqQT3E#e*|(2$xBKz&H`7jgnycISC8*B%>%X1;4hW
zKqyEB_em<=k4nSKpVQEwkq*DqbbQ{D4xdly2p*Jy3E~V~++NGfElObq-?n2$&QW3p
zgf-D>yQ5_IqL!Q=SJN6vC3){Er_iH$L?-dH=4=pcO(yD5>P}799VolM4Gm8<rwdUA
zv{+q>ta4SUd$0nXpT{BFm!;@)aUAqcjYZX|5-8m-MiZwPrv1c-$`oN_jS%v`6`{tz
z2s76dLSLy6JGlk;bX$PHAl9Ca=HrZ4J}&I%Bg&hP-KTh15XpnlQ!bkQL-5o-44aNc
zAme%z6eh%A<l9)B9ukk)@d+3al87|3B&<J^1SUNh_g*GrPG$<KZlvIsA*t9{m5L3o
zQ(=QNyq=PV<kM;RNhuv_KI!0;)-rzFRA!B)9b-I6iJ925i4H0pr8y^RY1@fvx;Cbg
zPOmDbUzg?4?9Xvju{DS;`V&QFx>I_+0|l7dP(+|P1(E?d{iaF3MyQhg5CyWD$f1K*
zB$#$U0=;z-oSZGejIk1!j$zM7OK`zUf?#_IJd7kLR**pFT^U*%%FypZ8I~O=18U3A
zR$Yd3Ka}C3xD5W8WeAQc!;TSU7(B2H?UunP7LG>LhA{Y~Mc~-!NcgHpW0zeF(sg2S
z{!%Pd3ghszOgvUd;*r}Dk2~QBNZp!%Ao)Zz1t-FLK_Z@9O+<@E5{@xR2rf**_?1cM
zH)<2Jq$GuTS7OJE9jnCn$u!ZJ*GGu<*OK|3YMKyMN%phLsi{1Ve!Cw>R~7})^8Q5c
zM!A#aW(WFN-G;n<%&Esf1A6;flLm#bbFqs8X^ds{<ckEQA0_DjP6C-03FKc&Alod#
z!^aZTG)Yj|AVKwA3C`S-z~H6?f@>1oxy)j9QG(OwBp7i<f>S5i>*MS*M<kejSOTv@
z61?3n!R-q{c=lox7OI9J@kBU;UXd7A5QQ*aG%9>zaP~tCp3jZNeX}?mUl<2(xp?HK
z#Y266JSw^+z&tDgPZuV@{%QhFc1wh|S0XsMiP*j{5xw;4m{{LrX5rr6Oz9^@=E#sn
z;zu5#N~V@Zg;vw8*6Cz8jK!oqm+H82l;;{mvA2EbgS<NxXFJf;3)W<<Y)<a-22_7g
zlM1X<X@IH%aYH%O+Acv`XKH8jzb0w<M>aiS(<L_TV$)JKjb~FPn}XTo&Zh2cQf1SJ
zHVGQobec`0%Y*Q{a1;ubhN3+?9D-93=+_j9sh6T)zdRbYqhs*)RSYH+$D;0iEUHrC
z(Cb(nu9(N8IWHc^cgADY=XkWbCEzSK0dXr6P<bH%iC+?6lU2*e7AG+;di7=o%u-~$
zHa?)$PwGkjX$=ifuO{si(`ltqIc0ClrPp?GwEKM^IWF;`pc8I%xt{}-R#{W-4KvEV
zsZS>sYLaQI3ROLpCw*rQaX97pyOLRZAZC9GL9afc7<4HNW<4XYmDQCePElAO6OHA&
zqVXm&2Gx&Z@H8bBrk7)p;~Ix8GvknPKMv2?`7S*-9>X`rV`^hO1o{bRU=nbi#b?#x
z8YXmmBGdhNFNW){$T-<Ppz8d3dNr?xUYxI@H8ZBu(qj^e2+gIJw_?d~XCQ^;_>jp8
zH`3Q|Ag>Z@D%@j6&-Uq4)i_PEdaOd}SLEqnkG~g_`M(EXrbY-_EJj1;P$*0d!tv29
z0$NUyFjJ0#^3f=0WJY83r)cyQ$Dr+547LZyLcJyyn_FYi*RKN~$@w^>E5~CE>zn(r
z7#S^!$J*=~Mxi8;G3WMTe%DiE${se*S^s+S$*Q5DE30T^?sU3QC85Rra%s)ZSQ5+%
zq*TU-TE%Yk^jTk$kFlnomza_5N_{d<(WEJ-RcO}^c@oO}y_kIT90|jaAT-B>K%_nz
z!9}5HTpfnyjp0~V8G+%!k?8(85-~HQpr{iKiz(4iYKcZ(Pz=`m6a(!ytp2#iVzoFH
zK0n98=4mWq*!5pG-#Dx?+Q_K%PGDZ`wPmI*QD6@DZ6L?1hiR~J4Xuf;qNwiE=|g)N
z)vnJWo#C<ctbZWAdFM@wo!zK$U0>R7W=-s*M*-|Q)~UZH)htyZgX!|LX-m7$-(@$_
zb_X9D2LB6wIA|Dv7w-b0SQ3n!Ga+bwJ_^YXMq|{DP<U~}FitBRD$B#M)hq%ZrbnR1
z>j;#Lh=gW!Bmy5rqOwmEW@bfU;i@RCVdwPciqXiLvyr*#8p}L6XUoi+q`>r3Y@lbm
z57VYTtnE}+(F5Uh66};vx3F9)eH2TTCjx0wg%7>i<3_pm4y3%;nsn}(k<9~rGFYWa
zr&LvG{6~3uP8^E+!5^QF4#C2Eu23~|!rgj%w1xLZ@hlq*@aPF%pC0f?=#Ir_&9N@m
z4B4)xnBCn3Cbq^1@-@QrVnZxCWB`2w1Kb;@kE?I>Sbfk#qKqEC9FIqNVLEbFW#bmR
zE<HSukM}q8@rnJe>HD|<v#JX5p-&O+uVU9eN<z5D3$bghko{g1Vura0Q9&ZSD;7a-
zr3k+t5uv?Ngnexy+%OcQgP_Ct2clSKfczV4MEZUG?~<Gn`*wVd|BWQ44ths!XUIIW
z^XvbAMDO|*-|ylabv7W)S`(KP{mCF{u<I!D>#ETe=I<o<E&WDu<8Y$>sr;ss-w0_t
ze=n!ajT6Eg{_@Xb*F~!8%KzQ5d-I<fJG*c6Mh@ryXzadmF#g5ZRc-xoq7-vQ6X|fs
z{4rKaDaqV4mRs4y?dT$3c2(#kU*?LsNEin$+_3-K2ahGwz4PQc{V{mfQlX{HxGux%
z8oa+=|Ie9m6C9*%b2#4G9eEf1>2GDmNi9q|QXR(nq0>^bJ0xJ~Y+Wavd5=1U#(x~V
zE}~7x;N^<`AkK6WW4nko|C$(^E1J`(VNNEOZhsrg|EmN&>D#~2h?G0Tm(qC#`h`VF
z-}G04o_9x;{{TU+tIwS;<h&Yd-u$m5`p*&cx+Ll%=+*h%*=#nog3+20z&t$KlZhOt
z$ec87q8p7z>GPpF@|(Jt-e0Mtw8Q1(Ta-upSH@9dbP(;-Br=O|C)XPebSBt_{O6jJ
z-E0FII$VoBtyZO1lN4z9Ne&$!XM;z6y>ML722Yc%5FTZT=JggRF*ZlTYE$$YWdau)
zV|dsafg5g!3!4nE+Rgyx`}9!|t&dx+dQe%Thp$0;@N3n@vn{&F%hJUT8(pMVM`OqJ
zM7S<bh3WJRY|_Yrg>4Sf_vK<^OCD-pa@p-O9^%LHp{|#YmGkoPrY#?w7y;VX39#q6
z0Fqwp_Doy>Tq+Cj)6N1Y-z`8Grx5WcbQ$e%J7!I#DN`j<Va^Z#K)3GPBAc3eD(hB9
zMK0AereYc$IX#xPogYKRS0l;ukU#1E3>vBHMrQ>Mv}dLbZ5d@j!}l3d{{!0e)?JOh
zs43Btm|$cN3dcsXP$+7JV3#NmzaH~P!Sdmlp+6L73cTT_>4~3g2BGP!I|hGph2jq`
zh>RZq+o=AKO>;&~ffG)Q?}w(jj+nj10hjmm#esA7SgXh&DklKEfDm@uI~;48qcHq_
zEZSWY@WwO=D=sJFnkW_WoOB4|GEg>c3_h@1w$=PBn7d`;@}q3LoR|YkgIsu2=fYnn
z58U!R{Q5W#+nuK~cBRRTO^zFrSgy~U5Bx|I{Vvch^X-&lxsr0^=Fyd9Q>e6hENwE*
zCWnL=;`s;AxicWOV=lB`*`7+qSW(0=6Kc`WBjX+#6n{pUS{voadCvgQ4JY)wXpeEV
zy>Mlf6)u}tpnsSt7P^`sLEZ$@)lIPIsWA=|7(-=`5t6SNqUDJJtW^x)FkBy>*Xm(h
zFFiOO)x~U~F7CPOLM*F`L5Bixgx&T&Tp5d_=h!WyyD6A4Fas?unYeT>8%tetVd|KN
zA+PfgP{~E6CJ%!oJoIkhq0NPlLlQpT9p)qdD<3ue@{tsmk9`yJ(Y!Vv(F2z=JGdFl
zkHyZ6hNc?RKi~<;%sD}yw``#|8<tVQ#aVPXV*<&1%BQ%@B#Ik5ik__<N_|QPQpEax
z6x3)-f=?E-<GLX+LT%DL+>H**QKbB1G8FXC2*rN-aKEYz{~~ov^6rMSx!sTgHP~EI
zLrJR|<~>z|)j>7LSF2&qWHp?cpoa8HHRLZ=gK?c29<yb^7t~O3R}DeW)bQc08t!~i
zgYD;0c*=gFc@9p6fl>xwnT3_299TWd#rflTFwEdW`vDi-86KMF^04eG4}-h$;lwU+
z^ON|HjOQbB1-o6hmyfNN_*nCZ50&5eNKnj2=GFyF`EQ}jP34|U&m=kK!MHo5+;osU
zj;<%(gdb_*+8-!u@;Is=mqjxoqN%xPBppcyX`gqYms$37watnqC7aTh3Oy3UYS6-8
zm1)gQIkM)q`3!#F;#FAw#Y^pBt5-^0yB9A(2KP1Pu%uWX_6roCEK)>@vl6afP(nbo
zGP++@hL5KT9xPG8FK<<l<*f>%DXIuQp^Ba|-C*e24I|RJp?Ov}==6=nqxLj7sb|8c
zM>hJt&w-759^~2W{e=Q9W}oMxN{@%~kvs@QJUA}qVfJnw-e2S)`4JDBKJrkez{g-6
zKCI07D759H#E}mXyVYw{8^~ne>dv%mmSYZ1y-(Y1>PhcL4J8e(rp524)5)IYWR#Oj
z-KNIUJEcH+EA%FAi7Q#0v8T6=J?WLel+N(=Xr+k;EzVV@VPSH#F1gjG>G*q}=c^BT
z%{zb7t9bKcujOA~d;OI3n^&eztC!-B9K@^1;s+x+%)cjxmqK~8ype}yi~<htP{1Q)
zMH~rN#FeFrXnCLrIV&aP#VbL5z7i^qD`CS2CA6rd;(2r~RL-&?oU{38NiRUMxDZAC
zitu!E5l*TL@i<Hft4bl7cL@=4TZqY@gfP+;!N67o@({slgb4BxBIu-uFez7rHAN!q
zEfs;!3}qUoIx~wWn=%<yip;2ck7&L5QOaaCP-@eBI#N4{h9wo!j=~grdp(riBo8OG
z%pp`Utv?-GU`J!e^q_msj7if>mz-qPX;7sSDJ_zv8OmRM`V0GEb3}g}p5Y9;UQTe_
z*B1-U_rk0hR&aV{j*&Sg*sW}Yxf2Yqd6GU(p3uWWb3Hs1=wk2@9URlvflHb;PW+;U
zq`q3<uh+zo0h)MxN(0O1rD4LjT#SFu#SS*C<&Y;pih3bVM;2jkln`^QMCi?a&e{2j
z5wk^%!K%gB8c~eNi;A)DQZXLLmB8Am1bre)Fs!%)Ulx`?ZEFceo-D!D4-*;Vs%YkV
zh%>V;sT*^%_%4|Q?WVfPD@oKamn?sqLXW$Z($dK}^b^I>g_!~LX&Go+mJ3bZ*@sfy
zt?2Sv6N+rmrN%4j6quz<@}+V#>~5>i`E9;<6z7A_y@z0YfeYqT+Cw$j8Y{9bA@kM(
z=0z6JeP9lwOJ)dT{ch?8W6akug3Cq&^qr!Q^rd=OeODK+2I)drr-NlqI*>oC4Y@dN
z7_fe1#@kHX*_4B#&v__k{pJ^S0kqh#+4Z@F$j&UnQB5I~eiUN8k_c5}L@3`YLY=G_
zCw#={Uo3{{MmE1DM#C2|PMQ|ul1DKbM;9aK#8O6OST57N&4F<%l4ZDl*XfSQelp*&
zfsVCS(ZlBHG$&a?a~rbh_vUEYUq6zJw}Ku!x>C&pd$RkXCk^~!N&*{wdaj^JmnW;x
zpo#MIbx*s`>mXlfuJvT~oSk3&>~QI$72NMzVut}6CcbTfLGR2k;EV}QTs6WYGQ=k>
z15E3wk7r4G7_>(hQ>=CIX_*c#*y`ZcUTypqrH#RDTIjW23tn4NFzD16v<K!OV;>jW
zo$|4xNr3q}g>Y^uM8uLJME4f5zDbDLZ9=RJ5#iPX5msMhw*pvP7CDNsFH8)3p%|m*
ziLtd-jMjQK)O1}8<+62*+mJlw{3bi5<+&VFeW!t{uN<RkTAS(3=fyOqwvrlWmeWR2
z9tm3GDE(#-ovtMMc*LEyxHwYVQX9&@WKN~O7*OwQEiyWwN;?-SQ1mqpMa2!oq2_*=
z25S_4GQ+D-Q@sA!1kL)!*#6J}ZR{3E|9QGNvRVf>TC_1NMjO+9)q;Jl7L3%jaBYVs
zJ|=6TMN1R=k7~e5qJd&>4Ol8`!1hKO_E%-&UNaZ2xAQTYU2{!;%7%Rp7Qrb|2)VmL
zHY_T_`Pm|5UJ}7YQH%%9VhAF{7*ix>{j?aj){1d(uNaTdu$bHuW8o7q6yJ$q=d+Q?
z{w<xEJKKs0ULeDS?YKwQTkGk1N)3H|u!s(nOry&8rIh<Lo7(DQ$nRkQZD&Ko@g8ne
zccd?+z?w=XnbDzX`eYxdN$+>7(3%bMH1u=3&(K&m<mcMsSy^}ZPBOu>ea1LE-UuHO
z4UiYA2j%-Z@O-3=R6A|lo~ecL8d{jKN)xL|6aG&%;5J7CzFr#G`$iq2I&~b(S4V`W
zI!ZOv@wNXL$g*p#!n1tXrxc*a+d?d4LuR?#giwSCr8Oc%aKzZluzf5M;}^F7b!^|~
zHjA-So(;9JwldkV7;K;&dO^iF98-+T>BZQ<D~3YcUPgu=#|){rWHOwTnP#KMWbAW<
zo^xtQZC({Q@1IT+gUTtcB9Df<#FK|?FnyXxbg#mlBqto`nS~7<jxZ;ONCSGJuSFZi
zs#12U0tM7?sMI(f^%DcJ;sS#m@4TTkWe5iJamN+)0f;biLiKzHgpam|N{Af}jO~Tx
z_3XERr8OQdu|mp_9*}A8j%W8RAb;B&qhFc9+0YCTqfAk>$OHi|+3ojXNk|`)f!))x
zaf6+Qje7Cmv7C?BEOkP;MgV>$OB+xu#DU5}1aOM5C9MdX_7|a4Rfy-oZ0LTb5Yvte
zafO|anr%dQ!)~|9uyfC<$-9`l3Q^4FgZ9j<HWg+9w~^?FBQ#36jub{NCS#8ov~OYs
z*~sx|)`xg95C)S%8_^>Z4+<AK(u#{VRP))KN*)@}_L*8V`JpN)99E!qjtt$@3c$Nv
z1EH<q2KUuYSRm*Fo13->?{AHYh92;jv&4qc76_{~LufzNu3j1AU850Zd^LoOlOYBb
z8NlkKK2G=6$D>7h$THHy@2gpVJD6R^dL}@4l?{={vg`6h9@b~%<HdCW9PSig>Q9Au
z;a!9QtQ}3Z7Q(hf2wpwA-Pq2CUtL7dO%@?_st69OE)6-&`j<x{++(*I<MdfOvJ->%
z_-DrbQZ)0q$bs=#ugnMsG|~XWqvYaSM-?j<6DM;9t>stHc$S7!{V1NArUjE}cZSm4
zJxFb;Bb~c$Lo>dZ(|i__`?IvD;;|~7u2-NVGBR{!svjop^uRT5SDX!W1aqo4+-BN9
z=b{xF7F%M~4hu}WZwARgQ@pul4B>7gY&v8J|0V-8nHnH2Lm$fz>tTh19$Yu-Vx_At
z4xQBjH&+L5uP0#7=uFra<zaawANRZk*gv}f_2Ub%XmAlWpDDr#FCjdZ3UTha5OoeB
z1g49yalQ!Y2Sj+$z|KRgZ8eyQ;qNBK<p43ZB#GfB5JRc|I;OiyEaSbR52JQTo(X+%
zpX|Krsam_1-cP8eMouN2vo5EsQ*8JmDUL=nK_q196?_AC3K2U{^IdC}YGF>UCaj*^
z&?FB}RT{z4D^Bq_WU(w95sLoUr8tz`P9F^aFn5faIRNrpCyZI@0Q1LvP*KqvuS#s;
zw#OO`);)1xYY%*lv&0<R?nqX%V8fl}cwk|Mts_iPIN1a{uNq@lKV!U&4aW*M)^<WO
zusR?cjcvK;e~*i!3?CyR^YPMIfFB;Qdc)4mieC$`EUOR&Ckrvkv<RzN8(O!v2r5nN
z++-p|SO7aeiG`T3f}OjM3UPXG9kXL@4kI4Wml5qyVT7%ZX?)fRIwIOk+P9X_5!;z$
z?J%CAV)!&RDUnRig^;D)P+Ds+khc5xqhXV6$#by<If)HvrJ*)Cr*)(49*WdkAVYR`
z!(pOFaL^o#HQU`#5;_3FVkaDZ;D8tN?P2?q9b#Vh!rnAn4FBC4r?>P(@B}M}c|GuB
zz9srk?v7;}Eim?;InG;|!&792@Lx@#GRzd4PKCl7@vw7E#XQ?F@U73n%#}HqtD1*7
z)?B=j=fPwjyM7<d$Ii=qJoL`T^%eP8_aPq*J}f<JiU6073lPaEK*@jt_8X`Gy`~mm
IdrblU3vFTTqW}N^

literal 0
HcmV?d00001

diff --git a/src/LIB/RAD/ecrad-1.4.0_mnh/data/YHCO_ice_scattering_rrtm.nc b/src/LIB/RAD/ecrad-1.4.0_mnh/data/YHCO_ice_scattering_rrtm.nc
new file mode 100644
index 0000000000000000000000000000000000000000..3b3ad6fef7c98f08bc9d2811d298c2dc0bc1911b
GIT binary patch
literal 16067
zcmeHt2UHZ<)@~!BfaIKmfJ6}xQ9<}%7a&N^C^@I5yMgXTRA@1wV#0u88go|6sOabz
zF(=Gn9COw&C-kotdglIf?_2NwZ@szey|t{>yLW|Cr}p{I-n;6Y^Gyiy4^&k)QtqUv
z$dPT|I<^Bo?Ppv%E({#Kc40)2e{?66t)04mKPTU%<C5%_ziQ=9964zb4Gro4Z3XE+
zhcl~PaX6i2J6(*Wdvcw>7iRYqv-ml=xoKtPy*srR>Eh5{pyQCEpT^5_Pvh}<X<2-J
zX*T=cUcOwuy*x*?L;tC~UGoyYo3GkWPEJR9zLT7soLWZ#nU>d)TC8bl-eH^tM_-Yn
zFJVh5d~b`qa(i2T?y!A&$Njn{-*1rb%=s$@vSSrTX;FDj7mlK|C===z6T{JFPf4+n
zw#z>_EFe54I3hePhVcsx;HV5?Az^P&Xul!cFD&4DdEbumvH1mM7VIy>f|ny?`~xkr
z3yO1iWd(emMW#q-!Q%`6yd=&jEy&2^<~Un)Y_9GcsX?9C`qelpQiIaMbA(yhx$gc2
z#Rb_M4b~spn`rUfIQE(B)`M5M<)1MOXursj>3<i)$PRz_{}n@z{{J%!+r7JWah*{+
zhF)@V)>7y?v5qM7KN*I9J?_BN;$QfsveaYc9ol>A<p@PBX`S+N90RH5a17Mjwai01
z{jWXr`@`iqI;=ym0U~n<MjBxSybNJZ*?)HnHhg69DAUqMNt>yX-0r*Ttar0j{ps4F
z!7;HML$-v}>-AZY<{fDu%4(45?_hN#CF!|;JeCG|x2YbhBAB2MS#xD!){ng=-Bh^J
zV3X&lb6lhrL<a;0L<fX30cqbY@L;QyT9BO~%qYtd{%Oh|F~m{)Zy~^1_}v2b9-00l
z*GM17ZeNBaO8%8aIvqFs=PZ(0xSBt3h}5F?lx1z#Zf8)tlEpWBBy$b9|E+7NwdbO=
za3->pFUTns7O+ksO=8_yk#>_fr&OAOyxpAGq{OBZZ+GW@&aO`R*|~kZU7cOsoLws)
z?yY1C{@DIsKwuyf5&GR*zQ3UV^zDwo<GkoFQrflB?XTx$uUpx9+55UZ8U0z_z_$JV
zj*S0DJ+kw%4+&X4t#<9@W&i0NNlbS9L;pbF2LeA3_<_I=1b!g!1A!k1{6OFb0zVM=
zfxr(0ejxDgBXIFp6Vp^(!f<c7GCCJk7_XLx#5F%mE8W+T`Pey>WjT?GE>w}f7LUwF
zrIUqk9KBl`LgDcY<-PGD!`MF59NmjNFWJ+=Q`Yp^-hz7FFedAp`V_j@8ScyM;pSw6
z$?q-kF~c0c!4xSHW4Ngsq2PA|Y%(;!R-rz;-{?VoiXN&x^<ei#7g<fZI9siYD^cuz
z4_&PJsssJ2I&j~igGF<6aP3hH9vw)=l+PJ(5ai&(ojfdRF2G~kB0Oif@G&aJ<lV*C
z?Zv~w#XKax;bBQXKIBUIP+iK$^pkv?f5pcxZ2{U`1aJ%yV0WefnnMLx{bm+Z#)LE3
z&1THPHQ(su`wMiu%MN<ocQJijK9#O)9!0ail+aJxGAVFSJUx09Lc=W>GUs`b+Z=b=
zRN_nvuG!J}y;dYpF{fX*8quZodgS<Ai?ZFLF#19m23Ll__IeO@P7K8IN`G9fCm3G#
zMQswW#KZ?R%HAk39)MFpUhtpfiHuwQ5#rY$)$95p)uJDE%<GG7rhTz|xd&dkdBEgQ
zABY@N@m`dLxk0%w8jz1k%L<^ls}PEdif}KRi*%)8EU7PsmmUx2YI!LBz{AaCKD2l7
z@jycW4YvI@8w6PVR)81%N{~0S1Vgr#K>b+>8r`-syT^!`70-J!s>9TnsQM?gCAx*o
z-mWM9n+7V1uchgtk(B&|Ppv&ODd|W&?b8V(*X90n(8-(T5A~pJIWDy0sskk*w4nrT
zOY&<rq1H_X^uD|jdI=-ouQVJf$13n_%1~U0DMzxA7@ZyqG5T^DvL2S=sb(qUBTMjn
zp#XND`PiJnN7xx2$^&?C-B*lLC`Ri6F7AYIvFTC~IC(|rW*CXAZm~Fwcr06&fLA*c
z@pSrN1cwem@XH}MKRO91+R1pkA{leNQ&4t31sifwapzSk#tlou@Goh&I5ZvmUZkTa
zBLixuGVoiUOjtBzqP*#6rv6wdGcCF|bF;G=GcfWoU5hwEm1Zm^0~%<;j)|myMMA&4
z=h5TH3{q{2Befx+)bO66u@eVS(XBpIaIqJK$2gFlzYQI3upsLo6MFA&KyC}mu_T}j
zy5W2nPAS5<mw9Lzn~h$687Om0MHK5OheHM<Cn*7=R>$FkeJmcIjmGdfQE)Gh#7_kg
zSj!DZpQ<pdn;(ixheNRX^B}l;4#ElhNbs9t5wSc0U(XLl-7>bW{weTUnTpk~)8KEA
zfidQp*rA&RL0cA1U&u!5!W;}6oD0RSdDyxt50jCP(I@gz7+HX};{~{bLTp`Eh)Jot
z7|C%yv#80Hu^+0!T-AL@3u_P2j;(8nuRDjD_Dvwoo>g>gL@_0QP9yijv9#Yigo<7f
zIW%}t+c$Tr>)eZ+r`VIxL~EMz$eczmGNup9^l5`u6>fH`#A>w>a8MtP*;>P}#-IXi
z=0lO^T#iX`V)U3T!k(u>q(%!7dZ-L1s0_s|rAUq~McTa*tQl5<?b;=<+9<%f2my@P
z9Dvz;9QY85@tq^FW=Rw#Z;!^Cni$Bt$70UWSRg76;b-Gu7atGn`|(&`m4KN>iD=%I
zi2AU>*mG+zM)8JV#k(PJu1G@oyClT%lCk_ZRw7wkHMcW;@_EcYB^Rb9svFZM@jjV&
z9Hi5$){yCd*|f0h1R8UH1U*dT()rh^#9W9W_q0JY*NQ0hvL`(aaVLjlXR3N>M?UYY
zXl<%F9WXVf$@cnWc&id2$49_(`*7&27zW>Y6;Q1oipH7cSUg{hq`e{(e-)x{j1c)d
z%W$`68Fs8^@$e|clS3siNh*Op8}A`A1c>Y-z=f-P%o)u`YA-(G<ijwvG6M5%MPkDF
zD9o7}jb^79>}`s{3y)aTZH&cx*EooOiGzAlJWM~t<IRi&L^>sczb_GeBL}1L+F<O@
z9D)J&hTu|W5+bf8K@gq{U6-wl{Xi~rp~;yUnWoG*JibdO^bXL6L#xSWY(4eijwg##
z!|8l<5uN&$Ldox=sYgRFeJk^$ngN~^^h<AAcB>~vjIyKkGpy*9f;s7)G@_or>XF3-
z3C8~{!M>FeG)<FWLWKlrnG)E9NHEY{f+TYZ>Qy8d{kaNHo>XDZjVkz@slv<yRVZkx
z!jiRBa9vadm04BjIk5^;Bvr5%RpDb{6)t3j;+j)9ymTY*dSL{1e2jpxPb7wjBC&i~
zBt(}YF~uYb#T8K)^(YG3Ink`oL?gB!22D?6u$>o+pP$8ommi0#H{wti7LS^ycsw^w
zfU*1*rbeNVdCK%;cH}BCLA7`2QMdh+`eYS-JwA)%PuEaG%P{&dypYxnN+Fl{XyV8R
z)9t&y^l)>3av9K@EYo_@TzNa<>RVC3Br_TjZ$yb)Ju>^mV)9M`gVz$YzL22ysRX)@
zBq+Qm!KvRQm~&Htlh-7OxFUh_B?(@)O7QI*i`iKTQcp{;@}z`~a|zOqvE^GN2s|Re
z;X@LHACzEMQz*if!?ACDI9gOA&^s&wQ-?+1)B@IbevW|6;|LfyN5Xq@BqnmA@OxDh
z=E+54;FxHXcaFiHkuiw-5Q8U$vEW~h#W&wL99b3zPUm<$dbW|Vo}I@W*zd?lauk`j
z8n<cUoqZJlbtO3}&!S?78oIoqg2tB@(2mq(8Y79KOrIdK?&V9lO8x2Sd^eig;zY}Y
zw$y!UcT!R|qqkQKDf_K1n`<Qa`WN}uar~LpZ&-Q6$`w`)v$Bbm`K*j*MaW79E0L^#
z6=zn=S<z%go|U|Jp{TeShM@3pc-DpE^3UN=xE+pEg$S&1ihyl$1gu#9S@kXgE*X)i
zyBLW(kx}qD5rvSzXl&gRjlHfh__{a-ebizRD2|2w)mS|2ABS4IMyB`U9HwWp1M_-$
zS7!U~x9GI;K5C6yNfWAP(&ugCXwLPagroWNsU?ZNeu<=$M+Z{rK~R?!{b&o^Xh?|@
z&9Jhi82|2Md&rc0=Ni(*ow~%~jKtq1D^7+&C<?<bEn!&uDGYaQ!_hf596Lqf@SPG4
zAJ%`|?}TH%djx)(7l8}vk!YC6V$(GWZ%0SrOIsAC2&19=AQ~2NY%Kj8gI23pG>nPG
zmGNtsWivCFuf}#v?F9wqtoKbC^JWjdTfLmTEoYEJ?O3{TMogby<`NkUp&)()Ij02D
z1F;XC6!xVFwyt!K=Scf|_n_(7mX!J0gwC`YQ0h0GzZH|}Hz7D$7K)#?hT_7jP}Ez8
z;j?!b<Oi}kEDY7zVOTYX^`oa@*p?6u*R$c!kBY#UbF3dlN22SQNSp|c!q~l0c<2(1
z!waIJ*d+#=*nD%tW(C8)mCQ8sv|$>r$TRbjuF<aT%{2PjQgU5dM{5JeknchvO_a|e
zgRDdv`yh<$Kl#%WS8vJ<@*ov;7aE!CK&$)N(0ic;1*(})(?@+WcK&-Y=@c~xE*l3S
z=ffbZbqawt6M{`4Ay^O|f-x~6a4ZOc*}@R4WU;%N9g4KOp}1EVhOe)}P*)KSPFpyJ
zjflX+_Yt^T9Ep@)Bk>|63XNN$uyyQ8CWxQPjFRiY2y;6#ihZurs;u2)TeyrC7)_`1
ztH)3<i70z=4&A<-NOqy&bR{ByOh$Oqs9Fz-h;SjRl@2sxxDB1yWkE~gO-RBR(CPBO
z7n3_D0-)VH5S)rYY^x2#jk-XDP7B1sDS=R*5r{ME0#WiL5blu!VRvRA)};o)<aH3H
zjSohDlR=1GHwcIOhTzGb5ET1_!ggmUbe+Qxu^<febQ81qaxqgj&W(w>rOMoCdrVtm
zTIh<`1~N~WOM*kSbou^B5*-jw)|*Un8J<8pwujLIj{vf2^rpOL9^`e(g<ShN(kJ~M
z6q{;E$G)1-$L9uQKQ9>Sp9bLVHg_EIbH?*<d$b5_Fd(8kE*&z*$aAKsPB1}>i!s6j
zjWDdn5N!_(pq*xbL-+MzIb0vrhWaqtr-zlrdhoE+gUVH1)XmmK)u&j9_9ek=cqS@l
z<l#>LA_SZ%#sPLbW7Vq!K2!=*wK5D}Q-(SlA)>|$;dn)elO`gJiWi}Nya=ndityVt
z5u3Y3=xQ#8-T*N+#f#BS(Bb?8QLJNu{5xwzI==pQNls}_?W6I(k>u1)@8EVUndf(m
z{{KhxuJ7ynGR{#)1ya?PHed87gXETd`=wx63E47#2f=T_cZwT_llV{9FS`T>N$dH0
zIc*-CNamMM|J-*nQk6{qyM6cjf2!~7ztbBzod2V~`_94m7k$^TY4#W?=3R`W-68i!
zUv<sQ;AOMiN*T99M!u94=pbKm%Vi{tJ*O`@{_UN|lIdQ1^IZPuJZq`ZQf8bi@MN9$
z*W>>=Gj4QGY26%-uV%ZJ(VzZSW}Gy_q&?LhtZzFarLf%smd+;I>Cl$z6>k4y=gEjR
z?VZOh|ARQwL5!6VYyLGcmRnxe;bAUDXRm$l%m1qcJ?YcG(})z?&6m=7`Ul4*N}u#s
zf}U@Ck^cZePuAv+2XdYdH+}rCEc(w8^kf#v2zu*-uQ$G%Fp=?DpTaDU@5xN^=*C=a
zxKF3U4wCZOHPm&@Y+7}20xh{wMS%-=G;?qU8G6Oj-3_56jP$3Q*9Q<|>p=rFT&QfM
z11-q0p;H?zXnLFp<pdj$!eJ{YKeR#f=<cYTXND>#_BA`s0Kcx#gL%F#x-8Yf?+!Y+
zcU>C~mTIFwsEwUL+Sp~L4aav{7=Bs{)@!wpF@e>ET5t%{f{vROMi^+pu_XfDYZH(+
zH3=FTX&5M%iGiir7%?&zZ$0vnUS9z5=0fcMsR(oOxnQ)55xu+^bF6rHP{)JaYaXTt
z^Kr4BkLc^{E8AEAk5~Z)j~BpxrV?ZO#Em(=#f_P&W5)1(Ix|tYMyZKAY2CTSq|<vk
z^$D#Z`>&PcH%UN|9yugTN}`eV(WJF%5OvOHsI_?jU76-VZ!Wvgy77+WSl@%1KUmV^
zBf&8G7K<T`ahSI~3iF*q(7KXA<ADKKHnlIhI=SOPS2sk>azR*iFTDH78GEnvghBV7
zSXk<Wv=&Erb$7%ki32-ywnx_>dpzE3hZWuJ;2LR-E5Uxq*cFP>(~;2oDGFArqmbbq
z4Tq{2Xi^;ddndpza4>cZNJ5QD3c77bg}!Av_T*=v$AnDm8JC5-IoW7&$U*d#94L!&
zp(U4xL))e?zH+0Op^1S^vZ^jKd&Dz(dHN`Yt=&YEzbqu5@zaR<j3utwP)av0pb<-x
zDaAONW-15M5ko(+>E=oEH}@ummpv(TtsS+UwW7v2b9(4tOwSrz&|k?7%s5+oE;fQr
zr6xP%QiI!gH7tIqhM3Q4*tJIu2Kj18GgQO%o2odoOBJDORIzQ9DukO<F?XLTRL`qo
z+<jG4epH28H#M9#Qo~d~PdKNBVferpDD92INB6-<UXcR585uY`DI3Wfa-rNe9~DCi
zQ0`HPoSTJE%V4KGj}}4Rk_#@Ei{Bc#IP#2((YD1nlTwV%b?h|uFU81H-o#A&y^6WJ
ztq-#|UWJJ&c}M~AN9gxQ>#2XQxpZb(E%_WBNo%JF==h;b(uz)?)ahYVsTV-krhC)0
zRu2l=;X)bKj#Q-3gUmuL$>6C88C^4=8e4tnOwqur*9zEBA_un@a=2D5kE7E&qajEE
z=T9mi%dZPG8oS_~epiGJ>xz@DT@h@n2$U)!f4d?o-Y7!7rxG@$E8+T7CCE1`q3*U4
z#%v8hUAH)_6el6`c^dY-$i!ykz}hAc6PD$}VOIenMifH6EBi^yE5d@kMNsL?#TQ>L
z94fe&y`GDtn_T#HD~5x6F&qaMW5%#zT%J=58>eMVl1diiyVH)DpV@_(`1)76_+c-t
z{&^+6wwpyo6*bgWG>n>k3+bzV3hg(KrkVpmbabvSRi^f*nK#|2@uL&XnQcoimUkx?
zJu@1A-H^J!(517RANf3={?ccq=(P`qeDIkiZu9wajDz&-PDofG5Br}w<5IQ)GL*aE
z%bYGK(d~+$(X8C-it*lxm^)b!n=UHih@lczCMe;<G$r^SR>GI}O0eA(fy~Slv{q%}
z+4U?$cF)1*ow?X+laDQ)1*rU7fXSl@VfU~Qex60pA6W!uPZ0`V7a`7yi&<e@oG<0#
z%nU9PHgREgf{W(6T+I8-g?D`-^VY$hk-So51{B|=t6%q%)}qyP{^waV)NLG%^%K#O
z(^*ub7*8Kh45C+sepH{;pRT32(Kb^jx;v}~tsiVjIx9@*O(Z)l5UxYBXKK)BV^xZD
z*x}Rg_JU7a=zX8^DKC9A<=*=&ar*4zrTf)qNp}wPK6B7*og9<{J7M<ePVfqm$BA9?
z?6h@f_>Sm|LuWdp+DrkItH8#R0!r>H;I(-dtc~k}<!M7u9hV25LIDif_E(Q$+x80&
z-$DhbjxE6z<x+&sEJccZ8ThGX5U(mj_1!W&&=sN?LJY}fhdpD3kSt}V!}kht|FRJL
zmqI9Y5#dgC2-A04AhRgLhFQ;*XC@zSr4L10$neE{3VA<?-WOL>j$a9Nd7nv7+u~{Q
zkx&};j3KVo09xtRhvvL?raM9Qq+()C2}$NO?1K@Vc%Vni&9zB!i8H#C_CfwESG-tk
zkB8$eaeuur^0N)l?Vuix-_gZAKV5ik&_R=#4w`3b!%SNnC+BHlbT2I^pVEXxtcjJj
znsB_OftH0DD9Y4;r-KH>U(?ZaB^NCUnK*2f0}IbWeCWo*tEmF?W~cdwFDk>~IYLy#
zir{-kgmKYg<g6FN_p=ypP!8wHa_rnvj^MlH*rYxbtsX;BK4d7|hY!Wj<I|ZnC}XaR
zoS9?SJ28p6mno`f53PB&imv3&roEyG)H0-!_NH^G%l=e)JR^p9-3QUd=YG_<!jlH+
zxzlSCXS%!0j^Y+t(f%)HG-{g>>F?De|04sznZ)4uGY@w8#Q}R0t#Ph*cT``ofUcrB
zDpX9-YmhOf{$_~u(+1#j3=o;B52Z_bkWbTt<4au_jn{>MT@GD-Tn7V(>R_pp4i-Pu
zMw4p>p1e)Pi6zN6qn?G9wtSpqV<1A6k88&SP%th<q7FOtf4U5u3?U3I3bCxO2xBIS
zu;#1?iRxmE3J~LvNQ{`pVti{6Bl?*bTh+?Zb>e1btYZaZveAhd)<=m+@_S73M~+e3
z<c;LoHjlWSCllvQHEmy9O8HB&sq$MQ-8YD!ev&}ycHf7_8T6yCFI{O<j1$#6+EQMA
zcS>qAC68N%WG@cF(Dh!hKjwrZ);-{_VGZ>MmWX#Z$FOD-q}?_`<y1q=nr478hxKv7
zLLXN3dYEsaha=6pxSyqqW?fy(IHiNQNjj(r(?LH?9X!0qe#1=7z~-)L=zlUDQwsBN
z>}(O5c|4><3lMa?1O~sBV$||7{NgFZo<<?&sfrLiM1<_QBDk@bXvv8&(NT<9QDWRH
z6GLH!7)6a@R2&r}C1p3WGKj~-9&}`;o>O7EXgni&c$CKUXrw^-dE}ftiO!Z+lmF8a
zdhC!*?#+o5@HU*r<OP!QQ6C!iq%U1)bS23*N1EKygA~44lHp`iIw3G5&pIzWoM(@C
z#vEJfOcD0l7zZXA!d9ZsPVMPo(I*{fd+A_yqc$!DYGcVOEr?cWVPK{f9+_%k(lt%^
z&)0-zh9-tNXkz3C4V*ltfnOJCz<OpThN-7xn|mH2cW}{p5+6&_OJKgQ6c)df;hb2A
zho6MdOcSwFN+RrjEW#mEF-(KSSj`h-;Z!kxT_=XY5i!o)6r=2o7-5R+Qj<|Rc6;n*
zj!wv9Tm_DdX{`#gZSG^5-*AM&%GQ&3djma)t);7rN7B;@0i~SIB=f8UdNMzZ=II2`
zgVEl!be9K})Vt8;HV3M`ZbNofmbCPs2?cLAAge8Xaa`91FLs*Z*HRN`!U&<a^bumH
zhwW!{VDU&BPyMw~cSH*<sahzI)552<npmHri9$V1yu6@+wevJ^DOUrlJv30Ptbs|t
zsbgxhI#O%0u%Ley7FHEPb2ksM1|=B7;?eY~3<u8$v7awOH#U!WdWtclN{kiT#IU|2
zMuSp0PB@g~<G^yf%_v7g1-k?}r5q0zmt#GfU(^qlWAm5;Oh!f?GpfEX({+ay^ZND+
z5>y{2<(7>U+%TU+M<>(r6{G0YzA{prm_tr2gX#A92r}tCkUXb@PVMYRtz+D%`ic{6
zTxLsQm%7u(Y%@w^j3_TV1^f52%Tdb$akU0WukD9U7QJC~*co14j!5C#p|)QSXz;Dk
z_@z5s&RAl{WeXH_wZOYbb8KC2hDQcw_<5=+iZo2Ip}_>>986%o-x#{=bo;0mF&N^I
z!Y+SjLVZOpCfq2%;RG%w2k|g^FCVL~3b3`Q1R43IXn9|XkHuy9cBTyR>{MF-PY90<
zLOi`M#7;vI<_;7gtVD!-c6oO7mu-w^O&KHj*q^Df)MRFxzM!sU$EoPVMtYMzpCSq-
zlOS*u9q&;_6-GJqP&AmXO^Tp-j{>RO7gV3#k5uj42$fEBA;gyQM|Y>y>^I7mPlhx+
zDF~zAdqMAaceZbOV##w`G?-XJ+`|&bn#}O-q6y@9#yB$A2=${3(et7KrUV$^9GhG8
z;`EXKTn`)S^sw1m4_d$LV$2#{a0+yB&`}pHliBHYcIiE*AqPj;<?j5k#o+GeqcyDr
z&C^Q3u<gI2DnxpX5MS;I5$h_#X`YDXbBN$`N`#s>A`H?KBh_7utOzk$^TbeN{V8oa
zJLRaji8&Q6V(vBdWiCEaXHL{SBa5hGM8_IQQa+FB7EGdz4b>#ysgyKzv*|`vA~`M!
zr=V_uG*Re76)XEvki?a0e{-Y}yL!-yPnJ|R%arm*8d76j0D5V7Vcaw~OzPqYqpvow
zOX&{JDsw1Xo5C~67!4{$n5bchmq7;jyi^|o1AV+(q=)GKdboH)7qe=0F~OIm;(XSD
zf3ps1CG2~9AS>NsurDJWGyL<QbB6s+ZDsT8Edg@b_i#>j8JcZ{(Ay@&Aal0g+5QjN
zB7)o_5j0pFR<cvANEPG3a51v$#R%LW#)ku99JwUM^G9rciCoK6-WkfYnfGCi{-VY-
zSw1BvqocHH#s)GdnM>gZYiZfVk)$|HKr41+lDuyMX;p`jnPLEKtnsGk6CQMAvkPsp
zbfkOixBY-fOUig-Lf;-3(6^7Vu&NG)*BU?YCiI7ET5nVpu;aXb4oEoO1JmuTp>N+E
zedk&rIL#ctW||>yhAHOXF@b2139Q%{c^qhr*o#Kk%rn9nH6t9{W{6{%hWJt&gh9{Z
z@v46+mTb;~`;t7A%_v07C3YM$kB8?k_!#+IfaonH2#H{q;ZKy}Xzw!2oWN2Q&Xpll
zQ;4`AA@nPR*u6{$%M(Jpd?Cbo4H2eXX=ZfZRWNdV59WxL2BWs+8P(oCN<8OA`sLMJ
zTH13G4RWofygnsV5tT(}k0(&`w=lXfEP&>|^rq=O`;v*8D}@eqq_V^wr2dm7)xeaB
zdl`~WL@)%C{ZMkWKdigCV}D{VNYWgk-^~v40XDdyYlR1SmMGD;fK4}Z7`vKb#&A<i
zyJdoe!6p#=W{hJZV>s#<<KZqNH0B$DYhnc73x-(La}YkAjzt?g-ccElj@vC+Xq}pi
zvX*=#Oe;jz%Occqim~iUF+NV<VS^DLPp0#c{DBWImZmXvt^i}M3$RPC1m%NDu%D$S
Ml`biPXG;nG3sUp?M*si-

literal 0
HcmV?d00001

diff --git a/src/LIB/RAD/ecrad-1.4.0_mnh/data/YPLA_ice_scattering_rrtm.nc b/src/LIB/RAD/ecrad-1.4.0_mnh/data/YPLA_ice_scattering_rrtm.nc
new file mode 100644
index 0000000000000000000000000000000000000000..92baaf1e9b331dfc351a01b1841f2bd813bd07d3
GIT binary patch
literal 16067
zcmeHt2UHZ<)@}omo19Tef<zHSQAwY@i>RUqD5B(?Izu-mTE~Q8bj&$t#~cP7#fWLl
zIXmV!<}sk6|7oFT=AV1tdiQ_p-Mijf%Ub>Q>8d()YM<}yy}Rn1Z%k}tw26_Uk+y*W
zPqTe#+V)2@9tmhV5i|V9`8{JJ6SetUwB5cQ)6r`>rMcwqrI9vIOI^j>T>X7pSN+Z7
zO=~Q8yk?podd=0pw3>ab%zhc<it`H!GUdu3?Zzs4K8+Qcc6oN0qTKeGBC#klS1gw0
zvG2zEN}a~~Jd>vK&mugqFy*WHCY`jjtkuW0wY0QMn<{8Z9c#747Ou`s#<}q940v`b
zR!jG5Uv!Kb`?9Oa_Q_4x+qC$4flf2t-#O6i^L^D-b$EI_19erdyC5lvXT=Vwxl#8k
zGA=%<M^aqR9+^p;pnDXL-;YtkF3@efpobtn>TCTjP4$zDisdfslkFnPFXtkoUGj>B
z`679-SmcsZQSKrVm;ZT6ku59E&MC|HcWK%J+w;^0X|wiCd3?1&nLYB$bMp$?M-~f<
z^LXaW9~!&p^3^zYPqwud4cYvU97Z*s<Y~%(mBWN4fB65E!{E;UGY%WQ``n^Qj*T3K
zYH7KvsnceTsPR7yhkxI1;?(7z_@$BBV|9EQN9)C4150(AI$At?b;;w|n>Cgi57qu(
zd+67j>+r0ZL$CnRI0UDTuwqemdA|I=y9Em$O+0E!^{v$1<fk?It{L-g*3_S_-90WT
zndiW2sJ-5fL7jKhfv9OgQ+_3@DJiLs{dHd*<i-<&SN7y$yJ@<s3A0Y@oO<Jzse?_2
zXT}RqTaXwP9hDf>gNw@iYC$k-Qf)zAc6qiuzx+>A{)i!-!G8+@X5m*0*fpB+k6fd^
z9os$)?5p}`7Fjo4@Sn3tW8o@sOh2_njVa6BtI^Kb#zGU{Y)|7FTK`+uFm235nLW4!
znOKr9D=%hFp-y6M4C-N0<jd3<D7>8?OG+%Agtzza<R9o)lvmIpJkUR|oqynvhr5Tc
zir=>X7Z`})dUpTnEnokj|Mcyqz~en{GEzOX>h15xHGj97<C_2L#$@zoeS6Qw>zgwE
zAMI$4YaS#t?VR&!tgrdbYD!|7{cq(r2EH-yje&0rd}H7n1K$|<#=tiQzA^BPfo}|a
zW8fPD|2hL#6Zdk(-Gp54IqkTYetd4j)Z26<aS!#=T19p)Gie#d(&_<&X_sLUIkrrp
zQ*LoIvR612@`7meL2nu>aV6(R4rJcSiURm1lr>hLI;ga%)$)%5zg|8#*}gTl>w03#
zMRz23c173cPEaN|VRbJj^!v~f)=n*9X43-8Lz?5_N=LMea)jiu1MYw4fVajD*fr4}
zHyrJeyTuOq{q5joWC!;hwn)9!6N7)t#G$?U2=Oh3@rY9R83`e`5MkT_5iDcHNLeLD
z$nRn-4wT?Rz65<IOW<20LDzE<WIU0;;EM!S7E;`GmBK1e3f~ARwBw}EFWAEAhV|zr
z{%+58Eq+TUTV9|k*J??}eg#zyo<alP4JVx|68d#xCduCPBC88g^kP&9mDjhW#l@{D
ztC<tI?6;*jxj8jEV?-mC=~CEX9-T>jFIau7Hw<!mLTny~7dN8O<+o@=PKt(3eI#;M
z5$pzZLCD7l+zg7qvi;#`hz$ofJq$jZL(!%#1h&sQLznA}L%TX*Ue8WAe!C<1!j5?N
zTQIIsUkpjj#o%W}czmN2yUqx)uD2N5ewLu_r4+s2$`E!|j<Ms)(XUkn<}a_nNEZbj
zO;Vu6YXy!%iBW@;_^?Ka`b$bgYgOWnM<ujkDiLIRn0q41;x@}VbMA-jIBDi1%Bk5&
zBdV9votLxd>G1C<t8^$Wt1Bk%RtkO8=}Pw!BPi&45T&bp$hD&zwLa%aH?6H{#d#A7
zwKbq#*R^SE+$TY+{X<Y^G#Hz>L6{vk5XXP4M5wtEFLf1|kx_w$x#c*qOb*?pGW6dk
z#nEO`4C*UEhmB&~G7}@eQiQrkLQKvSqV!@JY?I3H@O&u-W|yL4MmK!0?}Imc(=nxA
z4wg*F$I?%Qh~bss?BY^rX$x^aK!ic<#c;Hd;K?lsep)QWn<N=RU&wH?LXM6P<fw@)
z$MgB+2!2`)X@?4=S5{!)u1(zTk^x-J$1twzk2;+7uNTO^`$4Mexq+s{%%CxrqbMM9
z0NJ-Ip-N>c-9{X}i4LR4)OO^5)SL1axKjLm2l~0*lJ-90liyH%>NZrHCe8mSSkPJt
z?-O#IVLp55SP8;B3b84u07F9au>4~dt`%e;zDqjp<)-0Ov;K%KNQIMiU!-`Y;9dVd
z$Q#fb_eLh8c777No=${|O(LeH^+MR<1Z;ZH6)iiZV9@alm~GBQ)r>-%>s<nQ?=p0I
zBSh;mF>;qnP(DYB?~7%q^p@lPW#%{Ha_B!S$HcA`7_*=PE1y)LV}JtHLIox-Rba;j
z1zs8a$h{3vaKC=<%*8Cw;wnYAY4q?~I&pgmIqsfHbycIu&`?QnCv$05crvBkjiv^_
zP;xyPKvx!e(uh&c6xqw30*+eHuJgu(d_DSDtVPNN^@4{lhk)ll1ar3x!pFk{aiqmS
zbTS`+PWF{ZEm1=FlL9B-R6r}L0uz^&V~A-v&JLENHyg)KtRERqq<AWqB4102LyXrI
zUJ_XE7sEAMj40h$<W(l2Y*9ai^v^`mlpK6C&&QVVLfrmPgf7J;C|_HOgXhc8<CYLT
zkBTsNiWoH!5+t9MAR<MI1IMLU+ewDuvt;=7r3_x29Ptz62<%<Q^)625>YMv;21j(c
zRR&kd|L`s<YPo{WYfYn9ZlkH_u99M|<&o-CGOg(zL-LWK)Ym<bHfwuP$USEoIoqCk
z=vq=AQ$Fd<*P~AJwCLX3dci^CD*QZiFlgl<44ys^2BQXGL`@}JuPd?BLWy>13f$XT
zfmOB@nEYKic7K%P2IFLZGdZSCl_A<!hOvjFuqc*dshJeZSgz<;EWtiU2@Zc3i!<Gm
zkT$R%4xG-wWYZk<Sd@p~b{C+Nq6lx^7UM?uQg~IAAxtjBimoD@)D|ObsTkLTBrIo3
zVCN;pqaUQ;UPzG?BSXRr8M@t)VV-jxmp?y`>#^96E6{1iop8EEwcBb*d)P{HsGCNP
z>qnED?f{xREuWgL>rLN3i=j#Nq2x6qklN;X(Uf2p^1f$JUAtP+hIl?rxui$Y=d`GI
z<_Ce+6&2d-Qh~opg_sE{lnqp&E?Wg@4;2=6QbFXdf`Nq!pFUMV|8W(1UaZ22x+)yp
zR0Y#TjL(TxI8s#wUr`lqrB~tC_$ru1R3Y2H3h^Cdk-a$)m)%m)zc3wgaTZ=*%t7s)
zJd}<uzyy8~$}@`LR9S*?g{9aVT!#J6%5Y|i5JE2ztk;V$-&>3yW{Kg=lVD(m1naj-
za9&S}1*f)ge&NO3om?+2e7i39vgkV1?A%2Y;#W}H;nQf~@X@4fTS*HA`Q#DVn-ZtS
z(8QlZDY}0kJ!$7fseBh|z001Cds@-~e?Hmm)T1@qwP;eC4}#6_RT%$91;1Ae&sF&N
zM1_kFRS@4(!RodO_im{0;HnDNmsRL<L50<48Gol$P@YiX?oky|j;L_{pbFV_D%9^)
zVfJnny40%h@MjgqtcZnsO(I^}r^44M9g|yTp>{+LVrJ)IXi5Rf4izHyV-fc2lwiV>
z5(I27g`Ttwezroi`%#F5b|NIM72(I0VoVz$M*BNrYz~pYbi4%H*4J>?dKYkC!o0Y1
zj=G%W%{6-0e;2L%xSXsarqP+m(KPRvl2$#=qaQvb)2lu)beqM_I=?`Y8+g%U7CTSp
z*i*2fC4I8w)3C*Q)UZg4Mr(Z#Ec~p3?H3gi87dgQXPD2hnqdpW4u+o@HZaU)P%-2)
zbYlo)uxH>gJYwyfXV}NEkzqbVX;>V(UQWb#qg3!-r(w#gObn{Xf%&*RO!6*3n_-1G
zzp4mV*A^pxS_vNKmEu>YGT0p~L(6m_zT6XHK#mAH=R|lKDn{WPF;=}6V`G#Ax2<cq
z7WH{tO_nD&YeF;bcGxv~URO(FCoQMKms6>O(P%m^Q<D9{JnFP9nPMGd=tPfDiue>j
zIX66M<Q`|*Il!KNd}Kk7-x$-1(RySyMvE4ltrys_{uD5%xmd=a=3zI(VTKb7XBjRr
zTxGbyaFgK<!##%k3=bI|F+65?!tfizQ-)^@#S{ly<|kK<_Qj*EX=r^e6O*|d^b5<w
zfm8Wd)w&Q{l8W#?y%^<jCE%KuqW#5E#11WktAh}Gmk4pjRfJ<xMbP>z!s-k$zSM~E
z+DHQBuUomX*?C;NwI>&u(u|w)`&F7Wv6gJ|mecZIrqZd0qiAoilD3uQQR(<(8u~h#
zmidNKpX&j1t<IAo);P0gH+$N7&4Q|*7?ZzBj|L3WqMrNf1!>IosSN8G_Ane}Q1ft^
z;X1=DhPw=!{%iWK>GyMX?ghh3hF48^%}O=CoQp-Ge<FGX_eEuwG<1#0MERm@T%VVV
z2Vwc>FsA_4CkhdLqX<(@7Nhg>5{xV;MVeI^g4dTJKTL?zJB4@`D1ya25xRd7;YF$#
zw(G<geq{@1<eSTRm$&AAY0%-`4!ueb0&7WSvYZ|-n@SJPjiR_O3X1i~qbog=$@4%o
z^)?Kp-FpHkXq6{9{NPNR;_a!`J`4JE#+Y6e>Cp<I7Ij`&FBr*eSMyQFaEd|g6Ka3B
z&7kH%?E~uZ*KqLH{%iXFJ3IfT32#}c#ybYTqcJFK(F?U<DL5U`A2D18jxNnY{n{K{
z%E-g1WBC|iSO~dm5o#QZvG#p23U-&^H)$!xSe9Yok}_m@3Nd_v5dBO<Fq4V!<b(+3
z9_(G<{zmSoDw~^L<-xh9>2TXtU!g6VcG8NI%V>aPD$O}Gigt%6X@7AZ4e6ImzbuWW
z@s~o#enJ2_5AdXmInH#)&7R)RvY^vzjp<sn9*ymxMPWttf;r6gT?{7~)V`sP4Ru_o
zV?x7)#uwCnpyol%frf{_a-inn{a>&#I11_C_kic7WSrRE7n`2;$8ALh?#;=<n!+6X
z^fDLcqx0cYQGnAdw=57BA)PD6I`*79Xmtsua4e2al%lwM83ykz!(%TY0>%naeV4^l
ztE=45+M(R9d)jb$GxfQZtL~Ed$b;0(YAw}1n@RfP$I#|kgQ)Ue5fxgblA?Jx%32da
z3pNLlz}bhE?sFxpEJwP$(~64MnNUD;19E<?O+WSfBsihVo;g2sL9ttVbT#tFwjI8Z
zZ1YCh_11WQt`&}^x?;ovXPD<W;h3N$9E@Aw;%-Mg8RCF<S@sZ0?XY~AE!>{kAUf0r
zwd1U@`=%Az!wRQGz42;CKQ!#fgTuoTEK!Q!wNwJbb}|f2k|Q^u9A|En!>hCcv#wR(
zbw>sAhAFUQr-HSoz)Kq?RGpO&rz$a4p~ROTl=x+(5)W#XXcXx1{vs4>S|I<b8d1Mr
z|EnUWOuxpr@qeSpsZrih*0f|^)b#fMKazKSUEkMej+z=!7Y&(nEB;iF-1Kfd6{o48
zS>|sN_|5w&apUp&{=@l80dcYFcK*GZwqRZYxA((8j-5tSr78c_v3vQS8awZ=@<tx-
z|7h&KYB2uE*v;BFV}zP>JxBF$X#Fu(`Z?L6Jf^MGXgf6GOHG9)@nxA(Bf_XXdC~XZ
zK6p$__c~k@@W<e}tBqDG<1~S%8N9#m|Id|i!`rCa=JC2%G?p6q)4x?2r;af7NHqrQ
zo2E!9X|#aJ*)%_zO4qfzxBohL8lg?&;FT%=5Y99SV>QB>e=dwIQ%-8~umH!?SHF(s
z|5bsW`tDz4L<WuKtK~eM<C6QT@AP+pUYEuy{{exXrq4|`<UJeI^6@{j=szdW(^#Yt
z=&g;rww|9qh3or4#NAhga;KNsa*O@yX-WTcq_?1k7W}k`9A^AL(ktK5R)}fAgftqu
zqdVER1~sGy(E&?uGBI?aj+XY+Qs08i4j9pPUtL<(oJXhMy%jtax?*c9Z>&7zjk?`l
zkcE2Uz)%m^DO-Una>cAlXQXdziOy^aHfO#gCPX^m7j1j|e9;!S>TEFKxHZh)SfNF@
z72YnjL?3TUv_D{hZgLA8x9Nt98&gnkoQVxBa`Ccl9v(O3<JgZyxYw){24~CA)mwy+
zKrzOdN^tU+1P!H97(A0=ONtCp>ttxjrX&^Z<WT0zvFryqx^0zX&N(?|e)obN2KaM1
zY5v^uOnZ($OOG>1`kgK}T&E6~kI<3w8gi>#PJ#>5Y0{i(dOKf9J1=C>lg-^J`Ev*j
zQ~A<4eJ5I3YDv?E7|;PD9*wMdDd<1so<Q1?!<oz^6pl*9+!?(fj!nRvfeAP{IRVFf
z5->#26LT-cV^DoJ?AsWJJr*(W@{Png69Jwb3P){82>gsYqIPc(#x4)QieqiC0KWJ^
z*BgoLOz}Y_z{z!qF!tz+v1gJnB0d=n>r#QZOiUV>3+wO#NIZ*>SyPNEyHd1>EyJof
zAzrzQ5PL_2p<~79=)|U`7fIl1BZb>wDW+VM!vDhrZs&K^+>PnUT<f8J+^rj?+)2rE
zGP!Y>7KW~-=+Bc#U{giKce1E@Mpv>}*qM%R_oL@`+~{LMbMi~Kq1tM5(v$Nk`@8|U
z4AG??QQEY)&1XSKgD*a&1i{IvEzCZ&L7S2`xaHuBw?>`_z2b_`mMzdD$rf|`EMUYl
z!KVg8r28A7SgDUA%k^=4hyi?d8DdeOF$S#Tqm7*@&h|A!V39eRy9U5{XiqqvNW;FW
zbfm>(VDZKrv}4mqYf_6bq`m~6Y?`FFix8viL|D65gru%wwA(JmJSz!|*%XQ8TnR*{
zBxuth!PQpGf1=oQNv;%~Joa$L;oozDqMn=|&xJE{X`q%aCu!BfwRAgvI(@%;B(+VD
z(*(P0x^*a#mN$zbhi73__%?`+9`dD`f$sFkp(WXs+fvsMa|$jpqTk*&qec6^2%4X?
z!R$G9Sk%@Q?b_ME<){S~rkP;3r6HWP^>DPH8NBS9!EsMBoO-JZu0S8d&lq5yjS<?m
zFh=!FW30>P<LG8Ss?YP`dx4Lhd-w>O$A^&~;3?~i(muJUzLSli*YZ&Cq!5Nrijmf;
z6#Wg#5WcMpS3-s8yjX~HuZ4KqL4<>4B0Qcgg5N$79x-3>*Ala-Ix#x17;A_Wqq>h6
z1}C<2X7&R)3so@JOV5-$+w}pJb*!U1X{(6}XHdb_YFadCAay7#B(KaAI<hQ|%zq9i
zXF(8+c<)Vp*16K7P)ExD&Wc79nUL`%eRAKSO?uj&1aXoF5yKCDh}dGT4Tr(a5L%&&
z+Xi}=yipGuhv}ms-vBXP4dLfzgtEIvupe)XB@TQXo5;t*mwYVbOdzi|LE(NAJb%mJ
zZi?$YOi^253LTF)tWC*4=;j=3_>ha}?F9&%ScGXCi{bKJ2{w3?Lb-&!`s$UzrB4~2
z&M3pH6U=WK%HZfBL}H8(#`!|%tAw~WO$eXWLd3CG<3-2VE6$!$E>tUs3ye17qBcCE
zb+3-msCio`!hI2W@1I0ZJC3C7@(Rjdl1puqlc|q>45>0gNoWyBD;{`~c#AU~NV2Ey
zD=lc%EMw|tsYihgJi0mMy}-w5U4&cw`G_GMev8oV^f5vv(#HB(x~RUQ4>O@5UXC%s
z{zPNEykQKBE__Ib@$u_wKDO-R<N9$voKEqPbBvE0b$kTw;A7ebKHe?kquE?OE|nx<
z)9`#;7+QehF(nx7B!r8Z2(M1D*Zd)37+Fh@w?cv*!BW)yBt>~E87v0M@NlOLXI{zh
z&Q*?w@p3G#l;ikfIocePW8Lp^ENM}Wg@Uf!m8z~>Mo&L(*-C4!e{*flVb)!G-R>~0
zaM?hEm(C$BWGoG{SCVCQ1_jx6rPi%G(-mJ|I<eB3LfY8UtiGnSnEA)kLTx%d<D+2F
z($@lq6~75W&U<6?>-HEjsRMRL1mkdIFbZzA$BD>xXnn;Wrl)+6mC_m`R<b3*4=oT|
z#g_Rptr4`(99la~Frtq!tR@;ll4^hhi}i6cSs(3F^})^7NB@L=n310g-u4XS7ZxFh
zCxQ=~f<I&-$7eR}Kk-xroN5%P5h^i_En#I2ti;vpl}HX4fZgQ-5V3Xu-d`DjCk+Fz
z*=8WNyA4EI>w$=4rSis9ZupjA+^RX<xb@`8{W!pYOPzm{a{afH`Tg0{O?w1Ad!A3N
zEaFKP9ZFWa{OR3m5BjpvnI0K8r=(;%T5WDkKRz?3h#e-hFV%>4?9ihZH??U!7mY2^
z(dbx1n4i!EqmG9n;c*AtU*V7IyIP~*h7+o4?J!Ko5@q(LxNpRV-%MlVCGfGj)&$K5
znZY950!vR>BFN4f##%P0uC{^idK=g|q(P<IAMdr((Y$jZ#<wa%UlTEIc9-JuuQD{V
zEXPRg3as5<f$7-_OufU<Lx~&9m2iBaL}Hste9Wpu!}v;E*;I)=XDhMlMI{oA2EcgV
z9*#RbfeWum=1MQSa7S0XrJ0-dQ~L)CNPj^!-JUF`JsY#gFC&?_x!vgb;7F3)2&XIF
z9ch$D0DU>*MZ1$+$?t*#rH`>B)l_3Lx6mbPm(PO3)8f#ko-l4_IQq>9hGS|QgwJb*
zYjKWjIl&6)1!nN=&;Co|&EOtl0VG>tr?U<Eq}w9%9m6YIY|OUBrYsxmG_^*H@s@CZ
zXAYxMGaNjVf$(?z!DVFQ@sMKF6bf<NM1nr`Qg~*|F-2O=7UbA+(d`P<q$`lNTY-5d
zN(_xxV&Qj63|z0oiStUFd8<UaRVDQN*s`Rc67flu2$7xVTIx>V7T${G#?@JHqeR!~
z%Yx0c?C>l)GG{Ewy@t`j1?4p3ToDb+NT+j~d(q(QQFPxcg!cXJPd`m|r@Nh+)0NMb
z6m-j&st@VXuCY9lYP}acQ~>?7g7HS-j|(eXp=&&w2AOV#H#>}QsooH84jE%gi76f(
zG6&bs5-A@nA*r!M@Gc81k28nQEK|5W<l|<6F+O%M!h$?Qc)c>ft<wg0xHTKIx2D5=
zLm@s%g>Y3$uuvt#h$nJveO!(@wybx@LxF&C3aq)Lz$i9Fk{_)EU!+9n1SO8HP{M78
z619hwkep_?poG_DCEi_8;&|E>ZuyFF+~kX0IPWJ$9IZT0hD$e6ukg7vv)@E|o-m5S
zy9}X;l_J``Et97D^rAp;biOc{2E6p4o*!N4>>zsz?QTJRI~db%^}4huoh|8Pd=TV+
z3B!pef$*vF!r2im@Zz}zT(26VKG_I(WQ<{BOyPaf9J^dBk)3Uch%FYlm~IZg^`_YJ
zhL4sKV;qe(LedyRj0iMDtdAkiq!?nrd_#;&%0rz`Hd5D?;CvK|89f;u?3JU_yK+4I
zg~d#k0tu%S_}N{Ft+`5En52Z^W+kTnqQu;*EPfs-k^F{DC9zy{#-tMO9V#)6Et`)B
zsKm3QZ@JmJ6S+_PKyL1$H?-6C5M5rqmG&yuP@m6>Y4Z8$^rEhsmPJ-jvQ-+zZjGVV
zrJZPwv<=~V54t$enJhXv(i9gP>hs8qj#L}bK6_o7I_ityp)42PefvO>6OEIjJ7WEw
z0E}G|h}hfhaeq%R!c&8>|3VPDTxpBKC|~R~^n~M$R`{ahhEdV3=<n|ecY!N=7IlSU
zs|&8(aYn>`C*(+5;_a|hJQF4%j7=}rrWeBUUMaTKiqL{BRV}lZVS$?*uV2X#H?JJ0
zeJk)}Nd@K^D!^l0tzM<TYv#mPoRxUplSv8&Dp9jgi5dHqxOqp3B{|Qyx&Fhr-xFOq
z>5c}vH}Nioh|W-S!w%vtUqZvHz9+p!m6RHkO(8uJX;J@JTGF8ljepyb+UT~WRjWN{
z^!OHZnX{%_9wrpkS)YPqwaMA=lc3~gB6yB~Y^*Pij_|<XPhNPb<%dN*eId8;M3)8^
z1VlH-wVigb*lUZ`xptVDVvpxr?D5Rg9=%W4Vvw#4Y{RUe$Da54Tr-C;dxouvGzWGj
zW5>)iq<R*j#V=*>VbT&Uw&eT5T#kt!<QTu99NXE_bH{xZ=;NtC4}}6zn-#EOx%RT5
z5@USW6eEkJ=qx2Zvgt<II3;$?Q)2I0B^2ZCbBmr<aMl-{xV7Fo+)JHj^l;;O+C6t0
zB@5<}m0mSHy)36I({pHyM+$99kEg|s9Ca}1L?@s7QNji{+WXLfe*MXk4qo8XZx{7x
zXLlV+x!fQ)?2&+BokGwm%?qj>tzmQB2hN9mvA&ZRI$dzZ6KQj(s_jrc-xe7~b|{aw
zM=OCn9K7t1m1u)b30BBoYJp!So8xC|bIg2chR=7+&|#+;ruR$4=zAHMHM0n+6KqO(
zwFKkSWGJ5`#~XP$X1G_N%cctUT%y37K??L{Ikx(D1-3EYxWFlKEQjTkAxg}fro{f$
zN@TOC*SzCOY-IVwDdaMDXJZ!EV!ktX_>&3e{@q8a_d84L$FCue)X7wAID+o%sw9u?
zW#r#1lhTDfDNDBtHSBIjS*c#sxxj^7UptVEt2Ny%H>EQd3}{(_4%sdFB$(Vc4c!c4
zpfCu-%DYTvBZ$D{b)7NyW+3#LuiQN4g_(W5(WRXa2Gw}u?lw;hedLZgp>Eig?84p|
zozZcp6HdViF5kC=&Wsiau56C<R7cD#iG-|OU(_|v#k+)JTzFoF8YVAEzA8a*a~ZB$
z%b~g;NB6RF6y7PvqsR(;o?d~4*DFx(sDLa%fr%;wTOL%vndP0TR|*_qau_XNB~m<&
zaPLY>xVf`Ca*J)OIG37-bW*a13~40=InJWos`1nzZz%OzTtf3b`%?0xI4XJ)PN%;F
z5&xqP#sA<&Bab*zah(<2A8kTEe$*$cH`=s+?kB+qCV_ahGz|07!f<P6XIM-PLic!o
zZ13WS^e8{PF!Y0XqYup2c;T1Bt>HDq1A-cN6uG$LrnnV0PjrLFR993la6y+N&Iq@5
zMs~RqR$gg|y+flg_NP>6_s)S?WFef_l_1|(2+<Z1+MXAqVwnW{Bc&`yNO90l2JKNY
nq@0nVqm>-n<K*yV&#&Jtl4I^6_5YL{0Y>b3#-kjI!^-hr*$}7L

literal 0
HcmV?d00001

diff --git a/src/LIB/RAD/ecrad-1.4.0_mnh/ifs/radiation_setup.F90 b/src/LIB/RAD/ecrad-1.4.0_mnh/ifs/radiation_setup.F90
index 0677d3f94..69d2d44fc 100644
--- a/src/LIB/RAD/ecrad-1.4.0_mnh/ifs/radiation_setup.F90
+++ b/src/LIB/RAD/ecrad-1.4.0_mnh/ifs/radiation_setup.F90
@@ -40,7 +40,7 @@ MODULE RADIATION_SETUP
   USE radiation_config, ONLY : config_type, &
        &                       ISolverMcICA, ISolverSpartacus, &
        &                       ILiquidModelSlingo, ILiquidModelSOCRATES, &
-       &                       IIceModelFu, IIceModelBaran, &
+       &                       IIceModelFu, IIceModelBaran, IIceModelShapes, &
        &                       IOverlapExponential
   USE MODD_PARAM_ECRAD_n , ONLY : rad_config
   
@@ -161,6 +161,8 @@ CONTAINS
       rad_config%i_ice_model = IIceModelFu
     ELSEIF (NICEOPT == 4) THEN
       rad_config%i_ice_model = IIceModelBaran
+    ELSEIF (NICEOPT == 7) THEN
+      rad_config%i_ice_model = IIceModelShapes
     ELSE
       WRITE(NULERR,'(a,i0)') 'Unavailable ice optics model in modular radiation scheme: NICEOPT=', &
            &  NICEOPT
diff --git a/src/LIB/RAD/ecrad-1.4.0_mnh/radiation/radiation_cloud_optics.F90 b/src/LIB/RAD/ecrad-1.4.0_mnh/radiation/radiation_cloud_optics.F90
index 0fa20eb82..f636fd604 100644
--- a/src/LIB/RAD/ecrad-1.4.0_mnh/radiation/radiation_cloud_optics.F90
+++ b/src/LIB/RAD/ecrad-1.4.0_mnh/radiation/radiation_cloud_optics.F90
@@ -18,6 +18,7 @@
 module radiation_cloud_optics
 
   implicit none
+
   public
 
 contains
@@ -37,7 +38,7 @@ contains
     use radiation_io,     only : nulerr, radiation_abort
     use radiation_config, only : config_type, IIceModelFu, IIceModelBaran, &
          &                       IIceModelBaran2016, IIceModelBaran2017, &
-         &                       IIceModelYi, &
+         &                       IIceModelYi,IIceModelShapes, &
          &                       ILiquidModelSOCRATES, ILiquidModelSlingo
     use radiation_cloud_optics_data, only  : cloud_optics_type
     use radiation_ice_optics_fu, only    : NIceOpticsCoeffsFuSW, &
@@ -48,6 +49,8 @@ contains
          &                                 NIceOpticsGeneralCoeffsBaran2017
     use radiation_ice_optics_yi, only    : NIceOpticsCoeffsYiSW, &
          &                                 NIceOpticsCoeffsYiLW
+    use radiation_ice_optics_shapes, only    : NIceOpticsCoeffsShapesSW, &
+         &                                 NIceOpticsCoeffsShapesLW
     use radiation_liquid_optics_socrates, only : NLiqOpticsCoeffsSOCRATES
     use radiation_liquid_optics_slingo, only : NLiqOpticsCoeffsSlingoSW, &
          &                                     NLiqOpticsCoeffsLindnerLiLW
@@ -188,6 +191,23 @@ contains
              &  ') does not match number expected (', NIceOpticsCoeffsYiSW,')'
         call radiation_abort()
       end if
+    else if (config%i_ice_model == IIceModelShapes) then
+      if (size(config%cloud_optics%ice_coeff_lw, 2) &
+           &  /= NIceOpticsCoeffsShapesLW) then
+        write(nulerr,'(a,i0,a,i0,a,i0,a)') &
+             &  '*** Error: number of LW ice-particle optical coefficients (', &
+             &  size(config%cloud_optics%ice_coeff_lw, 2), &
+             &  ') does not match number expected (', NIceOpticsCoeffsShapesLW,')'
+        call radiation_abort()
+      end if
+      if (size(config%cloud_optics%ice_coeff_sw, 2) &
+           &  /= NIceOpticsCoeffsShapesSW) then
+        write(nulerr,'(a,i0,a,i0,a,i0,a)') &
+             &  '*** Error: number of SW ice-particle optical coefficients (', &
+             &  size(config%cloud_optics%ice_coeff_sw, 2), &
+             &  ') does not match number expected (', NIceOpticsCoeffsShapesSW,')'
+        call radiation_abort()
+      end if
     end if
 
     if (lhook) call dr_hook('radiation_cloud_optics:setup_cloud_optics',1,hook_handle)
@@ -208,7 +228,7 @@ contains
     use radiation_io,     only : nulout, nulerr, radiation_abort
     use radiation_config, only : config_type, IIceModelFu, IIceModelBaran, &
          &                       IIceModelBaran2016, IIceModelBaran2017, &
-         &                       IIceModelYi, &
+         &                       IIceModelYi,IIceModelShapes, &
          &                       ILiquidModelSOCRATES, ILiquidModelSlingo
     use radiation_thermodynamics, only    : thermodynamics_type
     use radiation_cloud, only             : cloud_type
@@ -221,6 +241,8 @@ contains
     use radiation_ice_optics_baran2017, only  : calc_ice_optics_baran2017
     use radiation_ice_optics_yi, only     : calc_ice_optics_yi_sw, &
          &                                  calc_ice_optics_yi_lw
+    use radiation_ice_optics_shapes, only     : calc_ice_optics_shapes_sw, &
+         &                                  calc_ice_optics_shapes_lw
     use radiation_liquid_optics_socrates, only:calc_liq_optics_socrates
     use radiation_liquid_optics_slingo, only:calc_liq_optics_slingo, &
          &                                   calc_liq_optics_lindner_li
@@ -344,10 +366,12 @@ contains
               call radiation_abort()
             end if
 
+            ! Delta-Eddington scaling in the shortwave only
             if (.not. config%do_sw_delta_scaling_with_gases) then
-              ! Delta-Eddington scaling in the shortwave only
               call delta_eddington_scat_od(od_sw_liq, scat_od_sw_liq, g_sw_liq)
             end if
+            !call delta_eddington_scat_od(od_lw_liq, scat_od_lw_liq, g_lw_liq)
+
           else
             ! Liquid not present: set properties to zero
             od_lw_liq = 0.0_jprb
@@ -427,6 +451,17 @@ contains
                    &  od_lw_ice, scat_od_lw_ice, g_lw_ice)
               ! Compute shortwave properties
               call calc_ice_optics_yi_sw(config%n_bands_sw, &
+                   &  config%cloud_optics%ice_coeff_sw, &
+                   &  iwp_in_cloud, cloud%re_ice(jcol,jlev), &
+                   &  od_sw_ice, scat_od_sw_ice, g_sw_ice)        
+            else if (config%i_ice_model == IIceModelShapes) then
+              ! Compute longwave properties
+              call calc_ice_optics_shapes_lw(config%n_bands_lw, &
+                   &  config%cloud_optics%ice_coeff_lw, &
+                   &  iwp_in_cloud, cloud%re_ice(jcol,jlev), &
+                   &  od_lw_ice, scat_od_lw_ice, g_lw_ice)
+              ! Compute shortwave properties
+              call calc_ice_optics_shapes_sw(config%n_bands_sw, &
                    &  config%cloud_optics%ice_coeff_sw, &
                    &  iwp_in_cloud, cloud%re_ice(jcol,jlev), &
                    &  od_sw_ice, scat_od_sw_ice, g_sw_ice)
@@ -436,14 +471,14 @@ contains
               call radiation_abort()
             end if
 
+            ! Delta-Eddington scaling in both longwave and shortwave
+            ! (assume that particles are larger than wavelength even
+            ! in longwave)
             if (.not. config%do_sw_delta_scaling_with_gases) then
-              ! Delta-Eddington scaling in both longwave and shortwave
-              ! (assume that particles are larger than wavelength even
-              ! in longwave)
               call delta_eddington_scat_od(od_sw_ice, scat_od_sw_ice, g_sw_ice)
             end if
-
             call delta_eddington_scat_od(od_lw_ice, scat_od_lw_ice, g_lw_ice)
+
           else
             ! Ice not present: set properties to zero
             od_lw_ice = 0.0_jprb
diff --git a/src/LIB/RAD/ecrad-1.4.0_mnh/radiation/radiation_config.F90 b/src/LIB/RAD/ecrad-1.4.0_mnh/radiation/radiation_config.F90
index 2a77f3a3e..0f641bb40 100644
--- a/src/LIB/RAD/ecrad-1.4.0_mnh/radiation/radiation_config.F90
+++ b/src/LIB/RAD/ecrad-1.4.0_mnh/radiation/radiation_config.F90
@@ -32,7 +32,7 @@
 !
 
 module radiation_config
-
+  USE MODD_PARAM_LIMA,          ONLY : CPRISTINE_ICE_LIMA
   use parkind1,                      only : jprb
 
   use radiation_cloud_optics_data,   only : cloud_optics_type
@@ -69,7 +69,6 @@ module radiation_config
        & IEntrapmentExplicit, & ! Estimate horiz migration dist, account for fractal clouds
        & IEntrapmentExplicitNonFractal, & ! As above but ignore fractal nature of clouds
        & IEntrapmentMaximum ! Complete horizontal homogenization within regions (old SPARTACUS assumption)
-
   end enum
   
   ! Names available in the radiation namelist for variable
@@ -110,14 +109,15 @@ module radiation_config
   enum, bind(c) 
      enumerator IIceModelMonochromatic, IIceModelFu, &
           &  IIceModelBaran, IIceModelBaran2016, IIceModelBaran2017,   &
-          &  IIceModelYi
+          &  IIceModelYi, IIceModelShapes
   end enum
-  character(len=*), parameter :: IceModelName(0:5) = (/ 'Monochromatic', &
+  character(len=*), parameter :: IceModelName(0:6) = (/ 'Monochromatic', &
        &                                                'Fu-IFS       ', &
        &                                                'Baran        ', &
        &                                                'Baran2016    ', &
        &                                                'Baran2017    ', &
-       &                                                'Yi           ' /)
+       &                                                'Yi           ', &
+       &                                                'Shapes       ' /)
 
   ! Cloud PDF distribution shapes
   enum, bind(c)
@@ -238,7 +238,7 @@ module radiation_config
     ! Codes describing the gas and cloud scattering models to use, the
     ! latter of which is currently not used
     integer :: i_gas_model = IGasModelIFSRRTMG
-    !     integer :: i_cloud_model
+    !     integer :: i_cloud_model    
 
     ! Optics if i_gas_model==IGasModelMonochromatic.
     ! The wavelength to use for the Planck function in metres. If this
@@ -274,10 +274,10 @@ module radiation_config
 
     ! User-defined monotonically increasing wavelength bounds (m)
     ! between input surface albedo/emissivity intervals. Implicitly
-    ! the first interval starts at zero and the last ends at infinity.
+    ! the first interval starts at zero and the last ends at
+    ! infinity. 
     real(jprb) :: sw_albedo_wavelength_bound(NMaxAlbedoIntervals-1) = -1.0_jprb
     real(jprb) :: lw_emiss_wavelength_bound( NMaxAlbedoIntervals-1)  = -1.0_jprb
-
     ! The index to the surface albedo/emissivity intervals for each of
     ! the wavelength bounds specified in sw_albedo_wavelength_bound
     ! and lw_emiss_wavelength_bound
@@ -419,8 +419,8 @@ module radiation_config
     ! then that will be used instead. If the user assigns one and it
     ! doesn't start with a '/' character then it will be prepended by
     ! the contents of directory_name.
-    character(len=511) :: ice_optics_override_file_name = ''
-    character(len=511) :: liq_optics_override_file_name = ''
+    character(len=511) :: ice_optics_override_file_name     = ''
+    character(len=511) :: liq_optics_override_file_name     = ''
     character(len=511) :: aerosol_optics_override_file_name = ''
 
     ! Optionally override the look-up table file for the cloud-water
@@ -428,8 +428,7 @@ module radiation_config
     character(len=511) :: cloud_pdf_override_file_name = ''
 
     ! Has "consolidate" been called?  
-    logical :: is_consolidated = .false.
-
+    logical :: is_consolidated = .false.    
     ! COMPUTED PARAMETERS
     ! Users of this library should not edit these parameters directly;
     ! they are set by the "consolidate" routine
@@ -441,6 +440,7 @@ module radiation_config
     real(jprb), allocatable, dimension(:) :: wavenumber1_lw
     real(jprb), allocatable, dimension(:) :: wavenumber2_lw
 
+
     ! If the nearest surface albedo/emissivity interval is to be used
     ! for each SW/LW band then the following arrays will be allocated
     ! to the length of the number of bands and contain the index to
@@ -502,7 +502,7 @@ module radiation_config
     character(len=511) :: ice_optics_file_name, &
          &                liq_optics_file_name, &
          &                aerosol_optics_file_name
-    
+
     ! McICA PDF look-up table file name
     character(len=511) :: cloud_pdf_file_name
 
@@ -601,7 +601,7 @@ contains
     integer :: i_aerosol_type_map(NMaxAerosolTypes) ! More than 256 is an error
 
     logical :: do_nearest_spectral_sw_albedo = .true.
-    logical :: do_nearest_spectral_lw_emiss  = .true.
+    logical :: do_nearest_spectral_lw_emiss = .true.
     real(jprb) :: sw_albedo_wavelength_bound(NMaxAlbedoIntervals-1)
     real(jprb) :: lw_emiss_wavelength_bound( NMaxAlbedoIntervals-1)
     integer :: i_sw_albedo_index(NMaxAlbedoIntervals)
@@ -1006,6 +1006,9 @@ contains
     else if (this%i_ice_model == IIceModelYi) then
       this%ice_optics_file_name &
            &   = trim(this%directory_name) // "/yi_ice_scattering_rrtm.nc"
+    else if (this%i_ice_model == IIceModelShapes) then
+      this%ice_optics_file_name &
+           &   = trim(this%directory_name) // "/" // CPRISTINE_ICE_LIMA // "_ice_scattering_rrtm.nc"
     end if
 
     ! Set cloud-water PDF look-up table file name
@@ -1326,7 +1329,7 @@ contains
 
     if (this%n_bands_sw <= 0) then
       write(nulerr,'(a)') '*** Error: get_sw_weights called before number of shortwave bands set'
-      call radiation_abort()      
+      call radiation_abort()
     end if
 
     ! Convert wavelength range (m) to wavenumber (cm-1)
@@ -1976,5 +1979,5 @@ contains
     
   end function out_of_bounds_3d
 
-
+    
 end module radiation_config
diff --git a/src/LIB/RAD/ecrad-1.4.0_mnh/radiation/radiation_ice_optics_shapes.F90 b/src/LIB/RAD/ecrad-1.4.0_mnh/radiation/radiation_ice_optics_shapes.F90
new file mode 100644
index 000000000..011f2666e
--- /dev/null
+++ b/src/LIB/RAD/ecrad-1.4.0_mnh/radiation/radiation_ice_optics_shapes.F90
@@ -0,0 +1,142 @@
+! radiation_ice_optics_shapes.F90 - 
+!
+! (C) Copyright 2017- ECMWF.
+!
+! This software is licensed under the terms of the Apache Licence Version 2.0
+! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
+!
+! In applying this licence, ECMWF does not waive the privileges and immunities
+! granted to it by virtue of its status as an intergovernmental organisation
+! nor does it submit to any jurisdiction.
+!
+! Authors: Marie Taufour
+! Email:   marie.taufour@aero.obs-mip.fr
+!
+! The reference for this ice optics parameterization is ....
+
+module radiation_ice_optics_shapes
+
+  implicit none
+  public
+
+  ! The number of ice coefficients depends on the parameterization
+  integer, parameter :: NIceOpticsCoeffsShapesSW  = 69
+  integer, parameter :: NIceOpticsCoeffsShapesLW  = 69
+
+  integer, parameter :: NSingleCoeffs = 23
+
+contains
+
+  !---------------------------------------------------------------------
+  ! Compute shortwave ice-particle scattering properties using Yi et
+  ! al. (2013) parameterization
+  subroutine calc_ice_optics_shapes_sw(nb, coeff, ice_wp, &
+       &  re, od, scat_od, g)
+
+    use parkind1, only : jprb, jpim
+    !use yomhook,  only : lhook, dr_hook
+
+    ! Number of bands
+    integer, intent(in)  :: nb
+    ! Coefficients read from a data file
+    real(jprb), intent(in) :: coeff(:,:)
+    ! Ice water path (kg m-2)
+    real(jprb), intent(in) :: ice_wp
+    ! Effective radius (m)
+    real(jprb), intent(in) :: re
+    ! Total optical depth, scattering optical depth and asymmetry factor
+    real(jprb), intent(out) :: od(nb), scat_od(nb), g(nb)
+
+    ! Yi's effective diameter (microns)
+    real(jprb) :: de_um
+    ! Ice water path in g m-2
+    real (jprb) :: iwp_gm_2
+    ! LUT temp variables
+    real(jprb) :: wts_1, wts_2
+    integer(jpim) :: lu_idx
+    real(kind=jprb), parameter    :: lu_scale  = 0.2_jprb
+    real(kind=jprb), parameter    :: lu_offset = 1.0_jprb
+    !real(jprb)  :: hook_handle
+
+    !if (lhook) call dr_hook('radiation_ice_optics:calc_ice_optics_yi_sw',0,hook_handle)
+
+    ! Convert to effective diameter using the relationship in the IFS
+    !de_um     = re * (1.0e6_jprb / 0.64952_jprb)
+    de_um     = re * 2.0e6_jprb
+
+    ! limit de_um to validity of LUT
+    de_um = max(de_um,10.0_jprb)
+    de_um = min(de_um,119.99_jprb) !avoid greater than or equal to 120 um
+
+    iwp_gm_2  = ice_wp * 1000.0_jprb
+
+    lu_idx = floor(de_um * lu_scale - lu_offset)
+    wts_2  = (de_um * lu_scale - lu_offset) - lu_idx
+    wts_1  = 1.0_jprb - wts_2
+    od     = 0.001_jprb * iwp_gm_2 * & 
+             & ( wts_1 * coeff(1:nb,lu_idx) + wts_2 * coeff(1:nb,lu_idx+1) )
+    scat_od = od * & 
+             & ( wts_1 * coeff(1:nb,lu_idx+NSingleCoeffs) + wts_2 * coeff(1:nb,lu_idx+NSingleCoeffs+1) )
+    g = wts_1 * coeff(1:nb,lu_idx+2*NSingleCoeffs) + wts_2 * coeff(1:nb,lu_idx+2*NSingleCoeffs+1)
+
+    !if (lhook) call dr_hook('radiation_ice_optics:calc_ice_optics_yi_sw',1,hook_handle)
+
+  end subroutine calc_ice_optics_shapes_sw
+
+
+  !---------------------------------------------------------------------
+  ! Compute longwave ice-particle scattering properties using .....
+  subroutine calc_ice_optics_shapes_lw(nb, coeff, ice_wp, &
+       &  re, od, scat_od, g)
+
+    use parkind1, only : jprb, jpim
+    !use yomhook,  only : lhook, dr_hook
+
+    ! Number of bands
+    integer, intent(in)  :: nb
+    ! Coefficients read from a data file
+    real(jprb), intent(in) :: coeff(:,:)
+    ! Ice water path (kg m-2)
+    real(jprb), intent(in) :: ice_wp
+    ! Effective radius (m)
+    real(jprb), intent(in) :: re
+    ! Total optical depth, scattering optical depth and asymmetry factor
+    real(jprb), intent(out) :: od(nb), scat_od(nb), g(nb)
+
+    ! Yi's effective diameter (microns)
+    real(jprb) :: de_um
+    ! Ice water path in g m-2
+    real (jprb) :: iwp_gm_2
+    ! LUT temp variables
+    real(jprb) :: wts_1, wts_2
+    integer(jpim) :: lu_idx
+    real(kind=jprb), parameter    :: lu_scale  = 0.2_jprb
+    real(kind=jprb), parameter    :: lu_offset = 1.0_jprb
+    !real(jprb)  :: hook_handle
+
+    !if (lhook) call dr_hook('radiation_ice_optics:calc_ice_optics_yi_sw',0,hook_handle)
+
+    ! Convert to effective diameter using the relationship in the IFS
+    !de_um     = re * (1.0e6_jprb / 0.64952_jprb)
+    de_um     = re * 2.0e6_jprb
+
+    ! limit de_um to validity of LUT
+    de_um = max(de_um,10.0_jprb)
+    de_um = min(de_um,119.99_jprb) !avoid greater than or equal to 120 um
+
+    iwp_gm_2  = ice_wp * 1000.0_jprb
+
+    lu_idx = floor(de_um * lu_scale - lu_offset)
+    wts_2  = (de_um * lu_scale - lu_offset) - lu_idx
+    wts_1  = 1.0_jprb - wts_2
+    od     = 0.001_jprb * iwp_gm_2 * & 
+             & ( wts_1 * coeff(1:nb,lu_idx) + wts_2 * coeff(1:nb,lu_idx+1) )
+    scat_od = od * & 
+             & ( wts_1 * coeff(1:nb,lu_idx+NSingleCoeffs) + wts_2 * coeff(1:nb,lu_idx+NSingleCoeffs+1) )
+    g = wts_1 * coeff(1:nb,lu_idx+2*NSingleCoeffs) + wts_2 * coeff(1:nb,lu_idx+2*NSingleCoeffs+1)
+
+     !if (lhook) call dr_hook('radiation_ice_optics:calc_ice_optics_yi_lw',1,hook_handle)
+
+  end subroutine calc_ice_optics_shapes_lw
+
+end module radiation_ice_optics_shapes
-- 
GitLab