Changeset 3263 for trunk/LMDZ.MARS


Ignore:
Timestamp:
Mar 12, 2024, 10:43:26 PM (8 months ago)
Author:
jliu
Message:

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

Location:
trunk/LMDZ.MARS
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.MARS/changelog.txt

    r3261 r3263  
    45494549Reversing -r 3250 (the ""bug"" I fixed was not a bug but my mystake)
    45504550
     4551== 12/03/2024 ==Jliu
     4552Synchronize all default non-orographic gravity wave parameters in callphys.def with the ones in scheme. The parameters are as follows:
     4553#Eliassen-Palm FLux(only if calllott_nonoro=.true.)
     4554nonoro_gwd_epflux_max=5.E-4
     4555# saturation parameter non-orographic gravity waves(only if calllott_nonoro=.true.)
     4556nonoro_gwd_sat=1.5
     4557# maximum of the phase speed of the gravity waves(only if calllott_nonoro=.true.)
     4558nonoro_gwd_cmax=50
     4559# value of the dissaption coefficiet(only if calllott_nonoro=.true.)
     4560nonoro_gwd_rdiss=0.07
     4561# value of the max wave number
     4562nonoro_gwd_kmax=1.E-4
     4563# value of the min wave number
     4564nonoro_gwd_kmin=7.E-6
     4565# value to control the launch altitude
     4566nonoro_gwd_xlaunch=0.6
     4567# effective factor
     4568nonoro_gwimixing_eff=0.1
     4569# effective factor
     4570nonoro_gwimixing_eff1=0.1
     4571# mixing vertical decrease rate
     4572nonoro_gwimixing_vdl=1.5
     4573
     4574

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

    r3144 r3263  
    199199     IF (firstcall) THEN
    200200        write(*,*) "nonoro_gwd_mix: non-oro GW-induced mixing is active!"
    201         epflux_max = 7.E-7 ! Mars' value !!
     201        epflux_max = 5.E-4 ! Mars' value !!
    202202        call getin_p("nonoro_gwd_epflux_max", epflux_max)
    203203        write(*,*) "NONORO_GWD_MIX: epflux_max=", epflux_max
    204         sat = 1. ! default gravity waves saturation value !!
     204        sat = 1.5 ! default gravity waves saturation value !!
    205205        call getin_p("nonoro_gwd_sat", sat)
    206206        write(*,*) "NONORO_GWD_MIX: sat=", sat     
    … …  
    208208!        call getin_p("nonoro_gwd_cmax", cmax)
    209209!        write(*,*) "NONORO_GWD_MIX: cmax=", cmax
    210         rdiss=0.01 ! default coefficient of dissipation !!
     210        rdiss=0.07 ! default coefficient of dissipation !!
    211211        call getin_p("nonoro_gwd_rdiss", rdiss)
    212212        write(*,*) "NONORO_GWD_MIX: rdiss=", rdiss
    213         kmax=7.E-4 ! default Max horizontal wavenumber !!
     213        kmax=1.E-4 ! default Max horizontal wavenumber !!
    214214        call getin_p("nonoro_gwd_kmax", kmax)
    215215        write(*,*) "NONORO_GWD_MIX: kmax=", kmax
    216         kmin=2.e-5 ! default Max horizontal wavenumber !!
     216        kmin=7.E-6 ! default Max horizontal wavenumber !!
    217217        call getin_p("nonoro_gwd_kmin", kmin)
    218218        write(*,*) "NONORO_GWD_MIX: kmin=", kmin
    219         xlaunch=0.4 ! default GW luanch altitude controller !!
     219        xlaunch=0.6 ! default GW luanch altitude controller !!
    220220        call getin_p("nonoro_gwd_xlaunch", xlaunch)
    221221        write(*,*) "NONORO_GWD_MIX: xlaunch=", xlaunch
    222         eff=5.E-3 ! Diffusion effective factor !!
     222        eff=0.1 ! Diffusion effective factor !!
    223223        call getin_p("nonoro_gwimixing_eff", eff)
    224224        write(*,*) "NONORO_GWD_MIX: eff=", eff
    225         eff1=1.0 ! Diffusion effective factor !!
     225        eff1=0.1 ! Diffusion effective factor !!
    226226        call getin_p("nonoro_gwimixing_eff1", eff1)
    227227        write(*,*) "NONORO_GWD_MIX: eff1=", eff1

  • trunk/LMDZ.MARS/libf/phymars/nonoro_gwd_ran_mod.F90

    r2934 r3263  
    168168     IF (firstcall) THEN
    169169        write(*,*) "nonoro_gwd_ran: FLott non-oro GW scheme is active!"
    170         epflux_max = 7.E-7 ! Mars' value !!
     170        epflux_max = 5.E-4 ! Mars' value !!
    171171        call getin_p("nonoro_gwd_epflux_max", epflux_max)
    172172        write(*,*) "nonoro_gwd_ran: epflux_max=", epflux_max
    173         sat = 1. ! default gravity waves saturation value !!
     173        sat = 1.5 ! default gravity waves saturation value !!
    174174        call getin_p("nonoro_gwd_sat", sat)
    175175        write(*,*) "nonoro_gwd_ran: sat=", sat     
    … …  
    177177!        call getin_p("nonoro_gwd_cmax", cmax)
    178178!        write(*,*) "nonoro_gwd_ran: cmax=", cmax
    179         rdiss=0.01 ! default coefficient of dissipation !!
     179        rdiss=0.07 ! default coefficient of dissipation !!
    180180        call getin_p("nonoro_gwd_rdiss", rdiss)
    181181        write(*,*) "nonoro_gwd_ran: rdiss=", rdiss
    182         kmax=7.E-4 ! default Max horizontal wavenumber !!
     182        kmax=1.E-4 ! default Max horizontal wavenumber !!
    183183        call getin_p("nonoro_gwd_kmax", kmax)
    184184        write(*,*) "nonoro_gwd_ran: kmax=", kmax
    185         kmin=2.e-5 ! default Max horizontal wavenumber !!
     185        kmin=7.E-6 ! default Max horizontal wavenumber !!
    186186        call getin_p("nonoro_gwd_kmin", kmin)
    187187        write(*,*) "nonoro_gwd_ran: kmin=", kmin
    188         xlaunch=0.4 ! default GW luanch altitude controller !!
     188        xlaunch=0.6 ! default GW luanch altitude controller !!
    189189        call getin_p("nonoro_gwd_xlaunch", xlaunch)
    190190        write(*,*) "nonoro_gwd_ran: xlaunch=", xlaunch
Note: See TracChangeset for help on using the changeset viewer.