Changeset 3935

Timestamp:

Oct 24, 2025, 4:39:28 PM (3 months ago)

Author:

jliu

Message:

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

Location:

trunk/LMDZ.PLUTO

Files:

• changelog.txt (modified) (1 diff)
• libf/phypluto/callkeys_mod.F90 (modified) (1 diff)
• libf/phypluto/inifis_mod.F90 (modified) (3 diffs)
• libf/phypluto/nonoro_gwd_mix_mod.F90 (added)
• libf/phypluto/nonoro_gwd_ran_mod.F90 (added)
• libf/phypluto/phyetat0_mod.F90 (modified) (1 diff)
• libf/phypluto/physiq_mod.F90 (modified) (4 diffs)

trunk/LMDZ.PLUTO/changelog.txt

@@ -1879 +1879 @@
 flag to trigger outputing a diagsoil.nc file or not.
 
+== 24/10/2025 == Jliu
+Add non-orographic gravity waves and induced turbulence to the model. 
+The model is compileable with the two new shcemes, however, carefull 
+parameter selection is needed since the current parameters are for Mars.

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

@@ -6 +6 @@
       logical,save :: calladj,calltherm,n2cond,callsoil,calllott
 !$OMP THREADPRIVATE(calladj,calltherm,n2cond,callsoil,calllott)
+      logical,save :: calllott_nonoro
+!$OMP THREADPRIVATE(calllott_nonoro)      
       logical,save :: callconduct,callmolvis,callmoldiff
 !$OMP THREADPRIVATE(callconduct,callmolvis,callmoldiff)

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

@@ -23 +23 @@
                           mugaz, pi, avocado
   use planetwide_mod, only: planetwide_sumval
+  use nonoro_gwd_mix_mod, only: calljliu_gwimix
   use callkeys_mod
+  use nonoro_gwd_ran_mod, only: ini_nonoro_gwd_ran, &
+                                end_nonoro_gwd_ran
+  use nonoro_gwd_mix_mod, only: ini_nonoro_gwd_mix, &
+                                end_nonoro_gwd_mix
+                                
   use surfdat_h
   use wstats_mod, only: callstats
@@ -373 +379 @@
      if (is_master) write(*,*) trim(rname)//": callsoil = ",callsoil
 
+     if (is_master) write(*,*)trim(rname)//&
+       ": call orographic gravity waves ?"
      calllott=.false. ! default value
      call getin_p("calllott",calllott)
      if (is_master) write(*,*)trim(rname)//": calllott = ",calllott
+
+     if (is_master) write(*,*)trim(rname)//&
+       ": call  non-orographic gravity waves ?"
+     calllott_nonoro=.false. ! default value
+     call getin_p("calllott_nonoro",calllott_nonoro)
+     if (is_master) write(*,*)trim(rname)//&
+       ": calllott_nonoro = ",calllott_nonoro
+
+     if (is_master) write(*,*)trim(rname)//&       
+       ": call jliu's non-orogrphic GW-induced turbulence"
+     calljliu_gwimix=.false. ! default value
+     call getin_p("calljliu_gwimix",calljliu_gwimix)
+     if (is_master) write(*,*)trim(rname)//&    
+       ": calljliu_gwimix = ",calljliu_gwimix
 
      if (is_master) write(*,*)trim(rname)//&
@@ -1500 +1522 @@
   call ini_comsoil_h(ngrid)
 
+  ! allocate arrays in "nonoro_gwd_ran_mod"
+  call end_nonoro_gwd_ran
+  call ini_nonoro_gwd_ran(ngrid,nlayer)
+
+  ! allocate arrays in "nonoro_gwd_mix_mod"
+  call end_nonoro_gwd_mix
+  call ini_nonoro_gwd_mix(ngrid,nlayer,nq)
+
   END SUBROUTINE inifis
 

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

@@ -27 +27 @@
                      get_field, get_var, inquire_field, &
                      inquire_dimension, inquire_dimension_length
+  use nonoro_gwd_ran_mod,  only: du_nonoro_gwd, dv_nonoro_gwd
+  use nonoro_gwd_mix_mod,  only: du_eddymix_gwd, dv_eddymix_gwd, de_eddymix_rto, &
+                               df_eddymix_flx, dh_eddymix_gwd, dq_eddymix_gwd             
   use callkeys_mod, only: surfalbedo,surfemis, callsoil
   use mod_phys_lmdz_para, only : is_master

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 :