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

!     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

!     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     AreaBcTag
      INTEGER     HeffBcTag
      INTEGER     HsnowBcTag
      INTEGER     HsaltBcTag
      INTEGER     UiceBcTag
      INTEGER     ViceBcTag
      INTEGER     AreaTag
      INTEGER     HeffTag
      INTEGER     HsnowTag
      INTEGER     HsaltTag
      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 ( 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 ( 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	! AreaBcTag
23	! open boundary ice area (fractional: 0-1) at the tracer points
24	! 2(Nx+Ny)-4 Real8 array
25	! i=[1:Nx Nx*ones(1,Ny-1) (Nx-1):-1:1 ones(1,Ny-2)];
26	! j=[ones(1,Nx) 2:Ny Ny*ones(1,Nx-1) (Ny-1):-1:2 ];
27
28	! HeffBcTag
29	! open boundary effective ice thickness (m) at the tracer points
30	! actual ice thickness for single-class ice is: Heff / Area
31	! 2(Nx+Ny)-4 Real8 array
32	! i=[1:Nx Nx*ones(1,Ny-1) (Nx-1):-1:1 ones(1,Ny-2)];
33	! j=[ones(1,Nx) 2:Ny Ny*ones(1,Nx-1) (Ny-1):-1:2 ];
34
35	! HsnowBcTag
36	! open boundary effective snow thickness (m) at the tracer points
37	! actual snow thickness is: Hsnow / Area
38	! 2(Nx+Ny)-4 Real8 array
39	! i=[1:Nx Nx*ones(1,Ny-1) (Nx-1):-1:1 ones(1,Ny-2)];
40	! j=[ones(1,Nx) 2:Ny Ny*ones(1,Nx-1) (Ny-1):-1:2 ];
41
42	! HsaltBcTag
43	! open boundary effective sea ice salinity (g/m2) at the tracer points
44	! actual salinity in g/kg is: Hsalt / Heff / rhoIce
45	! where rhoice = 910 kg / m^3
46	! 2(Nx+Ny)-4 Real8 array
47	! i=[1:Nx Nx*ones(1,Ny-1) (Nx-1):-1:1 ones(1,Ny-2)];
48	! j=[ones(1,Nx) 2:Ny Ny*ones(1,Nx-1) (Ny-1):-1:2 ];
49
50	! UiceBcTag
51	! open boundary zonal velocity (m/s)
52	! specified at the inside edge of the outer grid cells
53	! southwest C-grid locations for normal component and at the
54	! southwest C-grid locations for tangential
55	! >0 from West to East
56	! 2(Nx+Ny)-6 Real8 array
57	! i=[2:Nx Nxones(1,Ny-1) (Nx-1):-1:2 2ones(1,Ny-2)];
58	! j=[ones(1,Nx-1) 2:Ny Ny*ones(1,Nx-2) (Ny-1):-1:2 ];
59
60	! ViceBcTag
61	! open boundary meridional velocity (m/s)
62	! specified at the inside edge of the outer grid cells
63	! southwest C-grid locations for normal component and at the
64	! southwest C-grid locations for tangential
65	! >0 from South to North
66	! 2(Nx+Ny)-6 Real8 array
67	! i=[1:Nx Nx*ones(1,Ny-2) (Nx-1):-1:1 ones(1,Ny-3)];
68	! j=[2ones(1,Nx) 3:Ny Nyones(1,Nx-1) (Ny-1):-1:3 ];
69
70	! AreaTag
71	! ice area (fractional: 0-1) at the tracer points
72	! Nx * Ny Real*8 array
73
74	! HeffTag
75	! effective ice thickness (mean thickness in the grid box in m)
76	! at the tracer points
77	! actual ice thickness for single-class ice is: Heff / Area
78	! Nx * Ny Real*8 array
79
80	! HsnowTag
81	! effective snow thickness (mean thickness in the grid box in m)
82	! at the tracer points
83	! actual snow thickness is: Hsnow / Area
84	! Nx * Ny Real*8 array
85
86	! HsaltTag
87	! effective sea ice salinity (g/m2) at the tracer points
88	! actual salinity in g/kg is: Hsalt / Heff / rhoIce
89	! where rhoice = 910 kg / m^3
90	! Nx * Ny Real*8 array
91
92	! UwindTag
93	! 10-m u-wind at the tracer points (m/s)
94	! >0 from West to East
95	! Nx * Ny Real*8 array
96
97	! VwindTag
98	! 10-m v-wind at the tracer points (m/s)
99	! >0 from South to North
100	! Nx * Ny Real*8 array
101
102	! LwDownTag
103	! downward longwave at the tracer points (W/m2)
104	! > 0 for increase in theta (ocean warming)
105	! Nx * Ny Real*8 array
106
107	! SwDownTag
108	! downward shortwave at the tracer points (W/m2)
109	! > 0 for increase in theta (ocean warming)
110	! Nx * Ny Real*8 array
111
112	! AtempTag
113	! 2-m air temperature at the tracer points (deg K)
114	! Nx * Ny Real*8 array
115
116	! AqhTag
117	! 2-m specific humidity at the tracer points (kg/kg)
118	! Nx * Ny Real*8 array
119
120	! PrecipTag
121	! precipitation at the tracer points (m/s)
122	! > 0 for decrease in salt (ocean salinity)
123	! Nx * Ny Real*8 array
124
125	! SstTag
126	! sea surface temperature at the tracer points (deg C)
127	! Nx * Ny Real*8 array
128
129	! SssTag
130	! sea surface salinity at the tracer points (g/kg)
131	! Nx * Ny Real*8 array
132
133	! UvelTag
134	! sea surface zonal velocity (m/s)
135	! at southwest c-grid locations
136	! Nx * Ny Real*8 array
137
138	! VvelTag
139	! sea surface meridional velocity (m/s)
140	! at southwest c-grid locations
141	! Nx * Ny Real*8 array
142
143	! UstressTag
144	! zonal stress under ice at southwest c-grid locations (N/m^2)
145	! Ustress is computed for ice-covered area only
146	! > 0 increases ocean surface velocity in the U direction
147	! Nx * Ny Real*8 array
148
149	! VstressTag
150	! meridional stress under ice at southwest c-grid locations (N/m^2)
151	! Vstress is computed for ice-covered area only
152	! > 0 increases ocean surface velocity in the V direction
153	! Nx * Ny Real*8 array
154
155	! SwResidTag
156	! residual shortwave under ice at tracer point (W/m2)
157	! SwResid is computed for ice-covered area only
158	! > 0 for increase in theta (ocean warming)
159	! Nx * Ny Real*8 array
160
161	! HeatFluxTag
162	! heat flux (less shortwave) at tracer point (W/m2)
163	! HeatFlux is computed for ice-covered area only
164	! > 0 for increase in theta (ocean warming)
165	! Nx * Ny Real*8 array
166
167	! WaterFluxTag
168	! freshwater flux at tracer point ( kg/m^2/s )
169	! WaterFlux is computed for ice-covered area only
170	! > 0 for decrease in ocean salinity
171	! Nx * Ny Real*8 array
172
173	! SaltFluxTag
174	! salt flux at tracer point (g/m^2/s)
175	! SaltFlux is computed for ice-covered area only
176	! > 0 for increase in ocean salinity
177	! Nx * Ny Real*8 array
178
179	INTEGER TimeIntervalTag
180	INTEGER OceanTimeTag
181	INTEGER IceTimeTag
182	INTEGER OceanGridsizeTag
183	INTEGER AreaBcTag
184	INTEGER HeffBcTag
185	INTEGER HsnowBcTag
186	INTEGER HsaltBcTag
187	INTEGER UiceBcTag
188	INTEGER ViceBcTag
189	INTEGER AreaTag
190	INTEGER HeffTag
191	INTEGER HsnowTag
192	INTEGER HsaltTag
193	INTEGER UwindTag
194	INTEGER VwindTag
195	INTEGER LwDownTag
196	INTEGER SwDownTag
197	INTEGER AtempTag
198	INTEGER AqhTag
199	INTEGER PrecipTag
200	INTEGER SstTag
201	INTEGER SssTag
202	INTEGER UvelTag
203	INTEGER VvelTag
204	INTEGER UstressTag
205	INTEGER VstressTag
206	INTEGER SwResidTag
207	INTEGER HeatFluxTag
208	INTEGER WaterFluxTag
209	INTEGER SaltFluxTag
210
211	PARAMETER ( TimeIntervalTag = 10001000 )
212	PARAMETER ( OceanTimeTag = 10001001 )
213	PARAMETER ( IceTimeTag = 10001002 )
214	PARAMETER ( OceanGridsizeTag = 10002001 )
215	PARAMETER ( AreaBcTag = 10003001 )
216	PARAMETER ( HeffBcTag = 10003002 )
217	PARAMETER ( HsnowBcTag = 10003003 )
218	PARAMETER ( HsaltBcTag = 10003004 )
219	PARAMETER ( UiceBcTag = 10003005 )
220	PARAMETER ( ViceBcTag = 10003006 )
221	PARAMETER ( AreaTag = 10004001 )
222	PARAMETER ( HeffTag = 10004002 )
223	PARAMETER ( HsnowTag = 10004003 )
224	PARAMETER ( HsaltTag = 10004004 )
225	PARAMETER ( UwindTag = 10005001 )
226	PARAMETER ( VwindTag = 10005002 )
227	PARAMETER ( LwDownTag = 10005101 )
228	PARAMETER ( SwDownTag = 10005102 )
229	PARAMETER ( AtempTag = 10005103 )
230	PARAMETER ( AqhTag = 10005104 )
231	PARAMETER ( PrecipTag = 10005105 )
232	PARAMETER ( SstTag = 10005201 )
233	PARAMETER ( SssTag = 10005202 )
234	PARAMETER ( UvelTag = 10005203 )
235	PARAMETER ( VvelTag = 10005204 )
236	PARAMETER ( UstressTag = 10006001 )
237	PARAMETER ( VstressTag = 10006002 )
238	PARAMETER ( SwResidTag = 20006101 )
239	PARAMETER ( HeatFluxTag = 20006102 )
240	PARAMETER ( WaterFluxTag = 20006103 )
241	PARAMETER ( SaltFluxTag = 20006104 )