! ! $Header$ ! SUBROUTINE rotat (klevel, x, y, rot ) ! ! Auteur : P.Le Van !************************************************************** !. calcule le rotationnel ! a tous les niveaux d'1 vecteur de comp. x et y .. ! x et y etant des composantes covariantes ... !******************************************************************** ! klevel, x et y sont des arguments d'entree pour le s-prog ! rot est un argument de sortie pour le s-prog ! USE comgeom_mod_h USE dimensions_mod, ONLY: iim, jjm, llm, ndm USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, & ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm IMPLICIT NONE ! ! ! ..... variables en arguments ...... ! INTEGER :: klevel REAL :: rot( ip1jm,klevel ) REAL :: x( ip1jmp1,klevel ), y( ip1jm,klevel ) ! ! ... variables locales ... ! INTEGER :: l, ij ! ! DO l = 1,klevel ! DO ij = 1, ip1jm - 1 rot( ij,l ) = y( ij+1 , l ) - y( ij,l ) + & x(ij +iip1, l ) - x( ij,l ) ENDDO ! ! .... correction pour rot( iip1,j,l) .... ! .... rot(iip1,j,l)= rot(1,j,l) ... !DIR$ IVDEP DO ij = iip1, ip1jm, iip1 rot( ij,l ) = rot( ij -iim,l ) ENDDO ! END DO !cc CALL filtreg( rot, jjm, klevel, 2, 2, .FALSE., 1 ) DO l = 1, klevel DO ij = 1, ip1jm rot(ij,l) = rot(ij,l) * unsairez(ij) ENDDO ENDDO ! ! RETURN END SUBROUTINE rotat