| 1 |
MITgcm, 36km resolution, Arctic Ocean regional simulation ‘How to’ guide |
| 2 |
|
| 3 |
Dr Jonny Williams, School of Geographical Sciences, Bristol University, UK |
| 4 |
jonny.williams@bristol.ac.uk |
| 5 |
|
| 6 |
1. cd to your MITgcm central directory and check out the code and (some of the) required |
| 7 |
input files |
| 8 |
a. cvs co -d cs_36km_tutorial MITgcm_contrib/arctic/cs_36km/ |
| 9 |
|
| 10 |
2. cd to this new directory and make a ‘build’ directory |
| 11 |
|
| 12 |
3. Personally I recommend outputting NetCDF files. You will need to do the following if you |
| 13 |
want to do this |
| 14 |
a. add a line with “mnc” on it to cs_36km_tutorial/code/packages.conf |
| 15 |
b. add a line with “useMNC=.TRUE.,” to cs_36km_tutorial/input/data.pkg |
| 16 |
c. add a “data.mnc” file which you can copy from the input directory of any of the |
| 17 |
verification experiments that use NetCDF output |
| 18 |
|
| 19 |
4. Compile in the usual way (system dependent) |
| 20 |
|
| 21 |
5. The following files are all required in the run directory |
| 22 |
a. the executable you get when compiling “mitgcmuv” |
| 23 |
b. ftp://ecco2.jpl.nasa.gov/data2/data/runoff/runoff-360x180x12.bin |
| 24 |
c. ftp://ecco2.jpl.nasa.gov/data2/data/atmos/jra25/*1979* |
| 25 |
i. note that the following files may need their “.txt” suffixes removed |
| 26 |
1. jra25_v10m_1979 |
| 27 |
2. jra25_u10m_1979 |
| 28 |
3. jra25_spfh2m_1979 |
| 29 |
4. jra25_temp2m_1979 |
| 30 |
d. Also, the run requires the following files which are not included in the |
| 31 |
ftp://ecco2.jpl.nasa.gov/data2/data/atmos/jra25/ directory (I have simply created |
| 32 |
cymbolic links to the 1979 equivalents). |
| 33 |
i. jra25_dlw_1978 |
| 34 |
ii. jra25_dsw_1978 |
| 35 |
iii. jra25_rain_1978 |
| 36 |
e. ftp://ecco2.jpl.nasa.gov/data1/arctic/run_template2/* |
| 37 |
|
| 38 |
6. There are some duplicate file names in |
| 39 |
ftp://ecco2.jpl.nasa.gov/data1/arctic/run_template2/ and cs_36km_tutorial/input and it is |
| 40 |
crucial that those in the latter are used. I have done this by copying them to my run |
| 41 |
directory (platform dependent) second to automatically overwrite the other ones. |
| 42 |
|
| 43 |
7. The model also requires several files which need to be renamed and so I have created |
| 44 |
symbolic links to them. |
| 45 |
a. ln -s OBWv_arctic_210x192.bin OBWv_arctic_210x192_19792014m.bin |
| 46 |
b. ln -s OBWu_arctic_210x192.bin OBWu_arctic_210x192_19792014m.bin |
| 47 |
c. ln -s OBWt_arctic_210x192.stable OBWt_arctic_210x192_19792014m.stable |
| 48 |
d. ln -s OBWt_arctic_210x192.bin OBWt_arctic_210x192_19792014m.bin |
| 49 |
e. ln -s OBWs_arctic_210x192.stable OBWs_arctic_210x192_19792014m.stable |
| 50 |
f. ln -s OBWs_arctic_210x192.bin OBWs_arctic_210x192_19792014m.bin |
| 51 |
g. ln -s OBNv_arctic_210x192.bin OBNv_arctic_210x192_19792014m.bin |
| 52 |
h. ln -s OBNu_arctic_210x192.bin OBNu_arctic_210x192_19792014m.bin |
| 53 |
i. ln -s OBNt_arctic_210x192.stable OBNt_arctic_210x192_19792014m.stable |
| 54 |
j. ln -s OBNt_arctic_210x192.bin OBNt_arctic_210x192_19792014m.bin |
| 55 |
k. ln -s OBNs_arctic_210x192.stable OBNs_arctic_210x192_19792014m.stable |
| 56 |
l. ln -s OBNs_arctic_210x192.bin OBNs_arctic_210x192_19792014m.bin |
| 57 |
m. ln -s OBEv_arctic_210x192.bin OBEv_arctic_210x192_19792014m.bin |
| 58 |
n. ln -s OBEu_arctic_210x192.bin OBEu_arctic_210x192_19792014m.bin |
| 59 |
o. ln -s OBEt_arctic_210x192.stable OBEt_arctic_210x192_19792014m.stable |
| 60 |
p. ln -s OBEt_arctic_210x192.bin OBEt_arctic_210x192_19792014m.bin |
| 61 |
q. ln -s OBEs_arctic_210x192.stable OBEs_arctic_210x192_19792014m.stable |
| 62 |
r. ln -s OBEs_arctic_210x192.bin OBEs_arctic_210x192_19792014m.bin |
| 63 |
|
| 64 |
8. I suggest changing the following in the “data” file just to get you going |
| 65 |
a. endtime=86400., |
| 66 |
i. i.e. 1 day in seconds |
| 67 |
b. dumpFreq = 2400., |
| 68 |
i. i.e. outputting data on every timestep, i.e. every 40 minutes in |
| 69 |
seconds |
| 70 |
|
| 71 |
9. Submit the job (platform dependent) |
| 72 |
|
| 73 |
10. Note that the NetCDF output from this simulation will consist of one file for each processor |
| 74 |
(80 by default here, nPx × nPy in the cs_36km_tutorial/code/SIZE.h file). |
| 75 |
|
| 76 |
Acknowledgments |
| 77 |
Getting this simulation to run would not have been possible without Fanny Monteiro in the author’s |
| 78 |
home institution, members of the MITgcm-support online community and in particular, Dimitris |
| 79 |
Menemenlis, Renske Gelderloos and Patrick Heimbach. |
| 80 |
|
| 81 |
======================= |
| 82 |
|
| 83 |
Directory contents for |
| 84 |
ftp://ecco2.jpl.nasa.gov/data1/arctic/run_template2/ |
| 85 |
|
| 86 |
Matlab script "mk_run_template2.m" was used to generate many of the model |
| 87 |
input files described below. |
| 88 |
|
| 89 |
Grid information files needed when OLD_GRID_IO option is defined in |
| 90 |
CPP_OPTIONS.h. See GRID.h for a digram of various distances. |
| 91 |
LONC.bin longitude east of cell center |
| 92 |
LATC.bin latitude north of cell center |
| 93 |
LONG.bin longitude east of southwest corner of cell |
| 94 |
LATG.bin latitude north of southwest corner of cell |
| 95 |
DYF.bin meridional distance in m between V-points |
| 96 |
DXF.bin zonal distance in m between U-points |
| 97 |
DYU.bin meridional distance in m between U-points |
| 98 |
DXV.bin zonal distance in m between V-points |
| 99 |
DYC.bin meridional distance in m between tracer points |
| 100 |
DXC.bin zonal distance in m between tracer points |
| 101 |
DYG.bin meridional distance in m between cell corners |
| 102 |
DXG.bin zonal distance in m between cell corners |
| 103 |
RAZ.bin vertical face area in m^2 for vorticity points |
| 104 |
RAW.bin vertical face area in m^2 for u cells |
| 105 |
RAS.bin vertical face area in m^2 for v cells |
| 106 |
RA.bin vertical face area in m^2 for tracer cells |
| 107 |
|
| 108 |
Model bathymetry, initial, and surface boundary condition files. |
| 109 |
ETOPO2_210x192_arctic model bathymetry (m) |
| 110 |
WGHC_S_210x192x50_arctic initial salinity (g/kg) |
| 111 |
WGHC_T_210x192x50_arctic initial potential temperature (deg C) |
| 112 |
AREA_210x192_arctic.cube81 initial ice concentration (fractional >=0, <=1) |
| 113 |
HEFF_210x192_arctic.cube81 initial effective sea ice thickness (m) |
| 114 |
HSALT_210x192_arctic.cube81 initial effective sea ice salinity (g/m^2) |
| 115 |
HSNOW_210x192_arctic.cube81 initial effective snow thickness (m) |
| 116 |
|
| 117 |
Open boundary condition files. |
| 118 |
OBNs_arctic_210x192.stable North open boundary conditions, salinity (g/kg) |
| 119 |
OBNt_arctic_210x192.stable North open boundary conditions, temperature (deg C) |
| 120 |
OBNu_arctic_210x192.bin North open boundary conditions, U-velocity (m/s) |
| 121 |
OBNv_arctic_210x192.bin North open boundary conditions, V-velocity (m/s) |
| 122 |
OBEs_arctic_210x192.stable South open boundary conditions, salinity (g/kg) |
| 123 |
OBEt_arctic_210x192.stable South open boundary conditions, temperature (deg C) |
| 124 |
OBEu_arctic_210x192.bin South open boundary conditions, U-velocity (m/s) |
| 125 |
OBEv_arctic_210x192.bin South open boundary conditions, V-velocity (m/s) |
| 126 |
OBWs_arctic_210x192.stable West open boundary conditions, salinity (g/kg) |
| 127 |
OBWt_arctic_210x192.stable West open boundary conditions, temperature (deg C) |
| 128 |
OBWu_arctic_210x192.bin West open boundary conditions, U-velocity (m/s) |
| 129 |
OBWv_arctic_210x192.bin West open boundary conditions, V-velocity (m/s) |
| 130 |
|
| 131 |
Runtime parameter files. |
| 132 |
eedata |
| 133 |
data |
| 134 |
data.pkg |
| 135 |
data.cal |
| 136 |
data.exf |
| 137 |
data.gmredi |
| 138 |
data.seaice |
| 139 |
data.obcs |
| 140 |
data.kpp |
| 141 |
data.salt_plume |
| 142 |
data.diagnostics |
Parent Directory
| 
Revision Log
| 
Revision Graph