| 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 |