1 |
|
|
! |
2 |
|
|
! $Id: minmaxqfi.F90 2311 2015-06-25 07:45:24Z emillour $ |
3 |
|
|
! |
4 |
|
1152 |
SUBROUTINE minmaxqfi(zq,qmin,qmax,comment) |
5 |
|
|
USE dimphy |
6 |
|
|
USE print_control_mod, ONLY: prt_level |
7 |
|
|
IMPLICIT NONE |
8 |
|
|
|
9 |
|
|
! Entrees |
10 |
|
|
REAL,DIMENSION(klon,klev), INTENT(IN) :: zq |
11 |
|
|
REAL,INTENT(IN) :: qmin,qmax |
12 |
|
|
CHARACTER(LEN=*),INTENT(IN) :: comment |
13 |
|
|
|
14 |
|
|
! Local |
15 |
|
2304 |
INTEGER,DIMENSION(klon) :: jadrs |
16 |
|
|
INTEGER :: i, jbad, k |
17 |
|
|
|
18 |
✓✓ |
46080 |
DO k = 1, klev |
19 |
|
|
jbad = 0 |
20 |
✓✓ |
44703360 |
DO i = 1, klon |
21 |
✓✗✗✓
|
44703360 |
IF (zq(i,k).GT.qmax .OR. zq(i,k).LT.qmin) THEN |
22 |
|
|
jbad = jbad + 1 |
23 |
|
|
jadrs(jbad) = i |
24 |
|
|
ENDIF |
25 |
|
|
ENDDO |
26 |
✗✓ |
46080 |
IF (jbad.GT.0) THEN |
27 |
|
|
WRITE(*,*)comment, "k = ", k, "jbad > 0" |
28 |
|
|
if (prt_level >= 1) then |
29 |
|
|
DO i = 1, jbad |
30 |
|
|
WRITE(*,*) "i,q=", jadrs(i), zq(jadrs(i),k) |
31 |
|
|
ENDDO |
32 |
|
|
end if |
33 |
|
|
ENDIF |
34 |
|
|
ENDDO |
35 |
|
|
|
36 |
|
1152 |
END SUBROUTINE minmaxqfi |