[MITgcm-support] some questions

[Dimitris Menemenlis]

Hi Liuyu,

> 1 Does the thsice package include sea ice dynamic section?  

Not yet.

> 2 How to link the global ocean to the seaice package?

There are instructions how to do this in:
http://mitgcm.org/cgi-bin/viewcvs.cgi/MITgcm_contrib/high_res_cube/README_ice
I am reproducing relevant section below:

cvs co MITgcm_contrib/high_res_cube/README_ice
cvs co MITgcm_contrib/high_res_cube/code-mods
cvs co MITgcm_contrib/high_res_cube/input
cvs co MITgcm_contrib/high_res_cube/results
cvs co MITgcm_code
cvs co MITgcm/verification/global_ocean.cs32x15
cd MITgcm/verification/global_ocean.cs32x15/code
\cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
\cp ../../../utils/exch2/code-mods/s12t_16x32/* .
cd ../build
\rm *
../../../tools/genmake2 -mods=../code
make depend
make
cd ..
mkdir run
cd run
\rm *
\cp ../input/* .
\cp ../../../../MITgcm_contrib/high_res_cube/input/* .
../build/mitgcmuv >& output.txt

comparison output is in:
../../../../MITgcm_contrib/high_res_cube/results/output.txt

> 3 How to change the resolution of the global ocean model?

Changing vertical resolution is easy.  You change Nr in the SIZE.h header and
tRef, sRef, and delR in the runtime "data" file.  You also need to generate new
initial temperature and salinity files, hydrogThetaFile and hydrogSaltFile in
the runtime "data" file. that have the correct size.

Changing horizontal resolution for a lat/long grid is similarly easy but
changing resolution for a cubed-sphere is trickier as you need to generate new
definition files for the grid.  What resolution do you want to run at?
http://mitgcm.org/cgi-bin/viewcvs.cgi/MITgcm_contrib/high_res_cube/README_ice
also contains an example of a cubed-sphere configuration with approximately
18-km horizontal grid spacing.  Some results form integrations at that
resolution are described here http://ecco2.org/ and the input and output from
these integrations is freely available http://ecco2.org/products/

> 4 How to prepare needed input fields? 

Input fields can be direct access binary files on an arbitrary lat/long grid,
which is then described in the runtime data.exf file, for example,
http://mitgcm.org/cgi-bin/viewcvs.cgi/MITgcm_contrib/high_res_cube/input/data.exf

Dimitris

-- 
Dimitris Menemenlis <menemenlis at jpl.nasa.gov>
Jet Propulsion Lab, California Institute of Technology
MS 300-323, 4800 Oak Grove Dr, Pasadena CA 91109-8099
tel: 818-354-1656;  fax: 818-393-6720