Changeset 2476 for LMDZ5/trunk

Timestamp:

Mar 23, 2016, 10:48:19 PM (9 years ago)

Author:

oboucher

Message:

Introduction of a mass fixer for water in physiq
int q+ql dz at the end = int q+ql dz at the beginning + evap*pdtphys - precip*pdtphys

File:

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

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -937 +937 @@
     !albedo SB <<<
 
+    !--OB variables for mass fixer 
+    logical, parameter :: mass_fixer=.true.
+    real qql1(klon),qql2(klon),zdz,corrqql
+
     ! Ehouarn: set value of jjmp1 since it is no longer a "fixed parameter"
     jjmp1=nbp_lat
@@ -1588 +1592 @@
        ENDDO
     ENDDO
+    !
+    !--OB mass fixer 
+    IF (mass_fixer) THEN
+    !--store initial water burden
+    qql1(:)=0.0
+    DO i = 1, klon
+      DO k = 1, klev
+        zdz=(paprs(i,k)-paprs(i,k+1))/RG
+        qql1(i)=qql1(i)+(q_seri(i,k)+ql_seri(i,k))*zdz
+      ENDDO
+    ENDDO
+    ENDIF
+    !--fin mass fixer
+
     tke0(:,:)=pbl_tke(:,:,is_ave)
     !CR:Nombre de traceurs de l'eau: nqo
@@ -4528 +4546 @@
     ENDIF
 
+    !--OB mass fixer 
+    !--profile is corrected to force mass conservation of water
+    IF (mass_fixer) THEN
+    qql2(:)=0.0
+    DO i = 1, klon
+      DO k = 1, klev
+        zdz=(paprs(i,k)-paprs(i,k+1))/RG
+        qql2(i)=qql2(i)+(q_seri(i,k)+ql_seri(i,k))*zdz
+      ENDDO
+    ENDDO
+    DO i = 1, klon
+      !--compute ratio of what q+ql should be with conservation to what it is
+      corrqql=(qql1(i)+(evap(i)-rain_fall(i)-snow_fall(i))*pdtphys)/qql2(i)
+      DO k = 1, klev
+        q_seri(i,k) =q_seri(i,k)*corrqql
+        ql_seri(i,k)=ql_seri(i,k)*corrqql
+      ENDDO
+    ENDDO
+    ENDIF
+    !--fin mass fixer
 
     ! 22.03.04 END