55 |
|
|
56 |
! Check decreasing order for "rlatimod": |
! Check decreasing order for "rlatimod": |
57 |
DO i = 2, jjm |
DO i = 2, jjm |
58 |
IF (rlatimod(i) >= rlatimod(i-1)) stop & |
IF (rlatimod(i) >= rlatimod(i-1)) then |
59 |
'"inter_barxy": "rlatimod" should be strictly decreasing' |
print *, '"inter_barxy": "rlatimod" should be strictly decreasing' |
60 |
|
stop 1 |
61 |
|
end IF |
62 |
ENDDO |
ENDDO |
63 |
|
|
64 |
yjmod(:jjm) = ord_coordm(rlatimod) |
yjmod(:jjm) = ord_coordm(rlatimod) |
65 |
IF (jmods == jjm + 1) THEN |
IF (jmods == jjm + 1) THEN |
66 |
IF (90. - yjmod(jjm) < 0.01) stop & |
IF (90. - yjmod(jjm) < 0.01) then |
67 |
'"inter_barxy": with jmods = jjm + 1, yjmod(jjm) should be < 90.' |
print *, '"inter_barxy": with jmods = jjm + 1, ' & |
68 |
|
// 'yjmod(jjm) should be < 90.' |
69 |
|
stop 1 |
70 |
|
end IF |
71 |
ELSE |
ELSE |
72 |
! jmods = jjm |
! jmods = jjm |
73 |
IF (ABS(yjmod(jjm) - 90.) > 0.01) stop & |
IF (ABS(yjmod(jjm) - 90.) > 0.01) then |
74 |
'"inter_barxy": with jmods = jjm, yjmod(jjm) should be 90.' |
print *, '"inter_barxy": with jmods = jjm, yjmod(jjm) should be 90.' |
75 |
|
stop 1 |
76 |
|
end IF |
77 |
ENDIF |
ENDIF |
78 |
|
|
79 |
if (jmods == jjm + 1) yjmod(jjm + 1) = 90. |
if (jmods == jjm + 1) yjmod(jjm + 1) = 90. |
404 |
! Check monotonicity: |
! Check monotonicity: |
405 |
decrois = xi(2) < xi(1) |
decrois = xi(2) < xi(1) |
406 |
DO i = 3, nmax |
DO i = 3, nmax |
407 |
IF (decrois .neqv. xi(i) < xi(i-1)) stop & |
IF (decrois .neqv. xi(i) < xi(i-1)) then |
408 |
'"ord_coord": latitudes are not monotonic' |
print *, '"ord_coord": latitudes are not monotonic' |
409 |
|
stop 1 |
410 |
|
end IF |
411 |
ENDDO |
ENDDO |
412 |
|
|
413 |
IF (abs(xi(1)) < pi) then |
IF (abs(xi(1)) < pi) then |
422 |
print *, "ord_coord" |
print *, "ord_coord" |
423 |
PRINT *, '"xi" should contain the latitudes of the boundaries of ' & |
PRINT *, '"xi" should contain the latitudes of the boundaries of ' & |
424 |
// 'grid cells, not the centers of grid cells.' |
// 'grid cells, not the centers of grid cells.' |
425 |
STOP |
STOP 1 |
426 |
ENDIF |
ENDIF |
427 |
|
|
428 |
IF (decrois) xo(:nmax) = xo(nmax:1:- 1) |
IF (decrois) xo(:nmax) = xo(nmax:1:- 1) |