Changeset 3393 for trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90

Timestamp:

Jul 4, 2024, 4:14:10 PM (5 months ago)

Author:

jliu

Message:

Update of non-orographic gravity waves mixing scheme. 1)mixing in potential
temperature are added. 2)all mixing in q,u,theta now are implemented by AR-1
algorithm. Tests (MY29,MY32-35) runs show that this implementation has limited
impact to the temperature/tides fields.

File:

• trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90 (modified) (3 diffs)

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

@@ -24 +24 @@
 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 !dr_depflux_gwd
+                               df_eddymix_flx, dh_eddymix_gwd, dq_eddymix_gwd
 use compute_dtau_mod,    only: dtau
 use dust_rad_adjust_mod, only: dust_rad_adjust_prev,dust_rad_adjust_next
@@ -623 +623 @@
 
 if (startphy_file) then
+   call get_field("dh_eddymix_gwd",dh_eddymix_gwd,found,indextime)
+   if (.not.found) then
+      write(*,*) "phyetat0: <dh_eddymix_gwd> not in file"
+      dh_eddymix_gwd(:,:)=0.
+   endif
+else
+dh_eddymix_gwd(:,:)=0.
+endif ! if (startphy_file)
+write(*,*) "phyetat0: Memory of PT tendency due to non-orographic GW mixing"
+write(*,*) " <dh_eddymix_gwd> range:", &
+             minval(dh_eddymix_gwd), maxval(dh_eddymix_gwd)
+
+
+if (startphy_file) then
    call get_field("dv_nonoro_gwd",dv_nonoro_gwd,found,indextime)
    if (.not.found) then
@@ -648 +662 @@
              minval(dv_eddymix_gwd), maxval(dv_eddymix_gwd)
 
+if (startphy_file) then
+   call get_field("dq_eddymix_gwd",dq_eddymix_gwd,found,indextime)
+   if (.not.found) then
+      write(*,*) "phyetat0: <dq_eddymix_gwd> not in file"
+      dq_eddymix_gwd(:,:,:)=0.
+   endif
+else ! ! if (startphy_file)
+dq_eddymix_gwd(:,:,:)=0.
+endif ! if (startphy_file)
+write(*,*) "phyetat0: Memory of tracers tendency due to non-orographic GW mixing"
+write(*,*) " <dq_eddymix_gwd> range:", &
+             minval(dq_eddymix_gwd), maxval(dq_eddymix_gwd)
+
+     
 !if (startphy_file) then
 !   call get_field("dr_depflux_gwd",dr_depflux_gwd,found,indextime)