1 |
Tutorial Example: "P coordinate Global Ocean" |
2 |
(Global Ocean Simulation at 4o Resolution in Pressure Coordinates) |
3 |
================================================================== |
4 |
(formerly "global_ocean_pressure") |
5 |
|
6 |
|
7 |
Configure and compile the code: |
8 |
cd build |
9 |
../../../tools/genmake2 -mods ../code [-of my_platform_optionFile] |
10 |
make depend |
11 |
make |
12 |
|
13 |
To run: |
14 |
cd build |
15 |
ln -s ../input/* . |
16 |
./mitgcmuv > output.txt |
17 |
|
18 |
There is comparison output in the directory: |
19 |
results/output.txt |
20 |
|
21 |
Comments: |
22 |
o the set up is similar to that of global_ocean.90x40x15 |
23 |
o the code directory contains calc_phi_hyd.F, where the potential is computed |
24 |
according to the more natural finite volume discretization. Finite difference |
25 |
discretization is energy conserving, but the representation of the "fixed" |
26 |
surface (interface ocean-atmosphere) is less consistent. |
27 |
o the code directory also contains dynamics.F which calls |
28 |
remove_mean_rl.F, a generic routine, to remove the mean from the |
29 |
diagnostic variable phiHydLow (sea surface height/gravity in pressure |
30 |
coordinates) |
31 |
|
32 |
changes: 07 Feb. 2003 (jmc): |
33 |
o find difficult to maintain the local version of dynamics.F up to date. |
34 |
therefore, has been remove from the code directory. |
35 |
One can recover the same version (but up to date) simply |
36 |
by activating the commented lines [between lines Cml( and Cml) ], |
37 |
at the end of the standard version of dynamics.F |
38 |
o finite volume form of calc_phi_hyd.F is now a standard option. |
39 |
only needs to set integr_GeoPot=1 in file "data" to select this form. |