Changeset 5007
- Timestamp:
- Jul 4, 2024, 9:29:24 AM (6 months ago)
- Location:
- LMDZ6/trunk
- Files:
-
- 11 edited
-
DefLists/field_def_lmdz.xml (modified) (3 diffs)
-
libf/phylmd/lmdz_lscp.F90 (modified) (23 diffs)
-
libf/phylmd/lmdz_lscp_ini.F90 (modified) (14 diffs)
-
libf/phylmd/lmdz_lscp_tools.F90 (modified) (4 diffs)
-
libf/phylmd/phys_local_var_mod.F90 (modified) (3 diffs)
-
libf/phylmd/phys_output_ctrlout_mod.F90 (modified) (1 diff)
-
libf/phylmd/phys_output_write_mod.F90 (modified) (3 diffs)
-
libf/phylmd/physiq_mod.F90 (modified) (5 diffs)
-
libf/phylmdiso/phys_local_var_mod.F90 (modified) (3 diffs)
-
libf/phylmdiso/phys_output_ctrlout_mod.F90 (modified) (1 diff)
-
libf/phylmdiso/physiq_mod.F90 (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/DefLists/field_def_lmdz.xml
r4999 r5007 634 634 <field id="tke" long_name="TKE" unit="m2/s2" /> 635 635 <field id="tke_dissip" long_name="TKE DISSIPATION" unit="m2/s3" /> 636 <field id="tke_buoy" long_name="TKE Buoyancy term" unit="m2/s3" /> 637 <field id="tke_shear" long_name="TKE Shear term" unit="m2/s3" /> 636 638 <field id="tke_ter" long_name="Max Turb. Kinetic Energy ter" unit="m2/s2" /> 637 639 <field id="tke_lic" long_name="Max Turb. Kinetic Energy lic" unit="m2/s2" /> … … 642 644 <field id="tke_max_oce" long_name="Max Turb. Kinetic Energy oce" unit="m2/s2" operation="maximum"/> 643 645 <field id="tke_max_sic" long_name="Max Turb. Kinetic Energy sic" unit="m2/s2" operation="maximum"/> 644 <field id="l_mix_ter" long_name="PBL mixing length ter" unit="m" /> 646 647 <field id="qrain_lsc" long_name="LS specific rain content" unit="kg/kg" /> 648 <field id="qsnow_lsc" long_name="LS specific snow content" unit="kg/kg" /> 649 <field id="dqreva" long_name="LS rain tendency due to evaporation" unit="kg/kg/s" /> 650 <field id="dqssub" long_name="LS snow tendency due to sublimation" unit="kg/kg/s" /> 651 <field id="dqrauto" long_name="LS rain tendency due to autoconversion" unit="kg/kg/s" /> 652 <field id="dqrcol" long_name="LS rain tendency due to collection" unit="kg/kg/s" /> 653 <field id="dqrmelt" long_name="LS rain tendency due to melting" unit="kg/kg/s" /> 654 <field id="dqrfreez" long_name="LS rain tendency due to freezing" unit="kg/kg/s" /> 655 <field id="dqsauto" long_name="LS snow tendency due to autoconversion" unit="kg/kg/s" /> 656 <field id="dqsagg" long_name="LS snow tendency due to aggregation" unit="kg/kg/s" /> 657 <field id="dqsrim" long_name="LS snow tendency due to riming" unit="kg/kg/s" /> 658 <field id="dqsmelt" long_name="LS snow tendency due to melting" unit="kg/kg/s" /> 659 <field id="dqsfreez" long_name="LS snow tendency due to freezing" unit="kg/kg/s" /> 660 661 <field id="l_mix_ter" long_name="PBL mixing length ter" unit="m" /> 645 662 <field id="l_mix_lic" long_name="PBL mixing length lic" unit="m" /> 646 663 <field id="l_mix_oce" long_name="PBL mixing length oce" unit="m" /> … … 687 704 <field id="zfull" long_name="Altitude of full pressure levels" unit="m" /> 688 705 <field id="zhalf" long_name="Altitude of half pressure levels" unit="m" /> 689 <field id="rneb" long_name="Cloud fraction" unit="-" /> 706 <field id="rneb" long_name="Cloud fraction" unit="-" /> 707 <field id="cldfraliq" long_name="Fraction of liquid cloud" unit="-" /> 708 <field id="sigma2_icefracturb" long_name="Variance of the diagnostic supersaturation distribution (icefrac_turb)" unit="-" /> 709 <field id="mean_icefracturb" long_name="Mean of the diagnostic supersaturation distribution (icefrac_turb)" unit="-" /> 690 710 <field id="rnebcon" long_name="Convective Cloud Fraction" unit="-" /> 691 711 <field id="rnebls" long_name="LS Cloud fraction" unit="-" /> -
LMDZ6/trunk/libf/phylmd/lmdz_lscp.F90
r4915 r5007 7 7 !++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 8 8 SUBROUTINE lscp(klon,klev,dtime,missing_val, & 9 paprs,pplay,temp,qt, ptconv,ratqs,&9 paprs,pplay,temp,qt,qice_save,ptconv,ratqs, & 10 10 d_t, d_q, d_ql, d_qi, rneb, rneblsvol, rneb_seri, & 11 pfraclr,pfracld, & 11 pfraclr, pfracld, & 12 cldfraliq, sigma2_icefracturb,mean_icefracturb, & 12 13 radocond, radicefrac, rain, snow, & 13 14 frac_impa, frac_nucl, beta, & 14 prfl, psfl, rhcl, qta, fraca, &15 tv, pspsk, tla, thl, iflag_cld_th, &15 prfl, psfl, rhcl, qta, fraca, & 16 tv, pspsk, tla, thl, iflag_cld_th, & 16 17 iflag_ice_thermo, ok_ice_sursat, qsatl, qsats, & 17 distcltop, temp_cltop,&18 distcltop, temp_cltop, tke, tke_dissip, & 18 19 qclr, qcld, qss, qvc, rnebclr, rnebss, gamma_ss, & 19 20 Tcontr, qcontr, qcontr2, fcontrN, fcontrP, & … … 96 97 ! USE de modules contenant des fonctions. 97 98 USE lmdz_cloudth, ONLY : cloudth, cloudth_v3, cloudth_v6, cloudth_mpc 98 USE lmdz_lscp_tools, ONLY : calc_qsat_ecmwf, icefrac_lscp, calc_gammasat 99 USE lmdz_lscp_tools, ONLY : calc_qsat_ecmwf, calc_gammasat 100 USE lmdz_lscp_tools, ONLY : icefrac_lscp, icefrac_lscp_turb 99 101 USE lmdz_lscp_tools, ONLY : fallice_velocity, distance_to_cloud_top 100 102 USE ice_sursat_mod, ONLY : ice_sursat … … 111 113 USE lmdz_lscp_ini, ONLY : RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT, RD, RG 112 114 USE lmdz_lscp_ini, ONLY : ok_poprecip 113 115 USE lmdz_lscp_ini, ONLY : iflag_icefrac 114 116 IMPLICIT NONE 115 117 … … 129 131 REAL, DIMENSION(klon,klev), INTENT(IN) :: temp ! temperature (K) 130 132 REAL, DIMENSION(klon,klev), INTENT(IN) :: qt ! total specific humidity (in vapor phase in input) [kg/kg] 133 REAL, DIMENSION(klon,klev), INTENT(IN) :: qice_save ! ice specific from previous time step [kg/kg] 131 134 INTEGER, INTENT(IN) :: iflag_cld_th ! flag that determines the distribution of convective clouds 132 135 INTEGER, INTENT(IN) :: iflag_ice_thermo! flag to activate the ice thermodynamics … … 138 141 !Inputs associated with thermal plumes 139 142 140 REAL, DIMENSION(klon,klev), INTENT(IN) :: tv ! virtual potential temperature [K] 141 REAL, DIMENSION(klon,klev), INTENT(IN) :: qta ! specific humidity within thermals [kg/kg] 142 REAL, DIMENSION(klon,klev), INTENT(IN) :: fraca ! fraction of thermals within the mesh [-] 143 REAL, DIMENSION(klon,klev), INTENT(IN) :: pspsk ! exner potential (p/100000)**(R/cp) 144 REAL, DIMENSION(klon,klev), INTENT(IN) :: tla ! liquid temperature within thermals [K] 143 REAL, DIMENSION(klon,klev), INTENT(IN) :: tv ! virtual potential temperature [K] 144 REAL, DIMENSION(klon,klev), INTENT(IN) :: qta ! specific humidity within thermals [kg/kg] 145 REAL, DIMENSION(klon,klev), INTENT(IN) :: fraca ! fraction of thermals within the mesh [-] 146 REAL, DIMENSION(klon,klev), INTENT(IN) :: pspsk ! exner potential (p/100000)**(R/cp) 147 REAL, DIMENSION(klon,klev), INTENT(IN) :: tla ! liquid temperature within thermals [K] 148 REAL, DIMENSION(klon,klev), INTENT(IN) :: tke !--turbulent kinetic energy [m2/s2] 149 REAL, DIMENSION(klon,klev), INTENT(IN) :: tke_dissip !--TKE dissipation [m2/s3] 145 150 146 151 ! INPUT/OUTPUT variables 147 152 !------------------------ 148 153 149 REAL, DIMENSION(klon,klev), INTENT(INOUT) :: thl ! liquid potential temperature [K]150 REAL, DIMENSION(klon,klev), INTENT(INOUT) :: ratqs ! function of pressure that sets the large-scale154 REAL, DIMENSION(klon,klev), INTENT(INOUT) :: thl ! liquid potential temperature [K] 155 REAL, DIMENSION(klon,klev), INTENT(INOUT) :: ratqs ! function of pressure that sets the large-scale 151 156 152 157 153 158 ! Input sursaturation en glace 154 REAL, DIMENSION(klon,klev), INTENT(INOUT):: rneb_seri ! fraction nuageuse en memoire159 REAL, DIMENSION(klon,klev), INTENT(INOUT):: rneb_seri ! fraction nuageuse en memoire 155 160 156 161 ! OUTPUT variables 157 162 !----------------- 158 163 159 REAL, DIMENSION(klon,klev), INTENT(OUT) :: d_t ! temperature increment [K] 160 REAL, DIMENSION(klon,klev), INTENT(OUT) :: d_q ! specific humidity increment [kg/kg] 161 REAL, DIMENSION(klon,klev), INTENT(OUT) :: d_ql ! liquid water increment [kg/kg] 162 REAL, DIMENSION(klon,klev), INTENT(OUT) :: d_qi ! cloud ice mass increment [kg/kg] 163 REAL, DIMENSION(klon,klev), INTENT(OUT) :: rneb ! cloud fraction [-] 164 REAL, DIMENSION(klon,klev), INTENT(OUT) :: rneblsvol ! cloud fraction per unit volume [-] 165 REAL, DIMENSION(klon,klev), INTENT(OUT) :: pfraclr ! precip fraction clear-sky part [-] 166 REAL, DIMENSION(klon,klev), INTENT(OUT) :: pfracld ! precip fraction cloudy part [-] 167 REAL, DIMENSION(klon,klev), INTENT(OUT) :: radocond ! condensed water used in the radiation scheme [kg/kg] 168 REAL, DIMENSION(klon,klev), INTENT(OUT) :: radicefrac ! ice fraction of condensed water for radiation scheme 169 REAL, DIMENSION(klon,klev), INTENT(OUT) :: rhcl ! clear-sky relative humidity [-] 170 REAL, DIMENSION(klon), INTENT(OUT) :: rain ! surface large-scale rainfall [kg/s/m2] 171 REAL, DIMENSION(klon), INTENT(OUT) :: snow ! surface large-scale snowfall [kg/s/m2] 172 REAL, DIMENSION(klon,klev), INTENT(OUT) :: qsatl ! saturation specific humidity wrt liquid [kg/kg] 173 REAL, DIMENSION(klon,klev), INTENT(OUT) :: qsats ! saturation specific humidity wrt ice [kg/kg] 174 REAL, DIMENSION(klon,klev+1), INTENT(OUT) :: prfl ! large-scale rainfall flux in the column [kg/s/m2] 175 REAL, DIMENSION(klon,klev+1), INTENT(OUT) :: psfl ! large-scale snowfall flux in the column [kg/s/m2] 176 REAL, DIMENSION(klon,klev), INTENT(OUT) :: distcltop ! distance to cloud top [m] 177 REAL, DIMENSION(klon,klev), INTENT(OUT) :: temp_cltop ! temperature of cloud top [K] 178 REAL, DIMENSION(klon,klev), INTENT(OUT) :: beta ! conversion rate of condensed water 179 180 ! fraction of aerosol scavenging through impaction and nucleation (for on-line) 164 REAL, DIMENSION(klon,klev), INTENT(OUT) :: d_t ! temperature increment [K] 165 REAL, DIMENSION(klon,klev), INTENT(OUT) :: d_q ! specific humidity increment [kg/kg] 166 REAL, DIMENSION(klon,klev), INTENT(OUT) :: d_ql ! liquid water increment [kg/kg] 167 REAL, DIMENSION(klon,klev), INTENT(OUT) :: d_qi ! cloud ice mass increment [kg/kg] 168 REAL, DIMENSION(klon,klev), INTENT(OUT) :: rneb ! cloud fraction [-] 169 REAL, DIMENSION(klon,klev), INTENT(OUT) :: rneblsvol ! cloud fraction per unit volume [-] 170 REAL, DIMENSION(klon,klev), INTENT(OUT) :: pfraclr ! precip fraction clear-sky part [-] 171 REAL, DIMENSION(klon,klev), INTENT(OUT) :: pfracld ! precip fraction cloudy part [-] 172 REAL, DIMENSION(klon,klev), INTENT(OUT) :: cldfraliq ! liquid fraction of cloud [-] 173 REAL, DIMENSION(klon,klev), INTENT(OUT) :: sigma2_icefracturb ! Variance of the diagnostic supersaturation distribution (icefrac_turb) [-] 174 REAL, DIMENSION(klon,klev), INTENT(OUT) :: mean_icefracturb ! Mean of the diagnostic supersaturation distribution (icefrac_turb) [-] 175 REAL, DIMENSION(klon,klev), INTENT(OUT) :: radocond ! condensed water used in the radiation scheme [kg/kg] 176 REAL, DIMENSION(klon,klev), INTENT(OUT) :: radicefrac ! ice fraction of condensed water for radiation scheme 177 REAL, DIMENSION(klon,klev), INTENT(OUT) :: rhcl ! clear-sky relative humidity [-] 178 REAL, DIMENSION(klon), INTENT(OUT) :: rain ! surface large-scale rainfall [kg/s/m2] 179 REAL, DIMENSION(klon), INTENT(OUT) :: snow ! surface large-scale snowfall [kg/s/m2] 180 REAL, DIMENSION(klon,klev), INTENT(OUT) :: qsatl ! saturation specific humidity wrt liquid [kg/kg] 181 REAL, DIMENSION(klon,klev), INTENT(OUT) :: qsats ! saturation specific humidity wrt ice [kg/kg] 182 REAL, DIMENSION(klon,klev+1), INTENT(OUT) :: prfl ! large-scale rainfall flux in the column [kg/s/m2] 183 REAL, DIMENSION(klon,klev+1), INTENT(OUT) :: psfl ! large-scale snowfall flux in the column [kg/s/m2] 184 REAL, DIMENSION(klon,klev), INTENT(OUT) :: distcltop ! distance to cloud top [m] 185 REAL, DIMENSION(klon,klev), INTENT(OUT) :: temp_cltop ! temperature of cloud top [K] 186 REAL, DIMENSION(klon,klev), INTENT(OUT) :: beta ! conversion rate of condensed water 187 188 ! fraction of aerosol scavenging through impaction and nucleation (for on-line) 181 189 182 REAL, DIMENSION(klon,klev), INTENT(OUT) :: frac_impa ! scavenging fraction due tu impaction [-]183 REAL, DIMENSION(klon,klev), INTENT(OUT) :: frac_nucl ! scavenging fraction due tu nucleation [-]190 REAL, DIMENSION(klon,klev), INTENT(OUT) :: frac_impa ! scavenging fraction due tu impaction [-] 191 REAL, DIMENSION(klon,klev), INTENT(OUT) :: frac_nucl ! scavenging fraction due tu nucleation [-] 184 192 185 193 ! for supersaturation and contrails parameterisation 186 194 187 REAL, DIMENSION(klon,klev), INTENT(OUT) :: qclr ! specific total water content in clear sky region [kg/kg]188 REAL, DIMENSION(klon,klev), INTENT(OUT) :: qcld ! specific total water content in cloudy region [kg/kg]189 REAL, DIMENSION(klon,klev), INTENT(OUT) :: qss ! specific total water content in supersat region [kg/kg]190 REAL, DIMENSION(klon,klev), INTENT(OUT) :: qvc ! specific vapor content in clouds [kg/kg]191 REAL, DIMENSION(klon,klev), INTENT(OUT) :: rnebclr ! mesh fraction of clear sky [-]192 REAL, DIMENSION(klon,klev), INTENT(OUT) :: rnebss ! mesh fraction of ISSR [-]193 REAL, DIMENSION(klon,klev), INTENT(OUT) :: gamma_ss ! coefficient governing the ice nucleation RHi threshold [-]194 REAL, DIMENSION(klon,klev), INTENT(OUT) :: Tcontr ! threshold temperature for contrail formation [K]195 REAL, DIMENSION(klon,klev), INTENT(OUT) :: qcontr ! threshold humidity for contrail formation [kg/kg]196 REAL, DIMENSION(klon,klev), INTENT(OUT) :: qcontr2 ! // (2nd expression more consistent with LMDZ expression of q)197 REAL, DIMENSION(klon,klev), INTENT(OUT) :: fcontrN ! fraction of grid favourable to non-persistent contrails198 REAL, DIMENSION(klon,klev), INTENT(OUT) :: fcontrP ! fraction of grid favourable to persistent contrails199 REAL, DIMENSION(klon,klev), INTENT(OUT) :: cloudth_sth ! mean saturation deficit in thermals200 REAL, DIMENSION(klon,klev), INTENT(OUT) :: cloudth_senv ! mean saturation deficit in environment201 REAL, DIMENSION(klon,klev), INTENT(OUT) :: cloudth_sigmath ! std of saturation deficit in thermals202 REAL, DIMENSION(klon,klev), INTENT(OUT) :: cloudth_sigmaenv ! std of saturation deficit in environment195 REAL, DIMENSION(klon,klev), INTENT(OUT) :: qclr ! specific total water content in clear sky region [kg/kg] 196 REAL, DIMENSION(klon,klev), INTENT(OUT) :: qcld ! specific total water content in cloudy region [kg/kg] 197 REAL, DIMENSION(klon,klev), INTENT(OUT) :: qss ! specific total water content in supersat region [kg/kg] 198 REAL, DIMENSION(klon,klev), INTENT(OUT) :: qvc ! specific vapor content in clouds [kg/kg] 199 REAL, DIMENSION(klon,klev), INTENT(OUT) :: rnebclr ! mesh fraction of clear sky [-] 200 REAL, DIMENSION(klon,klev), INTENT(OUT) :: rnebss ! mesh fraction of ISSR [-] 201 REAL, DIMENSION(klon,klev), INTENT(OUT) :: gamma_ss ! coefficient governing the ice nucleation RHi threshold [-] 202 REAL, DIMENSION(klon,klev), INTENT(OUT) :: Tcontr ! threshold temperature for contrail formation [K] 203 REAL, DIMENSION(klon,klev), INTENT(OUT) :: qcontr ! threshold humidity for contrail formation [kg/kg] 204 REAL, DIMENSION(klon,klev), INTENT(OUT) :: qcontr2 ! // (2nd expression more consistent with LMDZ expression of q) 205 REAL, DIMENSION(klon,klev), INTENT(OUT) :: fcontrN ! fraction of grid favourable to non-persistent contrails 206 REAL, DIMENSION(klon,klev), INTENT(OUT) :: fcontrP ! fraction of grid favourable to persistent contrails 207 REAL, DIMENSION(klon,klev), INTENT(OUT) :: cloudth_sth ! mean saturation deficit in thermals 208 REAL, DIMENSION(klon,klev), INTENT(OUT) :: cloudth_senv ! mean saturation deficit in environment 209 REAL, DIMENSION(klon,klev), INTENT(OUT) :: cloudth_sigmath ! std of saturation deficit in thermals 210 REAL, DIMENSION(klon,klev), INTENT(OUT) :: cloudth_sigmaenv ! std of saturation deficit in environment 203 211 204 212 ! for POPRECIP … … 221 229 ! LOCAL VARIABLES: 222 230 !---------------- 223 224 REAL,DIMENSION(klon) :: qsl, qsi 231 REAL,DIMENSION(klon) :: qsl, qsi ! saturation threshold at current vertical level 225 232 REAL :: zct, zcl,zexpo 226 233 REAL, DIMENSION(klon,klev) :: ctot … … 229 236 REAL :: zdelta, zcor, zcvm5 230 237 REAL, DIMENSION(klon) :: zdqsdT_raw 231 REAL, DIMENSION(klon) :: gammasat,dgammasatdt ! coefficient to make cold condensation at the correct RH and derivative wrt T232 REAL, DIMENSION(klon) :: Tbef,qlbef,DT 238 REAL, DIMENSION(klon) :: gammasat,dgammasatdt ! coefficient to make cold condensation at the correct RH and derivative wrt T 239 REAL, DIMENSION(klon) :: Tbef,qlbef,DT ! temperature, humidity and temp. variation during lognormal iteration 233 240 REAL :: num,denom 234 241 REAL :: cste 235 REAL, DIMENSION(klon) :: zpdf_sig,zpdf_k,zpdf_delta 236 REAL, DIMENSION(klon) :: Zpdf_a,zpdf_b,zpdf_e1,zpdf_e2 242 REAL, DIMENSION(klon) :: zpdf_sig,zpdf_k,zpdf_delta ! lognormal parameters 243 REAL, DIMENSION(klon) :: Zpdf_a,zpdf_b,zpdf_e1,zpdf_e2 ! lognormal intermediate variables 237 244 REAL :: erf 238 245 REAL, DIMENSION(klon) :: zfice_th … … 251 258 REAL :: zmelt,zrain,zsnow,zprecip 252 259 REAL, DIMENSION(klon) :: dzfice 260 REAL, DIMENSION(klon) :: zfice_turb, dzfice_turb 253 261 REAL :: zsolid 254 262 REAL, DIMENSION(klon) :: qtot, qzero … … 281 289 REAL, DIMENSION(klon,klev) :: radocondi, radocondl 282 290 REAL :: effective_zneb 283 REAL, DIMENSION(klon) :: distcltop1D, temp_cltop1D284 291 REAL, DIMENSION(klon) :: zdistcltop, ztemp_cltop 292 REAL, DIMENSION(klon) :: zqliq, zqice, zqvapcl ! for icefrac_lscp_turb 285 293 286 294 INTEGER i, k, n, kk, iter … … 337 345 pfraclr(:,:)=0.0 338 346 pfracld(:,:)=0.0 347 cldfraliq(:,:)=0. 348 sigma2_icefracturb(:,:)=0. 349 mean_icefracturb(:,:)=0. 339 350 radocond(:,:) = 0.0 340 351 radicefrac(:,:) = 0.0 … … 346 357 zfice(:)=0.0 347 358 dzfice(:)=0.0 359 zfice_turb(:)=0.0 360 dzfice_turb(:)=0.0 348 361 zqprecl(:)=0.0 349 362 zqpreci(:)=0.0 … … 360 373 d_tot_zneb(:) = 0.0 361 374 qzero(:) = 0.0 362 distcltop1D(:)=0.0363 temp_cltop1D(:) = 0.0375 zdistcltop(:)=0.0 376 ztemp_cltop(:) = 0.0 364 377 ztupnew(:)=0.0 365 378 !--ice supersaturation … … 374 387 distcltop(:,:)=0. 375 388 temp_cltop(:,:)=0. 389 376 390 !-- poprecip 377 391 qraindiag(:,:)= 0. … … 393 407 394 408 395 396 409 !c_iso: variable initialisation for iso 397 410 … … 412 425 413 426 ! Initialisation temperature and specific humidity 427 ! temp(klon,klev) is not modified by the routine, instead all changes in temperature are made on zt 428 ! at the end of the klon loop, a temperature incremtent d_t due to all processes 429 ! (thermalization, evap/sub incoming precip, cloud formation, precipitation processes) is calculated 430 ! d_t = temperature tendency due to lscp 431 ! The temperature of the overlying layer is updated here because needed for thermalization 414 432 DO i = 1, klon 415 433 zt(i)=temp(i,k) … … 746 764 ELSEIF (iflag_cloudth_vert .EQ. 7) THEN 747 765 ! Updated version of Arnaud Jam (correction by E. Vignon) + adapted treatment 748 ! for boundary-layer mixed phase clouds following Vignon et al.766 ! for boundary-layer mixed phase clouds 749 767 CALL cloudth_mpc(klon,klev,k,mpc_bl_points,zt,zq,qta(:,k),fraca(:,k), & 750 768 pspsk(:,k),paprs(:,k+1),paprs(:,k),pplay(:,k), tla(:,k), & … … 766 784 767 785 ! lognormal 768 lognormale = .TRUE.786 lognormale(:) = .TRUE. 769 787 770 788 ELSEIF (iflag_cld_th .GE. 6) THEN 771 789 772 790 ! lognormal distribution when no thermals 773 lognormale = fraca(:,k) < min_frac_th_cld791 lognormale(:) = fraca(:,k) < min_frac_th_cld 774 792 775 793 ELSE 776 794 ! When iflag_cld_th=5, we always assume 777 795 ! bi-gaussian distribution 778 lognormale = .FALSE.796 lognormale(:) = .FALSE. 779 797 780 798 ENDIF … … 829 847 IF (iflag_t_glace.GE.4) THEN 830 848 ! For iflag_t_glace GE 4 the phase partition function dependends on temperature AND distance to cloud top 831 CALL distance_to_cloud_top(klon,klev,k,temp,pplay,paprs,rneb, distcltop1D,temp_cltop1D)849 CALL distance_to_cloud_top(klon,klev,k,temp,pplay,paprs,rneb,zdistcltop,ztemp_cltop) 832 850 ENDIF 833 CALL icefrac_lscp(klon, zt(:), iflag_ice_thermo, distcltop1D(:),temp_cltop1D(:),zfice(:),dzfice(:)) 851 852 CALL icefrac_lscp(klon, zt(:), iflag_ice_thermo, zdistcltop(:),ztemp_cltop(:),zfice(:),dzfice(:)) 834 853 835 854 DO i=1,klon !todoan : check if loop in i is needed … … 896 915 cste=RLSTT 897 916 ENDIF 898 917 918 ! LEA_R : check formule 899 919 qlbef(i)=max(0.,zqn(i)-zqs(i)) 900 920 num = -Tbef(i)+zt(i)+rneb(i,k)*((1-zfice(i))*RLVTT & … … 927 947 ! For iflag_t_glace GE 4 the phase partition function dependends on temperature AND distance to cloud top 928 948 IF (iflag_t_glace.GE.4) THEN 929 CALL distance_to_cloud_top(klon,klev,k,temp,pplay,paprs,rneb,distcltop1D,temp_cltop1D) 930 distcltop(:,k)=distcltop1D(:) 931 temp_cltop(:,k)=temp_cltop1D(:) 932 ENDIF 933 ! Partition function in stratiform clouds (will be overwritten in boundary-layer MPCs) 934 CALL icefrac_lscp(klon,zt,iflag_ice_thermo,distcltop1D,temp_cltop1D,zfice,dzfice) 935 949 CALL distance_to_cloud_top(klon,klev,k,temp,pplay,paprs,rneb,zdistcltop,ztemp_cltop) 950 distcltop(:,k)=zdistcltop(:) 951 temp_cltop(:,k)=ztemp_cltop(:) 952 ENDIF 953 954 ! Partition function depending on temperature 955 CALL icefrac_lscp(klon, Tbef, iflag_ice_thermo, zdistcltop, ztemp_cltop, zfice, dzfice) 956 957 ! Partition function depending on tke for non shallow-convective clouds 958 IF (iflag_icefrac .GE. 1) THEN 959 960 CALL icefrac_lscp_turb(klon, dtime, Tbef, pplay(:,k), paprs(:,k), paprs(:,k+1), qice_save(:,k), ziflcld, zqn, & 961 rneb(:,k), tke(:,k), tke_dissip(:,k), zqliq, zqvapcl, zqice, zfice_turb, dzfice_turb, cldfraliq(:,k),sigma2_icefracturb(:,k), mean_icefracturb(:,k)) 962 963 ENDIF 936 964 937 965 ! Water vapor update, Phase determination and subsequent latent heat exchange 938 966 DO i=1, klon 939 967 ! Overwrite phase partitioning in boundary layer mixed phase clouds when the 968 ! iflag_cloudth_vert=7 and specific param is activated 940 969 IF (mpc_bl_points(i,k) .GT. 0) THEN 941 942 970 zcond(i) = MAX(0.0,qincloud_mpc(i))*rneb(i,k) 943 971 ! following line is very strange and probably wrong … … 946 974 zq(i) = zq(i) - zcond(i) 947 975 zfice(i)=zfice_th(i) 948 949 976 ELSE 950 951 977 ! Checks on rneb, rhcl and zqn 952 978 IF (rneb(i,k) .LE. 0.0) THEN … … 964 990 ! following line is very strange and probably wrong: 965 991 rhcl(i,k)=(zqs(i)+zq(i))/2./zqs(i) 992 ! Overwrite partitioning for non shallow-convective clouds if iflag_icefrac>1 (icefrac turb param) 993 IF (iflag_icefrac .GE. 1) THEN 994 IF (lognormale(i)) THEN 995 zcond(i) = zqliq(i) + zqice(i) 996 zfice(i)=zfice_turb(i) 997 rhcl(i,k) = zqvapcl(i) * rneb(i,k) + (zq(i) - zqn(i)) * (1.-rneb(i,k)) 998 ENDIF 999 ENDIF 966 1000 ENDIF 967 1001 … … 1306 1340 znebprecipcld(i)=0.0 1307 1341 ENDIF 1308 1342 !IF ( ((1-zfice(i))*zoliq(i) .GT. 0.) .AND. (zt(i) .LE. 233.15) ) THEN 1343 !print*,'WARNING LEA OLIQ A <-40°C ' 1344 !print*,'zt,Tbef,oliq,oice,cldfraliq,icefrac,rneb',zt(i),Tbef(i),(1-zfice(i))*zoliq(i),zfice(i)*zoliq(i),cldfraliq(i,k),zfice(i),rneb(i,k) 1345 !ENDIF 1309 1346 ENDDO 1310 1347 -
LMDZ6/trunk/libf/phylmd/lmdz_lscp_ini.F90
r4915 r5007 6 6 !-------------------- 7 7 8 REAL RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT, RD, RG, R PI9 !$OMP THREADPRIVATE(RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT, RD, RG, R PI)8 REAL RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT, RD, RG, RV, RPI 9 !$OMP THREADPRIVATE(RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT, RD, RG, RV, RPI) 10 10 11 11 REAL, SAVE, PROTECTED :: seuil_neb=0.001 ! cloud fraction threshold: a cloud can precipitate when exceeded … … 67 67 !$OMP THREADPRIVATE(iflag_t_glace) 68 68 69 INTEGER, SAVE, PROTECTED :: iflag_cloudth_vert=0 ! option for determining cloud fraction and content in convective boundary layers69 INTEGER, SAVE, PROTECTED :: iflag_cloudth_vert=0 ! option for determining cloud fraction and content in convective boundary layers 70 70 !$OMP THREADPRIVATE(iflag_cloudth_vert) 71 71 72 INTEGER, SAVE, PROTECTED :: iflag_gammasat=0 ! which threshold for homogeneous nucleation below -40oC72 INTEGER, SAVE, PROTECTED :: iflag_gammasat=0 ! which threshold for homogeneous nucleation below -40oC 73 73 !$OMP THREADPRIVATE(iflag_gammasat) 74 74 75 INTEGER, SAVE, PROTECTED :: iflag_rain_incloud_vol=0 ! use of volume cloud fraction for rain autoconversion75 INTEGER, SAVE, PROTECTED :: iflag_rain_incloud_vol=0 ! use of volume cloud fraction for rain autoconversion 76 76 !$OMP THREADPRIVATE(iflag_rain_incloud_vol) 77 77 78 INTEGER, SAVE, PROTECTED :: iflag_bergeron=0 ! bergeron effect for liquid precipitation treatment78 INTEGER, SAVE, PROTECTED :: iflag_bergeron=0 ! bergeron effect for liquid precipitation treatment 79 79 !$OMP THREADPRIVATE(iflag_bergeron) 80 80 81 INTEGER, SAVE, PROTECTED :: iflag_fisrtilp_qsat=0 ! qsat adjustment (iterative) during autoconversion81 INTEGER, SAVE, PROTECTED :: iflag_fisrtilp_qsat=0 ! qsat adjustment (iterative) during autoconversion 82 82 !$OMP THREADPRIVATE(iflag_fisrtilp_qsat) 83 83 84 INTEGER, SAVE, PROTECTED :: iflag_pdf=0 ! type of subgrid scale qtot pdf84 INTEGER, SAVE, PROTECTED :: iflag_pdf=0 ! type of subgrid scale qtot pdf 85 85 !$OMP THREADPRIVATE(iflag_pdf) 86 86 87 INTEGER, SAVE, PROTECTED :: iflag_autoconversion=0 ! autoconversion option 87 INTEGER, SAVE, PROTECTED :: iflag_icefrac=0 ! which phase partitioning function to use 88 !$OMP THREADPRIVATE(iflag_icefrac) 89 90 INTEGER, SAVE, PROTECTED :: iflag_autoconversion=0 ! autoconversion option 88 91 !$OMP THREADPRIVATE(iflag_autoconversion) 89 92 90 LOGICAL, SAVE, PROTECTED :: reevap_ice=.false. ! no liquid precip for T< threshold 93 94 LOGICAL, SAVE, PROTECTED :: reevap_ice=.false. ! no liquid precip for T< threshold 91 95 !$OMP THREADPRIVATE(reevap_ice) 92 96 93 REAL, SAVE, PROTECTED :: cld_lc_lsc=2.6e-4 ! liquid autoconversion coefficient, stratiform rain97 REAL, SAVE, PROTECTED :: cld_lc_lsc=2.6e-4 ! liquid autoconversion coefficient, stratiform rain 94 98 !$OMP THREADPRIVATE(cld_lc_lsc) 95 99 … … 118 122 !$OMP THREADPRIVATE(coef_eva) 119 123 120 REAL, SAVE, PROTECTED :: coef_sub ! tuning coefficient ice precip sublimation124 REAL, SAVE, PROTECTED :: coef_sub ! tuning coefficient ice precip sublimation 121 125 !$OMP THREADPRIVATE(coef_sub) 122 126 … … 124 128 !$OMP THREADPRIVATE(expo_eva) 125 129 126 REAL, SAVE, PROTECTED :: expo_sub ! tuning coefficient ice precip sublimation130 REAL, SAVE, PROTECTED :: expo_sub ! tuning coefficient ice precip sublimation 127 131 !$OMP THREADPRIVATE(expo_sub) 128 132 … … 158 162 !$OMP THREADPRIVATE(thresh_precip_frac) 159 163 164 REAL, SAVE, PROTECTED :: tau_mixenv=100000 ! Homogeneization time of mixed phase clouds [s] 165 !$OMP THREADPRIVATE(tau_mixenv) 166 167 REAL, SAVE, PROTECTED :: capa_crystal=1. ! Sursaturation of ice part in mixed phase clouds [-] 168 !$OMP THREADPRIVATE(capa_crystal) 169 170 REAL, SAVE, PROTECTED :: lmix_mpc=1000 ! Length of turbulent zones in Mixed Phase Clouds [m] 171 !$OMP THREADPRIVATE(lmix_mpc) 172 173 REAL, SAVE, PROTECTED :: naero5=0.5 ! Number concentration of aerosol larger than 0.5 microns [scm-3] 174 !$OMP THREADPRIVATE(naero5) 175 176 REAL, SAVE, PROTECTED :: gamma_snwretro = 0. ! Proportion of snow taken into account in ice retroaction in icefrac_turb [-] 177 !$OMP THREADPRIVATE(gamma_snwretro) 178 179 REAL, SAVE, PROTECTED :: gamma_taud = 1. ! Tuning coeff for tau_dissipturb [-] 180 !$OMP THREADPRIVATE(gamma_taud) 181 160 182 REAL, SAVE, PROTECTED :: gamma_col=1. ! A COMMENTER TODO [-] 161 183 !$OMP THREADPRIVATE(gamma_col) … … 167 189 !$OMP THREADPRIVATE(gamma_rim) 168 190 169 REAL, SAVE, PROTECTED :: rho_rain=1000. ! A COMMENTER TODO[kg/m3]191 REAL, SAVE, PROTECTED :: rho_rain=1000. ! Rain density [kg/m3] 170 192 !$OMP THREADPRIVATE(rho_rain) 171 193 172 REAL, SAVE, PROTECTED :: rho_ice=920. ! A COMMENTER TODO[kg/m3]194 REAL, SAVE, PROTECTED :: rho_ice=920. ! Ice density [kg/m3] 173 195 !$OMP THREADPRIVATE(rho_ice) 174 196 175 REAL, SAVE, PROTECTED :: r_rain=500.E-6 ! A COMMENTER TODO[m]197 REAL, SAVE, PROTECTED :: r_rain=500.E-6 ! Rain droplets radius for POPRECIP [m] 176 198 !$OMP THREADPRIVATE(r_rain) 177 199 178 REAL, SAVE, PROTECTED :: r_snow=1.E-3 ! A COMMENTER TODO[m]200 REAL, SAVE, PROTECTED :: r_snow=1.E-3 ! Ice crystals radius for POPRECIP [m] 179 201 !$OMP THREADPRIVATE(r_snow) 180 202 181 REAL, SAVE, PROTECTED :: tau_auto_snow_min=100. ! A COMMENTER TODO [s]203 REAL, SAVE, PROTECTED :: tau_auto_snow_min=100. ! A COMMENTER TODO [s] 182 204 !$OMP THREADPRIVATE(tau_auto_snow_min) 183 205 … … 188 210 !$OMP THREADPRIVATE(eps) 189 211 190 REAL, SAVE, PROTECTED :: gamma_melt=1. ! A COMMENTER TODO [-]212 REAL, SAVE, PROTECTED :: gamma_melt=1. ! A COMMENTER TODO [-] 191 213 !$OMP THREADPRIVATE(gamma_melt) 192 214 193 REAL, SAVE, PROTECTED :: alpha_freez=4. ! A COMMENTER TODO [-]215 REAL, SAVE, PROTECTED :: alpha_freez=4. ! A COMMENTER TODO [-] 194 216 !$OMP THREADPRIVATE(alpha_freez) 195 217 196 REAL, SAVE, PROTECTED :: beta_freez=0.1 ! A COMMENTER TODO [m-3.s-1]218 REAL, SAVE, PROTECTED :: beta_freez=0.1 ! A COMMENTER TODO [m-3.s-1] 197 219 !$OMP THREADPRIVATE(beta_freez) 198 220 199 REAL, SAVE, PROTECTED :: gamma_freez=1. ! A COMMENTER TODO [-]221 REAL, SAVE, PROTECTED :: gamma_freez=1. ! A COMMENTER TODO [-] 200 222 !$OMP THREADPRIVATE(gamma_freez) 201 223 202 REAL, SAVE, PROTECTED :: rain_fallspeed=4. ! A COMMENTER TODO [m/s]224 REAL, SAVE, PROTECTED :: rain_fallspeed=4. ! A COMMENTER TODO [m/s] 203 225 !$OMP THREADPRIVATE(rain_fallspeed) 204 226 205 REAL, SAVE, PROTECTED :: rain_fallspeed_clr ! A COMMENTER TODO [m/s]227 REAL, SAVE, PROTECTED :: rain_fallspeed_clr ! A COMMENTER TODO [m/s] 206 228 !$OMP THREADPRIVATE(rain_fallspeed_clr) 207 229 208 REAL, SAVE, PROTECTED :: rain_fallspeed_cld ! A COMMENTER TODO [m/s]230 REAL, SAVE, PROTECTED :: rain_fallspeed_cld ! A COMMENTER TODO [m/s] 209 231 !$OMP THREADPRIVATE(rain_fallspeed_cld) 210 232 211 REAL, SAVE, PROTECTED :: snow_fallspeed=1. ! A COMMENTER TODO [m/s]233 REAL, SAVE, PROTECTED :: snow_fallspeed=1. ! A COMMENTER TODO [m/s] 212 234 !$OMP THREADPRIVATE(snow_fallspeed) 213 235 214 REAL, SAVE, PROTECTED :: snow_fallspeed_clr ! A COMMENTER TODO [m/s]236 REAL, SAVE, PROTECTED :: snow_fallspeed_clr ! A COMMENTER TODO [m/s] 215 237 !$OMP THREADPRIVATE(snow_fallspeed_clr) 216 238 217 REAL, SAVE, PROTECTED :: snow_fallspeed_cld ! A COMMENTER TODO [m/s]239 REAL, SAVE, PROTECTED :: snow_fallspeed_cld ! A COMMENTER TODO [m/s] 218 240 !$OMP THREADPRIVATE(snow_fallspeed_cld) 219 241 !--End of the parameters for poprecip … … 226 248 227 249 SUBROUTINE lscp_ini(dtime,lunout_in,prt_level_in,ok_ice_sursat, iflag_ratqs, fl_cor_ebil_in, RCPD_in, RLSTT_in, RLVTT_in, RLMLT_in, & 228 RVTMP2_in, RTT_in,RD_in,RG_in,R PI_in)250 RVTMP2_in, RTT_in,RD_in,RG_in,RV_in,RPI_in) 229 251 230 252 … … 238 260 239 261 REAL, INTENT(IN) :: RCPD_in, RLSTT_in, RLVTT_in, RLMLT_in 240 REAL, INTENT(IN) :: RVTMP2_in, RTT_in, RD_in, RG_in, R PI_in262 REAL, INTENT(IN) :: RVTMP2_in, RTT_in, RD_in, RG_in, RV_in, RPI_in 241 263 character (len=20) :: modname='lscp_ini_mod' 242 264 character (len=80) :: abort_message … … 254 276 RLMLT=RLMLT_in 255 277 RTT=RTT_in 256 R G=RG_in278 RV=RV_in 257 279 RVTMP2=RVTMP2_in 258 280 RPI=RPI_in … … 275 297 CALL getin_p('iflag_fisrtilp_qsat',iflag_fisrtilp_qsat) 276 298 CALL getin_p('iflag_pdf',iflag_pdf) 299 CALL getin_p('iflag_icefrac',iflag_icefrac) 277 300 CALL getin_p('reevap_ice',reevap_ice) 278 301 CALL getin_p('cld_lc_lsc',cld_lc_lsc) … … 296 319 CALL getin_p('dist_liq',dist_liq) 297 320 CALL getin_p('tresh_cl',tresh_cl) 321 CALL getin_p('tau_mixenv',tau_mixenv) 322 CALL getin_p('capa_crystal',capa_crystal) 323 CALL getin_p('lmix_mpc',lmix_mpc) 324 CALL getin_p('naero5',naero5) 325 CALL getin_p('gamma_snwretro',gamma_snwretro) 326 CALL getin_p('gamma_taud',gamma_taud) 298 327 CALL getin_p('iflag_oldbug_fisrtilp',iflag_oldbug_fisrtilp) 299 328 CALL getin_p('ok_poprecip',ok_poprecip) … … 334 363 WRITE(lunout,*) 'lscp_ini, iflag_fisrtilp_qsat:', iflag_fisrtilp_qsat 335 364 WRITE(lunout,*) 'lscp_ini, iflag_pdf', iflag_pdf 365 WRITE(lunout,*) 'lscp_ini, iflag_icefrac', iflag_icefrac 336 366 WRITE(lunout,*) 'lscp_ini, reevap_ice', reevap_ice 337 367 WRITE(lunout,*) 'lscp_ini, cld_lc_lsc', cld_lc_lsc … … 352 382 WRITE(lunout,*) 'lscp_ini, dist_liq', dist_liq 353 383 WRITE(lunout,*) 'lscp_ini, tresh_cl', tresh_cl 384 WRITE(lunout,*) 'lscp_ini, tau_mixenv', tau_mixenv 385 WRITE(lunout,*) 'lscp_ini, capa_crystal', capa_crystal 386 WRITE(lunout,*) 'lscp_ini, lmix_mpc', lmix_mpc 387 WRITE(lunout,*) 'lscp_ini, naero5', naero5 388 WRITE(lunout,*) 'lscp_ini, gamma_snwretro', gamma_snwretro 389 WRITE(lunout,*) 'lscp_ini, gamma_taud', gamma_taud 354 390 WRITE(lunout,*) 'lscp_ini, iflag_oldbug_fisrtilp', iflag_oldbug_fisrtilp 355 391 WRITE(lunout,*) 'lscp_ini, fl_cor_ebil', fl_cor_ebil -
LMDZ6/trunk/libf/phylmd/lmdz_lscp_tools.F90
r4818 r5007 136 136 CHARACTER (len = 80) :: abort_message 137 137 138 IF ((iflag_t_glace.LT.2) 138 IF ((iflag_t_glace.LT.2)) THEN !.OR. (iflag_t_glace.GT.6)) THEN 139 139 abort_message = 'lscp cannot be used if iflag_t_glace<2 or >6' 140 140 CALL abort_physic(modname,abort_message,1) … … 194 194 195 195 ! with CMIP6 function of temperature at cloud top 196 IF ( iflag_t_glace .EQ. 5) THEN196 IF ((iflag_t_glace .EQ. 5) .OR. (iflag_t_glace .EQ. 7)) THEN 197 197 liqfrac_tmp = (temp(i)-t_glace_min) / (t_glace_max-t_glace_min) 198 198 liqfrac_tmp = MIN(MAX(liqfrac_tmp,0.0),1.0) … … 232 232 ENDIF 233 233 ENDIF 234 234 235 235 236 236 ENDDO ! klon 237 238 237 RETURN 239 238 … … 241 240 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 242 241 242 SUBROUTINE ICEFRAC_LSCP_TURB(klon, dtime, temp, pplay, paprsdn, paprsup, qice_ini, snowcld, qtot_incl, cldfra, tke, tke_dissip, qliq, qvap_cld, qice, icefrac, dicefracdT, cldfraliq, sigma2_icefracturb, mean_icefracturb) 243 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 244 ! Compute the liquid, ice and vapour content (+ice fraction) based 245 ! on turbulence (see Fields 2014, Furtado 2016) 246 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 247 248 249 USE lmdz_lscp_ini, ONLY : prt_level, lunout 250 USE lmdz_lscp_ini, ONLY : RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT, RD, RG, RV, RPI 251 USE lmdz_lscp_ini, ONLY : seuil_neb, temp_nowater 252 USE lmdz_lscp_ini, ONLY : tau_mixenv, lmix_mpc, naero5, gamma_snwretro, gamma_taud, capa_crystal 253 USE lmdz_lscp_ini, ONLY : eps 254 255 IMPLICIT NONE 256 257 INTEGER, INTENT(IN) :: klon !--number of horizontal grid points 258 REAL, INTENT(IN) :: dtime !--time step [s] 259 260 REAL, INTENT(IN), DIMENSION(klon) :: temp !--temperature 261 REAL, INTENT(IN), DIMENSION(klon) :: pplay !--pressure in the middle of the layer [Pa] 262 REAL, INTENT(IN), DIMENSION(klon) :: paprsdn !--pressure at the bottom interface of the layer [Pa] 263 REAL, INTENT(IN), DIMENSION(klon) :: paprsup !--pressure at the top interface of the layer [Pa] 264 REAL, INTENT(IN), DIMENSION(klon) :: qtot_incl !--specific total cloud water content, in-cloud content [kg/kg] 265 REAL, INTENT(IN), DIMENSION(klon) :: cldfra !--cloud fraction in gridbox [-] 266 REAL, INTENT(IN), DIMENSION(klon) :: tke !--turbulent kinetic energy [m2/s2] 267 REAL, INTENT(IN), DIMENSION(klon) :: tke_dissip !--TKE dissipation [m2/s3] 268 269 REAL, INTENT(IN), DIMENSION(klon) :: qice_ini !--initial specific ice content gridbox-mean [kg/kg] 270 REAL, INTENT(IN), DIMENSION(klon) :: snowcld 271 REAL, INTENT(OUT), DIMENSION(klon) :: qliq !--specific liquid content gridbox-mean [kg/kg] 272 REAL, INTENT(OUT), DIMENSION(klon) :: qvap_cld !--specific cloud vapor content, gridbox-mean [kg/kg] 273 REAL, INTENT(OUT), DIMENSION(klon) :: qice !--specific ice content gridbox-mean [kg/kg] 274 REAL, INTENT(OUT), DIMENSION(klon) :: icefrac !--fraction of ice in condensed water [-] 275 REAL, INTENT(OUT), DIMENSION(klon) :: dicefracdT 276 277 REAL, INTENT(OUT), DIMENSION(klon) :: cldfraliq !--fraction of cldfra which is liquid only 278 REAL, INTENT(OUT), DIMENSION(klon) :: sigma2_icefracturb !--Temporary 279 REAL, INTENT(OUT), DIMENSION(klon) :: mean_icefracturb !--Temporary 280 281 REAL, DIMENSION(klon) :: qzero, qsatl, dqsatl, qsati, dqsati !--specific humidity saturation values 282 INTEGER :: i 283 284 REAL :: qvap_incl, qice_incl, qliq_incl, qiceini_incl !--In-cloud specific quantities [kg/kg] 285 REAL :: qsnowcld_incl 286 !REAL :: capa_crystal !--Capacitance of ice crystals [-] 287 REAL :: water_vapor_diff !--Water-vapour diffusion coefficient in air [m2/s] (function of T&P) 288 REAL :: air_thermal_conduct !--Thermal conductivity of air [J/m/K/s] (function of T) 289 REAL :: C0 !--Lagrangian structure function [-] 290 REAL :: tau_mixingenv 291 REAL :: tau_dissipturb 292 REAL :: invtau_phaserelax 293 REAL :: sigma2_pdf, mean_pdf 294 REAL :: ai, bi, B0 295 REAL :: sursat_iceliq 296 REAL :: sursat_env 297 REAL :: liqfra_max 298 REAL :: sursat_iceext 299 REAL :: nb_crystals !--number concentration of ice crystals [#/m3] 300 REAL :: moment1_PSD !--1st moment of ice PSD 301 REAL :: N0_PSD, lambda_PSD !--parameters of the exponential PSD 302 303 REAL :: rho_ice !--ice density [kg/m3] 304 REAL :: cldfra1D 305 REAL :: deltaz, rho_air 306 REAL :: psati !--saturation vapor pressure wrt i [Pa] 307 308 C0 = 10. !--value assumed in Field2014 309 rho_ice = 950. 310 sursat_iceext = -0.1 311 !capa_crystal = 1. !r_ice 312 qzero(:) = 0. 313 cldfraliq(:) = 0. 314 icefrac(:) = 0. 315 dicefracdT(:) = 0. 316 317 sigma2_icefracturb(:) = 0. 318 mean_icefracturb(:) = 0. 319 320 !--wrt liquid water 321 CALL calc_qsat_ecmwf(klon,temp(:),qzero(:),pplay(:),RTT,1,.false.,qsatl(:),dqsatl(:)) 322 !--wrt ice 323 CALL calc_qsat_ecmwf(klon,temp(:),qzero(:),pplay(:),RTT,2,.false.,qsati(:),dqsati(:)) 324 325 326 DO i=1,klon 327 328 329 rho_air = pplay(i) / temp(i) / RD 330 !deltaz = ( paprsdn(i) - paprsup(i) ) / RG / rho_air(i) 331 ! because cldfra is intent in, but can be locally modified due to test 332 cldfra1D = cldfra(i) 333 IF (cldfra(i) .LE. 0.) THEN 334 qvap_cld(i) = 0. 335 qliq(i) = 0. 336 qice(i) = 0. 337 cldfraliq(i) = 0. 338 icefrac(i) = 0. 339 dicefracdT(i) = 0. 340 341 ! If there is a cloud 342 ELSE 343 IF (cldfra(i) .GE. 1.0) THEN 344 cldfra1D = 1.0 345 END IF 346 347 ! T>0°C, no ice allowed 348 IF ( temp(i) .GE. RTT ) THEN 349 qvap_cld(i) = qsatl(i) * cldfra1D 350 qliq(i) = MAX(0.0,qtot_incl(i)-qsatl(i)) * cldfra1D 351 qice(i) = 0. 352 cldfraliq(i) = 1. 353 icefrac(i) = 0. 354 dicefracdT(i) = 0. 355 356 ! T<-38°C, no liquid allowed 357 ELSE IF ( temp(i) .LE. temp_nowater) THEN 358 qvap_cld(i) = qsati(i) * cldfra1D 359 qliq(i) = 0. 360 qice(i) = MAX(0.0,qtot_incl(i)-qsati(i)) * cldfra1D 361 cldfraliq(i) = 0. 362 icefrac(i) = 1. 363 dicefracdT(i) = 0. 364 365 ! MPC temperature 366 ELSE 367 ! Not enough TKE 368 IF ( tke_dissip(i) .LE. eps ) THEN 369 qvap_cld(i) = qsati(i) * cldfra1D 370 qliq(i) = 0. 371 qice(i) = MAX(0.,qtot_incl(i)-qsati(i)) * cldfra1D 372 cldfraliq(i) = 0. 373 icefrac(i) = 1. 374 dicefracdT(i) = 0. 375 376 ! Enough TKE 377 ELSE 378 !--------------------------------------------------------- 379 !-- ICE SUPERSATURATION PDF 380 !--------------------------------------------------------- 381 !--If -38°C< T <0°C and there is enough turbulence, 382 !--we compute the cloud liquid properties with a Gaussian PDF 383 !--of ice supersaturation F(Si) (Field2014, Furtado2016). 384 !--Parameters of the PDF are function of turbulence and 385 !--microphysics/existing ice. 386 387 sursat_iceliq = qsatl(i)/qsati(i) - 1. 388 psati = qsati(i) * pplay(i) / (RD/RV) 389 390 !-------------- MICROPHYSICAL TERMS -------------- 391 !--We assume an exponential ice PSD whose parameters 392 !--are computed following Morrison&Gettelman 2008 393 !--Ice number density is assumed equals to INP density 394 !--which is a function of temperature (DeMott 2010) 395 !--bi and B0 are microphysical function characterizing 396 !--vapor/ice interactions 397 !--tau_phase_relax is the typical time of vapor deposition 398 !--onto ice crystals 399 400 qiceini_incl = qice_ini(i) / cldfra1D 401 qsnowcld_incl = snowcld(i) * RG * dtime / ( paprsdn(i) - paprsup(i) ) / cldfra1D 402 sursat_env = max(0., (qtot_incl(i) - qiceini_incl)/qsati(i) - 1.) 403 IF ( qiceini_incl .GT. eps ) THEN 404 nb_crystals = 1.e3 * 5.94e-5 * ( RTT - temp(i) )**3.33 * naero5**(0.0264*(RTT-temp(i))+0.0033) 405 lambda_PSD = ( (RPI*rho_ice*nb_crystals*24.) / (6.*(qiceini_incl + gamma_snwretro * qsnowcld_incl)) ) ** (1./3.) 406 N0_PSD = nb_crystals * lambda_PSD 407 moment1_PSD = N0_PSD/2./lambda_PSD**2 408 ELSE 409 moment1_PSD = 0. 410 ENDIF 411 412 !--Formulae for air thermal conductivity and water vapor diffusivity 413 !--comes respectively from Beard and Pruppacher (1971) 414 !--and Hall and Pruppacher (1976) 415 416 air_thermal_conduct = ( 5.69 + 0.017 * ( temp(i) - RTT ) ) * 1.e-3 * 4.184 417 water_vapor_diff = 2.11*1e-5 * ( temp(i) / RTT )**1.94 * ( 101325 / pplay(i) ) 418 419 bi = 1./((qsati(i)+qsatl(i))/2.) + RLSTT**2 / RCPD / RV / temp(i)**2 420 B0 = 4. * RPI * capa_crystal * 1. / ( RLSTT**2 / air_thermal_conduct / RV / temp(i)**2 & 421 + RV * temp(i) / psati / water_vapor_diff ) 422 423 invtau_phaserelax = (bi * B0 * moment1_PSD ) 424 425 ! Old way of estimating moment1 : spherical crystals + monodisperse PSD 426 ! nb_crystals = rho_air * qiceini_incl / ( 4. / 3. * RPI * r_ice**3. * rho_ice ) 427 ! moment1_PSD = nb_crystals * r_ice 428 429 !----------------- TURBULENT SOURCE/SINK TERMS ----------------- 430 !--Tau_mixingenv is the time needed to homogeneize the parcel 431 !--with its environment by turbulent diffusion over the parcel 432 !--length scale 433 !--if lmix_mpc <0, tau_mixigenv value is prescribed 434 !--else tau_mixigenv value is derived from tke_dissip and lmix_mpc 435 !--Tau_dissipturb is the time needed turbulence to decay due to 436 !--viscosity 437 438 ai = RG / RD / temp(i) * ( RD * RLSTT / RCPD / RV / temp(i) - 1. ) 439 IF ( lmix_mpc .GT. 0 ) THEN 440 tau_mixingenv = ( lmix_mpc**2. / tke_dissip(i) )**(1./3.) 441 ELSE 442 tau_mixingenv = tau_mixenv 443 ENDIF 444 445 tau_dissipturb = gamma_taud * 2. * 2./3. * tke(i) / tke_dissip(i) / C0 446 447 !--------------------- PDF COMPUTATIONS --------------------- 448 !--Formulae for sigma2_pdf (variance), mean of PDF in Furtado2016 449 !--cloud liquid fraction and in-cloud liquid content are given 450 !--by integrating resp. F(Si) and Si*F(Si) 451 !--Liquid is limited by the available water vapor trough a 452 !--maximal liquid fraction 453 454 liqfra_max = MAX(0., (MIN (1.,( qtot_incl(i) - qiceini_incl - qsati(i) * (1 + sursat_iceext ) ) / ( qsatl(i) - qsati(i) ) ) ) ) 455 sigma2_pdf = 1./2. * ( ai**2 ) * 2./3. * tke(i) * tau_dissipturb / ( invtau_phaserelax + 1./tau_mixingenv ) 456 mean_pdf = sursat_env * 1./tau_mixingenv / ( invtau_phaserelax + 1./tau_mixingenv ) 457 cldfraliq(i) = 0.5 * (1. - erf( ( sursat_iceliq - mean_pdf) / (SQRT(2.* sigma2_pdf) ) ) ) 458 !IF (cldfraliq(i) .GT. liqfra_max) THEN 459 ! cldfraliq(i) = liqfra_max 460 !ENDIF 461 462 qliq_incl = qsati(i) * SQRT(sigma2_pdf) / SQRT(2.*RPI) * EXP( -1.*(sursat_iceliq - mean_pdf)**2. / (2.*sigma2_pdf) ) & 463 - qsati(i) * cldfraliq(i) * (sursat_iceliq - mean_pdf ) 464 465 sigma2_icefracturb(i)= sigma2_pdf 466 mean_icefracturb(i) = mean_pdf 467 !------------ ICE AMOUNT AND WATER CONSERVATION ------------ 468 469 IF ( (qliq_incl .LE. eps) .OR. (cldfraliq(i) .LE. eps) ) THEN 470 qliq_incl = 0. 471 cldfraliq(i) = 0. 472 END IF 473 474 !--Choice for in-cloud vapor : 475 !--1.Weighted mean between qvap in MPC parts and in ice-only parts 476 !--2.Always at ice saturation 477 qvap_incl = MAX(qsati(i), ( 1. - cldfraliq(i) ) * (sursat_iceext + 1.) * qsati(i) + cldfraliq(i) * qsatl(i) ) 478 479 IF ( qvap_incl .GE. qtot_incl(i) ) THEN 480 qvap_incl = qsati(i) 481 qliq_incl = qtot_incl(i) - qvap_incl 482 qice_incl = 0. 483 484 ELSEIF ( (qvap_incl + qliq_incl) .GE. qtot_incl(i) ) THEN 485 qliq_incl = MAX(0.0,qtot_incl(i) - qvap_incl) 486 qice_incl = 0. 487 ELSE 488 qice_incl = qtot_incl(i) - qvap_incl - qliq_incl 489 END IF 490 491 qvap_cld(i) = qvap_incl * cldfra1D 492 qliq(i) = qliq_incl * cldfra1D 493 qice(i) = qice_incl * cldfra1D 494 icefrac(i) = qice(i) / ( qice(i) + qliq(i) ) 495 dicefracdT(i) = 0. 496 !print*,'MPC turb' 497 498 END IF ! Enough TKE 499 500 END IF ! ! MPC temperature 501 502 END IF ! cldfra 503 504 ENDDO ! klon 505 END SUBROUTINE ICEFRAC_LSCP_TURB 506 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 243 507 244 508 -
LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90
r4998 r5007 478 478 REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: pfraclr,pfracld 479 479 !$OMP THREADPRIVATE(pfraclr,pfracld) 480 REAL, SAVE, ALLOCATABLE :: cldfraliq(:,:) 481 !$OMP THREADPRIVATE(cldfraliq) 482 REAL, SAVE, ALLOCATABLE ::mean_icefracturb(:,:) 483 !$OMP THREADPRIVATE(mean_icefracturb) 484 REAL, SAVE, ALLOCATABLE :: sigma2_icefracturb(:,:) 485 !$OMP THREADPRIVATE(sigma2_icefracturb) 480 486 481 487 ! variables de sorties MM … … 957 963 ALLOCATE(pfraclr(klon,klev),pfracld(klon,klev)) 958 964 pfraclr(:,:)=0. ; pfracld(:,:)=0. ! because not always defined 965 ALLOCATE(cldfraliq(klon,klev)) 966 ALLOCATE(sigma2_icefracturb(klon,klev)) 967 ALLOCATE(mean_icefracturb(klon,klev)) 959 968 ALLOCATE(distcltop(klon,klev)) 960 969 ALLOCATE(temp_cltop(klon,klev)) … … 1283 1292 DEALLOCATE(rneb) 1284 1293 DEALLOCATE(pfraclr,pfracld) 1294 DEALLOCATE(cldfraliq) 1295 DEALLOCATE(sigma2_icefracturb) 1296 DEALLOCATE(mean_icefracturb) 1285 1297 DEALLOCATE (zxsnow,snowhgt,qsnow,to_ice,sissnow,runoff,albsol3_lic) 1286 1298 DEALLOCATE(distcltop) -
LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90
r4998 r5007 1559 1559 TYPE(ctrl_out), SAVE :: o_rneb = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11, 11/), & 1560 1560 'rneb', 'Cloud fraction', '-', (/ ('', i=1, 10) /)) 1561 TYPE(ctrl_out), SAVE :: o_cldfraliq = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), & 1562 'cldfraliq', 'Liquid fraction of the cloud', '-', (/ ('', i=1, 10) /)) 1563 TYPE(ctrl_out), SAVE :: o_sigma2_icefracturb = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), & 1564 'sigma2_icefracturb', 'Variance of the diagnostic supersaturation distribution (icefrac_turb) [-]', '-', (/ ('', i=1, 10) /)) 1565 TYPE(ctrl_out), SAVE :: o_mean_icefracturb = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), & 1566 'mean_icefracturb', 'Mean of the diagnostic supersaturation distribution (icefrac_turb) [-]', '-', (/ ('', i=1, 10) /)) 1567 1561 1568 TYPE(ctrl_out), SAVE :: o_rnebjn = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11,11, 11/), & 1562 1569 'rnebjn', 'Cloud fraction in day', '-', (/ ('', i=1, 10) /)) -
LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90
r4998 r5007 140 140 o_zfull, o_zhalf, o_rneb, o_rnebjn, o_rnebcon, & 141 141 o_rnebls, o_rneblsvol, o_rhum, o_rhl, o_rhi, o_ozone, o_ozone_light, & 142 o_pfraclr, o_pfracld, &142 o_pfraclr, o_pfracld, o_cldfraliq, o_sigma2_icefracturb, o_mean_icefracturb, & 143 143 o_qrainlsc, o_qsnowlsc, o_dqreva, o_dqrauto, o_dqrcol, o_dqrmelt, o_dqrfreez, & 144 144 o_dqssub, o_dqsauto, o_dqsagg, o_dqsrim, o_dqsmelt, o_dqsfreez, & … … 363 363 ql_seri, qs_seri, qbs_seri, tr_seri, qbs_seri,& 364 364 zphi, u_seri, v_seri, omega, cldfra, & 365 rneb, rnebjn, rneblsvol, zx_rh, zx_rhl, zx_rhi, & 366 pfraclr, pfracld, & 365 rneb, rnebjn, rneblsvol, & 366 zx_rh, zx_rhl, zx_rhi, & 367 pfraclr, pfracld, cldfraliq, sigma2_icefracturb, mean_icefracturb, & 367 368 qraindiag, qsnowdiag, dqreva, dqssub, & 368 369 dqrauto,dqrcol,dqrmelt,dqrfreez, & … … 1999 2000 CALL histwrite_phy(o_pfraclr, pfraclr) 2000 2001 CALL histwrite_phy(o_pfracld, pfracld) 2002 CALL histwrite_phy(o_cldfraliq, cldfraliq) 2003 CALL histwrite_phy(o_sigma2_icefracturb, sigma2_icefracturb) 2004 CALL histwrite_phy(o_mean_icefracturb, mean_icefracturb) 2001 2005 IF (ok_poprecip) THEN 2002 2006 CALL histwrite_phy(o_qrainlsc, qraindiag) -
LMDZ6/trunk/libf/phylmd/physiq_mod.F90
r4976 r5007 332 332 ! 333 333 rneblsvol, & 334 pfraclr, pfracld,&335 distcltop, temp_cltop,&334 pfraclr, pfracld, cldfraliq, sigma2_icefracturb, mean_icefracturb, & 335 distcltop, temp_cltop, & 336 336 zqsatl, zqsats, & 337 337 qclr, qcld, qss, qvc, rnebclr, rnebss, gamma_ss, & … … 1863 1863 & RG,RD,RCPD,RKAPPA,RLVTT,RETV) 1864 1864 CALL ratqs_ini(klon,klev,iflag_thermals,lunout,nbsrf,is_lic,is_ter,RG,RV,RD,RCPD,RLSTT,RLVTT,RTT) 1865 CALL lscp_ini(pdtphys,lunout,prt_level,ok_ice_sursat,iflag_ratqs,fl_cor_ebil,RCPD,RLSTT,RLVTT,RLMLT,RVTMP2,RTT,RD,RG,R PI)1865 CALL lscp_ini(pdtphys,lunout,prt_level,ok_ice_sursat,iflag_ratqs,fl_cor_ebil,RCPD,RLSTT,RLVTT,RLMLT,RVTMP2,RTT,RD,RG,RV,RPI) 1866 1866 CALL blowing_snow_ini(RCPD, RLSTT, RLVTT, RLMLT, & 1867 1867 RVTMP2, RTT,RD,RG, RV, RPI) … … 1943 1943 ELSE IF (klon_glo==1) THEN 1944 1944 pbl_tke(:,:,is_ave) = 0. 1945 pbl_eps(:,:,is_ave) = 0. 1945 1946 DO nsrf=1,nbsrf 1946 1947 DO k = 1,klev+1 1947 1948 pbl_tke(:,k,is_ave) = pbl_tke(:,k,is_ave) & 1948 1949 +pctsrf(:,nsrf)*pbl_tke(:,k,nsrf) 1950 pbl_eps(:,k,is_ave) = pbl_eps(:,k,is_ave) & 1951 +pctsrf(:,nsrf)*pbl_eps(:,k,nsrf) 1949 1952 ENDDO 1950 1953 ENDDO … … 1952 1955 pbl_tke(:,:,is_ave) = 0. !ym missing init : maybe must be initialized in the same way that for klon_glo==1 ?? 1953 1956 !>jyg 1957 pbl_eps(:,:,is_ave) = 0. 1954 1958 ENDIF 1955 1959 !IM begin … … 3896 3900 3897 3901 CALL lscp(klon,klev,phys_tstep,missing_val,paprs,pplay, & 3898 t_seri, q_seri, ptconv,ratqs, &3902 t_seri, q_seri,qs_ancien,ptconv,ratqs, & 3899 3903 d_t_lsc, d_q_lsc, d_ql_lsc, d_qi_lsc, rneb, rneblsvol, rneb_seri, & 3900 pfraclr, pfracld,&3904 pfraclr, pfracld, cldfraliq, sigma2_icefracturb, mean_icefracturb, & 3901 3905 radocond, picefra, rain_lsc, snow_lsc, & 3902 3906 frac_impa, frac_nucl, beta_prec_fisrt, & 3903 3907 prfl, psfl, rhcl, & 3904 3908 zqasc, fraca,ztv,zpspsk,ztla,zthl,iflag_cld_th, & 3905 iflag_ice_thermo, ok_ice_sursat, zqsatl, zqsats, distcltop, temp_cltop, &3906 qclr, qcld, qss, qvc, rnebclr, rnebss, gamma_ss, &3909 iflag_ice_thermo, ok_ice_sursat, zqsatl, zqsats, distcltop, temp_cltop, & 3910 pbl_tke(:,:,is_ave), pbl_eps(:,:,is_ave), qclr, qcld, qss, qvc, rnebclr, rnebss, gamma_ss, & 3907 3911 Tcontr, qcontr, qcontr2, fcontrN, fcontrP , & 3908 3912 cloudth_sth,cloudth_senv,cloudth_sigmath,cloudth_sigmaenv, & -
LMDZ6/trunk/libf/phylmdiso/phys_local_var_mod.F90
r4982 r5007 608 608 REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: pfraclr,pfracld 609 609 !$OMP THREADPRIVATE(pfraclr,pfracld) 610 REAL, SAVE, ALLOCATABLE :: cldfraliq(:,:) 611 !$OMP THREADPRIVATE(cldfraliq) 612 REAL, SAVE, ALLOCATABLE ::mean_icefracturb(:,:) 613 !$OMP THREADPRIVATE(mean_icefracturb) 614 REAL, SAVE, ALLOCATABLE :: sigma2_icefracturb(:,:) 615 !$OMP THREADPRIVATE(sigma2_icefracturb) 610 616 611 617 ! variables de sorties MM … … 1156 1162 ALLOCATE(pfraclr(klon,klev),pfracld(klon,klev)) 1157 1163 pfraclr(:,:)=0. ; pfracld(:,:)=0. ! because not always defined 1164 ALLOCATE(cldfraliq(klon,klev)) 1165 ALLOCATE(sigma2_icefracturb(klon,klev)) 1166 ALLOCATE(mean_icefracturb(klon,klev)) 1158 1167 ALLOCATE(distcltop(klon,klev)) 1159 1168 ALLOCATE(temp_cltop(klon,klev)) … … 1543 1552 DEALLOCATE(rneb) 1544 1553 DEALLOCATE(pfraclr,pfracld) 1554 DEALLOCATE(cldfraliq) 1555 DEALLOCATE(sigma2_icefracturb) 1556 DEALLOCATE(mean_icefracturb) 1545 1557 DEALLOCATE (zxsnow,snowhgt,qsnow,to_ice,sissnow,runoff,albsol3_lic) 1546 1558 DEALLOCATE(distcltop) -
LMDZ6/trunk/libf/phylmdiso/phys_output_ctrlout_mod.F90
r4889 r5007 1553 1553 TYPE(ctrl_out), SAVE :: o_rneb = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11, 11/), & 1554 1554 'rneb', 'Cloud fraction', '-', (/ ('', i=1, 10) /)) 1555 TYPE(ctrl_out), SAVE :: o_cldfraliq = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), & 1556 'cldfraliq', 'Liquid fraction of the cloud', '-', (/ ('', i=1, 10) /)) 1557 TYPE(ctrl_out), SAVE :: o_sigma2_icefracturb = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), & 1558 'sigma2_icefracturb', 'Variance of the diagnostic supersaturation distribution (icefrac_turb) [-]', '-', (/ ('', i=1, 10) /)) 1559 TYPE(ctrl_out), SAVE :: o_mean_icefracturb = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), & 1560 'mean_icefracturb', 'Mean of the diagnostic supersaturation distribution (icefrac_turb) [-]', '-', (/ ('', i=1, 10) /)) 1561 1555 1562 TYPE(ctrl_out), SAVE :: o_rnebjn = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11,11, 11/), & 1556 1563 'rnebjn', 'Cloud fraction in day', '-', (/ ('', i=1, 10) /)) -
LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90
r4983 r5007 373 373 ! 374 374 rneblsvol, & 375 pfraclr, pfracld,&376 distcltop, temp_cltop,&375 pfraclr, pfracld, cldfraliq, sigma2_icefracturb, mean_icefracturb, & 376 distcltop, temp_cltop, & 377 377 zqsatl, zqsats, & 378 378 qclr, qcld, qss, qvc, rnebclr, rnebss, gamma_ss, & … … 2026 2026 & RG,RD,RCPD,RKAPPA,RLVTT,RETV) 2027 2027 CALL ratqs_ini(klon,klev,iflag_thermals,lunout,nbsrf,is_lic,is_ter,RG,RV,RD,RCPD,RLSTT,RLVTT,RTT) 2028 CALL lscp_ini(pdtphys,lunout,prt_level,ok_ice_sursat,iflag_ratqs,fl_cor_ebil,RCPD,RLSTT,RLVTT,RLMLT,RVTMP2,RTT,RD,RG,R PI)2028 CALL lscp_ini(pdtphys,lunout,prt_level,ok_ice_sursat,iflag_ratqs,fl_cor_ebil,RCPD,RLSTT,RLVTT,RLMLT,RVTMP2,RTT,RD,RG,RV,RPI) 2029 2029 CALL blowing_snow_ini(RCPD, RLSTT, RLVTT, RLMLT, & 2030 2030 RVTMP2, RTT,RD,RG, RV, RPI) … … 5076 5076 5077 5077 CALL lscp(klon,klev,phys_tstep,missing_val,paprs,pplay, & 5078 t_seri, q_seri, ptconv,ratqs, &5078 t_seri, q_seri,qs_ancien,ptconv,ratqs, & 5079 5079 d_t_lsc, d_q_lsc, d_ql_lsc, d_qi_lsc, rneb, rneblsvol, rneb_seri, & 5080 pfraclr, pfracld,&5080 pfraclr, pfracld, cldfraliq, sigma2_icefracturb, mean_icefracturb, & 5081 5081 radocond, picefra, rain_lsc, snow_lsc, & 5082 5082 frac_impa, frac_nucl, beta_prec_fisrt, & 5083 5083 prfl, psfl, rhcl, & 5084 5084 zqasc, fraca,ztv,zpspsk,ztla,zthl,iflag_cld_th, & 5085 iflag_ice_thermo, ok_ice_sursat, zqsatl, zqsats, distcltop, temp_cltop, &5086 qclr, qcld, qss, qvc, rnebclr, rnebss, gamma_ss, &5085 iflag_ice_thermo, ok_ice_sursat, zqsatl, zqsats, distcltop, temp_cltop, & 5086 pbl_tke(:,:,is_ave), pbl_eps(:,:,is_ave), qclr, qcld, qss, qvc, rnebclr, rnebss, gamma_ss, & 5087 5087 Tcontr, qcontr, qcontr2, fcontrN, fcontrP , & 5088 5088 cloudth_sth,cloudth_senv,cloudth_sigmath,cloudth_sigmaenv, &
Note: See TracChangeset
for help on using the changeset viewer.
