beaufort/code/CPL_MPMICE.h

!     ==========================================================
!     | CPL_MPMICE
!     | Declare MPI tags for each variable to be exchanged
!     | between MITgcm and MPMice
!     |----------------------------------------------------------
!     | Note: North/South/East/West and zonal/meridional below
!     | refer to grid, not geographical coordinates.
!     ==========================================================

!     TimeIntervalTag
!     time interval between communications (s): scalar Real*8

!     OceanTimeTag
!     ocean model time (s): scalar Real*8

!     IceTimeTag
!     ice model time (s): scalar Real*8

!     OceanGridsizeTag
!     ocean model Nx and Ny dimensions: 2 integers

!     xCtag
!     longitude East of center of grid cell
!     Nx * Ny Real*8 array

!     yCtag
!     latitude North of center of grid cell
!     Nx * Ny Real*8 array

!     xGtag
!     longitude East of SouthWest corner
!     Nx * Ny Real*8 array

!     yGtag
!     latitude North of SouthWest corner
!     Nx * Ny Real*8 array

!     dxGtag
!     distance in m between SouthWest and SouthEast corner
!     Nx * Ny Real*8 array

!     dyGtag
!     distance in m between SouthWest and NorthEast corner
!     Nx * Ny Real*8 array

!     aCStag
!     cosine(alpha) relative to geographic direction at grid cell center
!     Nx * Ny Real*8 array

!     aSNtag
!     sine(alpha) relative to geographic direction at grid cell center
!     Nx * Ny Real*8 array

!     hFacCtag
!     landmask of center of grid cell, 0 is land, >0 is ocean
!     Nx * Ny Real*8 array

!     AreaBcTag
!     open boundary ice area (fractional: 0-1) at the tracer points
!     2*(Nx+Ny)-4 Real*8 array
!     i=[1:Nx       Nx*ones(1,Ny-1) (Nx-1):-1:1     ones(1,Ny-2)];
!     j=[ones(1,Nx) 2:Ny            Ny*ones(1,Nx-1) (Ny-1):-1:2 ];

!     HeffBcTag
!     open boundary effective ice thickness (m) at the tracer points
!     actual ice thickness for single-class ice is: Heff / Area
!     2*(Nx+Ny)-4 Real*8 array
!     i=[1:Nx       Nx*ones(1,Ny-1) (Nx-1):-1:1     ones(1,Ny-2)];
!     j=[ones(1,Nx) 2:Ny            Ny*ones(1,Nx-1) (Ny-1):-1:2 ];

!     HsnowBcTag
!     open boundary effective snow thickness (m) at the tracer points
!     actual snow thickness is: Hsnow / Area 
!     2*(Nx+Ny)-4 Real*8 array
!     i=[1:Nx       Nx*ones(1,Ny-1) (Nx-1):-1:1     ones(1,Ny-2)];
!     j=[ones(1,Nx) 2:Ny            Ny*ones(1,Nx-1) (Ny-1):-1:2 ];

!     HsaltBcTag
!     open boundary effective sea ice salinity (g/m2) at the tracer points
!     actual salinity in g/kg is: Hsalt / Heff / rhoIce
!     where rhoice = 910 kg / m^3
!     2*(Nx+Ny)-4 Real*8 array
!     i=[1:Nx       Nx*ones(1,Ny-1) (Nx-1):-1:1     ones(1,Ny-2)];
!     j=[ones(1,Nx) 2:Ny            Ny*ones(1,Nx-1) (Ny-1):-1:2 ];

!     UiceBcTag
!     open boundary zonal velocity (m/s)
!     specified at the inside edge of the outer grid cells
!     southwest C-grid locations for normal component and at the
!     southwest C-grid locations for tangential
!     >0 from West to East
!     2*(Nx+Ny)-6 Real*8 array
!     i=[2:Nx         Nx*ones(1,Ny-1) (Nx-1):-1:2     2*ones(1,Ny-2)];
!     j=[ones(1,Nx-1) 2:Ny            Ny*ones(1,Nx-2) (Ny-1):-1:2   ];

!     ViceBcTag
!     open boundary meridional velocity (m/s)
!     specified at the inside edge of the outer grid cells
!     southwest C-grid locations for normal component and at the
!     southwest C-grid locations for tangential
!     >0 from South to North
!     2*(Nx+Ny)-6 Real*8 array
!     i=[1:Nx         Nx*ones(1,Ny-2) (Nx-1):-1:1     ones(1,Ny-3)];
!     j=[2*ones(1,Nx) 3:Ny            Ny*ones(1,Nx-1) (Ny-1):-1:3 ];

!     AreaTag
!     ice area (fractional: 0-1) at the tracer points
!     Nx * Ny Real*8 array

!     HeffTag
!     effective ice thickness (mean thickness in the grid box in m)
!     at the tracer points
!     actual ice thickness for single-class ice is: Heff / Area
!     Nx * Ny Real*8 array

!     HsnowTag
!     effective snow thickness (mean thickness in the grid box in m)
!     at the tracer points
!     actual snow thickness is: Hsnow / Area
!     Nx * Ny Real*8 array

!     HsaltTag
!     effective sea ice salinity (g/m2) at the tracer points
!     actual salinity in g/kg is: Hsalt / Heff / rhoIce
!     where rhoice = 910 kg / m^3
!     Nx * Ny Real*8 array

!     UiceTag
!     zonal ice velocity in m/s at South-West C-grid U point
!     >0 from West to East
!     Nx * Ny Real*8 array

!     ViceTag
!     meridional ice velocity in m/s at South-West C-grid V point
!     >0 from South to North
!     Nx * Ny Real*8 array

!     UwindTag
!     10-m u-wind at the tracer points (m/s)
!     >0 from West to East
!     Nx * Ny Real*8 array

!     VwindTag
!     10-m v-wind at the tracer points (m/s)
!     >0 from South to North
!     Nx * Ny Real*8 array

!     LwDownTag
!     downward longwave at the tracer points (W/m2)
!     > 0 for increase in theta (ocean warming)
!     Nx * Ny Real*8 array

!     SwDownTag
!     downward shortwave at the tracer points (W/m2)
!     > 0 for increase in theta (ocean warming)
!     Nx * Ny Real*8 array

!     AtempTag
!     2-m air temperature at the tracer points (deg K)
!     Nx * Ny Real*8 array

!     AqhTag
!     2-m specific humidity at the tracer points (kg/kg)
!     Nx * Ny Real*8 array

!     PrecipTag
!     precipitation at the tracer points (m/s)
!     > 0 for decrease in salt (ocean salinity)
!     Nx * Ny Real*8 array

!     SstTag
!     sea surface temperature at the tracer points (deg C)
!     Nx * Ny Real*8 array

!     SssTag
!     sea surface salinity at the tracer points (g/kg)
!     Nx * Ny Real*8 array

!     UvelTag
!     sea surface zonal velocity (m/s)
!     at southwest c-grid locations
!     Nx * Ny Real*8 array

!     VvelTag
!     sea surface meridional velocity (m/s)
!     at southwest c-grid locations
!     Nx * Ny Real*8 array

!     UstressTag
!     zonal stress under ice at southwest c-grid locations (N/m^2)
!     Ustress is computed for ice-covered area only
!     > 0 increases ocean surface velocity in the U direction
!     Nx * Ny Real*8 array

!     VstressTag
!     meridional stress under ice at southwest c-grid locations (N/m^2)
!     Vstress is computed for ice-covered area only
!     > 0 increases ocean surface velocity in the V direction
!     Nx * Ny Real*8 array

!     SwResidTag
!     residual shortwave under ice at tracer point (W/m2)
!     SwResid is computed for ice-covered area only
!     > 0 for increase in theta (ocean warming)
!     Nx * Ny Real*8 array

!     HeatFluxTag
!     heat flux (less shortwave) at tracer point (W/m2)
!     HeatFlux is computed for ice-covered area only
!     > 0 for increase in theta (ocean warming)
!     Nx * Ny Real*8 array

!     WaterFluxTag
!     freshwater flux at tracer point ( kg/m^2/s )
!     WaterFlux is computed for ice-covered area only
!     > 0 for decrease in ocean salinity
!     Nx * Ny Real*8 array

!     SaltFluxTag
!     salt flux at tracer point (g/m^2/s)
!     SaltFlux is computed for ice-covered area only
!     > 0 for increase in ocean salinity
!     Nx * Ny Real*8 array

      INTEGER     TimeIntervalTag
      INTEGER     OceanTimeTag
      INTEGER     IceTimeTag
      INTEGER     OceanGridsizeTag
      INTEGER     xCtag , yCtag
      INTEGER     xGtag , yGtag
      INTEGER     aCStag, aSNtag
      INTEGER     dxGtag, dyGtag
      INTEGER     hFacCtag
      INTEGER     AreaBcTag
      INTEGER     HeffBcTag
      INTEGER     HsnowBcTag
      INTEGER     HsaltBcTag
      INTEGER     UiceBcTag
      INTEGER     ViceBcTag
      INTEGER     AreaTag
      INTEGER     HeffTag
      INTEGER     HsnowTag
      INTEGER     HsaltTag
      INTEGER     UiceTag
      INTEGER     ViceTag
      INTEGER     UwindTag
      INTEGER     VwindTag
      INTEGER     LwDownTag
      INTEGER     SwDownTag
      INTEGER     AtempTag
      INTEGER     AqhTag
      INTEGER     PrecipTag
      INTEGER     SstTag
      INTEGER     SssTag
      INTEGER     UvelTag
      INTEGER     VvelTag
      INTEGER     UstressTag
      INTEGER     VstressTag
      INTEGER     SwResidTag
      INTEGER     HeatFluxTag
      INTEGER     WaterFluxTag
      INTEGER     SaltFluxTag

      PARAMETER ( TimeIntervalTag  = 10001000 )
      PARAMETER ( OceanTimeTag     = 10001001 )
      PARAMETER ( IceTimeTag       = 10001002 )
      PARAMETER ( OceanGridsizeTag = 10002001 )
      PARAMETER ( xCtag            = 10002002 )
      PARAMETER ( yCtag            = 10002003 )
      PARAMETER ( xGtag            = 10002004 )
      PARAMETER ( yGtag            = 10002005 )
      PARAMETER ( dxGtag           = 10002006 )
      PARAMETER ( dyGtag           = 10002007 )
      PARAMETER ( aCStag           = 10002008 )
      PARAMETER ( aSNtag           = 10002009 )
      PARAMETER ( hFacCtag         = 10002010 )
      PARAMETER ( AreaBcTag        = 10003001 )
      PARAMETER ( HeffBcTag        = 10003002 )
      PARAMETER ( HsnowBcTag       = 10003003 )
      PARAMETER ( HsaltBcTag       = 10003004 )
      PARAMETER ( UiceBcTag        = 10003005 )
      PARAMETER ( ViceBcTag        = 10003006 )
      PARAMETER ( AreaTag          = 10004001 )
      PARAMETER ( HeffTag          = 10004002 )
      PARAMETER ( HsnowTag         = 10004003 )
      PARAMETER ( HsaltTag         = 10004004 )
      PARAMETER ( UiceTag          = 10004005 )
      PARAMETER ( ViceTag          = 10004006 )
      PARAMETER ( UwindTag         = 10005001 )
      PARAMETER ( VwindTag         = 10005002 )
      PARAMETER ( LwDownTag        = 10005101 )
      PARAMETER ( SwDownTag        = 10005102 )
      PARAMETER ( AtempTag         = 10005103 )
      PARAMETER ( AqhTag           = 10005104 )
      PARAMETER ( PrecipTag        = 10005105 )
      PARAMETER ( SstTag           = 10005201 )
      PARAMETER ( SssTag           = 10005202 )
      PARAMETER ( UvelTag          = 10005203 )
      PARAMETER ( VvelTag          = 10005204 )
      PARAMETER ( UstressTag       = 10006001 )
      PARAMETER ( VstressTag       = 10006002 )
      PARAMETER ( SwResidTag       = 20006101 )
      PARAMETER ( HeatFluxTag      = 20006102 )
      PARAMETER ( WaterFluxTag     = 20006103 )
      PARAMETER ( SaltFluxTag      = 20006104 )
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	! 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	! 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	! 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	! hFacCtag
55	! landmask of center of grid cell, 0 is land, >0 is ocean
56	! Nx * Ny Real*8 array
57
58	! AreaBcTag
59	! open boundary ice area (fractional: 0-1) at the tracer points
60	! 2(Nx+Ny)-4 Real8 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	! open boundary effective ice thickness (m) at the tracer points
66	! actual ice thickness for single-class ice is: Heff / Area
67	! 2(Nx+Ny)-4 Real8 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	! open boundary effective snow thickness (m) at the tracer points
73	! actual snow thickness is: Hsnow / Area
74	! 2(Nx+Ny)-4 Real8 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	! 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	! 2(Nx+Ny)-4 Real8 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 Real8 array
93	! i=[2:Nx Nxones(1,Ny-1) (Nx-1):-1:2 2ones(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 Real8 array
103	! i=[1:Nx Nx*ones(1,Ny-2) (Nx-1):-1:1 ones(1,Ny-3)];
104	! j=[2ones(1,Nx) 3:Ny Nyones(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	! actual ice thickness for single-class ice is: Heff / Area
114	! 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	! actual snow thickness is: Hsnow / Area
120	! Nx * Ny Real*8 array
121
122	! HsaltTag
123	! 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	! Nx * Ny Real*8 array
127
128	! 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	! 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	! SssTag
176	! sea surface salinity at the tracer points (g/kg)
177	! Nx * Ny Real*8 array
178
179	! 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	! Ustress is computed for ice-covered area only
192	! > 0 increases ocean surface velocity in the U direction
193	! Nx * Ny Real*8 array
194
195	! VstressTag
196	! meridional stress under ice at southwest c-grid locations (N/m^2)
197	! Vstress is computed for ice-covered area only
198	! > 0 increases ocean surface velocity in the V direction
199	! Nx * Ny Real*8 array
200
201	! SwResidTag
202	! residual shortwave under ice at tracer point (W/m2)
203	! SwResid is computed for ice-covered area only
204	! > 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	! HeatFlux is computed for ice-covered area only
210	! > 0 for increase in theta (ocean warming)
211	! Nx * Ny Real*8 array
212
213	! WaterFluxTag
214	! freshwater flux at tracer point ( kg/m^2/s )
215	! WaterFlux is computed for ice-covered area only
216	! > 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	! SaltFlux is computed for ice-covered area only
222	! > 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	INTEGER xCtag , yCtag
230	INTEGER xGtag , yGtag
231	INTEGER aCStag, aSNtag
232	INTEGER dxGtag, dyGtag
233	INTEGER hFacCtag
234	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	INTEGER UiceTag
245	INTEGER ViceTag
246	INTEGER UwindTag
247	INTEGER VwindTag
248	INTEGER LwDownTag
249	INTEGER SwDownTag
250	INTEGER AtempTag
251	INTEGER AqhTag
252	INTEGER PrecipTag
253	INTEGER SstTag
254	INTEGER SssTag
255	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	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	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	PARAMETER ( UiceTag = 10004005 )
288	PARAMETER ( ViceTag = 10004006 )
289	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	PARAMETER ( SssTag = 10005202 )
298	PARAMETER ( UvelTag = 10005203 )
299	PARAMETER ( VvelTag = 10005204 )
300	PARAMETER ( UstressTag = 10006001 )
301	PARAMETER ( VstressTag = 10006002 )
302	PARAMETER ( SwResidTag = 20006101 )
303	PARAMETER ( HeatFluxTag = 20006102 )
304	PARAMETER ( WaterFluxTag = 20006103 )
305	PARAMETER ( SaltFluxTag = 20006104 )