3 |
! 12:53:07 lmdzadmin Exp $ |
! 12:53:07 lmdzadmin Exp $ |
4 |
|
|
5 |
SUBROUTINE prather(q, w, masse, pbaru, pbarv, nt, dt) |
SUBROUTINE prather(q, w, masse, pbaru, pbarv, nt, dt) |
6 |
USE dimens_m |
|
|
USE paramet_m |
|
7 |
USE comconst |
USE comconst |
8 |
|
USE dimens_m |
9 |
USE disvert_m |
USE disvert_m |
10 |
USE comgeom |
USE dynetat0_m, only: rlonv, rlonu |
11 |
USE nr_util, ONLY: pi |
USE nr_util, ONLY: pi |
12 |
|
USE paramet_m |
13 |
|
|
14 |
IMPLICIT NONE |
IMPLICIT NONE |
15 |
|
|
16 |
! ======================================================================= |
! ======================================================================= |
30 |
|
|
31 |
! Arguments: |
! Arguments: |
32 |
! ---------- |
! ---------- |
33 |
INTEGER iq, nt |
INTEGER nt |
34 |
REAL, INTENT (IN) :: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm) |
REAL, INTENT (IN) :: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm) |
35 |
REAL masse(iip1, jjp1, llm) |
REAL masse(iip1, jjp1, llm) |
36 |
REAL q(iip1, jjp1, llm, 0:9) |
REAL q(iip1, jjp1, llm, 0:9) |
37 |
REAL w(ip1jmp1, llm) |
REAL w(ip1jmp1, llm) |
|
INTEGER ordre, ilim |
|
38 |
|
|
39 |
! Local: |
! Local: |
40 |
! ------ |
! ------ |
41 |
LOGICAL limit |
LOGICAL limit |
|
REAL zq(iip1, jjp1, llm) |
|
42 |
REAL sm(iip1, jjp1, llm) |
REAL sm(iip1, jjp1, llm) |
43 |
REAL s0(iip1, jjp1, llm), sx(iip1, jjp1, llm) |
REAL s0(iip1, jjp1, llm), sx(iip1, jjp1, llm) |
44 |
REAL sy(iip1, jjp1, llm), sz(iip1, jjp1, llm) |
REAL sy(iip1, jjp1, llm), sz(iip1, jjp1, llm) |
53 |
|
|
54 |
REAL sinlon(iip1), sinlondlon(iip1) |
REAL sinlon(iip1), sinlondlon(iip1) |
55 |
REAL coslon(iip1), coslondlon(iip1) |
REAL coslon(iip1), coslondlon(iip1) |
|
REAL qmin, qmax |
|
|
SAVE qmin, qmax |
|
56 |
SAVE sinlon, coslon, sinlondlon, coslondlon |
SAVE sinlon, coslon, sinlondlon, coslondlon |
57 |
REAL dyn1, dyn2, dys1, dys2, qpn, qps, dqzpn, dqzps |
REAL dyn1, dyn2, dys1, dys2, qpn, qps, dqzpn, dqzps |
58 |
REAL masn, mass |
REAL masn, mass |
66 |
|
|
67 |
|
|
68 |
DATA first/.TRUE./ |
DATA first/.TRUE./ |
|
DATA qmin, qmax/ -1.E33, 1.E33/ |
|
|
|
|
69 |
|
|
70 |
! ========================================================================== |
! ========================================================================== |
71 |
! ========================================================================== |
! ========================================================================== |