Changeset 3935 for trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90

Timestamp:

Oct 24, 2025, 4:39:28 PM (7 weeks ago)

Author:

jliu

Message:

#Non-orographic gravity waves and induced turbulence
#are added to the model. Further tunning are still needed.

File:

• trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90 (modified) (4 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90

@@ -47 +47 @@
       use calldrag_noro_mod, only: calldrag_noro
       use time_phylmdz_mod, only: daysec
+      use nonoro_gwd_ran_mod, only: nonoro_gwd_ran
+      use nonoro_gwd_mix_mod, only: nonoro_gwd_mix, calljliu_gwimix
+
       use callkeys_mod, only: albedo_spectral_mode, calladj, calldifv,        &
                               calllott, callrad, callsoil, nosurf,            &
                               callconduct,callmolvis,callmoldiff,             &
-                              corrk,                                          &
+                              corrk, calllott_nonoro,                         &
                               diurnal, enertest, fat1au,                      &
                               icetstep, intheat, iradia, kastprof,            &
@@ -67 +70 @@
                               global1d, szangle,                              &
                               callmufi, evol1d
+
       use check_fields_mod, only: check_physics_fields
       use conc_mod, only: rnew, cpnew, ini_conc_mod
@@ -476 +480 @@
       REAL d_u_hin(ngrid,nlayer), d_v_hin(ngrid,nlayer)
       REAL d_t_hin(ngrid,nlayer)
+      REAL d_u_mix(ngrid,nlayer), d_v_mix(ngrid,nlayer)
+      REAL d_t_mix(ngrid,nlayer), zdq_mix(ngrid,nlayer,nq)
 !  Diagnostics 2D of gw_nonoro
       REAL zustrhi(ngrid), zvstrhi(ngrid)
@@ -1384 +1390 @@
       endif ! end of 'calladj'
 
+! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+! IV.b Non-orographic gravity waves and induced turbulence
+! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+        IF (calllott_nonoro) THEN
+         CALL nonoro_gwd_ran(ngrid,nlayer,ptimestep, &
+                      cpnew,rnew,                    &
+                      zplay,                         &
+                      zmax_th,                       &
+                      pt, pu, pv,                    &
+                      pdt, pdu, pdv,                 &
+                      zustrhi,zvstrhi,               &
+                      d_t_hin, d_u_hin, d_v_hin)
+           IF (calljliu_gwimix) THEN
+            CALL nonoro_gwd_mix(ngrid,nlayer,ptimestep,  &
+                      nq,cpnew, rnew,                    &
+                      zplay,                             &
+                      zmax_th,                           &
+                      pt, pu, pv, pq, zh,                &
+                  !loss,  chemical reaction loss rates
+                      pdt, pdu, pdv, pdq, zdh,           &
+                  !  zustrhi,zvstrhi,
+                      zdq_mix, d_t_mix, d_u_mix, d_v_mix)
+           ENDIF
+  
+  !  Update tendencies
+           pdt(1:ngrid,1:nlayer)=pdt(1:ngrid,1:nlayer) +   &
+                                 d_t_hin(1:ngrid,1:nlayer)
+           pdu(1:ngrid,1:nlayer)=pdu(1:ngrid,1:nlayer) +   &
+                                 d_u_hin(1:ngrid,1:nlayer)
+           pdv(1:ngrid,1:nlayer)=pdv(1:ngrid,1:nlayer) +   &
+                                 d_v_hin(1:ngrid,1:nlayer)
+  !  Update tendencies of gw mixing
+          IF (calljliu_gwimix) THEN
+           pdt(1:ngrid,1:nlayer)=pdt(1:ngrid,1:nlayer) +   &
+                                 d_t_mix(1:ngrid,1:nlayer)
+           pdu(1:ngrid,1:nlayer)=pdu(1:ngrid,1:nlayer) +   &
+                                 d_u_mix(1:ngrid,1:nlayer)
+           pdv(1:ngrid,1:nlayer)=pdv(1:ngrid,1:nlayer) +   &
+                                 d_v_mix(1:ngrid,1:nlayer)
+           pdq(1:ngrid,1:nlayer,1:nq)=pdq(1:ngrid,1:nlayer,1:nq) + &
+                                      zdq_mix(1:ngrid,1:nlayer,1:nq)
+          ENDIF
+  
+  
+        ENDIF ! of IF (calllott_nonoro)
+      
 !-----------------------------------------------
 !   V. Nitrogen condensation-sublimation :