Changeset 2622 for trunk/LMDZ.VENUS/libf/phyvenus/nonoro_gwd_ran_mod.F90

Timestamp:

Feb 16, 2022, 4:26:05 PM (3 years ago)

Author:

slebonnois

Message:

SL: VENUS update (i) bug correction (2 bugs, phytrac and physiq), affected meam molec mass computations... (ii) updates for VCD 2.0 (iii) aeropacity: for latitudinal variations of the cloud distribution

File:

• trunk/LMDZ.VENUS/libf/phyvenus/nonoro_gwd_ran_mod.F90 (modified) (7 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/nonoro_gwd_ran_mod.F90

@@ -96 +96 @@
 ! generic : kmin=1/sqrt(DxDy) Dx and Dy horizontal grid
 
-!     REAL, parameter:: cmax = 61.   ! Max horizontal absolute phase velocity  TestGW6
+!----------------------------------------
+! VCD 1.1 tuning
+!     REAL, parameter:: cmax = 61.   ! Max horizontal absolute phase velocity
 !----------------------------------------
 ! VCD 2.0 tuning
-     REAL, parameter:: cmax = 301.   ! Max horizontal absolute phase velocity  TestGW6
-!----------------------------------------
-
+      REAL, parameter:: cmax = 61.   ! Max horizontal absolute phase velocity
+!----------------------------------------
+    
 ! generic : cmax=zonal wind at launch
      REAL, parameter:: cmin = 1.    ! Min horizontal absolute phase velocity
@@ -126 +128 @@
 
      ! 0.3.2 Parameters of waves dissipations
-     REAL, parameter:: sat   = 0.85   ! saturation parameter
-     REAL, parameter:: rdiss = 0.1    ! coefficient of dissipation
+!----------------------------------------
+! VCD 1.1 tuning
+!    REAL, parameter:: sat   = 0.85   ! saturation parameter
+!    REAL, parameter:: rdiss = 0.1    ! coefficient of dissipation
+!----------------------------------------
+! VCD 2.0 tuning
+     REAL, parameter:: sat   = 0.6    ! saturation parameter
+     REAL, parameter:: rdiss = 8.e-4  ! coefficient of dissipation
+!----------------------------------------
      REAL, parameter:: zoisec = 1.e-8 ! security for intrinsic freguency
 
      ! 0.3.3 Background flow at 1/2 levels and vertical coordinate
      REAL H0bis(ngrid, nlayer)       ! characteristic height of the atmosphere
+     REAL min_k(ngrid)               ! min(kstar,kmin)
      REAL, save::  H0                ! characteristic height of the atmosphere
-     REAL, parameter:: pr = 5e5      ! Reference pressure [Pa]   ! VENUS: cloud layer
+     REAL, save::  MDR               ! characteristic mass density  
+!----------------------------------------
+! VCD 1.1 tuning
+!    REAL, parameter:: pr = 5e5      ! Reference pressure [Pa]   ! VENUS: cloud layer
+!----------------------------------------
+! VCD 2.0 tuning
+     REAL, parameter:: pr = 5e4      ! Reference pressure [Pa]   ! VENUS: cloud layer
+!----------------------------------------
      REAL, parameter:: tr = 300.     ! Reference temperature [K] ! VENUS: cloud layer
      REAL zh(ngrid, nlayer + 1)      ! Log-pressure altitude (constant H0) 
@@ -156 +173 @@
         ! Characteristic vertical scale height
         H0 = RD * tr / RG
+        ! Characteristic mass density at launch altitude
+        MDR = pr / ( RD * tr )
         ! Control
         if (deltat .LT. dtime) THEN
@@ -214 +233 @@
      uh(:, nlayer + 1) = uu(:, nlayer)
      vh(:, nlayer + 1) = vv(:, nlayer)
+
+     ! TN+GG April/June 2020
+     ! "Individual waves are not supposed
+     ! to occupy the entire domain, but
+     ! only a faction of it" Lott+2012
+     ! minimum value of k between kmin and kstar
+     DO ii = 1, ngrid
+        kstar = RPI / sqrt(cell_area(ii)) 
+        min_k(ii) = max(kmin,kstar)
+     end DO
 
 ! 3. WAVES CHARACTERISTICS CHOSEN RANDOMLY
@@ -235 +264 @@
            ! horizontal wavenumber amplitude
 !           zk(jw, ii) = kmin + (kmax - kmin) * ran_num_2
-           ! TN+GG April/June 2020
-           ! "Individual waves are not supposed
-           ! to occupy the entire domain, but
-           ! only a faction of it" Lott+2012
-           kstar = RPI / sqrt(cell_area(ii))
-           zk(jw, ii) = max(kmin,kstar) + (kmax - max(kmin,kstar)) * ran_num_2
+           zk(jw, ii) = min_k(ii) + (kmax - min_k(ii)) * ran_num_2
 
            ! horizontal phase speed
@@ -320 +344 @@
            ! Intrinsic frequency
            intr_freq_p(jw, :) = zo(jw, :) - zk(jw, :) * cos(zp(jw, :)) * uh(:, ll + 1)    &
-                      - zk(jw, :) * sin(zp(jw, :)) * vh(:, ll + 1)
+                                          - zk(jw, :) * sin(zp(jw, :)) * vh(:, ll + 1)
            ! Vertical velocity in flux formulation
            wwp(jw, :) = min(                                                              &
@@ -336 +360 @@
                       * abs(intr_freq_p(jw, :))**3 / bv(:, ll + 1)                        &
                       * exp(-zh(:, ll + 1) / H0)                                          &
-                      * sat**2 * kmin**2 / zk(jw, :)**4)
+!! Correction for VCD 2.0
+!                      * sat**2 * kmin**2 / zk(jw, :)**4)
+                      * sat**2 * min_k(:)**2 * MDR / zk(jw, :)**4)                       
         end DO