--- trunk/dyn3d/covcont.f 2013/11/15 18:45:49 76 +++ trunk/dyn3d/covcont.f90 2014/03/05 14:38:41 81 @@ -1,43 +1,44 @@ -! -! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/covcont.F,v 1.1.1.1 2004/05/19 12:53:07 lmdzadmin Exp $ -! - SUBROUTINE covcont (klevel,ucov, vcov, ucont, vcont ) - use dimens_m - use paramet_m - use comgeom - IMPLICIT NONE - -c======================================================================= -c -c Auteur: P. Le Van -c ------- -c -c Objet: -c ------ -c -c ********************************************************************* -c calcul des compos. contravariantes a partir des comp.covariantes -c ******************************************************************** -c -c======================================================================= - - - INTEGER klevel - REAL, intent(in):: ucov( ip1jmp1,klevel ), vcov( ip1jm,klevel ) - REAL ucont( ip1jmp1,klevel ), vcont( ip1jm,klevel ) - INTEGER l,ij - - - DO 10 l = 1,klevel - - DO 2 ij = iip2, ip1jm - ucont( ij,l ) = ucov( ij,l ) * unscu2( ij ) - 2 CONTINUE - - DO 4 ij = 1,ip1jm - vcont( ij,l ) = vcov( ij,l ) * unscv2( ij ) - 4 CONTINUE - - 10 CONTINUE - RETURN - END + +! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/covcont.F,v 1.1.1.1 2004/05/19 +! 12:53:07 lmdzadmin Exp $ + +SUBROUTINE covcont(klevel, ucov, vcov, ucont, vcont) + USE dimens_m + USE paramet_m + USE comgeom + IMPLICIT NONE + + ! ======================================================================= + + ! Auteur: P. Le Van + ! ------- + + ! Objet: + ! ------ + + ! ********************************************************************* + ! calcul des compos. contravariantes a partir des comp.covariantes + ! ******************************************************************** + + ! ======================================================================= + + + INTEGER klevel + REAL, INTENT (IN) :: ucov(ip1jmp1, klevel), vcov(ip1jm, klevel) + REAL ucont(ip1jmp1, klevel), vcont(ip1jm, klevel) + INTEGER l, ij + + + DO l = 1, klevel + + DO ij = iip2, ip1jm + ucont(ij, l) = ucov(ij, l)*unscu2(ij) + END DO + + DO ij = 1, ip1jm + vcont(ij, l) = vcov(ij, l)*unscv2(ij) + END DO + + END DO + RETURN +END SUBROUTINE covcont