pkg/gchem/gchem_surfmean.F

7f498da5b5 Oliv*0001 #include "GCHEM_OPTIONS.h"
daab022f42 Step*0002
                0003 C !INTERFACE: ==========================================================
7f498da5b5 Oliv*0004       SUBROUTINE GCHEM_SURFMEAN(
35050126c6 Oliv*0005      I                    tracer,
2471b84c83 Jean*0006      O                    gsm_tracer,
                0007      I                    myThid )
08536d17ba Step*0008
                0009 C !DESCRIPTION:
7f498da5b5 Oliv*0010 C  Calculate surface mean of tracer
08536d17ba Step*0011
                0012 C !USES: ===============================================================
daab022f42 Step*0013       IMPLICIT NONE
                0014 #include "SIZE.h"
                0015 #include "EEPARAMS.h"
                0016 #include "PARAMS.h"
f05c6d30c5 Jean*0017 #include "GRID.h"
08536d17ba Step*0018
daab022f42 Step*0019 C !INPUT PARAMETERS: ===================================================
                0020 C  myThid               :: thread number
29ad036528 Step*0021       INTEGER myThid
e32a77df02 Jean*0022       _RL tracer (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
08536d17ba Step*0023
                0024 C !OUTPUT PARAMETERS: ===================================================
f05c6d30c5 Jean*0025 C  gcm_tracer           :: surface mean of tracer
daab022f42 Step*0026       _RL gsm_tracer
                0027
08536d17ba Step*0028 C !LOCAL VARIABLES: ====================================================
daab022f42 Step*0029       INTEGER i,j,bi,bj
f05c6d30c5 Jean*0030       INTEGER ks
35050126c6 Oliv*0031       _RL sumTracerTile(nSx,nSy)
7ccec407f7 Davi*0032       _RL sumTracer
08536d17ba Step*0033 CEOP
daab022f42 Step*0034
35050126c6 Oliv*0035       IF ( usingPCoords ) THEN
f05c6d30c5 Jean*0036        ks = Nr
35050126c6 Oliv*0037       ELSE
f05c6d30c5 Jean*0038        ks = 1
35050126c6 Oliv*0039       ENDIF
                0040
daab022f42 Step*0041        DO bj=myByLo(myThid),myByHi(myThid)
                0042         DO bi=myBxLo(myThid),myBxHi(myThid)
f05c6d30c5 Jean*0043          sumTracerTile(bi,bj) = 0. _d 0
e32a77df02 Jean*0044          DO j=1,sNy
                0045           DO i=1,sNx
f05c6d30c5 Jean*0046             sumTracerTile(bi,bj) = sumTracerTile(bi,bj)
                0047      &        + tracer(i,j,ks,bi,bj)*rA(i,j,bi,bj)*maskInC(i,j,bi,bj)
daab022f42 Step*0048           ENDDO
                0049          ENDDO
                0050         ENDDO
e32a77df02 Jean*0051        ENDDO
35050126c6 Oliv*0052        CALL GLOBAL_SUM_TILE_RL(sumTracerTile,sumTracer,myThid)
f05c6d30c5 Jean*0053        gsm_tracer = sumTracer/globalArea
daab022f42 Step*0054
f05c6d30c5 Jean*0055       RETURN
                0056       END