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