/[MITgcm]/MITgcm/eesupp/src/global_max.F
ViewVC logotype

Contents of /MITgcm/eesupp/src/global_max.F

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

Revision 1.9 - (show annotations) (download)
Fri Sep 21 03:54:35 2001 UTC (24 years, 1 month ago) by cnh
Branch: MAIN
CVS Tags: checkpoint46n_post, checkpoint51k_post, checkpoint47e_post, checkpoint44e_post, checkpoint46l_post, checkpoint46g_pre, checkpoint47c_post, release1_p13_pre, checkpoint50c_post, checkpoint46f_post, checkpoint52d_pre, checkpoint48e_post, checkpoint50c_pre, checkpoint44f_post, checkpoint46b_post, checkpoint43a-release1mods, checkpoint51o_pre, ecco_c50_e32, ecco_c50_e33, ecco_c50_e30, ecco_c50_e31, release1_p13, checkpoint51l_post, checkpoint48i_post, checkpoint46l_pre, chkpt44d_post, checkpoint51, checkpoint50, checkpoint52, release1_p8, release1_p9, checkpoint50d_post, release1_p1, release1_p2, release1_p3, release1_p4, release1_p5, release1_p6, release1_p7, checkpoint50b_pre, checkpoint44e_pre, checkpoint51f_post, release1_b1, ecco_c51_e34d, ecco_c51_e34e, ecco_c51_e34f, ecco_c51_e34g, ecco_c51_e34a, ecco_c51_e34b, ecco_c51_e34c, checkpoint48b_post, checkpoint43, checkpoint51d_post, checkpoint48c_pre, checkpoint47d_pre, checkpoint51t_post, checkpoint51n_post, release1_chkpt44d_post, checkpoint51s_post, checkpoint47a_post, checkpoint48d_pre, checkpoint51j_post, checkpoint47i_post, release1_p11, checkpoint51n_pre, checkpoint47d_post, icebear5, icebear4, icebear3, icebear2, checkpoint46d_pre, checkpoint48d_post, release1-branch_tutorials, checkpoint48f_post, checkpoint45d_post, checkpoint52b_pre, checkpoint46j_pre, checkpoint51l_pre, chkpt44a_post, checkpoint44h_pre, checkpoint48h_post, checkpoint51q_post, ecco_c50_e29, checkpoint51b_pre, checkpoint46a_post, checkpoint47g_post, checkpoint52b_post, checkpoint52c_post, checkpoint46j_post, checkpoint51h_pre, checkpoint46k_post, ecco_c50_e28, chkpt44c_pre, checkpoint48a_post, checkpoint45a_post, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, ecco_c44_e19, ecco_c44_e18, ecco_c44_e17, ecco_c44_e16, release1_p12, release1_p10, release1_p16, release1_p17, release1_p14, release1_p15, checkpoint47j_post, ecco_c50_e33a, branch-exfmods-tag, checkpoint44g_post, branchpoint-genmake2, checkpoint46e_pre, checkpoint51r_post, checkpoint48c_post, checkpoint45b_post, checkpoint46b_pre, checkpoint51i_post, release1-branch-end, release1_final_v1, checkpoint51b_post, checkpoint51c_post, checkpoint46c_pre, checkpoint46, checkpoint47b_post, checkpoint44b_post, ecco_c51_e34, checkpoint46h_pre, checkpoint46m_post, checkpoint46a_pre, checkpoint50g_post, checkpoint45c_post, ecco_ice2, ecco_ice1, checkpoint44h_post, checkpoint46g_post, release1_p12_pre, checkpoint52a_pre, ecco_c44_e22, checkpoint50h_post, checkpoint50e_pre, checkpoint50i_post, ecco_c44_e25, checkpoint51i_pre, checkpoint47f_post, checkpoint50e_post, chkpt44a_pre, checkpoint46i_post, ecco_c44_e23, ecco_c44_e20, ecco_c44_e21, ecco_c44_e26, ecco_c44_e27, ecco_c44_e24, checkpoint46c_post, ecco-branch-mod1, ecco-branch-mod2, ecco-branch-mod3, ecco-branch-mod4, ecco-branch-mod5, branch-netcdf, checkpoint50d_pre, checkpoint46e_post, release1_beta1, checkpoint51e_post, checkpoint44b_pre, checkpoint42, checkpoint41, checkpoint47, checkpoint44, checkpoint45, checkpoint48, checkpoint49, checkpoint46h_post, checkpoint51o_post, checkpoint51f_pre, chkpt44c_post, checkpoint48g_post, checkpoint47h_post, checkpoint52a_post, checkpoint44f_pre, checkpoint51g_post, ecco_c52_e35, checkpoint46d_post, checkpoint50b_post, release1-branch_branchpoint, checkpoint51m_post, checkpoint51a_post, checkpoint51p_post, checkpoint51u_post
Branch point for: c24_e25_ice, branch-exfmods-curt, release1_final, release1-branch, branch-genmake2, release1, branch-nonh, tg2-branch, ecco-branch, release1_50yr, netcdf-sm0, icebear, checkpoint51n_branch, release1_coupled
Changes since 1.8: +56 -41 lines 
Starting to bring comments up to date and format comments
for document extraction of "prototypes".
1 C $Header: /u/gcmpack/models/MITgcmUV/eesupp/src/global_max.F,v 1.8 2001/05/29 14:01:36 adcroft Exp $
2 C $Name: $
3
4 C-- File global_max.F: Routines that perform global max reduction on an array
5 C of thread values.
6 C Contents
7 C o global_max_r4
8 C o global_max_r8
9 #include "CPP_EEOPTIONS.h"
10
11 CBOP
12
13 C !ROUTINE: GLOBAL_MAX_R4
14
15 C !INTERFACE:
16 SUBROUTINE GLOBAL_MAX_R4(
17 U maxPhi,
18 I myThid )
19 IMPLICIT NONE
20 C !DESCRIPTION:
21 C *==========================================================*
22 C | SUBROUTINE GLOBAL_MAX_R4
23 C | o Handle max for real*4 data.
24 C *==========================================================*
25 C | Perform max on array of one value per thread and then
26 C | max result of all the processes.
27 C | Notes
28 C | =====
29 C | Within a process only one thread does the max, each
30 C | thread is assumed to have already maxed its local data.
31 C | The same thread also does the inter-process max for
32 C | example with MPI and then writes the result into a shared
33 C | location. All threads wait until the max is avaiailable.
34 C *==========================================================*
35
36 C !USES:
37 C == Global data ==
38 #include "SIZE.h"
39 #include "EEPARAMS.h"
40 #include "EESUPPORT.h"
41 #include "GLOBAL_MAX.h"
42
43 C !INPUT/OUTPUT PARAMETERS:
44 C == Routine arguments ==
45 C maxPhi :: Result of max.
46 C myThid :: My thread id.
47 Real*4 maxPhi
48 INTEGER myThid
49
50 C !LOCAL VARIABLES:
51 C == Local variables ==
52 C I :: Loop counters
53 C mpiRC :: MPI return code
54 INTEGER I
55 Real*4 tmp
56 #ifdef ALLOW_USE_MPI
57 INTEGER mpiRC
58 #endif /* ALLOW_USE_MPI */
59 CEOP
60
61 CALL BAR2( myThid )
62 C-- write local max into array
63 phiGMRS(1,myThid) = maxPhi
64
65 C-- Can not start until everyone is ready
66 CALL BAR2( myThid )
67
68 C-- Max within the process first
69 _BEGIN_MASTER( myThid )
70 tmp = phiGMRS(1,1)
71 DO I=2,nThreads
72 tmp = MAX(tmp,phiGMRS(1,I))
73 ENDDO
74 maxPhi = tmp
75 #ifdef ALLOW_USE_MPI
76 #ifndef ALWAYS_USE_MPI
77 IF ( usingMPI ) THEN
78 #endif
79 CALL MPI_Allreduce(tmp,maxPhi,1,MPI_REAL,MPI_MAX,
80 & MPI_COMM_WORLD,mpiRC)
81 #ifndef ALWAYS_USE_MPI
82 ENDIF
83 #endif
84 #endif /* ALLOW_USE_MPI */
85 phiGMRS(1,1) = maxPhi
86 _END_MASTER( myThid )
87 C--
88 CALL BAR2( myThid )
89
90 C-- set result for every process
91 maxPhi = phiGMRS(1,1)
92 CALL BAR2( myThid )
93
94 RETURN
95 END
96
97 CBOP
98
99 C !ROUTINE: GLOBAL_MAX_R8
100
101 C !INTERFACE:
102 SUBROUTINE GLOBAL_MAX_R8(
103 O maxPhi,
104 I myThid )
105 IMPLICIT NONE
106 C !DESCRIPTION:
107 C *==========================================================*
108 C | SUBROUTINE GLOBAL_MAX_R8
109 C | o Handle max for real*8 data.
110 C *==========================================================*
111 C | Perform max on array of one value per thread and then
112 C | max result of all the processes.
113 C | Notes
114 C | =====
115 C | Within a process only one thread does the max, each
116 C | thread is assumed to have already maxed its local data.
117 C | The same thread also does the inter-process max for
118 C | example with MPI and then writes the result into a shared
119 C | location. All threads wait until the max is avaiailable.
120 C *==========================================================*
121
122 C !USES:
123 C === Global data ===
124 #include "SIZE.h"
125 #include "EEPARAMS.h"
126 #include "EESUPPORT.h"
127 #include "GLOBAL_MAX.h"
128
129 C !INPUT/OUTPUT PARAMETERS:
130 C === Routine arguments ===
131 C maxPhi :: Result of max.
132 C myThid :: My thread id.
133 Real*8 maxPhi
134 INTEGER myThid
135
136 C !LOCAL VARIABLES:
137 C === Local variables ===
138 C I :: Loop counters
139 C mpiRC :: MPI return code
140 INTEGER I
141 Real*8 tmp
142 #ifdef ALLOW_USE_MPI
143 INTEGER mpiRC
144 #endif /* ALLOW_USE_MPI */
145 CEOP
146
147 CALL BAR2( myThid )
148 C-- write local max into array
149 phiGMRL(1,myThid) = maxPhi
150
151 C-- Can not start until everyone is ready
152 CALL BAR2( myThid )
153
154 C-- Max within the process first
155 _BEGIN_MASTER( myThid )
156 tmp = phiGMRL(1,1)
157 DO I=2,nThreads
158 tmp = MAX(tmp,phiGMRL(1,I))
159 ENDDO
160 maxPhi = tmp
161 #ifdef ALLOW_USE_MPI
162 #ifndef ALWAYS_USE_MPI
163 IF ( usingMPI ) THEN
164 #endif
165 CALL MPI_Allreduce(tmp,maxPhi,1,MPI_DOUBLE_PRECISION,MPI_MAX,
166 & MPI_COMM_WORLD,mpiRC)
167 #ifndef ALWAYS_USE_MPI
168 ENDIF
169 #endif
170 #endif /* ALLOW_USE_MPI */
171 C-- Write solution to place where all threads can see it
172 phiGMRL(1,1) = maxPhi
173 _END_MASTER( myThid )
174
175 C-- Do not leave until we are sure that the max is done
176 CALL BAR2( myThid )
177
178 C-- set result for every process
179 maxPhi = phiGMRL(1,1)
180 CALL BAR2( myThid )
181
182 RETURN
183 END
  ViewVC Help
Powered by ViewVC 1.1.22