| 
                Last change
                  on this file since 3598 was
                  1300,
                  checked in by emillour, 11 years ago
           | 
        
        
          | 
               
Common dynamics: 
Some updates to keep up with LMDZ5 Earth model evolution (up to LMDZ5 rev 1955). 
Main change is the introduction of a "dyn3d_common" directory 
to store files common to dyn3d and dyn3dpar. 
See file "DOC/chantiers/commit_importants.log" for detailed list 
of changes. These changes do not change results on test cases. 
EM 
 
           | 
        
        | 
            File size:
            911 bytes
           | 
      
      
        
  | Line |   | 
|---|
| 1 | ! | 
|---|
| 2 | ! $Id: dump2d.F 1279 2009-12-10 09:02:56Z fairhead $ | 
|---|
| 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 | 
|---|
       
      
      Note: See 
TracBrowser
        for help on using the repository browser.