! ! $Header$ ! SUBROUTINE nxgraro2 (klevel,xcov, ycov, lr, grx, gry ) ! ! P.Le Van . ! *********************************************************** ! lr ! calcul de ( nxgrad (rot) ) du vect. v .... ! ! xcov et ycov etant les compos. covariantes de v ! *********************************************************** ! xcov , ycov et lr sont des arguments d'entree pour le s-prog ! grx et gry sont des arguments de sortie pour le s-prog ! ! USE comdissipn_mod_h USE dimensions_mod, ONLY: iim, jjm, llm, ndm USE paramet_mod_h IMPLICIT NONE ! ! ! ...... variables en arguments ....... ! INTEGER :: klevel REAL :: xcov( ip1jmp1,klevel ), ycov( ip1jm,klevel ) REAL :: grx( ip1jmp1,klevel ), gry( ip1jm,klevel ) ! ! ...... variables locales ........ ! REAL :: rot(ip1jm,llm) , signe, nugradrs INTEGER :: l,ij,iter,lr ! ........................................................ ! ! ! signe = (-1.)**lr nugradrs = signe * crot ! CALL SCOPY ( ip1jmp1* klevel, xcov, 1, grx, 1 ) CALL SCOPY ( ip1jm * klevel, ycov, 1, gry, 1 ) ! CALL rotatf ( klevel, grx, gry, rot ) ! CALL laplacien_rot ( klevel, rot, rot,grx,gry ) ! ! ..... Iteration de l'operateur laplacien_rotgam ..... ! DO iter = 1, lr -2 CALL laplacien_rotgam ( klevel, rot, rot ) ENDDO ! ! CALL filtreg( rot, jjm, klevel, 2,1, .FALSE.,1) CALL nxgrad ( klevel, rot, grx, gry ) ! DO l = 1, klevel DO ij = 1, ip1jm gry( ij,l ) = gry( ij,l ) * nugradrs ENDDO DO ij = 1, ip1jmp1 grx( ij,l ) = grx( ij,l ) * nugradrs ENDDO ENDDO ! RETURN END SUBROUTINE nxgraro2