Changeset 2149 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Jul 19, 2019, 11:40:36 AM (5 years ago)

Author:

emillour

Message:

Mars GCM:
Add F.Lott's non-orographic GW parametrization. Disabled by default for now, activated by setting calllott_nonoro=.true. in callphys.def
GG+EM

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• callkeys.h (modified) (2 diffs)
• conf_phys.F (modified) (1 diff)
• nonoro_gwd_ran_mod.F90 (added)
• physiq_mod.F (modified) (14 diffs)
• yoegwd.h (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/callkeys.h

@@ -15 +15 @@
      &   ,calltherm,callrichsl,callslope,tituscap,callyamada4,co2clouds &
      &   ,co2useh2o,meteo_flux,CLFvaryingCO2,spantCO2,CLFvarying        &
-     &   ,satindexco2,rdstorm
+     &   ,satindexco2,rdstorm,calllott_nonoro
      
       COMMON/callkeys_i/iradia,iaervar,iddist,ilwd,ilwb,ilwn,ncouche    &
@@ -25 +25 @@
      
       LOGICAL callrad,calldifv,calladj,callcond,callsoil,               &
-     &   season,diurnal,lwrite,calllott                                 &
+     &   season,diurnal,lwrite,calllott,calllott_nonoro                 &
      &   ,callstats,calleofdump                                         &
      &   ,callnirco2,callnlte,callthermos,callconduct,                  &

trunk/LMDZ.MARS/libf/phymars/conf_phys.F

@@ -280 +280 @@
          write(*,*)" calllott = ",calllott
 
+         write(*,*)"call Lott's non-oro GWs parameterisation ",
+     &             "scheme ?"
+         calllott_nonoro=.false. ! default value
+         call getin("calllott_nonoro",calllott_nonoro)
+         write(*,*)" calllott_nonoro = ",calllott_nonoro
+
 ! rocket dust storm injection scheme
          write(*,*)"call rocket dust storm and slope lifting", 

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -62 +62 @@
       use iono_h, only: allocate_param_iono
       use compute_dtau_mod, only: compute_dtau
+      use nonoro_gwd_ran_mod, only: nonoro_gwd_ran
 #ifdef MESOSCALE
       use comsoil_h, only: mlayer,layer
@@ -414 +415 @@
       real*8  varerr
 
+C  Non-oro GW drag & Calcul of Brunt-Vaisala freq. (BV2) 
+      REAL ztetalev(ngrid,nlayer)
+      real zdtetalev(ngrid,nlayer), zdzlev(ngrid,nlayer)
+      REAL bv2(ngrid,nlayer)    ! BV2 at zlev   
+c  Non-oro GW tendencies
+      REAL d_u_hin(ngrid,nlayer), d_v_hin(ngrid,nlayer)
+      REAL d_t_hin(ngrid,nlayer)
+c  Diagnostics 2D of gw_nonoro
+      REAL zustrhi(ngrid), zvstrhi(ngrid)
+
 c Variables for PBL
       REAL zz1(ngrid)
@@ -1270 +1281 @@
       ENDIF ! of IF(calladj)
 
+c-----------------------------------------------------
+c    8. Non orographic Gravity waves :
+c -------------------------------------------------
+
+      IF (calllott_nonoro) THEN
+
+         CALL nonoro_gwd_ran(ngrid,nlayer,ptimestep,zplay,
+     &               zmax_th,                      ! max altitude reached by thermals (m)
+     &               pt, pu, pv,
+     &               pdt, pdu, pdv,
+     &               zustrhi,zvstrhi,
+     &               d_t_hin, d_u_hin, d_v_hin)
+
+!  Update tendencies
+         pdt(1:ngrid,1:nlayer)=pdt(1:ngrid,1:nlayer)
+     &                         +d_t_hin(1:ngrid,1:nlayer)
+!        d_t_hin(:,:)= d_t_hin(:,:)/ptimestep ! K/s (for outputs?)
+         pdu(1:ngrid,1:nlayer)=pdu(1:ngrid,1:nlayer)
+     &                         +d_u_hin(1:ngrid,1:nlayer) 
+!        d_u_hin(:,:)= d_u_hin(:,:)/ptimestep ! (m/s)/s (for outputs?)
+         pdv(1:ngrid,1:nlayer)=pdv(1:ngrid,1:nlayer)
+     &                         +d_v_hin(1:ngrid,1:nlayer)
+!        d_v_hin(:,:)= d_v_hin(:,:)/ptimestep ! (m/s)/s (for outputs?)
+
+      ENDIF ! of IF (calllott_nonoro)
 
 c-----------------------------------------------------------------------
-c   8. Specific parameterizations for tracers 
+c   9. Specific parameterizations for tracers 
 c:   -----------------------------------------
 
       if (tracer) then 
 
-c   8a. Water and ice
+c   9a. Water and ice
 c     ---------------
 
@@ -1357 +1393 @@
          END IF  ! of IF (water)
 
-c   8a bis. CO2 clouds (CL & JA)
+c   9a bis. CO2 clouds (CL & JA)
 c        ---------------------------------------
 c        CO2 ice cloud condensation in the atmosphere
@@ -1479 +1515 @@
       END IF                    ! of IF (co2clouds)
 
-c   8b. Aerosol particles
+c   9b. Aerosol particles
 c     -------------------
 c        ----------
@@ -1571 +1607 @@
       ENDIF
 c
-c   8c. Chemical species
+c   9c. Chemical species
 c     ------------------
 
@@ -1632 +1668 @@
 #endif
 
-c   8d. Updates
+c   9d. Updates
 c     ---------
 
@@ -1650 +1686 @@
 #ifndef MESOSCALE
 c-----------------------------------------------------------------------
-c   9. THERMOSPHERE CALCULATION
+c   10. THERMOSPHERE CALCULATION
 c-----------------------------------------------------------------------
 
@@ -1675 +1711 @@
 #endif
 c-----------------------------------------------------------------------
-c   10. Carbon dioxide condensation-sublimation:
+c   11. Carbon dioxide condensation-sublimation:
 c     (should be the last atmospherical physical process to be computed)
 c   -------------------------------------------
@@ -1748 +1784 @@
 
 c-----------------------------------------------------------------------
-c   11. Surface  and sub-surface soil temperature
+c   12. Surface  and sub-surface soil temperature
 c-----------------------------------------------------------------------
 c
 c
-c   11.1 Increment Surface temperature:
+c   12.1 Increment Surface temperature:
 c   ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
@@ -1793 +1829 @@
 
 c
-c   11.2 Compute soil temperatures and subsurface heat flux:
+c   12.2 Compute soil temperatures and subsurface heat flux:
 c   ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
       IF (callsoil) THEN
@@ -1841 +1877 @@
 
 c-----------------------------------------------------------------------
-c  12. Write output files
+c  13. Write output files
 c  ----------------------
 
@@ -2539 +2575 @@
          call WRITEDIAGFI(ngrid,"v","Meridional wind","m.s-1",3,zv)
          call WRITEDIAGFI(ngrid,"w","Vertical wind","m.s-1",3,pw)
-         call WRITEDIAGFI(ngrid,"rho","density","none",3,rho)
+         call WRITEDIAGFI(ngrid,"rho","density","kg.m-3",3,rho)
 c        call WRITEDIAGFI(ngrid,"q2","q2","kg.m-3",3,q2)
 c        call WRITEDIAGFI(ngrid,'Teta','T potentielle','K',3,zh)
@@ -2878 +2914 @@
          end if  ! (tracer.and.(dustbin.ne.0))
 
+c        ----------------------------------------------------------
+c        GW non-oro outputs
+c        ----------------------------------------------------------
+
+         if(calllott_nonoro) then 
+           call WRITEDIAGFI(ngrid,"dugwno","GW non-oro dU","m/s2",
+     $           3,d_u_hin/ptimestep)
+           call WRITEDIAGFI(ngrid,"dvgwno","GW non-oro dV","m/s2",
+     $           3,d_v_hin/ptimestep)
+         endif                  !(calllott_nonoro)
 
 c        ----------------------------------------------------------

trunk/LMDZ.MARS/libf/phymars/yoegwd.h

@@ -1 @@
-C     -----------------------------------------------------------------
-C*    *COMMON* *YOEGWD* - PARAMETERS FOR GRAVITY WAVE DRAG CALCULATIONS
-C     -----------------------------------------------------------------
-C
-      real GFRCRIT,GRCRIT,GVCRIT,GKDRAG,GKDRAGL,GHMAX
-     *             , GRAHILO,GSIGCR,GSSEC, GTSEC,  GVSEC
-     *             ,GKWAKE
-
-      integer NKTOPG
-
-      COMMON/YOEGWD/ GFRCRIT,GRCRIT,GVCRIT,GKDRAG,GKDRAGL,GHMAX
-     *             , GRAHILO,GSIGCR,NKTOPG,GSSEC, GTSEC,  GVSEC
-     *             ,GKWAKE
-C
+!     -----------------------------------------------------------------
+!*    *COMMON* *YOEGWD* - PARAMETERS FOR GRAVITY WAVE DRAG CALCULATIONS
+!     -----------------------------------------------------------------
+!
+      real    ::  GFRCRIT,GRCRIT,GVCRIT,GKDRAG,GKDRAGL,GHMAX
+      real    ::   GRAHILO,GSIGCR,GSSEC, GTSEC,  GVSEC, GKWAKE
+      integer ::  NKTOPG
+      COMMON/YOEGWD/ GFRCRIT,GRCRIT,GVCRIT,GKDRAG,GKDRAGL,GHMAX         &
+     &             ,GRAHILO,GSIGCR,NKTOPG,GSSEC, GTSEC, GVSEC,GKWAKE
 
 
+