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

Timestamp:

Jan 29, 2026, 12:52:39 AM (2 months ago)

Author:

fhourdin

Message:

Travail pour la replayisation de la convection

Reunion de tous les anciens common devenus modules, dans lmdz_cv_ini.
Déplacement de presque toutes les routines d'initialisation dans lmdz_cv_ini.
Encapsulage de certains sous-programmes dans des modules.
Suppression de programmes inutilisés (cv3_crit et cv3_incp)
Reste :

• à sortir des routines d'initialisation "_pre" de cv_driver et

cva_driver

• à passer le variables argunement en intent(in/out/inout).

La convergence numérique a été testée pour
iflag_con=3/30/4
en 3D parallèle.
La compilation de la version isotopique fonctionne.

File:

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

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -102 +102 @@
     USE lmdz_reevap_ini, ONLY : reevap_ini
     USE lmdz_wake_ini, ONLY : wake_ini
+    USE lmdz_cv_ini, ONLY : cv_ini
+    USE lmdz_cv_ini, ONLY : epmax, coef_epmax_cape, cvl_comp_threshold, cvl_sig2feed, iflag_cvl_sigd, iflag_clw, ok_adj_ema
+
     USE lmdz_surf_wind_ini, ONLY : surf_wind_ini
     USE lmdz_surf_wind, ONLY : surf_wind
@@ -276 +279 @@
          JrNt,                             &
          dthmin, evap, snowerosion, icesub_ice, icemelt_ice, fder, plcl, plfc,   &
+
          prw, prlw, prsw, prbsw, water_budget,         &
          s_lcl, s_pblh, s_pblt, s_therm,   &
@@ -398 +402 @@
     USE yomcst_mod_h
     USE clesphys_mod_h
-    USE conema3_mod_h
+! Replayisation USE conema3_mod_h
     USE alpale_mod
     USE yoethf_mod_h
@@ -656 +660 @@
     !
     !    INTEGER izero
-    INTEGER k_upper_cv
+    INTEGER, SAVE :: k_upper_cv
+    !$OMP THREADPRIVATE(k_upper_cv)
+
     !------------------------------------------------------------------
     ! Compteur de l'occurence de cvpas=1
@@ -1388 +1394 @@
             read_climoz, &
             alp_offset)
+
        CALL init_etat0_limit_unstruct
        IF (.NOT. create_etat0_limit) CALL init_limit_read(days_elapsed)
@@ -1866 +1873 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
        CALL reevap_ini(fl_cor_ebil, iflag_ice_thermo, RCPD, RLSTT, RLVTT, RVTMP2)
+             !jyg<
+             ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+             ! Calculate the upmost level of deep convection loops: k_upper_cv
+             !  (near 22 km)
+             k_upper_cv = klev
+             !izero = klon/2+1/klon
+             !DO k = klev,1,-1
+             !   IF (pphi(izero,k) > 22.e4) k_upper_cv = k
+             !ENDDO
+             ! FH : nouveau calcul base sur un profil global sans quoi
+             ! le modele etait sensible au decoupage de domaines
+             DO k = klev,1,-1
+                IF (-7*log(presnivs(k)/presnivs(1)) > 25.) k_upper_cv = k
+             ENDDO
+             IF (prt_level .ge. 5) THEN
+                Print *, 'upmost level of deep convection loops: k_upper_cv = ', &
+                     k_upper_cv
+             ENDIF
+       CALL cv_ini(iflag_con,iflag_ice_thermo,rcpd,rcpv,rcw,rv,rcs,rd,rlvtt,rlstt,rg,retv,rtt,r,klev,k_upper_cv,phys_tstep,fl_cor_ebil,prt_level,lunout)
+
+
        CALL wake_ini(iflag_wake,rg,rd,rv,prt_level)
        CALL yamada_ini(klon,lunout,prt_level)
@@ -3268 +3296 @@
           ! (driver commun aux versions 3 et 4)
           !
+    print*,' lmdz_cv_ini, ONLY 4 :', epmax, coef_epmax_cape, cvl_comp_threshold, cvl_sig2feed, iflag_cvl_sigd, iflag_clw, ok_adj_ema
           IF (ok_cvl) THEN ! new driver for convectL
              !
-             !jyg<
-             ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-             ! Calculate the upmost level of deep convection loops: k_upper_cv
-             !  (near 22 km)
-             k_upper_cv = klev
-             !izero = klon/2+1/klon
-             !DO k = klev,1,-1
-             !   IF (pphi(izero,k) > 22.e4) k_upper_cv = k
-             !ENDDO
-             ! FH : nouveau calcul base sur un profil global sans quoi
-             ! le modele etait sensible au decoupage de domaines
-             DO k = klev,1,-1
-                IF (-7*log(presnivs(k)/presnivs(1)) > 25.) k_upper_cv = k
-             ENDDO
-             IF (prt_level .ge. 5) THEN
-                Print *, 'upmost level of deep convection loops: k_upper_cv = ', &
-                     k_upper_cv
-             ENDIF
              !
              !>jyg