!
! $Header$
!
      SUBROUTINE moyglo_pondaire(nhori, champ, aire,                                 &
       &            ok_msk, msk, moyglo)
!c
       USE dimphy
       IMPLICIT none
!c
!c ==================================================================
!c I. Musat, 07.2004
!c
!c Calcul moyenne globale ponderee par l'aire totale, avec ou sans masque
!c
!c moyenne = Somme_(champ* aire)/Somme_aire
!c
!c ==================================================================
!c
#include "dimensions.h"
!cym#include "dimphy.h"
       INTEGER i, nhori
       REAL champ(klon), aire(klon), msk(klon)
       LOGICAL ok_msk 
       REAL moyglo
!c
!c var locale
       REAL airetot
!c
!c      PRINT*,'moyglo_pondaire nhori',nhori
!c
       airetot=0.
       moyglo=0.
!c
       IF(ok_msk) THEN
        DO i=1, nhori
!c        IF(msk(i).EQ.1.) THEN 
         IF(msk(i).GT.0.) THEN 
!c
!c aire totale
          airetot=airetot+aire(i)*msk(i)
!c
!c ponderation par la masse
          moyglo=moyglo+champ(i)* aire(i)*msk(i) 
         ENDIF
        ENDDO
!c
       ELSE !ok_msk
        DO i=1, nhori
!c
!c aire totale
          airetot=airetot+aire(i)
!c
!c ponderation par la masse
          moyglo=moyglo+champ(i)* aire(i) 
        ENDDO 
!c
       ENDIF
!c 
!c moyenne ponderee par l'aire
       moyglo=moyglo/airetot
!c
       RETURN 
       END
!c
       SUBROUTINE moyglo_pondaima(nhori, nvert, champ,                               &
       & aire, pbord, moyglo)
       USE dimphy
       IMPLICIT none
!c ==================================================================
!c I. Musat, 07.2004
!c
!c Calcul moyenne globale ponderee par la masse d'air, divisee par l'aire
!c totale avec ou sans masque
!c
!c moyenne = Somme_(champ* masse_dair)/Somme_aire
!c
!c ==================================================================
#include "dimensions.h"
!cym#include "dimphy.h"
#include "YOMCST.h"
       INTEGER i, k, nhori, nvert
       REAL champ(klon,klev), aire(klon)
       REAL pbord(klon,klev+1)
       REAL moyglo
!c
!c var locale
       REAL airetot
!c
!c      PRINT*,'moyglo_pondaima RG, nhori, nvert',RG,nhori,nvert
!c
!c ponderation par la masse
       moyglo=0.
       DO k=1, nvert
       DO i=1, nhori
        moyglo=moyglo+                                                               &
       &  champ(i,k)*(pbord(i,k)-pbord(i,k+1))/RG*aire(i)
       ENDDO
       ENDDO
!c
!c aire totale
       airetot=0.
       DO i=1, nhori
        airetot=airetot+aire(i)
       ENDDO
!c
!c moyenne par mettre carre avec ponderation par la masse
       moyglo=moyglo/airetot
!c
      RETURN
      END
!c
       SUBROUTINE moyglo_pondmass(nhori, nvert, champ,                               &
       & aire, pbord, moyglo)
       USE dimphy
       IMPLICIT none
!c ==================================================================
!c I. Musat, 07.2004
!c
!c Calcul moyenne globale ponderee par la masse d'air, divisee par la
!c masse totale d'air, avec ou sans masque
!c
!c moyenne = Somme_(champ* masse_dair)/Somme_(masse_dair)
!c
!c ==================================================================
#include "dimensions.h"
!cym#include "dimphy.h"
#include "YOMCST.h"
       INTEGER i, k, nhori, nvert
       REAL champ(klon,klev), aire(klon)
       REAL pbord(klon,klev+1)
       REAL moyglo
!c
!c var locale
       REAL massetot
!c
!c      PRINT*,'moyglo_pondmass RG, nhori, nvert',RG,nhori,nvert
!c
!c ponderation par la masse
       moyglo=0.
       DO k=1, nvert
       DO i=1, nhori
        moyglo=moyglo+                                                               &
       &  champ(i,k)*(pbord(i,k)-pbord(i,k+1))/RG*aire(i)
       ENDDO
       ENDDO
!c
!c masse totale
       massetot=0.
       DO k=1, nvert
       DO i=1, nhori
        massetot=massetot+                                                           &
       &  (pbord(i,k)-pbord(i,k+1))/RG*aire(i)
       ENDDO
       ENDDO
!c
!c moyenne par mettre carre avec ponderation par la masse
       moyglo=moyglo/massetot
!c
      RETURN
      END
!c