source: LMDZ.3.3/trunk/libf/dyn3d/divgrad.F @ 2

Last change on this file since 2 was 2, checked in by lmdz, 25 years ago

Initial revision

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 1.3 KB
RevLine 
[2]1      SUBROUTINE divgrad (klevel,h, lh, divgra )
2      IMPLICIT NONE
3c
4c=======================================================================
5c
6c  Auteur :   P. Le Van
7c  ----------
8c
9c                              lh
10c      calcul de  (div( grad ))   de h  .....
11c      h  et lh  sont des arguments  d'entree pour le s-prog
12c      divgra     est  un argument  de sortie pour le s-prog
13c
14c=======================================================================
15c
16c   declarations:
17c   -------------
18c
19#include "dimensions.h"
20#include "paramet.h"
21#include "comgeom.h"
22#include "comdissipn.h"
23#include "logic.h"
24c
25      INTEGER klevel
26      REAL h( ip1jmp1,klevel ), divgra( ip1jmp1,klevel )
27c
28      REAL ghy(ip1jm,llm), ghx(ip1jmp1,llm)
29
30      INTEGER  l,ij,iter,lh
31c
32      EXTERNAL  filtreg
33      EXTERNAL  SCOPY, grad, covcont, diverg
34c
35c
36      CALL SCOPY ( ip1jmp1*klevel,h,1,divgra,1 )
37c
38      DO 10 iter = 1,lh
39
40      CALL filtreg ( divgra,jjp1,klevel,2,1,.true.,1  )
41
42      CALL    grad (klevel,divgra, ghx  , ghy          )
43      CALL  diverg (klevel,  ghx , ghy  , divgra       )
44
45      CALL filtreg ( divgra,jjp1,klevel,2,1,.true.,1)
46
47      DO 5 l = 1,klevel
48      DO 4  ij = 1, ip1jmp1
49      divgra( ij,l ) = - cdivh * divgra( ij,l )
50   4  CONTINUE
51   5  CONTINUE
52c
53  10  CONTINUE
54      RETURN
55      END
Note: See TracBrowser for help on using the repository browser.