1 |
dimitri |
1.1 |
! ========================================================== |
2 |
|
|
! | CPL_MPMICE |
3 |
|
|
! | Declare MPI tags for each variable to be exchanged |
4 |
|
|
! | between MITgcm and MPMice |
5 |
|
|
! |---------------------------------------------------------- |
6 |
|
|
! | Note: North/South/East/West and zonal/meridional below |
7 |
|
|
! | refer to grid, not geographical coordinates. |
8 |
|
|
! ========================================================== |
9 |
|
|
|
10 |
|
|
! TimeIntervalTag |
11 |
|
|
! time interval between communications (s): scalar Real*8 |
12 |
|
|
|
13 |
|
|
! OceanTimeTag |
14 |
|
|
! ocean model time (s): scalar Real*8 |
15 |
|
|
|
16 |
|
|
! IceTimeTag |
17 |
|
|
! ice model time (s): scalar Real*8 |
18 |
|
|
|
19 |
|
|
! OceanGridsizeTag |
20 |
|
|
! ocean model Nx and Ny dimensions: 2 integers |
21 |
|
|
|
22 |
dimitri |
1.6 |
! xCtag |
23 |
|
|
! longitude East of center of grid cell |
24 |
|
|
! Nx * Ny Real*8 array |
25 |
|
|
|
26 |
|
|
! yCtag |
27 |
|
|
! latitude North of center of grid cell |
28 |
|
|
! Nx * Ny Real*8 array |
29 |
|
|
|
30 |
dimitri |
1.5 |
! xGtag |
31 |
|
|
! longitude East of SouthWest corner |
32 |
|
|
! Nx * Ny Real*8 array |
33 |
|
|
|
34 |
|
|
! yGtag |
35 |
|
|
! latitude North of SouthWest corner |
36 |
|
|
! Nx * Ny Real*8 array |
37 |
|
|
|
38 |
|
|
! dxGtag |
39 |
|
|
! distance in m between SouthWest and SouthEast corner |
40 |
|
|
! Nx * Ny Real*8 array |
41 |
|
|
|
42 |
|
|
! dyGtag |
43 |
|
|
! distance in m between SouthWest and NorthEast corner |
44 |
|
|
! Nx * Ny Real*8 array |
45 |
|
|
|
46 |
dimitri |
1.6 |
! aCStag |
47 |
|
|
! cosine(alpha) relative to geographic direction at grid cell center |
48 |
|
|
! Nx * Ny Real*8 array |
49 |
|
|
|
50 |
|
|
! aSNtag |
51 |
|
|
! sine(alpha) relative to geographic direction at grid cell center |
52 |
|
|
! Nx * Ny Real*8 array |
53 |
|
|
|
54 |
dimitri |
1.5 |
! hFacCtag |
55 |
|
|
! landmask of center of grid cell, 0 is land, >0 is ocean |
56 |
|
|
! Nx * Ny Real*8 array |
57 |
|
|
|
58 |
dimitri |
1.1 |
! AreaBcTag |
59 |
|
|
! open boundary ice area (fractional: 0-1) at the tracer points |
60 |
|
|
! 2*(Nx+Ny)-4 Real*8 array |
61 |
|
|
! i=[1:Nx Nx*ones(1,Ny-1) (Nx-1):-1:1 ones(1,Ny-2)]; |
62 |
|
|
! j=[ones(1,Nx) 2:Ny Ny*ones(1,Nx-1) (Ny-1):-1:2 ]; |
63 |
|
|
|
64 |
|
|
! HeffBcTag |
65 |
dimitri |
1.2 |
! open boundary effective ice thickness (m) at the tracer points |
66 |
|
|
! actual ice thickness for single-class ice is: Heff / Area |
67 |
dimitri |
1.1 |
! 2*(Nx+Ny)-4 Real*8 array |
68 |
|
|
! i=[1:Nx Nx*ones(1,Ny-1) (Nx-1):-1:1 ones(1,Ny-2)]; |
69 |
|
|
! j=[ones(1,Nx) 2:Ny Ny*ones(1,Nx-1) (Ny-1):-1:2 ]; |
70 |
|
|
|
71 |
|
|
! HsnowBcTag |
72 |
dimitri |
1.2 |
! open boundary effective snow thickness (m) at the tracer points |
73 |
|
|
! actual snow thickness is: Hsnow / Area |
74 |
dimitri |
1.1 |
! 2*(Nx+Ny)-4 Real*8 array |
75 |
|
|
! i=[1:Nx Nx*ones(1,Ny-1) (Nx-1):-1:1 ones(1,Ny-2)]; |
76 |
|
|
! j=[ones(1,Nx) 2:Ny Ny*ones(1,Nx-1) (Ny-1):-1:2 ]; |
77 |
|
|
|
78 |
|
|
! HsaltBcTag |
79 |
dimitri |
1.2 |
! open boundary effective sea ice salinity (g/m2) at the tracer points |
80 |
|
|
! actual salinity in g/kg is: Hsalt / Heff / rhoIce |
81 |
|
|
! where rhoice = 910 kg / m^3 |
82 |
dimitri |
1.1 |
! 2*(Nx+Ny)-4 Real*8 array |
83 |
|
|
! i=[1:Nx Nx*ones(1,Ny-1) (Nx-1):-1:1 ones(1,Ny-2)]; |
84 |
|
|
! j=[ones(1,Nx) 2:Ny Ny*ones(1,Nx-1) (Ny-1):-1:2 ]; |
85 |
|
|
|
86 |
|
|
! UiceBcTag |
87 |
|
|
! open boundary zonal velocity (m/s) |
88 |
|
|
! specified at the inside edge of the outer grid cells |
89 |
|
|
! southwest C-grid locations for normal component and at the |
90 |
|
|
! southwest C-grid locations for tangential |
91 |
|
|
! >0 from West to East |
92 |
|
|
! 2*(Nx+Ny)-6 Real*8 array |
93 |
|
|
! i=[2:Nx Nx*ones(1,Ny-1) (Nx-1):-1:2 2*ones(1,Ny-2)]; |
94 |
|
|
! j=[ones(1,Nx-1) 2:Ny Ny*ones(1,Nx-2) (Ny-1):-1:2 ]; |
95 |
|
|
|
96 |
|
|
! ViceBcTag |
97 |
|
|
! open boundary meridional velocity (m/s) |
98 |
|
|
! specified at the inside edge of the outer grid cells |
99 |
|
|
! southwest C-grid locations for normal component and at the |
100 |
|
|
! southwest C-grid locations for tangential |
101 |
|
|
! >0 from South to North |
102 |
|
|
! 2*(Nx+Ny)-6 Real*8 array |
103 |
|
|
! i=[1:Nx Nx*ones(1,Ny-2) (Nx-1):-1:1 ones(1,Ny-3)]; |
104 |
|
|
! j=[2*ones(1,Nx) 3:Ny Ny*ones(1,Nx-1) (Ny-1):-1:3 ]; |
105 |
|
|
|
106 |
|
|
! AreaTag |
107 |
|
|
! ice area (fractional: 0-1) at the tracer points |
108 |
|
|
! Nx * Ny Real*8 array |
109 |
|
|
|
110 |
|
|
! HeffTag |
111 |
|
|
! effective ice thickness (mean thickness in the grid box in m) |
112 |
|
|
! at the tracer points |
113 |
dimitri |
1.2 |
! actual ice thickness for single-class ice is: Heff / Area |
114 |
dimitri |
1.1 |
! Nx * Ny Real*8 array |
115 |
|
|
|
116 |
|
|
! HsnowTag |
117 |
|
|
! effective snow thickness (mean thickness in the grid box in m) |
118 |
|
|
! at the tracer points |
119 |
dimitri |
1.2 |
! actual snow thickness is: Hsnow / Area |
120 |
dimitri |
1.1 |
! Nx * Ny Real*8 array |
121 |
|
|
|
122 |
|
|
! HsaltTag |
123 |
dimitri |
1.2 |
! effective sea ice salinity (g/m2) at the tracer points |
124 |
|
|
! actual salinity in g/kg is: Hsalt / Heff / rhoIce |
125 |
|
|
! where rhoice = 910 kg / m^3 |
126 |
dimitri |
1.1 |
! Nx * Ny Real*8 array |
127 |
|
|
|
128 |
dimitri |
1.4 |
! UiceTag |
129 |
|
|
! zonal ice velocity in m/s at South-West C-grid U point |
130 |
|
|
! >0 from West to East |
131 |
|
|
! Nx * Ny Real*8 array |
132 |
|
|
|
133 |
|
|
! ViceTag |
134 |
|
|
! meridional ice velocity in m/s at South-West C-grid V point |
135 |
|
|
! >0 from South to North |
136 |
|
|
! Nx * Ny Real*8 array |
137 |
|
|
|
138 |
dimitri |
1.1 |
! UwindTag |
139 |
|
|
! 10-m u-wind at the tracer points (m/s) |
140 |
|
|
! >0 from West to East |
141 |
|
|
! Nx * Ny Real*8 array |
142 |
|
|
|
143 |
|
|
! VwindTag |
144 |
|
|
! 10-m v-wind at the tracer points (m/s) |
145 |
|
|
! >0 from South to North |
146 |
|
|
! Nx * Ny Real*8 array |
147 |
|
|
|
148 |
|
|
! LwDownTag |
149 |
|
|
! downward longwave at the tracer points (W/m2) |
150 |
|
|
! > 0 for increase in theta (ocean warming) |
151 |
|
|
! Nx * Ny Real*8 array |
152 |
|
|
|
153 |
|
|
! SwDownTag |
154 |
|
|
! downward shortwave at the tracer points (W/m2) |
155 |
|
|
! > 0 for increase in theta (ocean warming) |
156 |
|
|
! Nx * Ny Real*8 array |
157 |
|
|
|
158 |
|
|
! AtempTag |
159 |
|
|
! 2-m air temperature at the tracer points (deg K) |
160 |
|
|
! Nx * Ny Real*8 array |
161 |
|
|
|
162 |
|
|
! AqhTag |
163 |
|
|
! 2-m specific humidity at the tracer points (kg/kg) |
164 |
|
|
! Nx * Ny Real*8 array |
165 |
|
|
|
166 |
|
|
! PrecipTag |
167 |
|
|
! precipitation at the tracer points (m/s) |
168 |
|
|
! > 0 for decrease in salt (ocean salinity) |
169 |
|
|
! Nx * Ny Real*8 array |
170 |
|
|
|
171 |
|
|
! SstTag |
172 |
|
|
! sea surface temperature at the tracer points (deg C) |
173 |
|
|
! Nx * Ny Real*8 array |
174 |
|
|
|
175 |
dimitri |
1.2 |
! SssTag |
176 |
|
|
! sea surface salinity at the tracer points (g/kg) |
177 |
|
|
! Nx * Ny Real*8 array |
178 |
|
|
|
179 |
dimitri |
1.1 |
! UvelTag |
180 |
|
|
! sea surface zonal velocity (m/s) |
181 |
|
|
! at southwest c-grid locations |
182 |
|
|
! Nx * Ny Real*8 array |
183 |
|
|
|
184 |
|
|
! VvelTag |
185 |
|
|
! sea surface meridional velocity (m/s) |
186 |
|
|
! at southwest c-grid locations |
187 |
|
|
! Nx * Ny Real*8 array |
188 |
|
|
|
189 |
|
|
! UstressTag |
190 |
|
|
! zonal stress under ice at southwest c-grid locations (N/m^2) |
191 |
dimitri |
1.2 |
! Ustress is computed for ice-covered area only |
192 |
|
|
! > 0 increases ocean surface velocity in the U direction |
193 |
dimitri |
1.1 |
! Nx * Ny Real*8 array |
194 |
|
|
|
195 |
|
|
! VstressTag |
196 |
|
|
! meridional stress under ice at southwest c-grid locations (N/m^2) |
197 |
dimitri |
1.2 |
! Vstress is computed for ice-covered area only |
198 |
|
|
! > 0 increases ocean surface velocity in the V direction |
199 |
dimitri |
1.1 |
! Nx * Ny Real*8 array |
200 |
|
|
|
201 |
|
|
! SwResidTag |
202 |
|
|
! residual shortwave under ice at tracer point (W/m2) |
203 |
dimitri |
1.2 |
! SwResid is computed for ice-covered area only |
204 |
dimitri |
1.1 |
! > 0 for increase in theta (ocean warming) |
205 |
|
|
! Nx * Ny Real*8 array |
206 |
|
|
|
207 |
|
|
! HeatFluxTag |
208 |
|
|
! heat flux (less shortwave) at tracer point (W/m2) |
209 |
dimitri |
1.2 |
! HeatFlux is computed for ice-covered area only |
210 |
dimitri |
1.1 |
! > 0 for increase in theta (ocean warming) |
211 |
|
|
! Nx * Ny Real*8 array |
212 |
|
|
|
213 |
|
|
! WaterFluxTag |
214 |
dimitri |
1.3 |
! freshwater flux at tracer point ( kg/m^2/s ) |
215 |
dimitri |
1.2 |
! WaterFlux is computed for ice-covered area only |
216 |
dimitri |
1.1 |
! > 0 for decrease in ocean salinity |
217 |
|
|
! Nx * Ny Real*8 array |
218 |
|
|
|
219 |
|
|
! SaltFluxTag |
220 |
|
|
! salt flux at tracer point (g/m^2/s) |
221 |
dimitri |
1.2 |
! SaltFlux is computed for ice-covered area only |
222 |
dimitri |
1.1 |
! > 0 for increase in ocean salinity |
223 |
|
|
! Nx * Ny Real*8 array |
224 |
|
|
|
225 |
|
|
INTEGER TimeIntervalTag |
226 |
|
|
INTEGER OceanTimeTag |
227 |
|
|
INTEGER IceTimeTag |
228 |
|
|
INTEGER OceanGridsizeTag |
229 |
dimitri |
1.6 |
INTEGER xCtag , yCtag |
230 |
|
|
INTEGER xGtag , yGtag |
231 |
|
|
INTEGER aCStag, aSNtag |
232 |
|
|
INTEGER dxGtag, dyGtag |
233 |
dimitri |
1.5 |
INTEGER hFacCtag |
234 |
dimitri |
1.1 |
INTEGER AreaBcTag |
235 |
|
|
INTEGER HeffBcTag |
236 |
|
|
INTEGER HsnowBcTag |
237 |
|
|
INTEGER HsaltBcTag |
238 |
|
|
INTEGER UiceBcTag |
239 |
|
|
INTEGER ViceBcTag |
240 |
|
|
INTEGER AreaTag |
241 |
|
|
INTEGER HeffTag |
242 |
|
|
INTEGER HsnowTag |
243 |
|
|
INTEGER HsaltTag |
244 |
dimitri |
1.4 |
INTEGER UiceTag |
245 |
|
|
INTEGER ViceTag |
246 |
dimitri |
1.1 |
INTEGER UwindTag |
247 |
|
|
INTEGER VwindTag |
248 |
|
|
INTEGER LwDownTag |
249 |
|
|
INTEGER SwDownTag |
250 |
|
|
INTEGER AtempTag |
251 |
|
|
INTEGER AqhTag |
252 |
|
|
INTEGER PrecipTag |
253 |
|
|
INTEGER SstTag |
254 |
dimitri |
1.2 |
INTEGER SssTag |
255 |
dimitri |
1.1 |
INTEGER UvelTag |
256 |
|
|
INTEGER VvelTag |
257 |
|
|
INTEGER UstressTag |
258 |
|
|
INTEGER VstressTag |
259 |
|
|
INTEGER SwResidTag |
260 |
|
|
INTEGER HeatFluxTag |
261 |
|
|
INTEGER WaterFluxTag |
262 |
|
|
INTEGER SaltFluxTag |
263 |
|
|
|
264 |
|
|
PARAMETER ( TimeIntervalTag = 10001000 ) |
265 |
|
|
PARAMETER ( OceanTimeTag = 10001001 ) |
266 |
|
|
PARAMETER ( IceTimeTag = 10001002 ) |
267 |
|
|
PARAMETER ( OceanGridsizeTag = 10002001 ) |
268 |
dimitri |
1.6 |
PARAMETER ( xCtag = 10002002 ) |
269 |
|
|
PARAMETER ( yCtag = 10002003 ) |
270 |
|
|
PARAMETER ( xGtag = 10002004 ) |
271 |
|
|
PARAMETER ( yGtag = 10002005 ) |
272 |
|
|
PARAMETER ( dxGtag = 10002006 ) |
273 |
|
|
PARAMETER ( dyGtag = 10002007 ) |
274 |
|
|
PARAMETER ( aCStag = 10002008 ) |
275 |
|
|
PARAMETER ( aSNtag = 10002009 ) |
276 |
|
|
PARAMETER ( hFacCtag = 10002010 ) |
277 |
dimitri |
1.1 |
PARAMETER ( AreaBcTag = 10003001 ) |
278 |
|
|
PARAMETER ( HeffBcTag = 10003002 ) |
279 |
|
|
PARAMETER ( HsnowBcTag = 10003003 ) |
280 |
|
|
PARAMETER ( HsaltBcTag = 10003004 ) |
281 |
|
|
PARAMETER ( UiceBcTag = 10003005 ) |
282 |
|
|
PARAMETER ( ViceBcTag = 10003006 ) |
283 |
|
|
PARAMETER ( AreaTag = 10004001 ) |
284 |
|
|
PARAMETER ( HeffTag = 10004002 ) |
285 |
|
|
PARAMETER ( HsnowTag = 10004003 ) |
286 |
|
|
PARAMETER ( HsaltTag = 10004004 ) |
287 |
dimitri |
1.4 |
PARAMETER ( UiceTag = 10004005 ) |
288 |
|
|
PARAMETER ( ViceTag = 10004006 ) |
289 |
dimitri |
1.1 |
PARAMETER ( UwindTag = 10005001 ) |
290 |
|
|
PARAMETER ( VwindTag = 10005002 ) |
291 |
|
|
PARAMETER ( LwDownTag = 10005101 ) |
292 |
|
|
PARAMETER ( SwDownTag = 10005102 ) |
293 |
|
|
PARAMETER ( AtempTag = 10005103 ) |
294 |
|
|
PARAMETER ( AqhTag = 10005104 ) |
295 |
|
|
PARAMETER ( PrecipTag = 10005105 ) |
296 |
|
|
PARAMETER ( SstTag = 10005201 ) |
297 |
dimitri |
1.2 |
PARAMETER ( SssTag = 10005202 ) |
298 |
|
|
PARAMETER ( UvelTag = 10005203 ) |
299 |
|
|
PARAMETER ( VvelTag = 10005204 ) |
300 |
dimitri |
1.1 |
PARAMETER ( UstressTag = 10006001 ) |
301 |
|
|
PARAMETER ( VstressTag = 10006002 ) |
302 |
|
|
PARAMETER ( SwResidTag = 20006101 ) |
303 |
|
|
PARAMETER ( HeatFluxTag = 20006102 ) |
304 |
|
|
PARAMETER ( WaterFluxTag = 20006103 ) |
305 |
|
|
PARAMETER ( SaltFluxTag = 20006104 ) |