Changeset 3144 for trunk/LMDZ.MARS/libf

Timestamp:

Nov 29, 2023, 4:04:35 PM (14 months ago)

Author:

jliu

Message:

Mars PCM
New scheme "nonoro_gwd_mix_mod" to simulate non-orographic gravity waves induced mixing. The scheme flag "calljliu_gwimix" has been put to "false" as default.
There are three tunable parameters that have been set up, as "nonoro_gwimixing_eff=0.1", "nonoro_gwimixing_eff1=0.1", "nonoro_gwimixing_vdl=1.5".
Remember to use full-layers model configuration (64x48x73) when you try to call this scheme because we need to evaluate the saturation altitude of each monochromatic wave.
More than ~75% of the waves saturated at altitude above 120 km.
Jliu

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• conf_phys.F (modified) (2 diffs)
• nonoro_gwd_mix_mod.F90 (added)
• phyetat0_mod.F90 (modified) (2 diffs)
• phys_state_var_init_mod.F90 (modified) (2 diffs)
• physiq_mod.F (modified) (4 diffs)

trunk/LMDZ.MARS/libf/phymars/conf_phys.F

@@ -41 +41 @@
       use microphys_h, only: mteta
       use comsoil_h, only: adsorption_soil, choice_ads
+      use nonoro_gwd_mix_mod, only: calljliu_gwimix
 
       IMPLICIT NONE
@@ -296 +297 @@
          call getin_p("calllott_nonoro",calllott_nonoro)
          write(*,*)" calllott_nonoro = ",calllott_nonoro
+
+         write(*,*)"call jliu's nonoro GWs-induced mixing ",
+     &             "scheme ?"
+         calljliu_gwimix=.false. ! default value
+         call getin_p("calljliu_gwimix",calljliu_gwimix)
+         write(*,*)" calljliu_gwimix = ",calljliu_gwimix
+
 
 ! rocket dust storm injection scheme

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

@@ -21 +21 @@
                      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 !dr_depflux_gwd
+
   use compute_dtau_mod, only: dtau
   use dust_rad_adjust_mod, only: dust_rad_adjust_prev,dust_rad_adjust_next
@@ -621 +624 @@
 write(*,*) " <dv_nonoro_gwd> range:", &
              minval(dv_nonoro_gwd), maxval(dv_nonoro_gwd)
+if (startphy_file) then 
+   call get_field("du_eddymix_gwd",du_eddymix_gwd,found,indextime)
+   if (.not.found) then
+      write(*,*) "phyetat0: <du_eddymix_gwd> not in file"
+      du_eddymix_gwd(:,:)=0.
+   endif
+else
+du_eddymix_gwd(:,:)=0.
+endif ! if (startphy_file)
+write(*,*) "phyetat0: Memory of zonal wind tendency due to non-orographic GW mixing"
+write(*,*) " <du_eddymix_gwd> range:", &
+             minval(du_eddymix_gwd), maxval(du_eddymix_gwd)
+ 
+
+if (startphy_file) then 
+   call get_field("dv_nonoro_gwd",dv_nonoro_gwd,found,indextime)
+   if (.not.found) then
+      write(*,*) "phyetat0: <dv_nonoro_gwd> not in file"
+      dv_nonoro_gwd(:,:)=0.
+   endif
+else ! ! if (startphy_file)
+dv_nonoro_gwd(:,:)=0.
+endif ! if (startphy_file)
+write(*,*) "phyetat0: Memory of meridional wind tendency due to non-orographic GW"
+write(*,*) " <dv_nonoro_gwd> range:", &
+             minval(dv_nonoro_gwd), maxval(dv_nonoro_gwd)
+
+if (startphy_file) then 
+   call get_field("dv_eddymix_gwd",dv_eddymix_gwd,found,indextime)
+   if (.not.found) then
+      write(*,*) "phyetat0: <dv_eddymix_gwd> not in file"
+      dv_eddymix_gwd(:,:)=0.
+   endif
+else ! ! if (startphy_file)
+dv_eddymix_gwd(:,:)=0.
+endif ! if (startphy_file)
+write(*,*) "phyetat0: Memory of meridional wind tendency due to non-orographic GW mixing"
+write(*,*) " <dv_eddymix_gwd> range:", &
+             minval(dv_eddymix_gwd), maxval(dv_eddymix_gwd)
+
+!if (startphy_file) then 
+!   call get_field("dr_depflux_gwd",dr_depflux_gwd,found,indextime)
+!   if (.not.found) then
+!      write(*,*) "phyetat0: <dr_depflux_gwd> not in file"
+!      dr_depflux_gwd(:,:,:)=0.
+!   endif
+!else ! ! if (startphy_file)
+!!dr_depflux_gwd(:,:,:)=0.
+!endif ! if (startphy_file)
+!write(*,*) "phyetat0: Memory of meridional wind tendency due to non-orographic GW mixing"
+!write(*,*) " <dr_depflux_gwd> range:", &
+!             minval(dr_depflux_gwd), maxval(dr_depflux_gwd)
+
+if (startphy_file) then 
+   call get_field("de_eddymix_rto",de_eddymix_rto,found,indextime)
+   if (.not.found) then
+      write(*,*) "phyetat0: <de_eddymix_rto> not in file"
+      de_eddymix_rto(:,:)=0.
+   endif
+else ! ! if (startphy_file)
+de_eddymix_rto(:,:)=0.
+endif ! if (startphy_file)
+write(*,*) "phyetat0: Memory of meridional wind tendency due to non-orographic GW mixing"
+write(*,*) " <de_eddymix_rto> range:", &
+             minval(de_eddymix_rto), maxval(de_eddymix_rto)
+
+if (startphy_file) then 
+   call get_field("df_eddymix_flx ",df_eddymix_flx ,found,indextime)
+   if (.not.found) then
+      write(*,*) "phyetat0: <df_eddymix_flx > not in file"
+      df_eddymix_flx (:,:)=0.
+   endif
+else ! ! if (startphy_file)
+df_eddymix_flx (:,:)=0.
+endif ! if (startphy_file)
+write(*,*) "phyetat0: Memory of meridional wind tendency due to non-orographic GW mixing"
+write(*,*) " <df_eddymix_flx > range:", &
+             minval(df_eddymix_flx ), maxval(df_eddymix_flx )
+
+
 
 ! tracer on surface

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

@@ -57 +57 @@
       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 dust_param_mod, only: ini_dust_param_mod, &
                                 end_dust_param_mod
@@ -170 +172 @@
       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)
+
       ! allocate arrays in "dust_param_mod"
       call end_dust_param_mod

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -90 +90 @@
       use compute_dtau_mod, only: compute_dtau
       use nonoro_gwd_ran_mod, only: nonoro_gwd_ran
+      use nonoro_gwd_mix_mod, only: nonoro_gwd_mix, calljliu_gwimix
       use check_fields_mod, only: check_physics_fields
       use surfini_mod, only: surfini
@@ -492 +493 @@
       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)
+
 c  Diagnostics 2D of gw_nonoro
       REAL zustrhi(ngrid), zvstrhi(ngrid)
@@ -1690 +1694 @@
      &               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, 
+                !loss,  chemical reaction loss rates
+     &               pdt, pdu, pdv, pdq,
+                !  zustrhi,zvstrhi,
+     &               zdq_mix, d_t_mix, d_u_mix, d_v_mix)
+         ENDIF
 
 !  Update tendencies
@@ -1698 +1713 @@
          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)