--- trunk/libf/dyn3d/geopot.f90 2008/02/27 13:16:39 3 +++ trunk/libf/dyn3d/geopot.f90 2010/04/06 17:52:58 32 @@ -1,18 +1,15 @@ SUBROUTINE geopot(ngrid,teta,pk,pks,phis,phi) - ! From libf/dyn3d/geopot.F,v 1.1.1.1 2004/05/19 + ! From libf/dyn3d/geopot.F, version 1.1.1.1 2004/05/19 + ! Author: P. Le Van + ! Objet : calcul du géopotentiel aux milieux des couches + ! L'integration se fait de bas en haut. USE dimens_m USE paramet_m IMPLICIT NONE - ! Auteur: P. Le Van - - ! Objet: - ! calcul du geopotentiel aux milieux des couches - ! l'integration se fait de bas en haut - ! Arguments: INTEGER, INTENT (IN):: ngrid REAL, INTENT (IN):: teta(ngrid,llm), pks(ngrid) @@ -25,12 +22,12 @@ ! ----------------------------------------------------------------------- - ! calcul de phi au niveau 1 pres du sol ..... + ! calcul de phi au niveau 1 pres du sol DO ij = 1, ngrid phi(ij,1) = phis(ij) + teta(ij,1)*(pks(ij)-pk(ij,1)) end DO - ! calcul de phi aux niveaux superieurs ....... + ! calcul de phi aux niveaux superieurs DO l = 2, llm DO ij = 1, ngrid phi(ij,l) = phi(ij,l-1) + 0.5 * (teta(ij,l) + teta(ij,l-1)) &