Changeset 4727 for LMDZ6/branches/LMDZ_ECRad/libf/phylmd/nuage.F90
- Timestamp:
- Oct 19, 2023, 4:02:57 PM (12 months ago)
- Location:
- LMDZ6/branches/LMDZ_ECRad
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/LMDZ_ECRad
- Property svn:mergeinfo changed
-
LMDZ6/branches/LMDZ_ECRad/libf/phylmd/nuage.F90
r3999 r4727 2 2 3 3 SUBROUTINE nuage(paprs, pplay, t, pqlwp,picefra, pclc, pcltau, pclemi, pch, pcl, pcm, & 4 pct, pctlwp, ok_aie, mass_solu_aero, mass_solu_aero_pi, bl95_b0, bl95_b1, &5 cldtaupi, re, fl)4 pct, pctlwp, ok_aie, mass_solu_aero, mass_solu_aero_pi, bl95_b0, bl95_b1, distcltop, & 5 temp_cltop, cldtaupi, re, fl) 6 6 USE dimphy 7 USE l scp_tools_mod, only: icefrac_lscp7 USE lmdz_lscp_tools, only: icefrac_lscp 8 8 USE icefrac_lsc_mod ! computes ice fraction (JBM 3/14) 9 USE lmdz_lscp_ini, only : iflag_t_glace 9 10 USE phys_local_var_mod, ONLY: ptconv 10 11 IMPLICIT NONE … … 49 50 50 51 REAL pct(klon), pctlwp(klon), pch(klon), pcl(klon), pcm(klon) 51 52 REAL distcltop(klon,klev) 53 REAL temp_cltop(klon,klev) 52 54 LOGICAL lo 53 55 … … 112 114 ! t_glace_max, exposant_glace) 113 115 IF (ok_new_lscp) THEN 114 CALL icefrac_lscp(klon,t(:,k), pplay(:,k)/paprs(:,1),zfice(:),dzfice(:))116 CALL icefrac_lscp(klon,t(:,k),iflag_ice_thermo,distcltop(:,k),temp_cltop(:,k),zfice(:),dzfice(:)) 115 117 ELSE 116 118 CALL icefrac_lsc(klon,t(:,k),pplay(:,k)/paprs(:,1),zfice(:)) … … 118 120 ENDIF 119 121 120 IF ((.NOT. ptconv(i,k)) .AND.ok_new_lscp .AND. ok_icefra_lscp) THEN121 122 123 124 125 zfice(i)=picefra(i,k)126 ENDDO122 IF (ok_new_lscp .AND. ok_icefra_lscp) THEN 123 ! EV: take the ice fraction directly from the lscp code 124 ! consistent only for non convective grid points 125 ! critical for mixed phase clouds 126 DO i=1,klon 127 IF (.NOT. ptconv(i,k)) THEN 128 zfice(i)=picefra(i,k) 127 129 ENDIF 130 ENDDO 131 ENDIF 128 132 129 133
Note: See TracChangeset
for help on using the changeset viewer.