Changeset 2682

Timestamp:

May 19, 2022, 2:29:25 PM (3 years ago)

Author:

emillour

Message:

Mars GCM:
Update Non-Orographic GW scheme: cmax is not longer fixed but evaluated
as the zonal wind at launch altitude.
JL

Location:

trunk/LMDZ.MARS

Files:

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

trunk/LMDZ.MARS/README

@@ -3711 +3711 @@
 + declare new variable "flux_1AU" for solar constant.
 
+== 19/05/2022 == JL
+Upate Non-Orographic GW scheme: cmax is not longer fixed but evaluated
+as the zonal wind at launch altitude.

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

@@ -107 +107 @@
 !$OMP THREADPRIVATE(kmin)
     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
+    REAL CPHA(ngrid)                      ! absolute PHASE VELOCITY frequency
     REAL ZK(NW, ngrid)             ! Horizontal wavenumber amplitude
     REAL ZP(NW, ngrid)             ! Horizontal wavenumber angle        
@@ -131 +130 @@
 !$OMP THREADPRIVATE(xlaunch)
     REAL, parameter:: zmaxth_top = 8000. ! Top of convective layer (approx. not used)
-
+    REAL cmax(ngrid,nlayer)             ! Max horizontal absolute phase velocity (at the maginitide of zonal wind u at the launch altitude)
 
     ! 0.3.2 PARAMETERS OF WAVES DISSIPATIONS
@@ -174 +173 @@
         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
+!        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 !!
         call getin_p("nonoro_gwd_rdiss", rdiss)
@@ -271 +270 @@
     UH(:, nlayer + 1) = UU(:, nlayer)
     VH(:, nlayer + 1) = VV(:, nlayer)
-
+    cmax(:,launch)=UU(:, launch)
     DO ii=1,ngrid
        KSTAR = PI/SQRT(cell_area(II))
@@ -300 +299 @@
 ! this computation allows to get a gaussian distribution centered on 0 (right ?)
 ! then cmin is not useful, and it favors small values.
-                CPHA = 0.
+                CPHA(:) = 0.
                 DO JJ = 1, NA
                     RAN_NUM_3=MOD(TT(II, JW+3*JJ)**2, 1.)
-                    CPHA = CPHA + 2.*CMAX*                         &
+                    CPHA(ii) = CPHA(ii) + 2.*CMAX(ii,launch)*      &
                            (RAN_NUM_3 -0.5)*                       &
                            SQRT(3.)/SQRT(NA*1.)
                 END DO
-                IF (CPHA.LT.0.)  THEN
-                   CPHA = -1.*CPHA
+                IF (CPHA(ii).LT.0.)  THEN
+                   CPHA(ii) = -1.*CPHA(ii)
                    ZP(JW,II) = ZP(JW,II) + PI
                 ENDIF
@@ -316 +315 @@
 
                 ! Intrinsic frequency
-                ZO(JW, II) = CPHA * ZK(JW, II)
+                ZO(JW, II) = CPHA(II) * ZK(JW, II)
                 ! Intrinsic frequency  is imposed
                 ZO(JW, II) = ZO(JW, II)                                            &