Changeset 2640 for trunk/LMDZ.MARS

Timestamp:

Mar 9, 2022, 5:02:26 PM (3 years ago)

Author:

emillour

Message:

Mars GCM:
Update for non-orographic Gravity Waves; some key parameters can now be
set via callphys.def
JL

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/nonoro_gwd_ran_mod.F90 (modified) (5 diffs)

trunk/LMDZ.MARS/README

@@ -3643 +3643 @@
 tracers.
 
-== 03/03/2022 == AB
+== 09/03/2022 == AB
 Update .def* files in LMDZ.MARS/deftank/
 *.def.GCM6 files are bound to change later but are ready to use.
 
-== 03/03/2022 == AB
+== 09/03/2022 == AB
 Following r2637, change default values for dust cycle parameters to make them
 consistent with callphys.def.GCM6
+
+== 09/03/2022 == JL
+Update for non-orographic Gravity Waves; some key parameters can now be
+set via callphys.def

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

@@ -104 +104 @@
     REAL, save :: kmax             ! Max horizontal wavenumber (lambda min,lambda=2*pi/kmax),kmax=N/u, u(=30~50) zonal wind velocity at launch altitude
 !$OMP THREADPRIVATE(kmax)
-    REAL, parameter:: kmin = 2.e-5 ! Min horizontal wavenumber (lambda max = 314 km,lambda=2*pi/kmin)
+    REAL, save :: kmin             ! Min horizontal wavenumber (lambda max = 314 km,lambda=2*pi/kmin)
+!$OMP THREADPRIVATE(kmim)
     REAL kstar                     ! Min horizontal wavenumber constrained by horizontal resolution
     REAL, save :: cmax             ! Max horizontal absolute phase velocity (at the maginitide of zonal wind u at the launch altitude)
 !$OMP THREADPRIVATE(cmax)
+    REAL :: max_k(ngrid)           ! max_k=max(kstar,kmin)
     REAL, parameter:: cmin = 1.    ! Min horizontal absolute phase velocity (not used)    
     REAL CPHA                      ! absolute PHASE VELOCITY frequency
@@ -126 +128 @@
 !$OMP THREADPRIVATE(epflux_max)
     INTEGER LAUNCH                       ! Launching altitude
-    REAL, parameter:: xlaunch = 0.4      ! Control the launching altitude by pressure
+    REAL, save :: xlaunch                ! Control the launching altitude by pressure
+!$OMP THREADPRIVATE(xlaunch)
     REAL, parameter:: zmaxth_top = 8000. ! Top of convective layer (approx. not used)
 
@@ -180 +183 @@
         call getin_p("nonoro_gwd_kmax", kmax)
         write(*,*) "nonoro_gwd_ran: kmax=", kmax
-
+        kmin=2.e-5 ! 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 !!
+        call getin_p("nonoro_gwd_xlaunch", xlaunch)
+        write(*,*) "nonoro_gwd_ran: xlaunch=", xlaunch
         ! Characteristic vertical scale height
         H0 = r * tr / g
@@ -264 +272 @@
     VH(:, nlayer + 1) = VV(:, nlayer)
 
+    DO ii=1,ngrid
+       KSTAR = PI/SQRT(cell_area(II))
+       MAX_K(II)=MAX(kmin,kstar)
+    ENDDO
     call writediagfi(ngrid,'nonoro_bv','nonoro_bv', 'm',3,BV(:,2:nlayer+1))
 
@@ -282 +294 @@
                 ! "Individual waves are not supposed to occupy the entire domain, but only a fraction of it" (Lott+2012)
                 ! ZK(JW, II) = KMIN + (KMAX - KMIN) *RAN_NUM_2
-                KSTAR = PI/SQRT(cell_area(II))
-                ZK(JW, II) = MAX(KMIN,KSTAR) + (KMAX - MAX(KMIN,KSTAR)) *RAN_NUM_2
+                KSTAR = PI/SQRT(cell_area(II))          
+                ZK(JW, II) = MAX_K(II) + (KMAX - MAX_K(II)) *RAN_NUM_2
                
                ! Horizontal phase speed