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

Timestamp:

Mar 11, 2026, 7:59:52 PM (4 weeks ago)

Author:

fhourdin

Message:

Debut separation init/compute dans physiq_mod

File:

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

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -4 +4 @@
 !#define IO_DEBUG
 MODULE physiq_mod
+
+! FH : commentaire a detruire avant mai 2026
+! ==========================================
+! 11 mars :
+! ---------
+!   remplacement de if (first) par debut
+!   élimination de la variable first
+!   Supression d'un print entre endif et if du meme debut
+! Deplacement des lignes 
+!    pmfd(:,:) = 0. 
+!    pen_u(:,:) = 0. 
+!    pen_d(:,:) = 0. 
+!    pde_d(:,:) = 0. 
+!    pde_u(:,:) = 0. 
+!    aam=0.
+!    d_t_adjwk(:,:)=0
+!    d_q_adjwk(:,:)=0
+!    alp_bl_conv(:)=0.
+!    torsfc=0.
+! APRES
+!   ****************     Fin  de   IF ( debut  )   ***************
+! Deplucation de 
+!       forall (k=1: nbp_lev) zmasse(:, k) = (paprs(:, k)-paprs(:, k+1)) / rg
+!     avant et apres la fin de debut
+
+
+
 
   IMPLICIT NONE
@@ -1075 +1102 @@
     ! Declaration des constantes et des fonctions thermodynamiques
     !
-    LOGICAL,SAVE :: first=.TRUE.
-    !$OMP THREADPRIVATE(first)
-
     ! VARIABLES RELATED TO OZONE CLIMATOLOGIES ; all are OpenMP shared
     ! Note that pressure vectors are in Pa and in stricly ascending order
@@ -1197 +1221 @@
     !======================================================================!
     if (debut) then                                                        !
+       print*,'DEBUT PHYSIQ_MOD TRAVAIL DEBUT 1'
        iflag_physiq=0
        call getin_p('iflag_physiq', iflag_physiq)                          !
@@ -1266 +1291 @@
          "physiq_mod paprs bad order")
 
-    IF (first) THEN
+    IF (debut) THEN
 
        CALL s2s_initialize     ! initialization of source to source tools
@@ -1366 +1391 @@
 
        itau_con=0
-       first=.FALSE.
-
-    ENDIF  ! first
-
-    !ym => necessaire pour iflag_con != 2
-    pmfd(:,:) = 0.
-    pen_u(:,:) = 0.
-    pen_d(:,:) = 0.
-    pde_d(:,:) = 0.
-    pde_u(:,:) = 0.
-    d_t_adjwk(:,:)=0
-    d_q_adjwk(:,:)=0
-
-    alp_bl_conv(:)=0.
-
-    forall (k=1: nbp_lev) zmasse(:, k) = (paprs(:, k)-paprs(:, k+1)) / rg
-
-
-    IF (debut) THEN
+
+
+       forall (k=1: nbp_lev) zmasse(:, k) = (paprs(:, k)-paprs(:, k+1)) / rg
+
+
        CALL suphel ! initialiser constantes et parametres phys.
        ! tau_gl : constante de rappel de la temperature a la surface de la glace - en
@@ -1465 +1476 @@
 #endif
 
-    ENDIF
-
-    IF (prt_level.ge.1) print *,'CONVERGENCE PHYSIQUE THERM 1 '
-
-    !======================================================================
-    ! Gestion calendrier : mise a jour du module phys_cal_mod
-    !
-    !     CALL phys_cal_update(jD_cur,jH_cur)
-
-    !
-    ! Si c'est le debut, il faut initialiser plusieurs choses
-    !          ********
-    !
-    IF (debut) THEN
        !rv CRinitialisation de wght_th et lalim_conv pour la
        !definition de la couche alimentation de la convection a partir
@@ -2245 +2242 @@
     ! Incrementer le compteur de la physique
     !
+    !ym => necessaire pour iflag_con != 2
+    pmfd(:,:) = 0.
+    pen_u(:,:) = 0.
+    pen_d(:,:) = 0.
+    pde_d(:,:) = 0.
+    pde_u(:,:) = 0.
+    d_t_adjwk(:,:)=0
+    d_q_adjwk(:,:)=0
+    alp_bl_conv(:)=0.
+
+    forall (k=1: nbp_lev) zmasse(:, k) = (paprs(:, k)-paprs(:, k+1)) / rg
+
     itap   = itap + 1
     IF (is_master .OR. prt_level > 9) THEN