Changeset 5189 for LMDZ6/trunk/libf/phylmd/physiq_mod.F90

Timestamp:

Sep 14, 2024, 2:33:38 PM (8 days ago)

Author:

evignon

Message:

correction in new cloud phase param. Ice feedback is now computed from ice concentration at the beginning of the physics (qs_ini)
Lea Raillard

File:

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

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -1263 +1263 @@
     !albedo SB <<<
 
+    !--Lea Raillard qs_ini 
+    REAL, dimension(klon,klev) :: qs_ini
+
     !--OB variables for mass fixer (hard coded for now)
     REAL qql1(klon),qql2(klon),corrqql
@@ -2449 +2452 @@
        ENDDO
     ENDDO
+    ! Lea Raillard qs_ini for cloud phase param.
+    qs_ini(:,:)=qs_seri(:,:)
     !
     !--OB water mass fixer 
@@ -3829 +3834 @@
 
     CALL lscp(klon,klev,phys_tstep,missing_val,paprs,pplay, &
-         t_seri, q_seri,qs_ancien,ptconv,ratqs, &
+         t_seri, q_seri,qs_ini,ptconv,ratqs, &
          d_t_lsc, d_q_lsc, d_ql_lsc, d_qi_lsc, rneb, rneblsvol, rneb_seri, & 
          pfraclr, pfracld, cldfraliq, sigma2_icefracturb, mean_icefracturb,  &