12/03/2024 ==Jliu
Synchronize all default non-orographic gravity wave parameters in callphys.def with the ones in scheme. The parameters are as follows:
#Eliassen-Palm FLux(only if calllott_nonoro=.true.)
nonoro_gwd_epflux_max=5.E-4
# saturation parameter non-orographic gravity waves(only if calllott_nonoro=.true.)
nonoro_gwd_sat=1.5
# maximum of the phase speed of the gravity waves(only if calllott_nonoro=.true.)
nonoro_gwd_cmax=50
# value of the dissaption coefficiet(only if calllott_nonoro=.true.)
nonoro_gwd_rdiss=0.07
# value of the max wave number
nonoro_gwd_kmax=1.E-4
# value of the min wave number
nonoro_gwd_kmin=7.E-6
# value to control the launch altitude
nonoro_gwd_xlaunch=0.6
# effective factor
nonoro_gwimixing_eff=0.1
# effective factor
nonoro_gwimixing_eff1=0.1
# mixing vertical decrease rate
nonoro_gwimixing_vdl=1.5
=========================================================
+== 12/03/2024 ==Jliu
+Synchronize all default non-orographic gravity wave parameters in callphys.def with the ones in scheme. The parameters are as follows:
+#Eliassen-Palm FLux(only if calllott_nonoro=.true.)
+nonoro_gwd_epflux_max=5.E-4
+# saturation parameter non-orographic gravity waves(only if calllott_nonoro=.true.)
+nonoro_gwd_sat=1.5
+# maximum of the phase speed of the gravity waves(only if calllott_nonoro=.true.)
+nonoro_gwd_cmax=50
+# value of the dissaption coefficiet(only if calllott_nonoro=.true.)
+nonoro_gwd_rdiss=0.07
+# value of the max wave number
+nonoro_gwd_kmax=1.E-4
+# value of the min wave number
+nonoro_gwd_kmin=7.E-6
+# value to control the launch altitude
+nonoro_gwd_xlaunch=0.6
+# effective factor
+nonoro_gwimixing_eff=0.1
+# effective factor
+nonoro_gwimixing_eff1=0.1
+# mixing vertical decrease rate
+nonoro_gwimixing_vdl=1.5
+
+
Index: libf/phymars/nonoro_gwd_mix_mod.F90
===================================================================
--- libf/phymars/nonoro_gwd_mix_mod.F90 (revision 3262)
+++ libf/phymars/nonoro_gwd_mix_mod.F90 (working copy)
@@ -198,31 +198,31 @@
!-----------------------------------------------------------------------------------------------------------------------
IF (firstcall) THEN
write(*,*) "nonoro_gwd_mix: non-oro GW-induced mixing is active!"
- epflux_max = 7.E-7 ! Mars' value !!
+ epflux_max = 5.E-4 ! Mars' value !!
call getin_p("nonoro_gwd_epflux_max", epflux_max)
write(*,*) "NONORO_GWD_MIX: epflux_max=", epflux_max
- sat = 1. ! default gravity waves saturation value !!
+ sat = 1.5 ! default gravity waves saturation value !!
call getin_p("nonoro_gwd_sat", sat)
write(*,*) "NONORO_GWD_MIX: sat=", sat
! cmax = 50. ! default gravity waves phase velocity value !!
! call getin_p("nonoro_gwd_cmax", cmax)
! write(*,*) "NONORO_GWD_MIX: cmax=", cmax
- rdiss=0.01 ! default coefficient of dissipation !!
+ rdiss=0.07 ! default coefficient of dissipation !!
call getin_p("nonoro_gwd_rdiss", rdiss)
write(*,*) "NONORO_GWD_MIX: rdiss=", rdiss
- kmax=7.E-4 ! default Max horizontal wavenumber !!
+ kmax=1.E-4 ! default Max horizontal wavenumber !!
call getin_p("nonoro_gwd_kmax", kmax)
write(*,*) "NONORO_GWD_MIX: kmax=", kmax
- kmin=2.e-5 ! default Max horizontal wavenumber !!
+ kmin=7.E-6 ! default Max horizontal wavenumber !!
call getin_p("nonoro_gwd_kmin", kmin)
write(*,*) "NONORO_GWD_MIX: kmin=", kmin
- xlaunch=0.4 ! default GW luanch altitude controller !!
+ xlaunch=0.6 ! default GW luanch altitude controller !!
call getin_p("nonoro_gwd_xlaunch", xlaunch)
write(*,*) "NONORO_GWD_MIX: xlaunch=", xlaunch
- eff=5.E-3 ! Diffusion effective factor !!
+ eff=0.1 ! Diffusion effective factor !!
call getin_p("nonoro_gwimixing_eff", eff)
write(*,*) "NONORO_GWD_MIX: eff=", eff
- eff1=1.0 ! Diffusion effective factor !!
+ eff1=0.1 ! Diffusion effective factor !!
call getin_p("nonoro_gwimixing_eff1", eff1)
write(*,*) "NONORO_GWD_MIX: eff1=", eff1
vdl=1.5 ! Diffusion effective factor !!
Index: libf/phymars/nonoro_gwd_ran_mod.F90
===================================================================
--- libf/phymars/nonoro_gwd_ran_mod.F90 (revision 3262)
+++ libf/phymars/nonoro_gwd_ran_mod.F90 (working copy)
@@ -167,25 +167,25 @@
!-----------------------------------------------------------------------------------------------------------------------
IF (firstcall) THEN
write(*,*) "nonoro_gwd_ran: FLott non-oro GW scheme is active!"
- epflux_max = 7.E-7 ! Mars' value !!
+ epflux_max = 5.E-4 ! Mars' value !!
call getin_p("nonoro_gwd_epflux_max", epflux_max)
write(*,*) "nonoro_gwd_ran: epflux_max=", epflux_max
- sat = 1. ! default gravity waves saturation value !!
+ sat = 1.5 ! default gravity waves saturation value !!
call getin_p("nonoro_gwd_sat", sat)
write(*,*) "nonoro_gwd_ran: sat=", sat
! cmax = 50. ! default gravity waves phase velocity value !!
! call getin_p("nonoro_gwd_cmax", cmax)
! write(*,*) "nonoro_gwd_ran: cmax=", cmax
- rdiss=0.01 ! default coefficient of dissipation !!
+ rdiss=0.07 ! default coefficient of dissipation !!
call getin_p("nonoro_gwd_rdiss", rdiss)
write(*,*) "nonoro_gwd_ran: rdiss=", rdiss
- kmax=7.E-4 ! default Max horizontal wavenumber !!
+ kmax=1.E-4 ! default Max horizontal wavenumber !!
call getin_p("nonoro_gwd_kmax", kmax)
write(*,*) "nonoro_gwd_ran: kmax=", kmax
- kmin=2.e-5 ! default Max horizontal wavenumber !!
+ kmin=7.E-6 ! default Max horizontal wavenumber !!
call getin_p("nonoro_gwd_kmin", kmin)
write(*,*) "nonoro_gwd_ran: kmin=", kmin
- xlaunch=0.4 ! default GW luanch altitude controller !!
+ xlaunch=0.6 ! default GW luanch altitude controller !!
call getin_p("nonoro_gwd_xlaunch", xlaunch)
write(*,*) "nonoro_gwd_ran: xlaunch=", xlaunch
! Characteristic vertical scale height
