| 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 |  | 
| 16 |  | 
| 17 | 1) use create_bc_grid.m (called by test_bc_grid.m) to generate bc grid files. | 
| 18 | create_bc_grid.m will generate files with names like | 
| 19 | bdry_grid_[snwe].mat for south, north, west, and east boundary. | 
| 20 | This script needs a directory path for MITgcm domain related files, | 
| 21 | specifically  YC.data, XC.data, hFacC.data, AngleCS.data, AngleSN.data. | 
| 22 |  | 
| 23 | 2) use test_tide_bc.m to generate tidal bc files which calls | 
| 24 | create_tide_bc_second.m | 
| 25 | or dummy tidal bc files (zero amplitude and phase, tide_bc_empty.m) | 
| 26 | for land boundary if there is such a need. | 
| 27 |  | 
| 28 | create_tide_bc_second.m will generate tidal boundary files with names | 
| 29 | like OB[SNWE]{am ph}.seaice_obcs. | 
| 30 | The tidal current and phase will be used in  obcs_add_tides.F. | 
| 31 |  | 
| 32 | 3) Modify MITgcm configuration | 
| 33 | 3.1) ALLOW_OBCS_TIDES in OBCS_OPTIONS.h | 
| 34 | 3.2) OBCS_PARAMS.h  (tidalComponents = ??, number of tidal constituent to be included.) | 
| 35 | 3.3) data.obcs as following | 
| 36 | The number of tidal constituents to be included and their period and | 
| 37 | order should be consistent with their order in tidal current amplitude | 
| 38 | and phase files,  OB?amFile and OB?phFile. | 
| 39 | useOBCStides = .TRUE. | 
| 40 | tidalperiod = ... | 
| 41 | OBSamFile =  ... | 
| 42 | OBNamFile =  ... | 
| 43 | OBWamFile = ... | 
| 44 | OBEamFile = ... | 
| 45 | OBSphFile = ... | 
| 46 | OBNphFile = ... | 
| 47 | OBWphFile = ... | 
| 48 | OBEphFile = ... | 
| 49 |  | 
| 50 |  | 
| 51 |  | 
| 52 | CAT Reference | 
| 53 | Padman, L., H. A. Fricker, R. Coleman, S. Howard, and S. Erofeeva, 2002: A new tidal model for the | 
| 54 | Antarctic ice shelves and seas, Ann. Glaciol., 34, 247-254. | 
| 55 |  | 
| 56 |  | 
| 57 | XC WANG 25 Mar  2013 | 
| 58 | XC WANG 07 Feb. 2013 | 
| 59 | XC WANG 20 Dec. 2012 | 
| 60 | xiaochun.wang@jpl.nasa.gov |