/[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.9 by dimitri, Sun Dec 14 02:06:15 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: December 11, 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  # 32*32*6*15, 1-cpu cube sphere  # 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 checkpoint52d_pre MITgcm_contrib/high_res_cube/README_ice  cvs co -r hrcube4a MITgcm_contrib/high_res_cube/README_ice
14  cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/code-mods  cvs co -r hrcube4a MITgcm_contrib/high_res_cube/code-mods
15  cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/matlab-grid-converter  cvs co -r hrcube4a MITgcm_contrib/high_res_cube/matlab-grid-converter
16  cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/input  cvs co -r hrcube4a MITgcm_contrib/high_res_cube/input
17  cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/results  cvs co -r hrcube4a MITgcm_contrib/high_res_cube/results
18  cvs co -r checkpoint52d_pre MITgcm_code  cvs co -r hrcube4a MITgcm_code
19  cvs co -r checkpoint52d_pre MITgcm/verification/global_ocean.cs32x15  cvs co -r hrcube4a MITgcm/verification/global_ocean.cs32x15
20    
21    
22    ##########################################################################
23    ##########################################################################
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
 \cp ../inp_thsice/* .  
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/* .
# Line 42  comparison output is in: Line 46  comparison output is in:
46    
47  to use matlab to look at the output  to use matlab to look at the output
48  matlab  matlab
49    
50  for fld={'Eta','VICE','UICE','HEFF','AREA','Qsw','Qnet','FV','FU','EmPmR'}  for fld={'Eta','VICE','UICE','HEFF','AREA','Qsw','Qnet','FV','FU','EmPmR'}
51      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]);
52      cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause      cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause
53  end  end
54    
55    heff=permute(readbin(['HEFF.0000000020.data'],[32,6,32]),[1 3 2]);
56    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  # 32*32*6*15, 2-cpu cube sphere on the altices with efc
69    
 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  
   
70  cd MITgcm/verification/global_ocean.cs32x15/code  cd MITgcm/verification/global_ocean.cs32x15/code
71  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
72  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s12t_16x32/* .  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s12t_16x32/* .
73  \cp SIZE.h_mpi SIZE.h  \cp SIZE.h_mpi SIZE.h
 \cp CPP_EEOPTIONS.h_mpi CPP_EEOPTIONS.h  
 \rm mdsio_readfield.F  
74    
75  cd ..  cd ../build
76  mkdir build  \rm *
77  cd build  ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi
 ../code/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi  
78  make depend  make depend
79    \rm getcon.F
80  make  make
81    
82  cd ..  cd ..
83  mkdir run  mkdir run
84  cd run  cd run
85  \cp ../inp_thsice/* .  \rm *
86  \cp ../input/* .  \cp ../input/* .
87  \cp ../../../../MITgcm_contrib/high_res_cube/matlab-grid-converter/tile00* .  \cp ../../../../MITgcm_contrib/high_res_cube/matlab-grid-converter/tile00* .
88  \cp ../../../../MITgcm_contrib/high_res_cube/input/* .  \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
# Line 89  mpirun -np 2 ../build/mitgcmuv Line 92  mpirun -np 2 ../build/mitgcmuv
92  ##########################################################################  ##########################################################################
93  ##########################################################################  ##########################################################################
94  # 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
95    # requires input files in sudirectories run_template and ncep_rgau
96  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  
97    
98  cd MITgcm/verification/global_ocean.cs32x15/code  cd MITgcm/verification/global_ocean.cs32x15/code
99  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
100  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* .  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* .
101  \cp SIZE.h_54 SIZE.h  \cp SIZE.h_54 SIZE.h
102  \cp CPP_EEOPTIONS.h_mpi CPP_EEOPTIONS.h  \cp packages.conf_hr packages.conf
 \rm mdsio_readfield.F  
103    
104  cd ..  cd ../build
105  mkdir build  \rm *
106  cd build  ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi
 ../code/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi  
107  make depend  make depend
108    \rm getcon.F
109  make  make
110    
111  cd ..  cd ..
112  mkdir run  mkdir run
113  cd run  cd run
114  \cp ../inp_thsice/* .  \rm *
115  \cp ../input/* .  \cp ../input/* .
116  \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 .  
117  qsub job_orion  qsub job_orion
118    
119    
120  ##########################################################################  ##########################################################################
121  ##########################################################################  ##########################################################################
122  # 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
123    # requires input files in sudirectories run_template and ncep_rgau
124    # presently residing under orion:/tmp1/dmenem/cube
125    
126  setenv CVSROOT :pserver:cvsanon@mitgcm.org:/u/u0/gcmpack  cd MITgcm/verification/global_ocean.cs32x15/code
127  cvs login ( CVS password: cvsanon )  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
128  cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/README_ice  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s205t_85x85/* .
129  cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/code-mods  \cp SIZE.h_41 SIZE.h
130  cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/matlab-grid-converter  \cp packages.conf_hr packages.conf
131  cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/input  
132  cvs co -r checkpoint52d_pre MITgcm_contrib/high_res_cube/results  cd ../build
133  cvs co -r checkpoint52d_pre MITgcm_code  \rm *
134  cvs co -r checkpoint52d_pre MITgcm/verification/global_ocean.cs32x15  ../../../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  cd MITgcm/verification/global_ocean.cs32x15/code
155  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
156  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* .  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* .
157  \cp SIZE.h_216 SIZE.h  \cp SIZE.h_216 SIZE.h
158  \cp CPP_EEOPTIONS.h_mpi CPP_EEOPTIONS.h  \cp packages.conf_hr packages.conf
 \rm mdsio_readfield.F  
159    
160  cd ..  cd ../build
161  mkdir build  \rm *
162  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  
163  make depend  make depend
164    \rm getcon.F
165  make  make
166    
167  cd ..  cd ..
168  mkdir run  mkdir run
169  cd run  cd run
170  \cp ../inp_thsice/* .  \rm *
171  \cp ../input/* .  \cp ../input/* .
172  \cp ../../../../MITgcm_contrib/high_res_cube/input/* .  \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
173  \cp data_hr data  qsub job_altix
174  \cp data.pkg_hr data.pkg  
175  \ln -sf ../../../../run_template/tile00* .  
176  \ln -sf ../../../../run_template/LEV* .  ##########################################################################
177  \ln -sf ../../../../run_template/BATHY_18360x85_filled.bin .  ##########################################################################
178  qsub job_orion  # 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.9  
changed lines
  Added in v.1.31
  ViewVC Help
Powered by ViewVC 1.1.22