Changeset 6132 for LMDZ6/trunk/libf

Timestamp:

Mar 26, 2026, 9:09:01 PM (6 days ago)

Author:

evignon

Message:

suppression de la param hines des sources, devenue totalement obsolete (apres verification avec Francois Lott)

Location:

LMDZ6/trunk/libf

Files:

• phylmd/ener_conserv_mod.f90 (modified) (2 diffs)
• phylmd/lmdz_call_gwd.F90 (modified) (6 diffs)
• phylmd/lmdz_gwd_hines.f90 (deleted)
• phylmd/lmdz_gwd_ini.f90 (modified) (4 diffs)
• phylmd/phyetat0_mod.f90 (modified) (2 diffs)
• phylmd/phys_local_var_mod.F90 (modified) (4 diffs)
• phylmd/phys_output_ctrlout_mod.F90 (modified) (3 diffs)
• phylmd/phys_output_var_mod.f90 (modified) (3 diffs)
• phylmd/phys_output_write_mod.F90 (modified) (6 diffs)
• phylmd/physiq_mod.F90 (modified) (2 diffs)
• phylmdiso/lmdz_gwd_hines.f90 (deleted)
• phylmdiso/phyetat0_mod.F90 (modified) (2 diffs)
• phylmdiso/phyredem.F90 (modified) (2 diffs)
• phylmdiso/physiq_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/ener_conserv_mod.f90

@@ -36 +36 @@
 USE phys_local_var_mod, ONLY : d_t_eva,d_t_lsc,d_q_eva,d_q_lsc
 USE phys_local_var_mod, ONLY : d_u_oro,d_v_oro,d_u_lif,d_v_lif
-USE phys_local_var_mod, ONLY : du_gwd_hines,dv_gwd_hines,dv_gwd_front,dv_gwd_rando
+USE phys_local_var_mod, ONLY : dv_gwd_front,dv_gwd_rando
 USE phys_state_var_mod, ONLY : du_gwd_front,du_gwd_rando
 USE phys_output_var_mod, ONLY : bils_ec,bils_ech,bils_tke,bils_kinetic,bils_enthalp,bils_latent,bils_diss
 USE add_phys_tend_mod, ONLY : fl_cor_ebil
 USE infotrac_phy, ONLY: nqtot
-USE lmdz_gwd_ini, ONLY : ok_gwd_rando, ok_hines
+USE lmdz_gwd_ini, ONLY : ok_gwd_rando
 
 USE yomcst_mod_h
@@ -152 +152 @@
       d_u(:,:)=d_u_vdf(:,:)+d_u_ajs(:,:)+d_u_con(:,:)+d_u_oro(:,:)+d_u_lif(:,:)
       d_v(:,:)=d_v_vdf(:,:)+d_v_ajs(:,:)+d_v_con(:,:)+d_v_oro(:,:)+d_v_lif(:,:) 
-      IF (ok_hines) THEN
-         d_u_vdf(:,:)=d_u_vdf(:,:)+du_gwd_hines(:,:)
-         d_v_vdf(:,:)=d_v_vdf(:,:)+dv_gwd_hines(:,:)
-      ENDIF
-      IF (.not. ok_hines .and. ok_gwd_rando) THEN
+      IF (ok_gwd_rando) THEN
          d_u_vdf(:,:)=d_u_vdf(:,:)+du_gwd_front(:,:)
          d_v_vdf(:,:)=d_v_vdf(:,:)+dv_gwd_front(:,:)

LMDZ6/trunk/libf/phylmd/lmdz_call_gwd.F90

@@ -26 +26 @@
                        d_u_oro, d_v_oro, d_t_oro, zustr_oro, zvstr_oro, &
                        d_u_lif, d_v_lif, d_t_lif, zustr_lif, zvstr_lif, &
-                       d_u_hines, d_v_hines, d_t_hines, zustr_hines, zvstr_hines, &
                        d_u_front, d_v_front, zustr_front, zvstr_front, &
                        d_u_precip, d_v_precip, zustr_precip, zvstr_precip, &
                        east_gwstress, west_gwstress, aam, torsfc)
 
-      USE lmdz_gwd_ini, ONLY: ok_strato, ok_orodr, ok_orolf, ok_hines, ok_gwd_rando
+      USE lmdz_gwd_ini, ONLY: ok_strato, ok_orodr, ok_orolf, ok_gwd_rando
       USE lmdz_gwd_ini, ONLY: addtkeoro, smallscales_tkeoro, alphatkeoro
       USE lmdz_gwd_ini, ONLY: zpmm_orodr_t, zstd_orodr_t, zpmm_orolf_t, nm_oro_t
@@ -38 +37 @@
       USE lmdz_gwd_ogwd, ONLY: drag_noro_strato, lift_noro_strato
       USE lmdz_gwd_ogwd_old, ONLY: drag_noro, lift_noro
-      USE lmdz_gwd_hines, ONLY: hines_gwd
       USE lmdz_gwd_front, ONLY: acama_gwd_rando
       USE lmdz_gwd_precip, ONLY: flott_gwd_rando
@@ -111 +109 @@
       REAL, DIMENSION(klon, klev), INTENT(OUT)  :: d_t_lif ! t increment due to sso lift [K]
 
-      REAL, DIMENSION(klon, klev), INTENT(OUT)  :: d_u_hines ! u increment due to Hines param [m/s]
-      REAL, DIMENSION(klon, klev), INTENT(OUT)  :: d_v_hines ! v increment due to Hines param [m/s]
-      REAL, DIMENSION(klon, klev), INTENT(OUT)  :: d_t_hines ! t increment due to Hines param [K]
-
       REAL, DIMENSION(klon), INTENT(OUT)  :: zustr_oro   ! vertically integrated zonal stress due to sso drag [kg.m/s/m2]
       REAL, DIMENSION(klon), INTENT(OUT)  :: zvstr_oro   ! vertically integrated meridional stress due to sso drag [kg.m/s/m2]
@@ -120 +114 @@
       REAL, DIMENSION(klon), INTENT(OUT)  :: zustr_lif   ! vertically integrated zonal stress due to sso lift [kg.m/s/m2]
       REAL, DIMENSION(klon), INTENT(OUT)  :: zvstr_lif   ! vertically integrated meridional stress due to sso lift [kg.m/s/m2]
-
-      REAL, DIMENSION(klon), INTENT(OUT)  :: zustr_hines   ! vertically integrated zonal stress due to Hines param [kg.m/s/m2]
-      REAL, DIMENSION(klon), INTENT(OUT)  :: zvstr_hines   ! vertically integrated meridional stress due to Hines param [kg.m/s/m2]
 
       REAL, DIMENSION(klon), INTENT(OUT)  :: zustr_front   ! vertically integrated zonal stress due to fronts gwd param [kg.m/s/m2]
@@ -168 +159 @@
       d_t_lif(:, :) = 0.
 
-      d_u_hines(:, :) = 0.
-      d_v_hines(:, :) = 0.
-      d_t_hines(:, :) = 0.
-
       ! DO NOT set tendencies associated with front and precip to 0
       ! as they are inout variables because
       ! the parameterizations considers some "process memory"
-
-      zustr_hines(:) = 0.
-      zvstr_hines(:) = 0.
 
       zustr_front(:) = 0.
@@ -322 +306 @@
 !--------------------------------------
 
-      ! old approach following Hines 1997
-      IF (ok_hines) THEN
-         ! this routine is now out of date. Do not use it except for specific purposes
-         CALL hines_gwd(klon, klev, phys_tstep, paprs, pplay, latitude_deg, t_seri, &
-                        u_seri, v_seri, zustr_hines, zvstr_hines, d_t_hines, &
-                        d_u_hines, d_v_hines)
-
-         !  Add tendencies
-         !-----------------------------------------------------------------------
-
-#ifdef ISO
-         STOP 'appel a add_phys_tend non traite dans call_gwd'
-#else
-         CALL add_phys_tend(d_u_hines, d_v_hines, d_t_hines, dq0, dql0, &
-                            dqi0, dqbs0, paprs, 'hin', abortphy, flag_inhib_tend, itap, 0)
-#endif
-
-         CALL prt_enerbil('hin', itap)
-
-         !  Local diagnostics
-         !-----------------------------------------------------------------------
-         DO k = 1, klev
-            zustr_hines(:) = zustr_hines(:) + d_u_hines(:, k)/phys_tstep &
-                             *(paprs(:, k) - paprs(:, k + 1))/rg
-            zvstr_hines(:) = zvstr_hines(:) + d_v_hines(:, k)/phys_tstep &
-                             *(paprs(:, k) - paprs(:, k + 1))/rg
-         END DO
-
-         east_gwstress(:, :) = 0.
-         west_gwstress(:, :) = 0.
-
-      END IF
-
       ! stochastic approach following De la Camara & Lott 2015
-      IF (.NOT. ok_hines .AND. ok_gwd_rando) THEN
+      IF (ok_gwd_rando) THEN
 
          CALL acama_gwd_rando(klon, klev, phys_tstep, pplay, presnivs, latitude_deg, t_seri, u_seri, &

LMDZ6/trunk/libf/phylmd/lmdz_gwd_ini.f90

@@ -78 +78 @@
 
    ! for non-orographic gravity wave drag
-   LOGICAL, SAVE, PROTECTED :: ok_hines = .false. ! activate non-orographic gravity wave drag from Hines
-   !$OMP THREADPRIVATE(ok_hines)
-
    LOGICAL, SAVE, PROTECTED :: ok_gwd_rando = .false.   ! activate non-orographic stochastic gravity wave drag params
    !$OMP THREADPRIVATE(ok_gwd_rando)
@@ -177 +174 @@
       CALL getin_p('ok_orodr', ok_orodr)
       CALL getin_p('ok_orolf', ok_orolf)
-      CALL getin_p('ok_hines', ok_hines)
       CALL getin_p('ok_gwd_rando', ok_gwd_rando)
       CALL getin_p('zpmm_orodr_t', zpmm_orodr_t)
@@ -200 +196 @@
       WRITE (lunout, *) 'gwd_ini, ok_orodr:', ok_orodr
       WRITE (lunout, *) 'gwd_ini, ok_orolf:', ok_orolf
-      WRITE (lunout, *) 'gwd_ini, ok_hines:', ok_hines
       WRITE (lunout, *) 'gwd_ini, ok_gwd_rando:', ok_gwd_rando
       WRITE (lunout, *) 'gwd_ini, zpmm_orodr_t:', zpmm_orodr_t
@@ -546 +541 @@
          END IF
 
-         IF (.NOT. ok_hines) THEN
-            IF (NW + 4*(NA - 1) + NA >= KLEV) THEN
-               abort_message = 'NW+3*NA>=KLEV Problem to generate waves associated with fronts'
-               CALL abort_physic(modname, abort_message, 1)
-            END IF
+         IF (NW + 4*(NA - 1) + NA >= KLEV) THEN
+             abort_message = 'NW+3*NA>=KLEV Problem to generate waves associated with fronts'
+             CALL abort_physic(modname, abort_message, 1)
          END IF
 

LMDZ6/trunk/libf/phylmd/phyetat0_mod.f90

@@ -66 +66 @@
   USE alpale_mod
   USE compbl_mod_h
-  USE lmdz_gwd_ini, ONLY: ok_hines, ok_gwd_rando, ok_orodr, ok_orolf
+  USE lmdz_gwd_ini, ONLY: ok_gwd_rando, ok_orodr, ok_orolf
 
   IMPLICIT NONE
@@ -671 +671 @@
   IF (ok_gwd_rando) found = &
        phyetat0_get(du_gwd_rando,"du_gwd_rando","du_gwd_rando",0.)
-  IF (.NOT. ok_hines .AND. ok_gwd_rando) found &
+  IF (ok_gwd_rando) found &
        = phyetat0_get(du_gwd_front,"du_gwd_front","du_gwd_front",0.)
 

LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -119 +119 @@
       !$OMP THREADPRIVATE(d_u_lif, d_v_lif)
 ! Tendances Ondes de G non oro (runs strato).
-      REAL, SAVE, ALLOCATABLE :: du_gwd_hines(:,:)
-      !$OMP THREADPRIVATE(du_gwd_hines)
-      REAL, SAVE, ALLOCATABLE :: dv_gwd_hines(:,:)
-      !$OMP THREADPRIVATE(dv_gwd_hines)
       REAL, SAVE, ALLOCATABLE :: dv_gwd_rando(:,:)
       !$OMP THREADPRIVATE(dv_gwd_rando)
@@ -131 +127 @@
       REAL, SAVE, ALLOCATABLE :: west_gwstress(:,:)
       !$OMP THREADPRIVATE(west_gwstress)
-      REAL, SAVE, ALLOCATABLE :: dt_gwd_hines(:,:)
-      !$OMP THREADPRIVATE(dt_gwd_hines)
 ! tendance due a l'oxydation du methane
       REAL, SAVE, ALLOCATABLE :: d_q_ch4(:,:)
@@ -981 +975 @@
       ALLOCATE(topsw0_aero(klon,naero_grp), solsw0_aero(klon,naero_grp))
       ALLOCATE(topswcf_aero(klon,3), solswcf_aero(klon,3))
-      ALLOCATE(du_gwd_hines(klon,klev),dv_gwd_hines(klon,klev))
       ALLOCATE(dv_gwd_rando(klon,klev),dv_gwd_front(klon,klev))
       ALLOCATE(east_gwstress(klon,klev),west_gwstress(klon,klev))
       east_gwstress(:,:)=0. !ym missing init
       west_gwstress(:,:)=0. !ym missing init
-      ALLOCATE(dt_gwd_hines(klon,klev))
       ALLOCATE(d_q_ch4(klon,klev))
       ALLOCATE(stratomask(klon,klev))
@@ -1503 +1495 @@
       DEALLOCATE(load_tmp9)
       DEALLOCATE(load_tmp10)
-      DEALLOCATE(du_gwd_hines,dv_gwd_hines,dt_gwd_hines)
       DEALLOCATE(d_q_ch4)
       DEALLOCATE(dv_gwd_rando,dv_gwd_front)

LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -2402 +2402 @@
   TYPE(ctrl_out), SAVE :: o_dvlif = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'dvlif', 'Orography dV', 'm/s2', (/ ('', i=1, 10) /))
-  TYPE(ctrl_out), SAVE :: o_du_gwd_hines = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
-    'du_gwd_hines', 'Hines GWD dU', 'm/s2', (/ ('', i=1, 10) /))
-  TYPE(ctrl_out), SAVE :: o_dv_gwd_hines = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
-    'dv_gwd_hines', 'Hines GWD dV', 'm/s2', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_du_gwd_front = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'du_gwd_front', 'Fronts GWD dU', 'm/s2', (/ ('', i=1, 10) /))
@@ -2418 +2414 @@
   TYPE(ctrl_out), SAVE :: o_dtlif = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'dtlif', 'Orography dT', 'K/s', (/ ('', i=1, 10) /))
-  TYPE(ctrl_out), SAVE :: o_dthin = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
-    'dthin', 'Hines GWD dT', 'K/s', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_dqch4 = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'dqch4', 'H2O due to CH4 oxidation & photolysis', '(kg/kg)/s', (/ ('', i=1, 10) /))
@@ -2429 +2423 @@
        = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), 'dv_gwd_rando', &
        "Random gravity waves dV/dt", "m/s2", (/ ('', i=1, 10) /))
-  type(ctrl_out), save:: o_ustr_gwd_hines &
-       = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), 'ustr_gwd_hines', &
-       "zonal wind stress Hines gravity waves", "Pa", (/ ('', i=1, 10) /))
-  type(ctrl_out), save:: o_vstr_gwd_hines &
-       = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), 'vstr_gwd_hines', &
-       "meridional wind stress Hines gravity waves", "Pa", (/ ('', i=1, 10) /))
   type(ctrl_out), save:: o_ustr_gwd_front &
        = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), 'ustr_gwd_front', &

LMDZ6/trunk/libf/phylmd/phys_output_var_mod.f90

@@ -142 +142 @@
   LOGICAL, SAVE :: vars_defined = .FALSE. ! ug PAS THREADPRIVATE ET C'EST NORMAL
 
-  REAL, allocatable:: zustr_gwd_hines(:), zvstr_gwd_hines(:) ! (klon)
   REAL, allocatable:: zustr_gwd_front(:), zvstr_gwd_front(:) ! (klon)
   REAL, allocatable:: zustr_gwd_rando(:), zvstr_gwd_rando(:) ! (klon)
-  !$OMP THREADPRIVATE(zustr_gwd_hines, zvstr_gwd_hines)
   !$OMP THREADPRIVATE(zustr_gwd_front, zvstr_gwd_front)
   !$OMP THREADPRIVATE(zustr_gwd_rando, zvstr_gwd_rando)
@@ -242 +240 @@
 
     ! stress from non orographic gravity wave drag params
-    allocate(zustr_gwd_hines(klon), zvstr_gwd_hines(klon))
     allocate(zustr_gwd_front(klon), zvstr_gwd_front(klon))
     allocate(zustr_gwd_rando(klon), zvstr_gwd_rando(klon))
@@ -295 +292 @@
     ! stress from non orographic gravity wave drag params
     
-    deallocate(zustr_gwd_hines, zvstr_gwd_hines)
     deallocate(zustr_gwd_front, zvstr_gwd_front)
     deallocate(zustr_gwd_rando, zvstr_gwd_rando)

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -33 +33 @@
     USE phystokenc_mod, ONLY: offline
     USE water_int_mod, ONLY: water_int
-    USE lmdz_gwd_ini, ONLY: ok_orodr, ok_orolf, ok_hines, ok_gwd_rando
+    USE lmdz_gwd_ini, ONLY: ok_orodr, ok_orolf, ok_gwd_rando
     USE conf_phys_m, ONLY : aerosol_couple, ok_ade, ok_aie, ok_volcan, &
                             flag_aerosol, flag_aerosol_strat, ok_cdnc
@@ -180 +180 @@
          o_duvdf, o_dvvdf, o_duoro, o_dvoro, &
          o_dtoro, o_dulif, o_dvlif, o_dtlif, &
-         o_du_gwd_hines, o_dv_gwd_hines, o_dthin, o_dqch4, o_rsu, &
+         o_dqch4, o_rsu, &
          o_du_gwd_front, o_dv_gwd_front, &
          o_east_gwstress, o_west_gwstress, &
@@ -202 +202 @@
          o_dtr_sat, o_dtr_uscav, o_dtr_wet_cv, o_dtr_wet, & 
          o_trac_cum, o_du_gwd_rando, o_dv_gwd_rando, &
-         o_ustr_gwd_hines,o_vstr_gwd_hines,o_ustr_gwd_rando,o_vstr_gwd_rando, &
+         o_ustr_gwd_rando,o_vstr_gwd_rando, &
          o_ustr_gwd_front,o_vstr_gwd_front, &
          o_sens_prec_liq_oce, o_sens_prec_liq_sic, &
@@ -439 +439 @@
          zw2, fraca, zmax_th, d_q_ajsb, d_t_ec, d_u_vdf, &
          d_v_vdf, d_u_oro, d_v_oro, d_t_oro, d_u_lif, &
-         d_v_lif, d_t_lif, du_gwd_hines, dv_gwd_hines, dt_gwd_hines, &
+         d_v_lif, d_t_lif,  &
          dv_gwd_rando, dv_gwd_front, &
          east_gwstress, west_gwstress, &
@@ -486 +486 @@
          bils_ec,bils_ech, bils_tke, bils_kinetic, bils_latent, bils_enthalp, &
          itau_con, nfiles, clef_files, nid_files, dryaod_diag, &
-         zustr_gwd_hines, zvstr_gwd_hines,zustr_gwd_rando, zvstr_gwd_rando, &
+         zustr_gwd_rando, zvstr_gwd_rando, &
          zustr_gwd_front, zvstr_gwd_front, sza_o,    &
          sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o, &
@@ -2785 +2785 @@
        ENDIF
 
-       IF (ok_hines) THEN
-          IF (vars_defined) zx_tmp_fi3d=du_gwd_hines/pdtphys
-          CALL histwrite_phy(o_du_gwd_hines, zx_tmp_fi3d)
-
-          IF (vars_defined) zx_tmp_fi3d= dv_gwd_hines/pdtphys          
-          CALL histwrite_phy(o_dv_gwd_hines, zx_tmp_fi3d)
-          
-          IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=dt_gwd_hines(1:klon,1:klev)/pdtphys
-          CALL histwrite_phy(o_dthin, zx_tmp_fi3d)
-          CALL histwrite_phy(o_ustr_gwd_hines, zustr_gwd_hines)
-          CALL histwrite_phy(o_vstr_gwd_hines, zvstr_gwd_hines)
-       ENDIF
-
-       IF (.not. ok_hines .and. ok_gwd_rando) THEN
+       IF (ok_gwd_rando) THEN
           IF (vars_defined)  zx_tmp_fi3d=du_gwd_front / pdtphys
           CALL histwrite_phy(o_du_gwd_front, zx_tmp_fi3d)

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -254 +254 @@
          d_t_ec, &
                                 !
-         du_gwd_hines,dv_gwd_hines,dt_gwd_hines, &
          dv_gwd_rando,dv_gwd_front, &
          east_gwstress,west_gwstress, &
@@ -4685 +4684 @@
                    d_u_oro, d_v_oro, d_t_oro, zustr_oro, zvstr_oro, &
                    d_u_lif, d_v_lif, d_t_lif, zustr_lif, zvstr_lif, &
-                   du_gwd_hines, dv_gwd_hines, dt_gwd_hines, zustr_gwd_hines, zvstr_gwd_hines, &
                    du_gwd_front, dv_gwd_front, zustr_gwd_front, zvstr_gwd_front, &
                    du_gwd_rando, dv_gwd_rando, zustr_gwd_rando, zvstr_gwd_rando, &

LMDZ6/trunk/libf/phylmdiso/phyetat0_mod.F90

@@ -78 +78 @@
   USE alpale_mod
   USE compbl_mod_h
-  USE lmdz_gwd_ini, ONLY: ok_hines, ok_gwd_rando, ok_orodr, ok_orolf
+  USE lmdz_gwd_ini, ONLY: ok_gwd_rando, ok_orodr, ok_orolf
 
   IMPLICIT NONE
@@ -675 +675 @@
   IF (ok_gwd_rando) found = &
        phyetat0_get(du_gwd_rando,"du_gwd_rando","du_gwd_rando",0.)
-  IF (.NOT. ok_hines .AND. ok_gwd_rando) found = &
+  IF (ok_gwd_rando) found = &
        phyetat0_get(du_gwd_front,"du_gwd_front","du_gwd_front",0.)
 

LMDZ6/trunk/libf/phylmdiso/phyredem.F90

@@ -13 +13 @@
   USE clesphys_mod_h
   USE dimphy, ONLY: klon, klev
-  USE lmdz_gwd_ini, ONLY: ok_gwd_rando, ok_hines, ok_orodr, ok_orolf 
+  USE lmdz_gwd_ini, ONLY: ok_gwd_rando, ok_orodr, ok_orolf 
   USE fonte_neige_mod,  ONLY : fonte_neige_final
   USE pbl_surface_mod,  ONLY : pbl_surface_final
@@ -461 +461 @@
          "tendency on zonal wind due to flott gravity waves", du_gwd_rando)
 
-    IF (.NOT. ok_hines .AND. ok_gwd_rando) CALL put_field(pass,"du_gwd_front", &
+    IF (ok_gwd_rando) CALL put_field(pass,"du_gwd_front", &
          "tendency on zonal wind due to acama gravity waves", du_gwd_front)
 

LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

@@ -264 +264 @@
        d_t_lif,d_u_lif,d_v_lif, &
        d_t_ec, &
-       !
-       ! du_gwd_hines,dv_gwd_hines,d_t_hin, & FH a enlever avant 2026/05 d_t_hin pas dans phylmdiso
-       du_gwd_hines,dv_gwd_hines,dt_gwd_hines, &
        dv_gwd_rando,dv_gwd_front, &
        east_gwstress,west_gwstress, &
@@ -6226 +6223 @@
                    d_u_oro, d_v_oro, d_t_oro, zustr_oro, zvstr_oro, &
                    d_u_lif, d_v_lif, d_t_lif, zustr_lif, zvstr_lif, &
-                   du_gwd_hines, dv_gwd_hines, dt_gwd_hines, zustr_gwd_hines, zvstr_gwd_hines, &
                    du_gwd_front, dv_gwd_front, zustr_gwd_front, zvstr_gwd_front, &
                    du_gwd_rando, dv_gwd_rando, zustr_gwd_rando, zvstr_gwd_rando, &