| 
                Last change
                  on this file since 1962 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:
            1.1 KB
           | 
      
      
        
  | Line |   | 
|---|
| 1 | ! | 
|---|
| 2 | ! $Header$ | 
|---|
| 3 | ! | 
|---|
| 4 |       SUBROUTINE laplacien_rotgam ( klevel, rotin, rotout ) | 
|---|
| 5 | c | 
|---|
| 6 | c     P. Le Van | 
|---|
| 7 | c | 
|---|
| 8 | c   ************************************************************ | 
|---|
| 9 | c   ... calcul de  (rotat x nxgrad)_gam  du rotationnel rotin .. | 
|---|
| 10 | c   ************************************************************ | 
|---|
| 11 | c     klevel et teta  sont des arguments  d'entree pour le s-prog | 
|---|
| 12 | c      divgra     est  un argument  de sortie pour le s-prog | 
|---|
| 13 | c | 
|---|
| 14 |       IMPLICIT NONE | 
|---|
| 15 | c | 
|---|
| 16 | #include "dimensions.h" | 
|---|
| 17 | #include "paramet.h" | 
|---|
| 18 | #include "comgeom.h" | 
|---|
| 19 |  | 
|---|
| 20 | c | 
|---|
| 21 | c    .............   variables  en  arguments    ........... | 
|---|
| 22 | c | 
|---|
| 23 |       INTEGER klevel | 
|---|
| 24 |       REAL rotin( ip1jm,klevel ), rotout( ip1jm,klevel ) | 
|---|
| 25 | c | 
|---|
| 26 | c   ............     variables   locales     ............... | 
|---|
| 27 | c | 
|---|
| 28 |       INTEGER l, ij | 
|---|
| 29 |       REAL ghy(ip1jm,llm), ghx(ip1jmp1,llm) | 
|---|
| 30 | c   ........................................................ | 
|---|
| 31 | c | 
|---|
| 32 | c | 
|---|
| 33 |  | 
|---|
| 34 |       CALL   nxgrad_gam ( klevel, rotin,   ghx ,   ghy  ) | 
|---|
| 35 |       CALL   rotat_nfil ( klevel, ghx  ,   ghy , rotout ) | 
|---|
| 36 | c | 
|---|
| 37 |       DO l = 1, klevel | 
|---|
| 38 |         DO ij = 1, ip1jm | 
|---|
| 39 |          rotout(ij,l) = rotout(ij,l) * unsairz_gam(ij) | 
|---|
| 40 |         ENDDO | 
|---|
| 41 |       ENDDO | 
|---|
| 42 |  | 
|---|
| 43 |       RETURN | 
|---|
| 44 |       END | 
|---|
       
      
      Note: See 
TracBrowser
        for help on using the repository browser.