Changeset 4869 for LMDZ6/trunk/libf/phylmd/lmdz_lscp.F90

Timestamp:

Mar 22, 2024, 2:31:51 PM (2 months ago)

Author:

evignon

Message:

identification d'un bug dans l'ajustement de la temperature
lors de la fonte dans lscp
Ajout d'un flag pour maintenir/enlever le probleme

File:

• LMDZ6/trunk/libf/phylmd/lmdz_lscp.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_lscp.F90

@@ -110 +110 @@
 USE lmdz_lscp_ini, ONLY : iflag_autoconversion, ffallv_con, ffallv_lsc
 USE lmdz_lscp_ini, ONLY : RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT, RD, RG
-USE lmdz_lscp_ini, ONLY : ok_poprecip
+USE lmdz_lscp_ini, ONLY : ok_poprecip, ok_bug_fonte_lscp
 
 IMPLICIT NONE
@@ -630 +630 @@
                     zrflcld(i)=zrflcld(i)+zmelt*ziflcld(i)
                     zrfl(i)=zrflclr(i)+zrflcld(i)
-
+                    IF (ok_bug_fonte_lscp) THEN
                     ziflclr(i)=ziflclr(i)*(1.-zmelt)
                     ziflcld(i)=ziflcld(i)*(1.-zmelt)
                     zifl(i)=ziflclr(i)+ziflcld(i)
-
+                    ENDIF
                 ELSE
                     zrfl(i)=zrfl(i)+zmelt*zifl(i)
-
+                    IF (ok_bug_fonte_lscp) THEN
                     zifl(i)=zifl(i)*(1.-zmelt)
+                    ENDIF
                 ENDIF
 
@@ -648 +649 @@
                 zt(i)=zt(i)-zifl(i)*zmelt*(RG*dtime)/(paprs(i,k)-paprs(i,k+1)) &
                 *RLMLT/RCPD/(1.0+RVTMP2*(zq(i)+zmqc(i)))
-                
+               
+                IF (.NOT. ok_bug_fonte_lscp) THEN
+                IF (iflag_evap_prec.GE.4) THEN
+                    ziflclr(i)=ziflclr(i)*(1.-zmelt)
+                    ziflcld(i)=ziflcld(i)*(1.-zmelt)
+                    zifl(i)=ziflclr(i)+ziflcld(i)
+                ELSE
+                    zifl(i)=zifl(i)*(1.-zmelt)
+                ENDIF
+                ENDIF
 
             ELSE