/[MITgcm]/MITgcm_contrib/high_res_cube/README_ice
ViewVC logotype

Contents of /MITgcm_contrib/high_res_cube/README_ice

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


Revision 1.7 - (show annotations) (download)
Wed Dec 10 19:40:13 2003 UTC (21 years, 7 months ago) by dimitri
Branch: MAIN
Changes since 1.6: +5 -6 lines
o changes to permit hi-res, cubed-sphere, configuration
  - added useSingleCpuIO capability to mdsio_readfield.F and exf_interp.F
  - added "#undef ALLOW_USE_MPI" support to eesupp/src/scatter_2d.F
  - added pkg/exf/exf_set_uv.F for on-the-fly interpolation and rotation
    of surface winds for the cube
  (A verification experiment, based on global_ocean.cs32x15, but with
   pkg/seaice turned on and with on-the-fly interpolation from the NCEP
   Gaussian grid is described in MITgcm_contrib/high_res_cube/README_ice,
   complete with example surface forcing files and matlab scripts to look
   at the output and compare it with that of global_ocean.cs32x15.)

1 Instructions for setting up a cube-sphere integration with sea-ice.
2 Date last tested: December 6, 2003.
3
4 setenv CVSROOT :pserver:cvsanon@mitgcm.org:/u/u0/gcmpack
5 cvs login ( CVS password: cvsanon )
6
7 cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/README_ice
8 cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/code-mods
9 cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/matlab-grid-converter
10 cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/input
11 cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/results
12
13 cvs co -r checkpoint52d_pre MITgcm_code
14 cvs co -r checkpoint52d_pre MITgcm/verification/global_ocean.cs32x15
15
16 cd MITgcm/verification/global_ocean.cs32x15/code
17 cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
18 cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s12t_16x32/* .
19 \rm mdsio_readfield.F
20
21 cd ..
22 mkdir build
23 cd build
24 ../code/genmake2 -mods=../code
25 make depend
26 make
27
28 cd ..
29 mkdir run
30 cd run
31 cp ../inp_thsice/* .
32 cp ../input/* .
33 cp ../../../../MITgcm_contrib/high_res_cube/matlab-grid-converter/tile00* .
34 cp ../../../../MITgcm_contrib/high_res_cube/input/* .
35 ../build/mitgcmuv >& output.txt
36
37 comparison output is in:
38 ../../../../MITgcm_contrib/high_res_cube/results/output.txt
39
40 to use matlab to look at the output
41 matlab
42 for fld={'Eta','VICE','UICE','HEFF','AREA','Qsw','Qnet','FV','FU','EmPmR'}
43 tmp=permute(readbin([fld{1} '.0000000020.data'],[32,6,32]),[1 3 2]);
44 cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause
45 end

  ViewVC Help
Powered by ViewVC 1.1.22