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

Diff of /MITgcm_contrib/high_res_cube/README_ice

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

revision 1.2 by dimitri, Sun Nov 23 20:11:39 2003 UTC revision 1.32 by dimitri, Wed Mar 10 04:49:51 2004 UTC
# Line 1  Line 1 
1  Instructions for setting up a cube-sphere integration with sea-ice.  Instructions for setting up a cube-sphere integration with sea-ice.
2  Date last tested: November 21, 2003.  Date last tested: March 9, 2004.
3    Code and input used for integrating cube5.
4    
5    
6    ##########################################################################
7    ##########################################################################
8    # getting the code from anonymous CVS server
9    
10  setenv CVSROOT :pserver:cvsanon@mitgcm.org:/u/u0/gcmpack  setenv CVSROOT :pserver:cvsanon@mitgcm.org:/u/u0/gcmpack
11  cvs login ( CVS password: cvsanon )  cvs login ( CVS password: cvsanon )
12  cvs co MITgcm_contrib/high_res_cube/code-mods  cvs co -r hrcube5 MITgcm_contrib/high_res_cube/README_ice
13  cvs co MITgcm_contrib/high_res_cube/matlab-grid-converter  cvs co -r hrcube5 MITgcm_contrib/high_res_cube/code-mods
14  cvs co MITgcm_contrib/high_res_cube/README_ice  cvs co -r hrcube5 MITgcm_contrib/high_res_cube/input
15  cvs co -r checkpoint52a_post MITgcm_code  cvs co -r hrcube5 MITgcm_contrib/high_res_cube/results
16  cvs co MITgcm/verification/global_ocean.cs32x15  cvs co -r hrcube5 MITgcm_code
17    cvs co -r hrcube5 MITgcm/verification/global_ocean.cs32x15
18    
19    
20    ##########################################################################
21    ##########################################################################
22    # 32*32*6*15, 1-cpu cube sphere
23    
24    cd MITgcm/verification/global_ocean.cs32x15/code
25    \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
26    \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s12t_16x32/* .
27    
28    cd ../build
29    ../../../tools/genmake2 -mods=../code
30    make depend
31    \rm getcon.F
32    make
33    
34    cd ..
35    mkdir run
36    cd run
37    \cp ../input/* .
38    \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
39    ../build/mitgcmuv >& output.txt
40    
41    comparison output is in:
42    ../../../../MITgcm_contrib/high_res_cube/results/output.txt
43    
44    to use matlab to look at the output
45    matlab
46    
47    for fld={'Eta','VICE','UICE','HEFF','AREA','Qsw','Qnet','FV','FU','EmPmR'}
48        tmp=permute(readbin([fld{1} '.0000000020.data'],[32,6,32]),[1 3 2]);
49        cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause
50    end
51    
52    heff=permute(readbin(['HEFF.0000000020.data'],[32,6,32]),[1 3 2]);
53    mask=0*heff; mask(find(heff>.001))=1;
54    uice=mask.*permute(readbin(['UICE.0000000020.data'],[32,6,32]),[1 3 2]);
55    vice=mask.*permute(readbin(['VICE.0000000020.data'],[32,6,32]),[1 3 2]);
56    clf
57    subplot(221), mypcolor(uice(:,:,3)'), colorbar, title('uice, tile 3')
58    subplot(222), mypcolor(uice(:,:,6)'), colorbar, title('uice, tile 6')
59    subplot(223), mypcolor(vice(:,:,3)'), colorbar, title('vice, tile 3')
60    subplot(224), mypcolor(vice(:,:,6)'), colorbar, title('vice, tile 6')
61    
62    
63    ##########################################################################
64    ##########################################################################
65    # 32*32*6*15, 2-cpu cube sphere on the altices with efc
66    
67    cd MITgcm/verification/global_ocean.cs32x15/code
68    \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
69    \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s12t_16x32/* .
70    \cp SIZE.h_mpi SIZE.h
71    
72    cd ../build
73    \rm *
74    ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi
75    make depend
76    \rm getcon.F
77    make
78    
79    cd ..
80    mkdir run
81    cd run
82    \rm *
83    \cp ../input/* .
84    \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
85    mpirun -np 2 ../build/mitgcmuv
86    
87    
88    ##########################################################################
89    ##########################################################################
90    # 510*510*6*50, 54-cpu cube sphere on the altices with efc
91    # requires input files in sudirectories run_template and ncep_rgau
92    # presently residing under orion:/tmp1/dmenem/cube
93    
94  cd MITgcm/verification/global_ocean.cs32x15/code  cd MITgcm/verification/global_ocean.cs32x15/code
95  cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
96  cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s12t_16x32/* .  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* .
97  \rm mdsio_readfield.F  \cp SIZE.h_54 SIZE.h
98    \cp packages.conf_hr packages.conf
99    
100    cd ../build
101    \rm *
102    ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi
103    make depend
104    \rm getcon.F
105    make
106    
107  cd ..  cd ..
108  mkdir build  mkdir run
109  cd build  cd run
110  ../code/genmake2 -mods=../code  \rm *
111    \cp ../input/* .
112    \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
113    qsub job_orion
114    
115    
116    ##########################################################################
117    ##########################################################################
118    # 510*510*6*50, 205-tile, 41-cpu cube sphere on the altices with efc
119    # requires input files in sudirectories run_template and ncep_rgau
120    # presently residing under orion:/tmp1/dmenem/cube
121    
122    cd MITgcm/verification/global_ocean.cs32x15/code
123    \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
124    \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s205t_85x85/* .
125    \cp SIZE.h_41 SIZE.h
126    \cp packages.conf_hr packages.conf
127    
128    cd ../build
129    \rm *
130    ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi
131  make depend  make depend
132    \rm getcon.F
133  make  make
134    
135  cd ..  cd ..
136  mkdir run  mkdir run
137  cd run  cd run
138  cp ../input/* .  \rm *
139  cp ../../../../MITgcm_contrib/high_res_cube/matlab-grid-converter/tile00* .  \cp ../input/* .
140  cp ../../../../MITgcm_contrib/high_res_cube/input/* .  \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
141  ../build/mitgcmuv >& output  qsub job_orion41
142    
143    
144    ##########################################################################
145    ##########################################################################
146    # 510*510*6*50, 216-cpu cube sphere on altix with efc
147    # requires input files in sudirectories run_template and ncep_rgau
148    # these can be obained from lou:/u/menemenl/cube
149    
150    cd MITgcm/verification/global_ocean.cs32x15/code
151    \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
152    \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* .
153    \cp SIZE.h_216 SIZE.h
154    \cp packages.conf_hr packages.conf
155    
156    cd ../build
157    \rm *
158    ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi_altix
159    make depend
160    \rm getcon.F
161    make
162    
163    cd ..
164    mkdir run
165    cd run
166    \rm *
167    \cp ../input/* .
168    \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
169    qsub job_altix
170    
171    
172    ##########################################################################
173    ##########################################################################
174    # matlab script for generating pChkptFreq
175    
176    nTimeSteps=26352;
177    for niter0=216:26352:(26352*11)
178     pChkptFreq=(niter0+nTimeSteps)*1200;
179     day=pChkptFreq/24/60/60;
180     disp(['niter0=' int2str(niter0) ',day=' int2str(day) ', ' ...
181     datestr(datenum(1992,1,1)+day) ', pChkptFreq=' int2str(pChkptFreq)])
182    end
183    
184    niter0 = 216   , day = 369 , 04-Jan-1993, pChkptFreq = 31881600
185    niter0 = 26568 , day = 735 , 05-Jan-1994, pChkptFreq = 63504000
186    niter0 = 52920 , day = 1101, 06-Jan-1995, pChkptFreq = 95126400
187    niter0 = 79272 , day = 1467, 07-Jan-1996, pChkptFreq = 126748800
188    niter0 = 105624, day = 1833, 07-Jan-1997, pChkptFreq = 158371200
189    niter0 = 131976, day = 2199, 08-Jan-1998, pChkptFreq = 189993600
190    niter0 = 158328, day = 2565, 09-Jan-1999, pChkptFreq = 221616000
191    niter0 = 184680, day = 2931, 10-Jan-2000, pChkptFreq = 253238400
192    niter0 = 211032, day = 3297, 10-Jan-2001, pChkptFreq = 284860800
193    niter0 = 237384, day = 3663, 11-Jan-2002, pChkptFreq = 316483200
194    niter0 = 263736, day = 4029, 12-Jan-2003, pChkptFreq = 348105600
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.32
  ViewVC Help
Powered by ViewVC 1.1.22