--- trunk/dyn3d/laplacien_rotgam.f 2013/11/15 18:45:49 76 +++ trunk/dyn3d/laplacien_rotgam.f90 2014/03/05 14:38:41 81 @@ -1,44 +1,45 @@ -! -! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/laplacien_rotgam.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $ -! - SUBROUTINE laplacien_rotgam ( klevel, rotin, rotout ) -c -c P. Le Van -c -c ************************************************************ -c ... calcul de (rotat x nxgrad)_gam du rotationnel rotin .. -c ************************************************************ -c klevel et teta sont des arguments d'entree pour le s-prog -c divgra est un argument de sortie pour le s-prog -c - use dimens_m - use paramet_m - use comgeom - IMPLICIT NONE -c - -c -c ............. variables en arguments ........... -c - INTEGER, intent(in):: klevel - REAL rotin( ip1jm,klevel ), rotout( ip1jm,klevel ) -c -c ............ variables locales ............... -c - INTEGER l, ij - REAL ghy(ip1jm,llm), ghx(ip1jmp1,llm) -c ........................................................ -c -c - - CALL nxgrad_gam ( klevel, rotin, ghx , ghy ) - CALL rotat_nfil ( klevel, ghx , ghy , rotout ) -c - DO l = 1, klevel - DO ij = 1, ip1jm - rotout(ij,l) = rotout(ij,l) * unsairz_gam(ij) - ENDDO - ENDDO - RETURN - END +! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/laplacien_rotgam.F,v 1.1.1.1 +! 2004/05/19 12:53:05 lmdzadmin Exp $ + +SUBROUTINE laplacien_rotgam(klevel, rotin, rotout) + + ! P. Le Van + + ! ************************************************************ + ! ... calcul de (rotat x nxgrad)_gam du rotationnel rotin .. + ! ************************************************************ + ! klevel et teta sont des arguments d'entree pour le s-prog + ! divgra est un argument de sortie pour le s-prog + + USE dimens_m + USE paramet_m + USE comgeom + IMPLICIT NONE + + + + ! ............. variables en arguments ........... + + INTEGER, INTENT (IN) :: klevel + REAL rotin(ip1jm, klevel), rotout(ip1jm, klevel) + + ! ............ variables locales ............... + + INTEGER l, ij + REAL ghy(ip1jm, llm), ghx(ip1jmp1, llm) + ! ........................................................ + + + + CALL nxgrad_gam(klevel, rotin, ghx, ghy) + CALL rotat_nfil(klevel, ghx, ghy, rotout) + + DO l = 1, klevel + DO ij = 1, ip1jm + rotout(ij, l) = rotout(ij, l)*unsairz_gam(ij) + END DO + END DO + + RETURN +END SUBROUTINE laplacien_rotgam