/[MITgcm]/MITgcm/pkg/dic/tracer_meanarea.F
ViewVC logotype

Diff of /MITgcm/pkg/dic/tracer_meanarea.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.6 by jmc, Wed Aug 9 01:32:54 2006 UTC revision 1.7 by jmc, Sun Aug 13 22:10:25 2006 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4  #include "DIC_OPTIONS.h"  #include "DIC_OPTIONS.h"
5  #include "GCHEM_OPTIONS.h"  c#include "GCHEM_OPTIONS.h"
6    
7  C !INTERFACE: ==========================================================  C !INTERFACE: ==========================================================
8        SUBROUTINE TRACER_MEANAREA(myThid,tracer, kLev,        SUBROUTINE TRACER_MEANAREA(myThid,tracer, kLev,
9       &                    gsm_tracer)       O                    gsm_tracer)
10    
11  C !DESCRIPTION:  C !DESCRIPTION:
12  C  Calculate surface means of tracer  C  Calculate surface means of tracer
# Line 12  C !USES: =============================== Line 15  C !USES: ===============================
15        IMPLICIT NONE        IMPLICIT NONE
16  #include "SIZE.h"  #include "SIZE.h"
17  #include "GRID.h"  #include "GRID.h"
18  #include "DYNVARS.h"  c#include "DYNVARS.h"
19  #include "EEPARAMS.h"  #include "EEPARAMS.h"
20  #include "PARAMS.h"  #include "PARAMS.h"
21    
# Line 20  C !INPUT PARAMETERS: =================== Line 23  C !INPUT PARAMETERS: ===================
23  C  myThid               :: thread number  C  myThid               :: thread number
24        INTEGER myThid        INTEGER myThid
25        INTEGER kLev        INTEGER kLev
26        _RL tracer (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy)        _RL tracer (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
27    
28  C !OUTPUT PARAMETERS: ===================================================  C !OUTPUT PARAMETERS: ===================================================
29  C  gcm_tracer               :: surface mean of tracer  C  gcm_tracer               :: surface mean of tracer
# Line 28  C  gcm_tracer               :: surface m Line 31  C  gcm_tracer               :: surface m
31    
32  C !LOCAL VARIABLES: ====================================================  C !LOCAL VARIABLES: ====================================================
33        INTEGER i,j,bi,bj        INTEGER i,j,bi,bj
34        _RL sumarea        _RL sumArea, sumTracer
35  CEOP  CEOP
36    
37         sumarea= 0. _d 0         sumArea= 0. _d 0
38         gsm_tracer= 0. _d 0         sumTracer= 0. _d 0
39    
40         DO bj=myByLo(myThid),myByHi(myThid)         DO bj=myByLo(myThid),myByHi(myThid)
41          DO bi=myBxLo(myThid),myBxHi(myThid)          DO bi=myBxLo(myThid),myBxHi(myThid)
42           DO i=1,sNx           DO j=1,sNy
43            DO j=1,sNy            DO i=1,sNx
44              sumarea=sumarea+              sumArea=sumArea+
45       &            rA(i,j,bi,bj)*maskC(i,j,kLev,bi,bj)       &            rA(i,j,bi,bj)*maskC(i,j,kLev,bi,bj)
46              gsm_tracer=gsm_tracer+              sumTracer=sumTracer+
47       &          tracer(i,j,kLev,bi,bj)*rA(i,j,bi,bj)*       &          tracer(i,j,kLev,bi,bj)*rA(i,j,bi,bj)*
48       &                        maskC(i,j,kLev,bi,bj)       &                        maskC(i,j,kLev,bi,bj)
49            ENDDO            ENDDO
           ENDDO  
50           ENDDO           ENDDO
51          ENDDO          ENDDO
52          _GLOBAL_SUM_R8(gsm_tracer,myThid)         ENDDO
53          _GLOBAL_SUM_R8(sumarea,myThid)         _GLOBAL_SUM_R8(sumTracer,myThid)
54          gsm_tracer=gsm_tracer/sumarea         _GLOBAL_SUM_R8(sumArea,myThid)
55           gsm_tracer=sumTracer/sumArea
56    
57         RETURN         RETURN
58         END         END
   
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22