1 |
Instructions for running optimized (A1) simulation from Nguyen et al (2011) |
2 |
http://ecco2.org/manuscripts/2011/NguyenJGR2011.pdf |
3 |
|
4 |
=============== |
5 |
|
6 |
DM 9-Jul-2015 |
7 |
Updated for checkpoint65m (2015/06/15) on pleiades |
8 |
|
9 |
-> get a copy of ftp://ecco2.jpl.nasa.gov/data6/arctic/run_template_cube81 |
10 |
also available on pfe:/nobackup/dmenemen/arctic18km/run_template_cube81 |
11 |
and harbor.mit.edu:/net/nares/raid8/ecco-shared/arctic18km/run_template_cube81 |
12 |
|
13 |
-> get a copy of ftp://ecco2.jpl.nasa.gov/data2/data/atmos/jra25 |
14 |
also available on pfe:/nobackup/hzhang1/forcing/jra25 |
15 |
|
16 |
cvs co MITgcm_contrib/arctic/cs_18km |
17 |
cvs co -r checkpoint65m MITgcm_code |
18 |
cd MITgcm |
19 |
mkdir build run |
20 |
cd build |
21 |
module purge |
22 |
module load comp-intel/2015.0.090 mpi-sgi/mpt.2.11r13 netcdf/4.0 |
23 |
../tools/genmake2 -of ../tools/build_options/linux_amd64_ifort+mpi_ice_nas -mods ../../MITgcm_contrib/arctic/cs_18km/code |
24 |
make depend |
25 |
make -j 16 |
26 |
|
27 |
cd ../run |
28 |
cp ../../MITgcm_contrib/arctic/cs_18km/input/* . |
29 |
rm data.obcs_* |
30 |
ln -sf ../../run_template_cube81/*.bin . |
31 |
ln -sf ../../run_template_cube81/BATHY_cube81_420x384_arctic . |
32 |
ln -sf ../../run_template_cube81/WOA05_* . |
33 |
ln -sf ../../run_template_cube81/*cube81 . |
34 |
ln -sf ../../run_template_cube81/+ob1979_2014_merge/* . |
35 |
ln -sf ../../jra25/jra25_*_199[1-9] . |
36 |
ln -sf ../../jra25/jra25_*_2* . |
37 |
rm jra25_rain_n* |
38 |
ln -sf ../build/mitgcmuv . |
39 |
qsub run_arctic_cs_18km_nas |
40 |
|
41 |
-> ftp://ecco2.jpl.nasa.gov/data1/arctic/arctic_18km/GRACErunoff/matlab/ |
42 |
ftp://ecco2.jpl.nasa.gov/data1/arctic/arctic_18km/GRACErunoff/matlab/lookat.m |
43 |
contains some routines for looking at this solution, in particular: |
44 |
ftp://ecco2.jpl.nasa.gov/data1/arctic/arctic_18km/GRACErunoff/figs/FIG9.jpg |
45 |
is a recreation of Figure 9 in Nguyen et al (2011): |
46 |
the results are similar but not identical, e.g., thinner sea ice near |
47 |
Canadian Archipelago for TBD reasons. |
48 |
|
49 |
=============== |
50 |
|
51 |
ATN 23-Dec-2013 |
52 |
|
53 |
SEAICE_OPTIONS.h_atn_May2011: Paper: Nguyen et al 2011 |
54 |
|
55 |
Cpp_OPTIONS.h_oldcode_before2012: code similar to what is used in Nguyen et al 2011, but compatible with codes around Sep2011 |
56 |
SEAICE_OPTIONS.h_Ian_Sep2011: Equivalence to above but updated to be compatible with Ian's code on that date |
57 |
|
58 |
After 2011, pkg/seaice has updated significantly with the merge of Ian's code into the main branch. |
59 |
These followings _phAug2012 files include appropriate flag adjustments + param name change for new code. |
60 |
They should reproduce almost identical results to Nguyen et al 2011 paper, except for whatever that has |
61 |
improved in ocean/sea-ice physics. An offline comparison done by ATN on Aug2012 showed the two runs Aug2012 |
62 |
and Nguyen et al 2011 papers are VERY close to each other, thus ATN considered the two solutions are equivalent. |
63 |
|
64 |
CPP_OPTIONS.h_phAug2012 |
65 |
SEAICE_OPTIONS.h_phAug2012 |
66 |
|
67 |
To run this set up, use binary inputs from: |
68 |
ftp://ecco2.jpl.nasa.gov/data6/arctic/run_template_cube81/ |
69 |
obcs for the years 1979--2014: +ob1979_2014_merge |
70 |
|
71 |
There are two bathymetries: |
72 |
1) original from Nguyen et al 2011: BATHY_cube81_420x384_arctic |
73 |
2) an updated version to blank out some CAA passages (set bathymetry to zero) |
74 |
because 18km can not resolve these passages resutling in excessive sea ice thickness |