Changeset 3754

Timestamp:

May 5, 2025, 5:23:50 PM (7 weeks ago)

Author:

afalco

Message:

Pluto: imported orographic gravity waves from Mars.
AF

Location:

trunk

Files:

• LMDZ.MARS/libf/phymars/sugwd.F90 (modified) (1 diff)
• LMDZ.PLUTO/libf/phypluto/calldrag_noro_mod.F90 (added)
• LMDZ.PLUTO/libf/phypluto/callkeys_mod.F90 (modified) (2 diffs)
• LMDZ.PLUTO/libf/phypluto/dimphy.F90 (modified) (2 diffs)
• LMDZ.PLUTO/libf/phypluto/drag_noro_mod.F90 (added)
• LMDZ.PLUTO/libf/phypluto/gwprofil_mod.F90 (added)
• LMDZ.PLUTO/libf/phypluto/gwstress_mod.F90 (added)
• LMDZ.PLUTO/libf/phypluto/inifis_mod.F90 (modified) (1 diff)
• LMDZ.PLUTO/libf/phypluto/orodrag_mod.F90 (added)
• LMDZ.PLUTO/libf/phypluto/orosetup.F90 (added)
• LMDZ.PLUTO/libf/phypluto/phys_state_var_mod.F90 (modified) (2 diffs)
• LMDZ.PLUTO/libf/phypluto/physiq_mod.F90 (modified) (4 diffs)
• LMDZ.PLUTO/libf/phypluto/sugwd.F90 (copied) (copied from trunk/LMDZ.MARS/libf/phymars/sugwd.F90) (4 diffs)
• LMDZ.PLUTO/libf/phypluto/yoegwd_h.F90 (added)

trunk/LMDZ.MARS/libf/phymars/sugwd.F90 ¶

@@ -1 +1 @@
       SUBROUTINE SUGWD(nlayer,sigtest)
-! ==============================================================================    
-!     Initialize common variables in yoegwd.h to control the orographic 
-!     graivty wave drag parameterization. That means, all the tunable parameters 
+! ==============================================================================
+!     Initialize common variables in yoegwd.h to control the orographic
+!     gravity wave drag parameterization. That means, all the tunable parameters
 !     for oro-GW scheme are in this subroutine.
 !     MARTIN MILLER             *ECMWF*               ORIGINAL : 90-01-01  

trunk/LMDZ.PLUTO/libf/phypluto/callkeys_mod.F90 ¶

@@ -4 +4 @@
       logical,save :: callrad,corrk,calldifv,UseTurbDiff
 !$OMP THREADPRIVATE(callrad,corrk,calldifv,UseTurbDiff)
-      logical,save :: calladj,calltherm,n2cond,callsoil
-!$OMP THREADPRIVATE(calladj,calltherm,n2cond,callsoil)
+      logical,save :: calladj,calltherm,n2cond,callsoil,calllott
+!$OMP THREADPRIVATE(calladj,calltherm,n2cond,callsoil,calllott)
       logical,save :: callconduct,callmolvis,callmoldiff
 !$OMP THREADPRIVATE(callconduct,callmolvis,callmoldiff)
@@ -199 +199 @@
       real,save    :: deltap    ! width of transition to alpha_top (Pa)
 !$OMP THREADPRIVATE(alpha_top,pref,deltap)
-      
+
 !! Microphysics-specific variables
       logical,save :: callmufi, call_haze_prod_pCH4

trunk/LMDZ.PLUTO/libf/phypluto/dimphy.F90 ¶

@@ -1 +1 @@
 MODULE dimphy
-  
+
   INTEGER,SAVE :: klon   ! number of atmospheric columns (for this OpenMP subgrid)
   INTEGER,SAVE :: klev   ! number of atmospheric layers, read by master
@@ -6 +6 @@
   INTEGER,SAVE :: klevm1 ! number of atmospheric layers-1, read by master
 !  INTEGER,SAVE :: kflev
+  integer,save :: ndomainsz !=(ngrid-1)/20 + 1
 
 !$OMP THREADPRIVATE(klon)
 
 CONTAINS
-  
+
   SUBROUTINE Init_dimphy(klon0,klev0)
   IMPLICIT NONE
-  
+
     INTEGER, INTENT(in) :: klon0
     INTEGER, INTENT(in) :: klev0
-    
+
     klon=klon0
-    
-!$OMP MASTER 
+
+!$OMP MASTER
     klev=klev0
     klevp1=klev+1
     klevm1=klev-1
 !    kflev=klev
-!$OMP END MASTER    
+!$OMP END MASTER
 !$OMP BARRIER
-    
+
   END SUBROUTINE Init_dimphy
 
-  
+
 END MODULE dimphy

trunk/LMDZ.PLUTO/libf/phypluto/inifis_mod.F90 ¶

@@ -355 +355 @@
      if (is_master) write(*,*) trim(rname)//": callsoil = ",callsoil
 
+     calllott=.true. ! default value
+     call getin_p("calllott",calllott)
+     write(*,*)" calllott = ",calllott
+
      if (is_master) write(*,*)trim(rname)//&
        ": Rad transfer is computed every iradia", &

trunk/LMDZ.PLUTO/libf/phypluto/phys_state_var_mod.F90 ¶

@@ -9 +9 @@
 !======================================================================
 ! Declaration des variables
-      USE dimphy, only : klon,klev
+      USE dimphy, only : klon,klev,ndomainsz
       USE comsoil_h, only : nsoilmx
       use comsaison_h, only: mu0, fract
@@ -117 +117 @@
 !rugoro(:) ! longueur de rugosite de l'OESM
 
+        ndomainsz=(klon-1)/20 + 1
         ALLOCATE(phisfi(klon))
         ALLOCATE(tsurf(klon))

trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90 ¶

@@ -44 +44 @@
                             obliquit, z0, adjust, tpal
       use comcstfi_mod, only: pi, g, rcp, r, rad, mugaz, cpp
+      use calldrag_noro_mod, only: calldrag_noro
       use time_phylmdz_mod, only: daysec
       use callkeys_mod, only: albedo_spectral_mode, calladj, calldifv,        &
-                              callrad, callsoil, nosurf,                      &
+                              calllott, callrad, callsoil, nosurf,            &
                               callconduct,callmolvis,callmoldiff,             &
                               corrk,                                          &
@@ -117 +118 @@
 !         II.2.b Option 2 : Atmosphere has no radiative effect.
 !
+!      II.3 Gravity wave and subgrid scale topography drag :
+!
 !      III. Vertical diffusion (turbulent mixing)
 !
@@ -385 +388 @@
       real zdhdif(ngrid,nlayer)                             ! Turbdiff/vdifc routines.
       real zdhadj(ngrid,nlayer)                             ! Convadj routine.
+      REAL zdtgw(ngrid,nlayer)                              ! Gravity waves (K/s)
+      REAL zdugw(ngrid,nlayer),zdvgw(ngrid,nlayer)          ! Gravity waves (m.s-2)
       REAL zdvc(ngrid,nlayer),zduc(ngrid,nlayer)            ! condense_n2 routine.
 
@@ -1143 +1148 @@
    endif
 
+
+!-----------------------------------------------------------------------
+! II.3 Gravity wave and subgrid scale topography drag :
+!    -------------------------------------------------
+
+      IF(calllott)THEN
+        CALL calldrag_noro(ngrid,nlayer,ptimestep, &
+                           zplay,zplev,pt,pu,pv,zdtgw,zdugw,zdvgw)
+
+        DO l=1,nlayer
+          DO ig=1,ngrid
+            pdv(ig,l)=pdv(ig,l)+zdvgw(ig,l)
+            pdu(ig,l)=pdu(ig,l)+zdugw(ig,l)
+            pdt(ig,l)=pdt(ig,l)+zdtgw(ig,l)
+          ENDDO
+        ENDDO
+      ENDIF
+
+
 !  --------------------------------------------
 !  III. Vertical diffusion (turbulent mixing) :

trunk/LMDZ.PLUTO/libf/phypluto/sugwd.F90 ¶

@@ -1 +1 @@
       SUBROUTINE SUGWD(nlayer,sigtest)
-! ==============================================================================    
-!     Initialize common variables in yoegwd.h to control the orographic 
-!     graivty wave drag parameterization. That means, all the tunable parameters 
+! ==============================================================================
+!     Initialize common variables in yoegwd.h to control the orographic
+!     gravity wave drag parameterization. That means, all the tunable parameters
 !     for oro-GW scheme are in this subroutine.
-!     MARTIN MILLER             *ECMWF*               ORIGINAL : 90-01-01  
+!     MARTIN MILLER             *ECMWF*               ORIGINAL : 90-01-01
 !     Update:    Jiandong Liu     2022/03/15          Rewirite into .F90 and
-!                                                     comment.   
+!                                                     comment.
 !     REFERENCE.
 !     ----------
@@ -20 +20 @@
       implicit none
 
-      ! 0.1 Inputs:         
+      ! 0.1 Inputs:
       integer,intent(in):: nlayer                ! Number of model levels
       REAL,intent(in):: sigtest(nlayer+1)        ! Vertical coordinate table
@@ -33 +33 @@
       integer jk
 
-!-------------------------------------------------------------------------------     
-! 1.   Set the values of the parameters 
-!-------------------------------------------------------------------------------   
+!-------------------------------------------------------------------------------
+! 1.   Set the values of the parameters
+!-------------------------------------------------------------------------------
 !     PRINT *,' Dans sugwd nlayer=',nlayer,' SIG=',sigtest
       GHMAX=10000.
-   
+
 !     old  ZSIGT=0.94
 !     old  ZPR=80000.
       ZSIGT=0.85      ! Sigmal levels
       ZPR=100000.     ! Surface (Reference) Pressure?
-      
+
       ! ! Condition to find NKTOPG layer, which NKTOPG is a condition to set
       ! 1*pvar and 2*pvar layers (OROSETUP)
@@ -49 +49 @@
          ZPM1R=0.5*ZPR*(sigtest(JK)+sigtest(JK+1))
          IF((ZPM1R/ZPR).GE.ZSIGT)THEN
-            NKTOPG=JK   
+            NKTOPG=JK
          ENDIF
       ENDDO
       WRITE(*,*) 'In sugwd NKTOPG=',NKTOPG
-    
-      GSIGCR=0.80  ! Sigmal levels to found the top of low level flow height (OROSETUP)  
+
+      GSIGCR=0.80  ! Sigmal levels to found the top of low level flow height (OROSETUP)
       GKDRAG= 0.1  ! used to be 0.1 for mcd Version 1 and 2 (before 10/2000) (OROSETUP)
-                      
+
       GFRCRIT=1.0
       GKWAKE=1.0   ! The G in equation (16)
-      GRCRIT=0.25  ! Critical value for Mean flow richardson number(OROSETUP)      
+      GRCRIT=0.25  ! Critical value for Mean flow richardson number(OROSETUP)
       GKDRAGL=4.*GKDRAG
       GRAHILO=1.
-      GVCRIT =0.0    
-!-------------------------------------------------------------------------------  
-! 2.    Set values of security parameters 
-!-------------------------------------------------------------------------------  
+      GVCRIT =0.0
+!-------------------------------------------------------------------------------
+! 2.    Set values of security parameters
+!-------------------------------------------------------------------------------
       GVSEC=0.10      ! Security values for For normal wind (pu^2+pv^2)^0.5(OROSETUP,GWSTRESS)
-      GSSEC=1.E-12    ! Security values for Brunt–Väisälä frequency N^2 (OROSETUP)   
+      GSSEC=1.E-12    ! Security values for Brunt–Väisälä frequency N^2 (OROSETUP)
       GTSEC=1.E-07    ! Security values for Sub-grid scale anisotropy(OROSETUP,GWSTRESS)