/[MITgcm]/MITgcm_contrib/MPMice/beaufort/CPL_README.txt
ViewVC logotype

Diff of /MITgcm_contrib/MPMice/beaufort/CPL_README.txt

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

revision 1.1 by dimitri, Sun May 31 13:25:42 2009 UTC revision 1.8 by dimitri, Wed Jul 13 23:13:31 2016 UTC
# Line 3  Initialization (myTime .EQ. startTime in Line 3  Initialization (myTime .EQ. startTime in
3  ===> Ocean Sends/ Ice Receives (fields that are sent only once)  ===> Ocean Sends/ Ice Receives (fields that are sent only once)
4  deltatimestep                1           Real*8  TimeIntervalTag  deltatimestep                1           Real*8  TimeIntervalTag
5  grid dimensions (Nx,Ny)      2           Integer OceanGridsizeTag  grid dimensions (Nx,Ny)      2           Integer OceanGridsizeTag
6  (grid locations - later?)  
7    xC                           Nx*Ny       Real*8  xCtag
8    yC                           Nx*Ny       Real*8  yCtag
9    xG                           Nx*Ny       Real*8  xGtag
10    yG                           Nx*Ny       Real*8  yGtag
11    dxG                          Nx*Ny       Real*8  dxGtag
12    dyG                          Nx*Ny       Real*8  dxGtag
13    AngleCS                      Nx*Ny       Real*8  aCStag
14    AngleSN                      Nx*Ny       Real*8  aSNtag
15    hFacC                        Nx*Ny       Real*8  hFacCtag
16  ice area                     Nx*Ny       Real*8  AreaTag  ice area                     Nx*Ny       Real*8  AreaTag
17  ice thickness                Nx*Ny       Real*8  HeffTag  ice thickness                Nx*Ny       Real*8  HeffTag
18  ice salinity                 Nx*Ny       Real*8  HsaltTag  ice salinity                 Nx*Ny       Real*8  HsaltTag
19  snow thickness               Nx*Ny       Real*8  HsnowTag  snow thickness               Nx*Ny       Real*8  HsnowTag
20    
 ===> Ocean Sends/ Ice Receives  
 ocean model time             1           Real*8  OceanTimeTag  
 boundary ice area            2*(Nx+Ny)-4 Real*8  AreaBcTag  
 boundary ice thickness       2*(Nx+Ny)-4 Real*8  HeffBcTag  
 boundary ice salinity        2*(Nx+Ny)-4 Real*8  HsaltBcTag  
 boundary snow thickness      2*(Nx+Ny)-4 Real*8  HsnowBcTag  
 boundary u ice               2*(Nx+Ny)-6 Real*8  UiceBcTag  
 boundary v ice               2*(Nx+Ny)-6 Real*8  ViceBcTag  
 u-wind velocity              Nx*Ny       Real*8  UwindTag  
 v-wind velocity              Nx*Ny       Real*8  VwindTag  
 downward longwave radiation  Nx*Ny       Real*8  LwDownTag  
 downward shortwave radiation Nx*Ny       Real*8  SwDownTag  
 air temperature              Nx*Ny       Real*8  AtempTag  
 humidity                     Nx*Ny       Real*8  AqhTag  
 precipitation                Nx*Ny       Real*8  PrecipTag  
 ocean surface temperature    Nx*Ny       Real*8  SstTag  
 surface u current            Nx*Ny       Real*8  UvelTag  
 surface v current            Nx*Ny       Real*8  VvelTag  
   
 ===> Ice Sends/Ocean Receives  
 ice model time               1           Real*8  IceTimeTag  
 ice area                     Nx*Ny       Real*8  AreaTag (initial ice area)  
 ice thickness                Nx*Ny       Real*8  HeffTag (initial ice thickness)  
 ice salinity                 Nx*Ny       Real*8  HsaltTag (initial ice salinity)  
 snow thickness               Nx*Ny       Real*8  HsnowTag (initial snow thickness)  
 u surface stress             Nx*Ny       Real*8  UstressTag (surface v current)  
 v surface stress             Nx*Ny       Real*8  VstressTag (array of 2)  
 residual shortwave           Nx*Ny       Real*8  SwResidTag (array of 1)  
 heat flux                    Nx*Ny       Real*8  HeatFluxTag (array of 1)  
 freshwater flux              Nx*Ny       Real*8  WaterFluxTag (array of 1)  
 salt flux                    Nx*Ny       Real*8  SaltFluxTag (array of 1)  
   
21  =====================================  =====================================
22    
23  Each deltatimestep:  Each deltatimestep:
# Line 61  air temperature              Nx*Ny Line 38  air temperature              Nx*Ny
38  humidity                     Nx*Ny       Real*8  AqhTag  humidity                     Nx*Ny       Real*8  AqhTag
39  precipitation                Nx*Ny       Real*8  PrecipTag  precipitation                Nx*Ny       Real*8  PrecipTag
40  ocean surface temperature    Nx*Ny       Real*8  SstTag  ocean surface temperature    Nx*Ny       Real*8  SstTag
41    ocean surface salinity       Nx*Ny       Real*8  SssTag
42  surface u current            Nx*Ny       Real*8  UvelTag  surface u current            Nx*Ny       Real*8  UvelTag
43  surface v current            Nx*Ny       Real*8  VvelTag  surface v current            Nx*Ny       Real*8  VvelTag
44    
# Line 70  ice area                     Nx*Ny Line 48  ice area                     Nx*Ny
48  ice thickness                Nx*Ny       Real*8  HeffTag  ice thickness                Nx*Ny       Real*8  HeffTag
49  ice salinity                 Nx*Ny       Real*8  HsaltTag  ice salinity                 Nx*Ny       Real*8  HsaltTag
50  snow thickness               Nx*Ny       Real*8  HsnowTag  snow thickness               Nx*Ny       Real*8  HsnowTag
51    u ice velocity               Nx*Ny       Real*8  UiceTag
52    v ice velocity               Nx*Ny       Real*8  ViceTag
53  u surface stress             Nx*Ny       Real*8  UstressTag  u surface stress             Nx*Ny       Real*8  UstressTag
54  v surface stress             Nx*Ny       Real*8  VstressTag  v surface stress             Nx*Ny       Real*8  VstressTag
55  residual shortwave           Nx*Ny       Real*8  SwResidTag  residual shortwave           Nx*Ny       Real*8  SwResidTag
# Line 89  DETAILS: Line 69  DETAILS:
69    
70  2. ocean sends  2. ocean sends
71    
72  deltatimestep  **** initial time step only *******************************
73    
74  ocean model time "myTime" (s)  deltatimestep                1           Real*8  TimeIntervalTag
75    
76    grid dimensions (Nx,Ny)      2           Integer OceanGridsizeTag
77    
78    Some grid information - see MITgcm/model/inc/GRID.h for details:
79    xC  :: longitude East of center of grid cell
80    yC  ::latitude North of center of grid cell
81    xG  :: longitude East of SouthWest corner
82    yG  :: latitude North of SouthWest corner
83    dxG :: distance in m between SouthWest and SouthEast corner
84    dyG :: distance in m between SouthWest and NorthEast corner
85    
86    aCS :: cosine(alpha) relative to geographic direction at grid cell center
87    aSN :: sine(alpha) relative to geographic direction at grid cell center
88    
89           alpha = angle of model uVel direction vs geographical East
90                 = angle of model vVel direction vs geographical North
91           (AngleCS*uVelc - AngleSN*vVelc, AngleSN*uVelc + AngleCN*vVelc)
92           rotates model velocity to geographical coordinates, where
93           (uVelc,vVelc) is model velocity vector at center of grid cell
94    
95  Then on tracer grid:  Then on tracer grid:
96    
97  grid information + landmask  TBD  hFacC :: landmask of center of grid cell, 0 is land, >0 is ocean
98    
99  initial ice area (fractional: 0-1) at the tracer points  initial ice area (fractional: 0-1) at the tracer points
100  C     AREA  - fractional ice-covered area in m^2/m^2  C     AREA  - fractional ice-covered area in m^2/m^2
# Line 110  C             at center of grid, i.e., t Line 109  C             at center of grid, i.e., t
109  C             note: for non-zero AREA, actual ice  C             note: for non-zero AREA, actual ice
110  C                thickness is HEFF / AREA  C                thickness is HEFF / AREA
111    
112    initial ice salinity at the tracer points
113    C     HSALT - effective sea ice salinity in g/m^2
114    C             at center of grid, i.e., tracer point
115    
116  initial snow thickness (mean thickness, m) at the tracer points  initial snow thickness (mean thickness, m) at the tracer points
117  C     HSNOW - effective snow thickness in m  C     HSNOW - effective snow thickness in m
118  C             at center of grid, i.e., tracer point  C             at center of grid, i.e., tracer point
119  C             note: for non-zero AREA, actual snow  C             note: for non-zero AREA, actual snow
120  C                thickness is HSNOW / AREA  C                thickness is HSNOW / AREA
121    
122  initial ice salinity at the tracer points  **** every time step **************************************
123  C     HSALT - effective sea ice salinity in g/m^2  
124  C             at center of grid, i.e., tracer point  ocean model time "myTime" (s)
125    
126    boundary ice area            2*(Nx+Ny)-4 Real*8  AreaBcTag
127    
128    boundary ice thickness       2*(Nx+Ny)-4 Real*8  HeffBcTag
129    
130    boundary ice salinity        2*(Nx+Ny)-4 Real*8  HsaltBcTag
131    
132    boundary snow thickness      2*(Nx+Ny)-4 Real*8  HsnowBcTag
133    
134  open boundary u-ice velocity  open boundary u-ice velocity
135  specified at the inside edge of the outer grid cells  specified at the inside edge of the outer grid cells
136  southwest C-grid locations for normal component and at the  southwest C-grid locations for normal component and at the
137  southwest C-grid locations for tangential  southwest C-grid locations for tangential
138  C     UICE  - zonal ice velocity in m/s at South-West B-grid  C     UICE  - zonal ice velocity in m/s at South-West C-grid U point
 C             (or C-grid #ifdef SEAICE_CGRID) U point  
139  C             >0 from West to East  C             >0 from West to East
140    
141  open boundary v-ice velocity  open boundary v-ice velocity
142  specified at the inside edge of the outer grid cells  specified at the inside edge of the outer grid cells
143  southwest C-grid locations for normal component and at the  southwest C-grid locations for normal component and at the
144  southwest C-grid locations for tangential  southwest C-grid locations for tangential
145  C     VICE  - meridional ice velocity in m/s at South-West B-grid  C     VICE  - meridional ice velocity in m/s at South-West C-grid V point
 C             (or C-grid #ifdef SEAICE_CGRID) V point  
146  C             >0 from South to North  C             >0 from South to North
 C             note: the South-West B-grid U and V points are on  
 C                the lower, left-hand corner of each grid cell  
147    
148  10-m u-wind at the tracer points  10-m u-wind at the tracer points
149  c     uwind     :: Surface (10-m) zonal wind velocity in m/s  c     uwind     :: Surface (10-m) zonal wind velocity in m/s
# Line 196  at southwest c-grid locations Line 203  at southwest c-grid locations
203  3. ice receives 2., computes initial ice conditions and sends on tracer grid  3. ice receives 2., computes initial ice conditions and sends on tracer grid
204    
205  ice model time (s)  ice model time (s)
206    
207  ice area at tracer point  ice area at tracer point
208    C     AREA  :: fractional ice-covered area in m^2/m^2
209    C              at center of grid, i.e., tracer point
210    C              0 is no cover, 1 is 100% cover
211    
212  ice thickness at tracer point  ice thickness at tracer point
213    C     HEFF  :: effective ice thickness in m
214    C              at center of grid, i.e., tracer point
215    C              note: for non-zero AREA,
216    C                    actual ice thickness is HEFF / AREA
217    
218  ice salinity at tracer point  ice salinity at tracer point
219    C     HSALT :: effective sea ice salinity in g/m^2
220    C              at center of grid, i.e., tracer point
221    
222  snow thickness at tracer point  snow thickness at tracer point
223  u-stress under ice at southwest b-grid locations  C     HSNOW :: effective snow thickness in m
224  v-stress under ice at southwest b-grid locations  C              at center of grid, i.e., tracer point
225    C              note: for non-zero AREA, actual snow thickness is HSNOW / AREA
226    
227    u ice velocity at southwest c-grid locations
228    C     UICE  - zonal ice velocity in m/s at South-West C-grid U point
229    C             >0 from West to East
230    
231    v ice velocity at southwest c-grid locations
232    C     VICE  - meridional ice velocity in m/s at South-West C-grid V point
233    C             >0 from South to North
234    
235    u-stress under ice at southwest c-grid locations
236    C     fu    :: Zonal surface wind stress in N/m^2
237    C              > 0 for increase in uVel, which is west to
238    C                  east for cartesian and spherical polar grids
239    C              Typical range: -0.5 < fu < 0.5
240    C              Southwest C-grid U point
241    
242    v-stress under ice at southwest c-grid locations
243    C     fv    :: Meridional surface wind stress in N/m^2
244    C              > 0 for increase in vVel, which is south to
245    C                  north for cartesian and spherical polar grids
246    C              Typical range: -0.5 < fv < 0.5
247    C              Southwest C-grid V point
248    
249  residual shortwave under ice at tracer point  residual shortwave under ice at tracer point
250    C     Qsw   :: Net upward shortwave radiation in W/m^2
251    C              Qsw = - ( downward - ice and snow absorption - reflected )
252    C              > 0 for decrease in theta (ocean cooling)
253    C              Typical range: -350 < Qsw < 0
254    C              Southwest C-grid tracer point
255    
256  heat flux (less shortwave) at tracer point  heat flux (less shortwave) at tracer point
257    C     HeatFlux
258    C           :: Upward surface heat flux (excluding shortwave) in W/m^2
259    C              HeatFlux = latent + sensible + net longwave
260    C              > 0 for decrease in theta (ocean cooling)
261    C              Typical range: -250 < HeatFlux < 950
262    C              Southwest C-grid tracer point
263    
264  freshwater flux at tracer point  freshwater flux at tracer point
265    C     EmPmR :: Net upward freshwater flux in kg/m2/s
266    C              EmPmR = Evaporation - precipitation - runoff
267    C              > 0 for increase in salt (ocean salinity)
268    C              Typical range: -1e-4 < EmPmR < 1e-4
269    C              Southwest C-grid tracer point
270    C           NOTE: for backward compatibility EmPmRfile is specified in
271    C                 m/s when using external_fields_load.F.  It is converted
272    C                 to kg/m2/s by multiplying by rhoConstFresh.
273    
274  salt flux at tracer point  salt flux at tracer point
275    C  saltFlux :: Net upward salt flux in psu.kg/m^2/s
276    C              flux of Salt taken out of the ocean per time unit (second).
277    C              Note: a) only used when salty sea-ice forms or melts.
278    C                    b) units: when salinity (unit= psu) is expressed
279    C                       in g/kg, saltFlux unit becomes g/m^2/s.
280    C              > 0 for decrease in SSS.
281    C              Southwest C-grid tracer point
282    
283  ===========================================================  ===========================================================
284    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22