pkg/dic/tracer_meanarea.F

#include "DIC_OPTIONS.h"
#include "GCHEM_OPTIONS.h"

C !INTERFACE: ==========================================================
      SUBROUTINE TRACER_MEANAREA(myThid,tracer, kLev,
     &                    gsm_tracer)

C !DESCRIPTION:
C  Calculate surface means of tracer

C !USES: ===============================================================
      IMPLICIT NONE
#include "SIZE.h"
#include "GRID.h"
#include "DYNVARS.h"
#include "EEPARAMS.h"
#include "PARAMS.h"

C !INPUT PARAMETERS: ===================================================
C  myThid               :: thread number
      INTEGER myThid
      INTEGER kLev
      _RL tracer (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy)

C !OUTPUT PARAMETERS: ===================================================
C  gcm_tracer               :: surface mean of tracer
      _RL gsm_tracer

C !LOCAL VARIABLES: ====================================================
      INTEGER i,j,bi,bj
      _RL sumarea
CEOP

       sumarea= 0. _d 0
       gsm_tracer= 0. _d 0

       DO bj=myByLo(myThid),myByHi(myThid)
        DO bi=myBxLo(myThid),myBxHi(myThid)
         DO i=1,sNx
          DO j=1,sNy
            sumarea=sumarea+
     &            rA(i,j,bi,bj)*maskC(i,j,kLev,bi,bj)
            gsm_tracer=gsm_tracer+
     &          tracer(i,j,kLev,bi,bj)*rA(i,j,bi,bj)*
     &                        maskC(i,j,kLev,bi,bj)
          ENDDO
          ENDDO
         ENDDO
        ENDDO
        _GLOBAL_SUM_R8(gsm_tracer,myThid)
        _GLOBAL_SUM_R8(sumarea,myThid)
        gsm_tracer=gsm_tracer/sumarea

       RETURN
       END

1	#include "DIC_OPTIONS.h"
2	#include "GCHEM_OPTIONS.h"
3
4	C !INTERFACE: ==========================================================
5	SUBROUTINE TRACER_MEANAREA(myThid,tracer, kLev,
6	& gsm_tracer)
7
8	C !DESCRIPTION:
9	C Calculate surface means of tracer
10
11	C !USES: ===============================================================
12	IMPLICIT NONE
13	#include "SIZE.h"
14	#include "GRID.h"
15	#include "DYNVARS.h"
16	#include "EEPARAMS.h"
17	#include "PARAMS.h"
18
19	C !INPUT PARAMETERS: ===================================================
20	C myThid :: thread number
21	INTEGER myThid
22	INTEGER kLev
23	_RL tracer (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy)
24
25	C !OUTPUT PARAMETERS: ===================================================
26	C gcm_tracer :: surface mean of tracer
27	_RL gsm_tracer
28
29	C !LOCAL VARIABLES: ====================================================
30	INTEGER i,j,bi,bj
31	_RL sumarea
32	CEOP
33
34	sumarea= 0. _d 0
35	gsm_tracer= 0. _d 0
36
37	DO bj=myByLo(myThid),myByHi(myThid)
38	DO bi=myBxLo(myThid),myBxHi(myThid)
39	DO i=1,sNx
40	DO j=1,sNy
41	sumarea=sumarea+
42	& rA(i,j,bi,bj)*maskC(i,j,kLev,bi,bj)
43	gsm_tracer=gsm_tracer+
44	& tracer(i,j,kLev,bi,bj)rA(i,j,bi,bj)
45	& maskC(i,j,kLev,bi,bj)
46	ENDDO
47	ENDDO
48	ENDDO
49	ENDDO
50	_GLOBAL_SUM_R8(gsm_tracer,myThid)
51	_GLOBAL_SUM_R8(sumarea,myThid)
52	gsm_tracer=gsm_tracer/sumarea
53
54	RETURN
55	END
56