| 1 |
The scripts in this directory are for adding tidal boundary condition |
| 2 |
to a regional configuration of MITgcm. |
| 3 |
|
| 4 |
The scripts are tested for Pine Island Model |
| 5 |
and also for verification/seaice_obcs using specified tidal boundary file. |
| 6 |
|
| 7 |
Steps to add tide boundary conditions to regional MITgcm |
| 8 |
0) Download the Circum-Antarctic Tide (CATS) Solution (Padman et al. 2002) |
| 9 |
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 |
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 |
|
| 17 |
|
| 18 |
1) use create_bc_grid.m (called by test_bc_grid.m) to generate bc grid files. |
| 19 |
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 |
specifically YC.data, XC.data, hFacC.data, AngleCS.data, AngleSN.data. |
| 23 |
|
| 24 |
2) use test_tide_bc.m to generate tidal bc files which calls |
| 25 |
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 |
|
| 29 |
create_tide_bc_second.m will generate tidal boundary files with names |
| 30 |
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 |
tidalperiod = ... |
| 42 |
OBSamFile = ... |
| 43 |
OBNamFile = ... |
| 44 |
OBWamFile = ... |
| 45 |
OBEamFile = ... |
| 46 |
OBSphFile = ... |
| 47 |
OBNphFile = ... |
| 48 |
OBWphFile = ... |
| 49 |
OBEphFile = ... |
| 50 |
|
| 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 |
XC WANG 02 Apr 2013 |
| 59 |
XC WANG 25 Mar 2013 |
| 60 |
XC WANG 07 Feb. 2013 |
| 61 |
XC WANG 20 Dec. 2012 |
| 62 |
xiaochun.wang@jpl.nasa.gov |