32 |
! Local |
! Local |
33 |
! --------- |
! --------- |
34 |
|
|
35 |
INTEGER ij, l, j, i, iju, ijq, indu(ip1jmp1), niju |
INTEGER ij, l |
|
INTEGER n0, iadvplus(ip1jmp1, llm), nl(llm) |
|
36 |
|
|
37 |
REAL q(ip1jmp1, llm) |
REAL q(ip1jmp1, llm) |
38 |
REAL dxq(ip1jmp1, llm) |
REAL dxq(ip1jmp1, llm) |
39 |
|
|
|
|
|
|
REAL new_m, zm |
|
40 |
REAL dxqu(ip1jmp1) |
REAL dxqu(ip1jmp1) |
41 |
REAL adxqu(ip1jmp1), dxqmax(ip1jmp1) |
REAL adxqu(ip1jmp1), dxqmax(ip1jmp1) |
42 |
|
|
43 |
LOGICAL extremum, first |
REAL ssum |
|
SAVE first |
|
|
|
|
|
REAL ssum, cvmgp, cvmgt |
|
44 |
INTEGER ismax, ismin |
INTEGER ismax, ismin |
45 |
EXTERNAL ssum, convflu, ismin, ismax |
EXTERNAL ssum, convflu, ismin, ismax |
46 |
|
|
|
DATA first/.TRUE./ |
|
|
|
|
|
|
|
47 |
DO l = 1, llm |
DO l = 1, llm |
48 |
DO ij = 1, ip1jmp1 |
DO ij = 1, ip1jmp1 |
49 |
q(ij, l) = s0(ij, l)/sm(ij, l) |
q(ij, l) = s0(ij, l)/sm(ij, l) |