Changeset 3767

Timestamp:

May 16, 2025, 3:35:51 PM (3 weeks ago)

Author:

jliu

Message:

implementation of GW mixing for tracers: tracers
vertical gradients are considered.

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/nonoro_gwd_mix_mod.F90 (modified) (5 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -4841 +4841 @@
 == 14/05/2025 == JBC
 Cleaning and improvement of robustness for "run0" and "run_month1" scripts with file checks and clearer errors.
+
+== 16/05/2025 == JLIU
+New implementation of non-orographic GWs mixing for tracers: now
+D_mix_tracer = D_eddy*tracers_density_gradient
+while keep the mixing coefficient unchanged for u and theta

trunk/LMDZ.MARS/libf/phymars/nonoro_gwd_mix_mod.F90

@@ -127 +127 @@
     REAL d_eddy_mix_tot(ngrid, nlayer+1)
     REAL d_eddy_mix(NW,ngrid)          ! Comprehensive Diffusion coefficients
+    REAL d_wave(NW,ngrid)              ! coefficients consider the tracers' gradients
     REAL u_eddy_mix_p(NW, ngrid)       ! Zonal Diffusion coefficients
     REAL v_eddy_mix_p(NW, ngrid)       ! Meridional Diffusion coefficients
@@ -141 +142 @@
     REAL pq_eddy_mix_tot(ngrid, nlayer+1,nq)
     REAL zq(ngrid,nlayer,nq) ! advected field nq
+    REAL zq_var(ngrid,nlayer,nq)
+    REAL dzq_var(ngrid,nlayer,nq)
+    REAL mdzq_var(nlayer,nq)
+    REAL zq_ave(nlayer,nq)
+    REAL dzq_ave(nlayer,nq)
+    REAL zq_ratio(nlayer,nq)
     REAL, save:: eff 
 !$OMP THREADPRIVATE(eff)
@@ -251 +258 @@
     zq(:,:,:)=pq(:,:,:)+dtime*pdq(:,:,:)
     hh(:,:)=pht(:,:)+dtime*pdht(:,:)
+
+!  tracer average and gradients for mixing
+    zq_ave(:,:)=0.
+    zq_ave(:,:) = SUM(zq(:,:,:), dim=1) / real(ngrid)
+    
+    zq_var(:,:,:)=0.
+    DO ii=1,ngrid
+      zq_var(ii,:,:)=zq(ii,:,:)-zq_ave(:,:)
+    endDO
+
+    dzq_var(:,:,:)=0.
+    mdzq_var(:,:) =0.
+    dzq_ave(:,:)  =0.
+    zq_ratio(:,:) =0.
+    DO LL=1,nlayer-1
+       dzq_var(:, LL, :) = (zq_var(:, LL+1, :) - zq_var(:, LL, :))**2.
+       mdzq_var(LL, :)   = SUM(dzq_var(:, LL, :), dim=1)/real(ngrid)
+       dzq_ave(LL, :)    = (zq_ave(LL+1, :) - zq_ave(LL, :))**2.
+       zq_ratio(LL,:)    = MAX(0., (mdzq_var(LL+1,:) + mdzq_var(LL,:)))&
+                          /MAX(1E-15, (dzq_ave(LL+1,:) + dzq_ave(LL,:)))
+      ! do qq=1,nq
+      !   print*, 'ratio=', zq_ratio(LL,QQ)
+      ! enddo
+    endDO
+    dzq_var(:,nlayer,:)    = dzq_var(:,nlayer-1, :)
+    mdzq_var(nlayer,:)     = mdzq_var(nlayer-1, :)
+    dzq_ave(nlayer,:)      = dzq_ave(nlayer-1, :)
 ! Compute the real mass density by rho=p/R(T)T
      DO ll=1,nlayer
@@ -623 +657 @@
     pq_eddy_mix_tot(:, :,:)=0.
     d_eddy_mix(:,:)=0.
+    d_wave(:, :) =0.
     d_eddy_mix_p_ll(nlayer,:,:)=d_eddy_mix_p_ll(nlayer-1,:,:)
     d_eddy_mix_tot(:, :) =0.
@@ -656 +691 @@
       ! print*, "this line is for tunning"
       !endif
+     ! d_wave(JW,:) = d_eddy_mix(JW, :)*(mdzq_var(LL,QQ)/dzq_ave(LL,QQ))**2.
       DO QQ=1,NQ
          DO JW=1,NW
-         pq_eddy_mix_p(JW, :, QQ) = d_eddy_mix(JW, :)* (zq(:, LL + 1,QQ)- zq(:, LL, QQ)) &
-                                   /(ZH(:, LL + 1)- ZH(:, LL))                           &
+         d_wave(JW, :) = d_eddy_mix(JW, :)*zq_ratio(LL,QQ)
+         pq_eddy_mix_p(JW, :, QQ) = d_wave(JW, :)* (zq(:, LL + 1,QQ)- zq(:, LL, QQ)) &
+                                   /(ZH(:, LL + 1)- ZH(:, LL))                       &
                                    *SIGN(1.,intr_freq_p(JW, :)) * SIN(ZP(JW, :))  
       !   pq_eddy_mix_tot(:, LL+1,QQ) = pq_eddy_mix_tot(:, LL+1,QQ)                   &