Changeset 4869

Timestamp:

Mar 22, 2024, 2:31:51 PM (7 weeks 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

Location:

LMDZ6/trunk/libf/phylmd

Files:

• lmdz_lscp.F90 (modified) (3 diffs)
• lmdz_lscp_ini.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

LMDZ6/trunk/libf/phylmd/lmdz_lscp_ini.F90

@@ -8 +8 @@
   REAL RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT, RD, RG, RPI
   !$OMP THREADPRIVATE(RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT, RD, RG, RPI)
+  
+  LOGICAL, SAVE, PROTECTED :: ok_bug_fonte_lscp=.true. ! maintains bug on temperature adjustment associated with snow melting
+  !$OMP THREADPRIVATE(ok_bug_fonte_lscp)
 
   REAL, SAVE, PROTECTED :: seuil_neb=0.001      ! cloud fraction threshold: a cloud really exists when exceeded
@@ -317 +320 @@
     CALL getin_p('snow_fallspeed_clr',snow_fallspeed_clr)
     CALL getin_p('snow_fallspeed_cld',snow_fallspeed_cld)
+    CALL getin_p('ok_bug_fonte_lscp',ok_bug_fonte_lscp)
 
 
@@ -367 +371 @@
     WRITE(lunout,*) 'lscp_ini, snow_fallspeed_clr:', snow_fallspeed_clr
     WRITE(lunout,*) 'lscp_ini, snow_fallspeed_cld:', snow_fallspeed_cld
+    WRITE(lunout,*) 'lscp_ini, ok_bug_fonte_lscp:', ok_bug_fonte_lscp