Changeset 4886

Timestamp:

Apr 1, 2024, 6:09:01 PM (5 weeks ago)

Author:

oboucher

Message:

Correcting a couple of bugs in the mass fixer for water vapour in the physics (the mass_fixer is switched to FALSE by default)

File:

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

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -2450 +2450 @@
              qs_seri(i,k) = qx(i,k,isol)
              IF (ok_ice_sursat) THEN
-             rneb_seri(i,k) = qx(i,k,irneb)
+               rneb_seri(i,k) = qx(i,k,irneb)
              ENDIF
              IF (ok_bs) THEN
-             qbs_seri(i,k)= qx(i,k,ibs)
+               qbs_seri(i,k)= qx(i,k,ibs)
              ENDIF
-
           ENDIF
-
-
        ENDDO
     ENDDO
@@ -5588 +5585 @@
     qql2(:)=0.0
     DO k = 1, klev
-      qql2(:)=qql2(:)+(q_seri(:,k)+ql_seri(:,k)+qs_seri(:,k))*zmasse(:,k)
+      qql2(:)=qql2(:)+(q_seri(:,k)+ql_seri(:,k)+qs_seri(:,k)+qbs_seri(:,k))*zmasse(:,k)
     ENDDO
 
@@ -5610 +5607 @@
         q_seri(i,k) =q_seri(i,k)*corrqql
         ql_seri(i,k)=ql_seri(i,k)*corrqql
+        qs_seri(i,k)=qs_seri(i,k)*corrqql
+        qbs_seri(i,k)=qbs_seri(i,k)*corrqql
       ENDDO
     ENDDO