/[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.7 by dimitri, Wed Dec 10 19:40:13 2003 UTC revision 1.33 by dimitri, Tue Mar 16 17:20:07 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: December 6, 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 -r hrcube5 MITgcm_contrib/high_res_cube/README_ice
13    cvs co -r hrcube5 MITgcm_contrib/high_res_cube/code-mods
14    cvs co -r hrcube5 MITgcm_contrib/high_res_cube/input
15    cvs co -r hrcube5 MITgcm_contrib/high_res_cube/results
16    cvs co -r hrcube5 MITgcm_code
17    cvs co -r hrcube5 MITgcm/verification/global_ocean.cs32x15
18    
 cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/README_ice  
 cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/code-mods  
 cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/matlab-grid-converter  
 cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/input  
 cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/results  
19    
20  cvs co -r checkpoint52d_pre MITgcm_code  ##########################################################################
21  cvs co -r checkpoint52d_pre MITgcm/verification/global_ocean.cs32x15  ##########################################################################
22    # 32*32*6*15, 1-cpu cube sphere
23    
24  cd MITgcm/verification/global_ocean.cs32x15/code  cd MITgcm/verification/global_ocean.cs32x15/code
25  cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
26  cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s12t_16x32/* .  \cp ../../../utils/exch2/code-mods/s12t_16x32/* .
 \rm mdsio_readfield.F  
27    
28  cd ..  cd ../build
29  mkdir build  ../../../tools/genmake2 -mods=../code
 cd build  
 ../code/genmake2 -mods=../code  
30  make depend  make depend
31    \rm getcon.F
32  make  make
33    
34  cd ..  cd ..
35  mkdir run  mkdir run
36  cd run  cd run
37  cp ../inp_thsice/* .  \cp ../input/* .
38  cp ../input/* .  \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
 cp ../../../../MITgcm_contrib/high_res_cube/matlab-grid-converter/tile00* .  
 cp ../../../../MITgcm_contrib/high_res_cube/input/* .  
39  ../build/mitgcmuv >& output.txt  ../build/mitgcmuv >& output.txt
40    
41  comparison output is in:  comparison output is in:
# Line 39  comparison output is in: Line 43  comparison output is in:
43    
44  to use matlab to look at the output  to use matlab to look at the output
45  matlab  matlab
46    
47  for fld={'Eta','VICE','UICE','HEFF','AREA','Qsw','Qnet','FV','FU','EmPmR'}  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]);      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      cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause
50  end  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 ../../../utils/exch2/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
95    \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
96    \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* .
97    \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 ..
108    mkdir run
109    cd run
110    \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
132    \rm getcon.F
133    make
134    
135    cd ..
136    mkdir run
137    cd run
138    \rm *
139    \cp ../input/* .
140    \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
141    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.7  
changed lines
  Added in v.1.33
  ViewVC Help
Powered by ViewVC 1.1.22