/[MITgcm]/MITgcm_contrib/eh3/llc/ecco-godae/README
ViewVC logotype

Contents of /MITgcm_contrib/eh3/llc/ecco-godae/README

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.5 - (show annotations) (download)
Wed Mar 28 18:29:45 2007 UTC (18 years, 4 months ago) by heimbach
Branch: MAIN
CVS Tags: HEAD
Changes since 1.4: +52 -4 lines
Add some description

1
2 LAT-LON-Polar-Cap setup
3 -----------------------
4 ed@eh3.com, cnh@mit.edu, heimbach@mit.edu
5
6 Setup code/, input/
7 -------------------
8 Note that the *.mitgrid files must be named using the following pattern
9 or they will not be read correctly by MITgcm:
10
11 cd run
12 for i in 1 2 3 4 5 6; do
13 ln -s ../input/llc_p90_${i}.mitgrid tile00${i}.mitgrid
14 done
15
16 Change topology:
17
18 cd MITgcm/utils/exch2/matlab-topology-generator/
19
20 modify matlab script driver.m
21 line
22 nr = 360; ng = 90; nb = 90; tnx=90;tny=90;
23 to (e.g.)
24 nr = 360; ng = 90; nb = 90; tnx=45;tny=45;
25 for compact format, set
26 mapIO=1
27 make sure blanklist.txt is correct
28 (or empty, if no tiles to be omitted)
29 then run driver.m
30
31 This generates routine and header
32 w2_e2setup.F
33 W2_EXCH2_TOPOLOGY.h
34 which need to go in code/
35
36 Accordingly, bathy.0??.0??.data need changing.
37
38 N.B. Domain decomposition nPy not active with exch2, i.e.
39 possible decomps are (for orig. with 17 tiles):
40 sNx = 90 30 90 30
41 sNy = 90 30 90 30
42 nSx = 17 17 1 1
43 nSy = 1 1 1 1
44 nPx = 1 9 17 153
45 nPy = 1 1 1 1
46
47 I've generated artificial tile #18 with all zeros
48 blanklist.txt needs to be generated by hand.
49
50 Setup code_15lev/, input_15lev/
51 -------------------------------
52
53 * Real bathymetry, 15 vertical levels,
54 * monthly climatological forcing via exf package, uses
55 4x4 deg. fields with on-the-fly interpolation
56 * Bathymetry is generated via script ... in ...
57 * Initial conditions for T, S are generated via script ... in ...
58
59 All input fields can be downloaded using script input_15lev/prepare_run
60 (which also contains link pattern to standard tile00${i}.mitgrid fields),
61 then link to run directory.
62
63 Setup code_50lev/, input_50lev/
64 -------------------------------
65
66 * Based on ETOPO2 bathy,
67 * 50 vertical level exactly matching high_res_cube (ECCO2), incl. partial cell params.
68 * forced by 6-hourly NCEP atmos. state, uses original reduced Gaussian grid
69 and exf interpolation on-the-fly
70 * initial conditions generated via scripts in
71 MITgcm_contrib/eh3/llc/ecco-godae/climatology/
72 resulting fields are
73 llc90x90x360_50lev_ptemp.bin, llc90x90x360_50lev_salt.bin
74 these fields are run are run with a diffusion/conv.adj. setup of MITgcm
75 for 10 days to remove some instabilities
76 * input fields can be obtained at
77 http://channel.mit.edu/datasets/MITgcm_input_files/MITgcm_contrib/eh3/llc/ecco-godae/data_50lev/
78
79 Required input to generate bathy & initial conditions:
80 (1)
81 (in directory input/ )
82 llc_p90_1.nc
83 llc_p90_2.nc
84 llc_p90_3.nc
85 llc_p90_4.nc
86 llc_p90_5.nc
87 llc_p90_6.nc
88 (2)
89 ETOPO2.raw.bin
90 (3)
91 directory CLIMannual/
92
93 Generation steps:
94
95 1) run script "gen_bathy"
96 Input:
97 -----
98 llc_p90_1.nc
99 llc_p90_2.nc
100 llc_p90_3.nc
101 llc_p90_4.nc
102 llc_p90_5.nc
103 llc_p90_6.nc
104 Output:
105 -----
106 llc_p90_bathy.f001.nc
107 llc_p90_bathy.f002.nc
108 llc_p90_bathy.f003.nc
109 llc_p90_bathy.f004.nc
110 llc_p90_bathy.f005.nc
111
112 2) run script "gen_initial_cond_50_lev"
113 Input:
114 -----
115 (a) above llc_p90_*.nc
116 (b) orig. climatology fields in CLIMannual/
117 Output:
118 -----
119
120 3) run script "t18_to_global_inputs.m"
121
122 4) run script "fill_shallow_regions_and_isolated_basins.m"

  ViewVC Help
Powered by ViewVC 1.1.22