Changeset 5803 for LMDZ6/trunk/libf/phylmd/calwake.f90

Timestamp:

Sep 5, 2025, 12:52:13 AM (5 months ago)

Author:

fhourdin

Message:

Work in progress on wake (JYG&FH)

File:

• LMDZ6/trunk/libf/phylmd/calwake.f90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/calwake.f90

@@ -1 +1 @@
 ! $Id$
-!$gpum horizontal klon nlon
 MODULE calwake_mod
   PRIVATE
@@ -58 +57 @@
   USE indice_sol_mod, ONLY: is_oce
   USE print_control_mod, ONLY: lunout, prt_level
-  USE lmdz_wake, ONLY : wake, wake2
+  USE lmdz_wake, ONLY : wake, wake2, wake3
   USE alpale_mod, ONLY: iflag_wake
   USE yomcst_mod_h
@@ -139 +138 @@
     print *, '-> calwake, wake_s, awake_s, wgen input ', wake_s(1), awake_s(1), wgen(1)
   ENDIF
+  print*,'NOUVEAU WAKE ================================='
 
   rdcp = 1./3.5
@@ -256 +256 @@
   ELSE IF (iflag_wake/10 == 2) THEN
     CALL wake2(klon,klev,znatsurf, p, ph, pi, dtime, &
+      te, qe, omgbe, &
+      dtdwn, dqdwn, amdwn, amup, dta, dqa, wgen, &
+      sigd0, Cin, &
+      dtw, dqw, sigmaw, asigmaw, wdens, awdens, &                      ! state variables
+      dth, hw, wape, fip, gfl, &
+      dtls, dqls, ktopw, omgbdth, dp_omgb, tx, qx, &
+      dtke, dqke, omg, dp_deltomg, spread, cstar, &
+      d_deltat_gw, &
+      d_deltatw, d_deltaqw, d_sigmaw, d_asigmaw, d_wdens, d_awdens)      ! tendencies
+
+  ELSE IF (iflag_wake/10 == 3) THEN
+    CALL wake3(klon,klev,znatsurf, p, ph, pi, dtime, &
       te, qe, omgbe, &
       dtdwn, dqdwn, amdwn, amup, dta, dqa, wgen, &