Changeset 4639 for LMDZ6/trunk/libf/phylmdiso

Timestamp:

Jul 24, 2023, 4:46:18 PM (18 months ago)

Author:

evignon

Message:

modification de Lea pour inclure une dependance de la phase a la temperature du sommet du nuage

Location:

LMDZ6/trunk/libf/phylmdiso

Files:

• phys_local_var_mod.F90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmdiso/phys_local_var_mod.F90

@@ -599 +599 @@
       REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: distcltop
 !$OMP THREADPRIVATE(distcltop)
+      REAL, SAVE, ALLOCATABLE :: temp_cltop(:,:) 
+!$OMP THREADPRIVATE(temp_cltop)
 
 ! variables de sorties MM
@@ -1090 +1092 @@
       pfraclr(:,:)=0. ; pfracld(:,:)=0. ! because not always defined
       ALLOCATE(distcltop(klon,klev))
+      ALLOCATE(temp_cltop(klon,klev))
 
       ALLOCATE (zxsnow(klon),snowhgt(klon),qsnow(klon),to_ice(klon))
@@ -1464 +1467 @@
       DEALLOCATE(pfraclr,pfracld)
       DEALLOCATE(distcltop)
+      DEALLOCATE(temp_cltop)
       DEALLOCATE (zxsnow,snowhgt,qsnow,to_ice,sissnow,runoff,albsol3_lic)
 #ifdef ISO

LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

@@ -334 +334 @@
        rneblsvol, &
        pfraclr,pfracld, &
-       distcltop, &
+       distcltop,temp_cltop, &
        zqsatl, zqsats, &
        qclr, qcld, qss, qvc, rnebclr, rnebss, gamma_ss, &
@@ -4806 +4806 @@
          prfl, psfl, rhcl,  &
          zqasc, fraca,ztv,zpspsk,ztla,zthl,iflag_cld_th, &
-         iflag_ice_thermo, ok_ice_sursat, zqsatl, zqsats, distcltop, &
+         iflag_ice_thermo, ok_ice_sursat, zqsatl, zqsats, distcltop, temp_cltop, &
          qclr, qcld, qss, qvc, rnebclr, rnebss, gamma_ss, &
          Tcontr, qcontr, qcontr2, fcontrN, fcontrP )
@@ -5574 +5574 @@
                flwp, fiwp, flwc, fiwc, &
                mass_solu_aero, mass_solu_aero_pi, &
-               cldtaupi, latitude_deg, distcltop, re, fl, ref_liq, ref_ice, &
+               cldtaupi, latitude_deg, distcltop, temp_cltop, re, fl, ref_liq, ref_ice, &
                ref_liq_pi, ref_ice_pi)
        ELSE
@@ -5582 +5582 @@
                ok_aie, &
                mass_solu_aero, mass_solu_aero_pi, &
-               bl95_b0, bl95_b1, distcltop, &
+               bl95_b0, bl95_b1, distcltop, temp_cltop, &
                cldtaupi, re, fl)
        ENDIF