1 |
dimitri |
1.1 |
The scripts in this directory are for adding tidal boundary condition |
2 |
|
|
to a regional configuration of MITgcm. |
3 |
|
|
|
4 |
dimitri |
1.2 |
The scripts are tested for Pine Island Model |
5 |
|
|
and also for verification/seaice_obcs using specified tidal boundary file. |
6 |
dimitri |
1.1 |
|
7 |
|
|
Steps to add tide boundary conditions to regional MITgcm |
8 |
dimitri |
1.2 |
0) Download the Circum-Antarctic Tide (CATS) Solution (Padman et al. 2002) |
9 |
dimitri |
1.3 |
from Earth and Space Research, http://www.esr.org/~padman/cats2008b/ |
10 |
|
|
Modify CATS directory in test_tide_bc.m |
11 |
|
|
Download the TMD2.03 software for extracting the tidal solution from |
12 |
|
|
http://www.esr.org. |
13 |
|
|
Modifications are needed in the original /TMD2.03/tmd_extract_HC.m in order |
14 |
xiaochun.wang |
1.4 |
to have the correct path of tidal solution as shown in the sample_tmd_extract_HC.m. |
15 |
|
|
Use CATS2008b not CAT2008a since CATS200a is not on north-south east-west coordinate |
16 |
dimitri |
1.1 |
|
17 |
dimitri |
1.3 |
|
18 |
|
|
1) use create_bc_grid.m (called by test_bc_grid.m) to generate bc grid files. |
19 |
dimitri |
1.1 |
create_bc_grid.m will generate files with names like |
20 |
|
|
bdry_grid_[snwe].mat for south, north, west, and east boundary. |
21 |
|
|
This script needs a directory path for MITgcm domain related files, |
22 |
dimitri |
1.3 |
specifically YC.data, XC.data, hFacC.data, AngleCS.data, AngleSN.data. |
23 |
dimitri |
1.1 |
|
24 |
dimitri |
1.2 |
2) use test_tide_bc.m to generate tidal bc files which calls |
25 |
dimitri |
1.3 |
create_tide_bc_second.m |
26 |
|
|
or dummy tidal bc files (zero amplitude and phase, tide_bc_empty.m) |
27 |
|
|
for land boundary if there is such a need. |
28 |
dimitri |
1.1 |
|
29 |
dimitri |
1.3 |
create_tide_bc_second.m will generate tidal boundary files with names |
30 |
dimitri |
1.2 |
like OB[SNWE]{am ph}.seaice_obcs. |
31 |
|
|
The tidal current and phase will be used in obcs_add_tides.F. |
32 |
|
|
|
33 |
|
|
3) Modify MITgcm configuration |
34 |
|
|
3.1) ALLOW_OBCS_TIDES in OBCS_OPTIONS.h |
35 |
|
|
3.2) OBCS_PARAMS.h (tidalComponents = ??, number of tidal constituent to be included.) |
36 |
|
|
3.3) data.obcs as following |
37 |
|
|
The number of tidal constituents to be included and their period and |
38 |
|
|
order should be consistent with their order in tidal current amplitude |
39 |
|
|
and phase files, OB?amFile and OB?phFile. |
40 |
|
|
useOBCStides = .TRUE. |
41 |
dimitri |
1.1 |
tidalperiod = ... |
42 |
|
|
OBSamFile = ... |
43 |
|
|
OBNamFile = ... |
44 |
dimitri |
1.2 |
OBWamFile = ... |
45 |
|
|
OBEamFile = ... |
46 |
dimitri |
1.1 |
OBSphFile = ... |
47 |
|
|
OBNphFile = ... |
48 |
|
|
OBWphFile = ... |
49 |
dimitri |
1.2 |
OBEphFile = ... |
50 |
dimitri |
1.1 |
|
51 |
|
|
|
52 |
|
|
|
53 |
|
|
CAT Reference |
54 |
|
|
Padman, L., H. A. Fricker, R. Coleman, S. Howard, and S. Erofeeva, 2002: A new tidal model for the |
55 |
|
|
Antarctic ice shelves and seas, Ann. Glaciol., 34, 247-254. |
56 |
|
|
|
57 |
|
|
|
58 |
xiaochun.wang |
1.4 |
XC WANG 02 Apr 2013 |
59 |
dimitri |
1.2 |
XC WANG 25 Mar 2013 |
60 |
dimitri |
1.1 |
XC WANG 07 Feb. 2013 |
61 |
|
|
XC WANG 20 Dec. 2012 |
62 |
|
|
xiaochun.wang@jpl.nasa.gov |