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

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

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

Revision 1.5 - (hide annotations) (download)
Thu Aug 5 21:50:26 2004 UTC (20 years, 1 month ago) by stephd
Branch: MAIN
CVS Tags: checkpoint57t_post, checkpoint58l_post, checkpoint57o_post, checkpoint58e_post, checkpoint57v_post, checkpoint57f_post, checkpoint55c_post, checkpoint57s_post, checkpoint57k_post, checkpoint55d_pre, checkpoint57d_post, checkpoint57g_post, checkpoint57b_post, checkpoint57c_pre, checkpoint55j_post, checkpoint56b_post, checkpoint57i_post, checkpoint57y_post, checkpoint58g_post, checkpoint57x_post, checkpoint57m_post, checkpoint55h_post, checkpoint57g_pre, checkpoint58h_post, checkpoint57e_post, checkpoint56c_post, checkpoint54e_post, checkpoint58j_post, checkpoint55b_post, checkpoint57h_post, checkpoint57y_pre, checkpoint55, checkpoint57f_pre, checkpoint57a_post, checkpoint54f_post, checkpoint55g_post, checkpoint55f_post, checkpoint57r_post, checkpoint58, checkpoint57a_pre, checkpoint55i_post, checkpoint57, checkpoint56, eckpoint57e_pre, checkpoint57h_done, checkpoint58f_post, checkpoint57n_post, checkpoint58d_post, checkpoint58c_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint58a_post, checkpoint58i_post, checkpoint57q_post, checkpoint57z_post, checkpoint57c_post, checkpoint55e_post, checkpoint58k_post, checkpoint55a_post, checkpoint57j_post, checkpoint58b_post, checkpoint57h_pre, checkpoint58m_post, checkpoint57l_post, checkpoint56a_post, checkpoint55d_post
Changes since 1.4: +2 -0 lines 
o correct bug in mean calculation
1 edhill 1.3 #include "DIC_OPTIONS.h"
2 stephd 1.1 #include "GCHEM_OPTIONS.h"
3    
4     C !INTERFACE: ==========================================================
5 stephd 1.2 SUBROUTINE TRACER_MEANAREA(myThid,tracer, kLev,
6 stephd 1.1 & gsm_tracer)
7 stephd 1.4
8     C !DESCRIPTION:
9     C Calculate surface means of tracer
10    
11     C !USES: ===============================================================
12 stephd 1.1 IMPLICIT NONE
13     #include "SIZE.h"
14     #include "GRID.h"
15     #include "DYNVARS.h"
16     #include "EEPARAMS.h"
17     #include "PARAMS.h"
18 stephd 1.4
19 stephd 1.1 C !INPUT PARAMETERS: ===================================================
20     C myThid :: thread number
21 stephd 1.2 INTEGER myThid
22 stephd 1.1 INTEGER kLev
23     _RL tracer (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy)
24 stephd 1.4
25     C !OUTPUT PARAMETERS: ===================================================
26     C gcm_tracer :: surface mean of tracer
27 stephd 1.1 _RL gsm_tracer
28    
29 stephd 1.4 C !LOCAL VARIABLES: ====================================================
30 stephd 1.1 INTEGER i,j,bi,bj
31     _RL sumarea
32 stephd 1.4 CEOP
33 stephd 1.1
34     _BEGIN_MASTER(myThid)
35    
36     sumarea= 0. _d 0
37     gsm_tracer= 0. _d 0
38    
39     DO bj=myByLo(myThid),myByHi(myThid)
40     DO bi=myBxLo(myThid),myBxHi(myThid)
41     DO i=1,sNx
42     DO j=1,sNy
43     sumarea=sumarea+
44     & rA(i,j,bi,bj)*maskC(i,j,kLev,bi,bj)
45     gsm_tracer=gsm_tracer+
46     & tracer(i,j,kLev,bi,bj)*rA(i,j,bi,bj)*
47     & maskC(i,j,kLev,bi,bj)
48     ENDDO
49     ENDDO
50     ENDDO
51     ENDDO
52 stephd 1.5 _GLOBAL_SUM_R8(gsm_tracer,myThid)
53     _GLOBAL_SUM_R8(sumarea,myThid)
54 stephd 1.1 gsm_tracer=gsm_tracer/sumarea
55    
56     _END_MASTER(myThid)
57    
58     C-- Everyone else must wait
59     _BARRIER
60    
61    
62     RETURN
63     END
64    
  ViewVC Help
Powered by ViewVC 1.1.22