4 |
|
|
5 |
contains |
contains |
6 |
|
|
7 |
SUBROUTINE gwprofil(nlon, nlev, kgwd, kdx, ktest, kkcrith, kcrit, paphm1, & |
SUBROUTINE gwprofil(nlon, nlev, ktest, kkcrith, kcrit, paphm1, & |
8 |
prho, pstab, pvph, pri, ptau, pdmod, psig, pvar) |
prho, pstab, pvph, pri, ptau, pdmod, psig, pvar) |
9 |
|
|
10 |
! Method. The stress profile for gravity waves is computed as |
! Method. The stress profile for gravity waves is computed as |
27 |
! 0.1 ARGUMENTS |
! 0.1 ARGUMENTS |
28 |
|
|
29 |
INTEGER nlon, nlev |
INTEGER nlon, nlev |
30 |
INTEGER kkcrith(nlon), kcrit(nlon), kdx(nlon), ktest(nlon) |
INTEGER kkcrith(nlon), kcrit(nlon), ktest(nlon) |
31 |
|
|
32 |
REAL paphm1(nlon, nlev+1), pstab(nlon, nlev+1), prho(nlon, nlev+1), & |
REAL paphm1(nlon, nlev+1), pstab(nlon, nlev+1), prho(nlon, nlev+1), & |
33 |
pvph(nlon, nlev+1), pri(nlon, nlev+1), ptau(nlon, nlev+1) |
pvph(nlon, nlev+1), pri(nlon, nlev+1), ptau(nlon, nlev+1) |
38 |
|
|
39 |
! 0.2 LOCAL ARRAYS |
! 0.2 LOCAL ARRAYS |
40 |
|
|
41 |
INTEGER ilevh, kgwd, jl, jk |
INTEGER ilevh, jl, jk |
42 |
REAL zsqr, zalfa, zriw, zdel, zb, zalpha, zdz2n |
REAL zsqr, zalfa, zriw, zdel, zb, zalpha, zdz2n |
43 |
REAL zdelp, zdelpt |
REAL zdelp, zdelpt |
44 |
REAL zdz2(klon, klev), znorm(klon), zoro(klon) |
REAL zdz2(klon, klev), znorm(klon), zoro(klon) |