1 |
|
|
! |
2 |
|
|
! $Id: dump2d.F 1952 2014-01-28 13:05:47Z lguez $ |
3 |
|
|
! |
4 |
|
|
SUBROUTINE dump2d(im,jm,z,nom_z) |
5 |
|
|
IMPLICIT NONE |
6 |
|
|
INTEGER im,jm |
7 |
|
|
REAL z(im,jm) |
8 |
|
|
CHARACTER (len=*) :: nom_z |
9 |
|
|
|
10 |
|
|
INTEGER i,j,imin,illm,jmin,jllm |
11 |
|
|
REAL zmin,zllm |
12 |
|
|
|
13 |
|
|
WRITE(*,*) "dump2d: ",trim(nom_z) |
14 |
|
|
|
15 |
|
|
zmin=z(1,1) |
16 |
|
|
zllm=z(1,1) |
17 |
|
|
imin=1 |
18 |
|
|
illm=1 |
19 |
|
|
jmin=1 |
20 |
|
|
jllm=1 |
21 |
|
|
|
22 |
|
|
DO j=1,jm |
23 |
|
|
DO i=1,im |
24 |
|
|
IF(z(i,j).GT.zllm) THEN |
25 |
|
|
illm=i |
26 |
|
|
jllm=j |
27 |
|
|
zllm=z(i,j) |
28 |
|
|
ENDIF |
29 |
|
|
IF(z(i,j).LT.zmin) THEN |
30 |
|
|
imin=i |
31 |
|
|
jmin=j |
32 |
|
|
zmin=z(i,j) |
33 |
|
|
ENDIF |
34 |
|
|
ENDDO |
35 |
|
|
ENDDO |
36 |
|
|
|
37 |
|
|
PRINT*,'MIN: ',zmin |
38 |
|
|
PRINT*,'MAX: ',zllm |
39 |
|
|
|
40 |
|
|
IF(zllm.GT.zmin) THEN |
41 |
|
|
DO j=1,jm |
42 |
|
|
WRITE(*,'(600i1)') (NINT(10.*(z(i,j)-zmin)/(zllm-zmin)),i=1,im) |
43 |
|
|
ENDDO |
44 |
|
|
ENDIF |
45 |
|
|
RETURN |
46 |
|
|
END |