phylmd/Orography/sugwd.f90

    SUBROUTINE sugwd(nlon,nlev,paprs,pplay)

!**** *SUGWD* INITIALIZE COMMON YOEGWD CONTROLLING GRAVITY WAVE DRAG

!     PURPOSE.
!     --------
!           INITIALIZE YOEGWD, THE COMMON THAT CONTROLS THE
!           GRAVITY WAVE DRAG PARAMETRIZATION.

!**   INTERFACE.
!     ----------
!        CALL *SUGWD* FROM *SUPHEC*
!              -----        ------

!        EXPLICIT ARGUMENTS :
!        --------------------
!        PSIG        : VERTICAL COORDINATE TABLE
!        NLEV        : NUMBER OF MODEL LEVELS

!        IMPLICIT ARGUMENTS :
!        --------------------
!        COMMON YOEGWD

!     METHOD.
!     -------
!        SEE DOCUMENTATION

!     EXTERNALS.
!     ----------
!        NONE

!     REFERENCE.
!     ----------
!        ECMWF Research Department documentation of the IFS

!     AUTHOR.
!     -------
!        MARTIN MILLER             *ECMWF*

!     MODIFICATIONS.
!     --------------
!        ORIGINAL : 90-01-01
!     ------------------------------------------------------------------
      USE yoegwd
      IMPLICIT NONE

!     -----------------------------------------------------------------
!      ----------------------------------------------------------------

      INTEGER nlon, nlev, jk
      REAL, INTENT (IN) :: paprs(nlon,nlev+1)
      REAL, INTENT (IN) :: pplay(nlon,nlev)
      REAL zpr, zstra, zsigt, zpm1r

!*       1.    SET THE VALUES OF THE PARAMETERS
!              --------------------------------

100   CONTINUE

      PRINT *, ' DANS SUGWD NLEV=', nlev
      ghmax = 10000.

      zpr = 100000.
      zstra = 0.1
      zsigt = 0.94
!old  ZPR=80000.
!old  ZSIGT=0.85

      DO 110 jk = 1, nlev
        zpm1r = pplay(nlon/2,jk)/paprs(nlon/2,1)
        IF (zpm1r>=zsigt) THEN
          nktopg = jk
        END IF
        zpm1r = pplay(nlon/2,jk)/paprs(nlon/2,1)
        IF (zpm1r>=zstra) THEN
          nstra = jk
        END IF
110   CONTINUE

!  inversion car dans orodrag on compte les niveaux a l'envers
      nktopg = nlev - nktopg + 1
      nstra = nlev - nstra
      PRINT *, ' DANS SUGWD nktopg=', nktopg
      PRINT *, ' DANS SUGWD nstra=', nstra

      gsigcr = 0.80

      gkdrag = 0.2
      grahilo = 1.
      grcrit = 0.01
      gfrcrit = 1.0
      gkwake = 0.50

      gklift = 0.50
      gvcrit = 0.0


!      ----------------------------------------------------------------

!*       2.    SET VALUES OF SECURITY PARAMETERS
!              ---------------------------------

200   CONTINUE

      gvsec = 0.10
      gssec = 1.E-12

      gtsec = 1.E-07

!      ----------------------------------------------------------------

      RETURN
    END
1	guez	23	SUBROUTINE sugwd(nlon,nlev,paprs,pplay)
2
3			!**** SUGWD INITIALIZE COMMON YOEGWD CONTROLLING GRAVITY WAVE DRAG
4
5			! PURPOSE.
6			! --------
7			! INITIALIZE YOEGWD, THE COMMON THAT CONTROLS THE
8			! GRAVITY WAVE DRAG PARAMETRIZATION.
9
10			!** INTERFACE.
11			! ----------
12			! CALL SUGWD FROM SUPHEC
13			! ----- ------
14
15			! EXPLICIT ARGUMENTS :
16			! --------------------
17			! PSIG : VERTICAL COORDINATE TABLE
18			! NLEV : NUMBER OF MODEL LEVELS
19
20			! IMPLICIT ARGUMENTS :
21			! --------------------
22			! COMMON YOEGWD
23
24			! METHOD.
25			! -------
26			! SEE DOCUMENTATION
27
28			! EXTERNALS.
29			! ----------
30			! NONE
31
32			! REFERENCE.
33			! ----------
34			! ECMWF Research Department documentation of the IFS
35
36			! AUTHOR.
37			! -------
38			! MARTIN MILLER ECMWF
39
40			! MODIFICATIONS.
41			! --------------
42			! ORIGINAL : 90-01-01
43			! ------------------------------------------------------------------
44			USE yoegwd
45			IMPLICIT NONE
46
47			! -----------------------------------------------------------------
48			! ----------------------------------------------------------------
49
50			INTEGER nlon, nlev, jk
51			REAL, INTENT (IN) :: paprs(nlon,nlev+1)
52			REAL, INTENT (IN) :: pplay(nlon,nlev)
53			REAL zpr, zstra, zsigt, zpm1r
54
55			!* 1. SET THE VALUES OF THE PARAMETERS
56			! --------------------------------
57
58			100 CONTINUE
59
60			PRINT *, ' DANS SUGWD NLEV=', nlev
61			ghmax = 10000.
62
63			zpr = 100000.
64			zstra = 0.1
65			zsigt = 0.94
66			!old ZPR=80000.
67			!old ZSIGT=0.85
68
69			DO 110 jk = 1, nlev
70			zpm1r = pplay(nlon/2,jk)/paprs(nlon/2,1)
71			IF (zpm1r>=zsigt) THEN
72			nktopg = jk
73			END IF
74			zpm1r = pplay(nlon/2,jk)/paprs(nlon/2,1)
75			IF (zpm1r>=zstra) THEN
76			nstra = jk
77			END IF
78			110 CONTINUE
79
80			! inversion car dans orodrag on compte les niveaux a l'envers
81			nktopg = nlev - nktopg + 1
82			nstra = nlev - nstra
83			PRINT *, ' DANS SUGWD nktopg=', nktopg
84			PRINT *, ' DANS SUGWD nstra=', nstra
85
86			gsigcr = 0.80
87
88			gkdrag = 0.2
89			grahilo = 1.
90			grcrit = 0.01
91			gfrcrit = 1.0
92			gkwake = 0.50
93
94			gklift = 0.50
95			gvcrit = 0.0
96
97
98			! ----------------------------------------------------------------
99
100			!* 2. SET VALUES OF SECURITY PARAMETERS
101			! ---------------------------------
102
103			200 CONTINUE
104
105			gvsec = 0.10
106			gssec = 1.E-12
107
108			gtsec = 1.E-07
109
110			! ----------------------------------------------------------------
111
112			RETURN
113			END