/[MITgcm]/MITgcm_contrib/darwin2/pkg/darwin/tracer_meanarea.F
ViewVC logotype

Diff of /MITgcm_contrib/darwin2/pkg/darwin/tracer_meanarea.F

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

revision 1.1 by jahn, Wed Apr 13 18:56:24 2011 UTC revision 1.2 by jahn, Sat Nov 9 21:25:30 2013 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4  #include "DARWIN_OPTIONS.h"  #include "DARWIN_OPTIONS.h"
5    
6  C !INTERFACE: ==========================================================  C !INTERFACE: ==========================================================
7        SUBROUTINE TRACER_MEANAREA(        SUBROUTINE TRACER_MEANAREA(
8       I                    tracer, kLev,       I                    tracer,
9       O                    gsm_tracer,       O                    gsm_tracer,
10       I                    myThid )       I                    myThid )
11    
# Line 12  C  Calculate surface means of tracer Line 15  C  Calculate surface means of tracer
15  C !USES: ===============================================================  C !USES: ===============================================================
16        IMPLICIT NONE        IMPLICIT NONE
17  #include "SIZE.h"  #include "SIZE.h"
 #include "GRID.h"  
18  #include "EEPARAMS.h"  #include "EEPARAMS.h"
19  #include "PARAMS.h"  #include "PARAMS.h"
20    #include "GRID.h"
21    
22  C !INPUT PARAMETERS: ===================================================  C !INPUT PARAMETERS: ===================================================
23  C  myThid               :: thread number  C  myThid               :: thread number
24        INTEGER myThid        INTEGER myThid
       INTEGER kLev  
25        _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)
26    
27  C !OUTPUT PARAMETERS: ===================================================  C !OUTPUT PARAMETERS: ===================================================
28  C  gcm_tracer               :: surface mean of tracer  C  gcm_tracer           :: surface mean of tracer
29        _RL gsm_tracer        _RL gsm_tracer
30    
31  C !LOCAL VARIABLES: ====================================================  C !LOCAL VARIABLES: ====================================================
32        INTEGER i,j,bi,bj        INTEGER i,j,bi,bj
33          INTEGER ks
34          _RL sumTracerTile(nSx,nSy)
35        _RL sumTracer        _RL sumTracer
36  CEOP  CEOP
37    
38         sumTracer= 0. _d 0        IF ( usingPCoords ) THEN
39           ks = Nr
40          ELSE
41           ks = 1
42          ENDIF
43    
44         DO bj=myByLo(myThid),myByHi(myThid)         DO bj=myByLo(myThid),myByHi(myThid)
45          DO bi=myBxLo(myThid),myBxHi(myThid)          DO bi=myBxLo(myThid),myBxHi(myThid)
46             sumTracerTile(bi,bj) = 0. _d 0
47           DO j=1,sNy           DO j=1,sNy
48            DO i=1,sNx            DO i=1,sNx
49              sumTracer=sumTracer+              sumTracerTile(bi,bj) = sumTracerTile(bi,bj)
50       &          tracer(i,j,kLev,bi,bj)*rA(i,j,bi,bj)*       &        + tracer(i,j,ks,bi,bj)*rA(i,j,bi,bj)*maskInC(i,j,bi,bj)
      &                        maskC(i,j,kLev,bi,bj)  
51            ENDDO            ENDDO
52           ENDDO           ENDDO
53          ENDDO          ENDDO
54         ENDDO         ENDDO
55         _GLOBAL_SUM_RL(sumTracer,myThid)         CALL GLOBAL_SUM_TILE_RL(sumTracerTile,sumTracer,myThid)
56         gsm_tracer=sumTracer/globalArea         gsm_tracer = sumTracer/globalArea
57    
58         RETURN        RETURN
59         END        END
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22