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

!     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     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 ( 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	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			! 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	dimitri	1.2	! open boundary effective ice thickness (m) at the tracer points
30			! actual ice thickness for single-class ice is: Heff / Area
31	dimitri	1.1	! 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	dimitri	1.2	! open boundary effective snow thickness (m) at the tracer points
37			! actual snow thickness is: Hsnow / Area
38	dimitri	1.1	! 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	dimitri	1.2	! 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	dimitri	1.1	! 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	dimitri	1.2	! actual ice thickness for single-class ice is: Heff / Area
78	dimitri	1.1	! 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	dimitri	1.2	! actual snow thickness is: Hsnow / Area
84	dimitri	1.1	! Nx * Ny Real*8 array
85
86			! HsaltTag
87	dimitri	1.2	! 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	dimitri	1.1	! Nx * Ny Real*8 array
91
92	dimitri	1.4	! UiceTag
93			! zonal ice velocity in m/s at South-West C-grid U point
94			! >0 from West to East
95			! Nx * Ny Real*8 array
96
97			! ViceTag
98			! meridional ice velocity in m/s at South-West C-grid V point
99			! >0 from South to North
100			! Nx * Ny Real*8 array
101
102	dimitri	1.1	! UwindTag
103			! 10-m u-wind at the tracer points (m/s)
104			! >0 from West to East
105			! Nx * Ny Real*8 array
106
107			! VwindTag
108			! 10-m v-wind at the tracer points (m/s)
109			! >0 from South to North
110			! Nx * Ny Real*8 array
111
112			! LwDownTag
113			! downward longwave at the tracer points (W/m2)
114			! > 0 for increase in theta (ocean warming)
115			! Nx * Ny Real*8 array
116
117			! SwDownTag
118			! downward shortwave at the tracer points (W/m2)
119			! > 0 for increase in theta (ocean warming)
120			! Nx * Ny Real*8 array
121
122			! AtempTag
123			! 2-m air temperature at the tracer points (deg K)
124			! Nx * Ny Real*8 array
125
126			! AqhTag
127			! 2-m specific humidity at the tracer points (kg/kg)
128			! Nx * Ny Real*8 array
129
130			! PrecipTag
131			! precipitation at the tracer points (m/s)
132			! > 0 for decrease in salt (ocean salinity)
133			! Nx * Ny Real*8 array
134
135			! SstTag
136			! sea surface temperature at the tracer points (deg C)
137			! Nx * Ny Real*8 array
138
139	dimitri	1.2	! SssTag
140			! sea surface salinity at the tracer points (g/kg)
141			! Nx * Ny Real*8 array
142
143	dimitri	1.1	! UvelTag
144			! sea surface zonal velocity (m/s)
145			! at southwest c-grid locations
146			! Nx * Ny Real*8 array
147
148			! VvelTag
149			! sea surface meridional velocity (m/s)
150			! at southwest c-grid locations
151			! Nx * Ny Real*8 array
152
153			! UstressTag
154			! zonal stress under ice at southwest c-grid locations (N/m^2)
155	dimitri	1.2	! Ustress is computed for ice-covered area only
156			! > 0 increases ocean surface velocity in the U direction
157	dimitri	1.1	! Nx * Ny Real*8 array
158
159			! VstressTag
160			! meridional stress under ice at southwest c-grid locations (N/m^2)
161	dimitri	1.2	! Vstress is computed for ice-covered area only
162			! > 0 increases ocean surface velocity in the V direction
163	dimitri	1.1	! Nx * Ny Real*8 array
164
165			! SwResidTag
166			! residual shortwave under ice at tracer point (W/m2)
167	dimitri	1.2	! SwResid is computed for ice-covered area only
168	dimitri	1.1	! > 0 for increase in theta (ocean warming)
169			! Nx * Ny Real*8 array
170
171			! HeatFluxTag
172			! heat flux (less shortwave) at tracer point (W/m2)
173	dimitri	1.2	! HeatFlux is computed for ice-covered area only
174	dimitri	1.1	! > 0 for increase in theta (ocean warming)
175			! Nx * Ny Real*8 array
176
177			! WaterFluxTag
178	dimitri	1.3	! freshwater flux at tracer point ( kg/m^2/s )
179	dimitri	1.2	! WaterFlux is computed for ice-covered area only
180	dimitri	1.1	! > 0 for decrease in ocean salinity
181			! Nx * Ny Real*8 array
182
183			! SaltFluxTag
184			! salt flux at tracer point (g/m^2/s)
185	dimitri	1.2	! SaltFlux is computed for ice-covered area only
186	dimitri	1.1	! > 0 for increase in ocean salinity
187			! Nx * Ny Real*8 array
188
189			INTEGER TimeIntervalTag
190			INTEGER OceanTimeTag
191			INTEGER IceTimeTag
192			INTEGER OceanGridsizeTag
193			INTEGER AreaBcTag
194			INTEGER HeffBcTag
195			INTEGER HsnowBcTag
196			INTEGER HsaltBcTag
197			INTEGER UiceBcTag
198			INTEGER ViceBcTag
199			INTEGER AreaTag
200			INTEGER HeffTag
201			INTEGER HsnowTag
202			INTEGER HsaltTag
203	dimitri	1.4	INTEGER UiceTag
204			INTEGER ViceTag
205	dimitri	1.1	INTEGER UwindTag
206			INTEGER VwindTag
207			INTEGER LwDownTag
208			INTEGER SwDownTag
209			INTEGER AtempTag
210			INTEGER AqhTag
211			INTEGER PrecipTag
212			INTEGER SstTag
213	dimitri	1.2	INTEGER SssTag
214	dimitri	1.1	INTEGER UvelTag
215			INTEGER VvelTag
216			INTEGER UstressTag
217			INTEGER VstressTag
218			INTEGER SwResidTag
219			INTEGER HeatFluxTag
220			INTEGER WaterFluxTag
221			INTEGER SaltFluxTag
222
223			PARAMETER ( TimeIntervalTag = 10001000 )
224			PARAMETER ( OceanTimeTag = 10001001 )
225			PARAMETER ( IceTimeTag = 10001002 )
226			PARAMETER ( OceanGridsizeTag = 10002001 )
227			PARAMETER ( AreaBcTag = 10003001 )
228			PARAMETER ( HeffBcTag = 10003002 )
229			PARAMETER ( HsnowBcTag = 10003003 )
230			PARAMETER ( HsaltBcTag = 10003004 )
231			PARAMETER ( UiceBcTag = 10003005 )
232			PARAMETER ( ViceBcTag = 10003006 )
233			PARAMETER ( AreaTag = 10004001 )
234			PARAMETER ( HeffTag = 10004002 )
235			PARAMETER ( HsnowTag = 10004003 )
236			PARAMETER ( HsaltTag = 10004004 )
237	dimitri	1.4	PARAMETER ( UiceTag = 10004005 )
238			PARAMETER ( ViceTag = 10004006 )
239	dimitri	1.1	PARAMETER ( UwindTag = 10005001 )
240			PARAMETER ( VwindTag = 10005002 )
241			PARAMETER ( LwDownTag = 10005101 )
242			PARAMETER ( SwDownTag = 10005102 )
243			PARAMETER ( AtempTag = 10005103 )
244			PARAMETER ( AqhTag = 10005104 )
245			PARAMETER ( PrecipTag = 10005105 )
246			PARAMETER ( SstTag = 10005201 )
247	dimitri	1.2	PARAMETER ( SssTag = 10005202 )
248			PARAMETER ( UvelTag = 10005203 )
249			PARAMETER ( VvelTag = 10005204 )
250	dimitri	1.1	PARAMETER ( UstressTag = 10006001 )
251			PARAMETER ( VstressTag = 10006002 )
252			PARAMETER ( SwResidTag = 20006101 )
253			PARAMETER ( HeatFluxTag = 20006102 )
254			PARAMETER ( WaterFluxTag = 20006103 )
255			PARAMETER ( SaltFluxTag = 20006104 )