1 |
Southern-Ocean box with Biochemistry, using Open-Boundary Conditions |
2 |
(pkg/obcs) at Northern, Eastern and Western edges of the domain. |
3 |
====================================================================== |
4 |
|
5 |
This experiment illustrates and tests the use of package DIC with OBCS. |
6 |
|
7 |
The configuration (e.g., resolution), model parameters and forcing are |
8 |
almost identical to tutorial_global_oce_biogeo expect that the horizontal |
9 |
domain is limited to a sub-domain around Drake passage with open-boundary |
10 |
conditions coming from the last year of a 2 yrs global simulation and |
11 |
initial conditions taken from t=1.yr of this same simulation. |
12 |
This enable to compare directly the results of the first year simulation |
13 |
of this regional set-up with the 2nd year results of the global set-up run. |
14 |
|
15 |
Configure and compile the code: |
16 |
cd build |
17 |
../../../tools/genmake2 -mods ../code [-of my_platform_optionFile] |
18 |
make depend |
19 |
make |
20 |
cd .. |
21 |
|
22 |
To run: |
23 |
cd run |
24 |
ln -s ../input/* . |
25 |
ln -s ../build/mitgcmuv . |
26 |
./mitgcmuv > output.txt |
27 |
cd .. |
28 |
|
29 |
There is comparison output in the directory: |
30 |
results/output.txt |
31 |
|
32 |
---------------------------------------------------------------------- |
33 |
To generate inital and open boundary conditions: |
34 |
a) The global set-up (using executable from: tutorial_global_oce_biogeo/build) |
35 |
was run for 2 years using model-parameter from inp_global/. |
36 |
The only differences with the ones used in tutorial_global_oce_biogeo are: |
37 |
- CD-Scheme is turned off (useCDscheme=F), since it is not implemented for OBCS ; |
38 |
- as a consequence, horizontal viscosity (viscAh) is increased from 2.E5 to 3.E5; |
39 |
- convective-adjustment diffusivity (ivdc_kappa) is reduced from 100. to 10.m^2/s; |
40 |
- implicit vertical viscosity is turned off (not needed); |
41 |
- free-surface exact conservation is turned on (exactConserv=T); |
42 |
- use the Avective-Form of GM (in data.gmredi: GM_AdvForm=T) instead of Skew-Flux; |
43 |
|
44 |
to repeat this global simutlation, from so_box_biogeo/ |
45 |
mkdir run_glob ; cd run_glob |
46 |
ln -s ../inp_global/* . |
47 |
../inp_global/prepare_run |
48 |
ln -s ../../tutorial_global_oce_biogeo/build/mitgcmuv . |
49 |
./mitgcmuv > output.glob |
50 |
cd .. |
51 |
|
52 |
b) Then regional model forcing files, initial conditions and open-boundary |
53 |
conditions files are generated using the matlab script: "mk_box_input.m" |
54 |
cd inp_global/ |
55 |
matlab |
56 |
>> mk_box_input |