/[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.11 by dimitri, Sun Dec 14 05:16:46 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 11, 2003.  Date last tested: March 9, 2004.
3    Code and input used for integrating cube5.
4    
5    
6  ##########################################################################  ##########################################################################
7  ##########################################################################  ##########################################################################
8  # 32*32*6*15, 1-cpu cube sphere  # 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 checkpoint52d_pre MITgcm_contrib/high_res_cube/README_ice  cvs co -r hrcube5 MITgcm_contrib/high_res_cube/README_ice
13  cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/code-mods  cvs co -r hrcube5 MITgcm_contrib/high_res_cube/code-mods
14  cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/matlab-grid-converter  cvs co -r hrcube5 MITgcm_contrib/high_res_cube/input
15  cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/input  cvs co -r hrcube5 MITgcm_contrib/high_res_cube/results
16  cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/results  cvs co -r hrcube5 MITgcm_code
17  cvs co -r checkpoint52d_pre MITgcm_code  cvs co -r hrcube5 MITgcm/verification/global_ocean.cs32x15
18  cvs co -r checkpoint52d_pre MITgcm/verification/global_ocean.cs32x15  
19    
20    ##########################################################################
21    ##########################################################################
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 packages.conf  
27    
28  cd ..  cd ../build
 mkdir build  
 cd build  
29  ../../../tools/genmake2 -mods=../code  ../../../tools/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 ../input/* .  \cp ../input/* .
 \cp ../../../../MITgcm_contrib/high_res_cube/matlab-grid-converter/tile00* .  
38  \cp ../../../../MITgcm_contrib/high_res_cube/input/* .  \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
39  ../build/mitgcmuv >& output.txt  ../build/mitgcmuv >& output.txt
40    
# Line 41  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  # 32*32*6*15, 2-cpu cube sphere on the altices with efc
66    
 setenv CVSROOT :pserver:cvsanon@mitgcm.org:/u/u0/gcmpack  
 cvs login ( CVS password: cvsanon )  
 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  
 cvs co -r checkpoint52d_pre MITgcm_code  
 cvs co -r checkpoint52d_pre MITgcm/verification/global_ocean.cs32x15  
   
67  cd MITgcm/verification/global_ocean.cs32x15/code  cd MITgcm/verification/global_ocean.cs32x15/code
68  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
69  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s12t_16x32/* .  \cp ../../../utils/exch2/code-mods/s12t_16x32/* .
70  \cp SIZE.h_mpi SIZE.h  \cp SIZE.h_mpi SIZE.h
 \cp CPP_EEOPTIONS.h_mpi CPP_EEOPTIONS.h  
 \rm mdsio_readfield.F  
71    
72  cd ..  cd ../build
73  mkdir build  \rm *
 cd build  
74  ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi  ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi
75  make depend  make depend
76    \rm getcon.F
77  make  make
78    
79  cd ..  cd ..
80  mkdir run  mkdir run
81  cd run  cd run
82    \rm *
83  \cp ../input/* .  \cp ../input/* .
 \cp ../../../../MITgcm_contrib/high_res_cube/matlab-grid-converter/tile00* .  
84  \cp ../../../../MITgcm_contrib/high_res_cube/input/* .  \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
85  mpirun -np 2 ../build/mitgcmuv  mpirun -np 2 ../build/mitgcmuv
86    
# Line 87  mpirun -np 2 ../build/mitgcmuv Line 88  mpirun -np 2 ../build/mitgcmuv
88  ##########################################################################  ##########################################################################
89  ##########################################################################  ##########################################################################
90  # 510*510*6*50, 54-cpu cube sphere on the altices with efc  # 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  setenv CVSROOT :pserver:cvsanon@mitgcm.org:/u/u0/gcmpack  # presently residing under orion:/tmp1/dmenem/cube
 cvs login ( CVS password: cvsanon )  
 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  
 cvs co -r checkpoint52d_pre MITgcm_code  
 cvs co -r checkpoint52d_pre MITgcm/verification/global_ocean.cs32x15  
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/s216t_85x85/* .  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* .
97  \cp SIZE.h_54 SIZE.h  \cp SIZE.h_54 SIZE.h
 \cp CPP_EEOPTIONS.h_mpi CPP_EEOPTIONS.h  
98  \cp packages.conf_hr packages.conf  \cp packages.conf_hr packages.conf
 \rm mdsio_readfield.F  
99    
100  cd ..  cd ../build
101  mkdir build  \rm *
 cd build  
102  ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi  ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi
103  make depend  make depend
104    \rm getcon.F
105  make  make
106    
107  cd ..  cd ..
108  mkdir run  mkdir run
109  cd run  cd run
110    \rm *
111  \cp ../input/* .  \cp ../input/* .
112  \cp ../../../../MITgcm_contrib/high_res_cube/input/* .  \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
 \cp data_hr data  
 \ln -sf ../../../../run_template/tile00* .  
 \ln -sf ../../../../run_template/LEV* .  
 \ln -sf ../../../../run_template/BATHY_18360x85_filled.bin .  
113  qsub job_orion  qsub job_orion
114    
115    
116  ##########################################################################  ##########################################################################
117  ##########################################################################  ##########################################################################
118  # 510*510*6*50, 216-cpu cube sphere on the altices with efc  # 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  setenv CVSROOT :pserver:cvsanon@mitgcm.org:/u/u0/gcmpack  cd MITgcm/verification/global_ocean.cs32x15/code
123  cvs login ( CVS password: cvsanon )  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
124  cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/README_ice  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s205t_85x85/* .
125  cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/code-mods  \cp SIZE.h_41 SIZE.h
126  cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/matlab-grid-converter  \cp packages.conf_hr packages.conf
127  cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/input  
128  cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/results  cd ../build
129  cvs co -r checkpoint52d_pre MITgcm_code  \rm *
130  cvs co -r checkpoint52d_pre MITgcm/verification/global_ocean.cs32x15  ../../../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  cd MITgcm/verification/global_ocean.cs32x15/code
151  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
152  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* .  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* .
153  \cp SIZE.h_216 SIZE.h  \cp SIZE.h_216 SIZE.h
154  \cp CPP_EEOPTIONS.h_mpi CPP_EEOPTIONS.h  \cp packages.conf_hr packages.conf
 \rm mdsio_readfield.F  
155    
156  cd ..  cd ../build
157  mkdir build  \rm *
158  cd build  ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi_altix
 ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi  
159  make depend  make depend
160    \rm getcon.F
161  make  make
162    
163  cd ..  cd ..
164  mkdir run  mkdir run
165  cd run  cd run
166    \rm *
167  \cp ../input/* .  \cp ../input/* .
168  \cp ../../../../MITgcm_contrib/high_res_cube/input/* .  \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
169  \cp data_hr data  qsub job_altix
170  \cp data.pkg_hr data.pkg  
171  \ln -sf ../../../../run_template/tile00* .  
172  \ln -sf ../../../../run_template/LEV* .  ##########################################################################
173  \ln -sf ../../../../run_template/BATHY_18360x85_filled.bin .  ##########################################################################
174  qsub job_orion  # 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.11  
changed lines
  Added in v.1.33
  ViewVC Help
Powered by ViewVC 1.1.22