Changeset 5189 for LMDZ6/trunk/libf/phylmdiso

Timestamp:

Sep 14, 2024, 2:33:38 PM (3 months 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/phylmdiso/physiq_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

@@ -1368 +1368 @@
     !albedo SB <<<
 
+    !--Lea Raillard qs_ini 
+    REAL, dimension(klon,klev) :: qs_ini
+
     !--OB variables for mass fixer (hard coded for now)
     LOGICAL, PARAMETER :: mass_fixer=.FALSE.
@@ -2620 +2623 @@
        ENDDO
     ENDDO
+    
+    ! Lea Raillard qs_ini for cloud phase param.
+    qs_ini(:,:)=qs_seri(:,:)
 
     ! C Risi: dispatcher les isotopes dans les xt_seri
@@ -5062 +5068 @@
 
     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,  &