*=*=*=*= gradiv.html =*=*=*=*
SUBROUTINE gradiv(klevel, xcont, ycont, ld, gdx, gdy ) c c Auteur : P. Le Van c c ********************************************************** c c ld c calcul de (grad (div) ) du vect. v .... c c xcont et ycont etant les composant.contravariantes de v c ********************************************************** c xcont , ycont et ld sont des arguments d'entree pour le s-prog c gdx et gdy sont des arguments de sortie pour le s-prog c c IMPLICIT NONE c #include "dimensions.h" #include "paramet.h" #include "comdissipn.h" #include "logic.h" INTEGER klevel c REAL xcont( ip1jmp1,klevel ), ycont( ip1jm,klevel ) REAL gdx( ip1jmp1,klevel ), gdy( ip1jm,klevel ) REAL div(ip1jmp1,llm) INTEGER l,ij,iter,ld c EXTERNAL SCOPY, diverg, covcont, grad EXTERNAL filtreg c c CALL SCOPY(ip1jmp1*klevel,xcont,1,gdx,1 ) CALL SCOPY(ip1jm*klevel,ycont,1,gdy,1 ) c DO 10 iter = 1,ld c CALL diverg(klevel, gdx , gdy, div ) CALL filtreg( div, jjp1, klevel, 2,1, .true.,2) CALL grad( klevel, div, gdx, gdy ) IF ( iter. lt . ld ) THEN CALL covcont (klevel,gdx, gdy, gdx, gdy ) END IF c DO 5 l = 1, klevel DO 3 ij = 1, ip1jmp1 gdx( ij,l ) = - gdx( ij,l ) * cdivu 3 CONTINUE DO 4 ij = 1, ip1jm gdy( ij,l ) = - gdy( ij,l ) * cdivu 4 CONTINUE 5 CONTINUE c 10 CONTINUE RETURN END