1 |
guez |
3 |
|
2 |
guez |
81 |
! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/dump2d.F,v 1.1.1.1 2004/05/19 |
3 |
|
|
! 12:53:05 lmdzadmin Exp $ |
4 |
guez |
3 |
|
5 |
guez |
81 |
SUBROUTINE dump2d(im, jm, z, nom_z) |
6 |
|
|
IMPLICIT NONE |
7 |
|
|
INTEGER im, jm |
8 |
|
|
REAL z(im, jm) |
9 |
|
|
CHARACTER *80 nom_z |
10 |
guez |
3 |
|
11 |
guez |
81 |
INTEGER i, j, imin, illm, jmin, jllm |
12 |
|
|
REAL zmin, zllm |
13 |
guez |
3 |
|
14 |
guez |
81 |
PRINT *, nom_z |
15 |
guez |
3 |
|
16 |
guez |
81 |
zmin = z(1, 1) |
17 |
|
|
zllm = z(1, 1) |
18 |
|
|
imin = 1 |
19 |
|
|
illm = 1 |
20 |
|
|
jmin = 1 |
21 |
|
|
jllm = 1 |
22 |
guez |
3 |
|
23 |
guez |
81 |
DO j = 1, jm |
24 |
|
|
DO i = 1, im |
25 |
|
|
IF (z(i,j)>zllm) THEN |
26 |
|
|
illm = i |
27 |
|
|
jllm = j |
28 |
|
|
zllm = z(i, j) |
29 |
|
|
END IF |
30 |
|
|
IF (z(i,j)<zmin) THEN |
31 |
|
|
imin = i |
32 |
|
|
jmin = j |
33 |
|
|
zmin = z(i, j) |
34 |
|
|
END IF |
35 |
|
|
END DO |
36 |
|
|
END DO |
37 |
|
|
|
38 |
|
|
PRINT *, 'MIN: ', zmin |
39 |
|
|
PRINT *, 'MAX: ', zllm |
40 |
|
|
|
41 |
|
|
IF (zllm>zmin) THEN |
42 |
|
|
DO j = 1, jm |
43 |
|
|
WRITE (*, '(72i1)')(nint(10.*(z(i,j)-zmin)/(zllm-zmin)), i=1, im) |
44 |
|
|
END DO |
45 |
|
|
END IF |
46 |
|
|
RETURN |
47 |
|
|
END SUBROUTINE dump2d |