/[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.8 by dimitri, Fri Dec 12 10:20:05 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: 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 ../../../../MITgcm_contrib/high_res_cube/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
 \cp ../inp_thsice/* .  
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 42  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 ../../../../MITgcm_contrib/high_res_cube/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 *
74  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  
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  \cp ../inp_thsice/* .  \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    
87    
88  ##########################################################################  ##########################################################################
89  ##########################################################################  ##########################################################################
90  # 32*32*6*15, 2-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
98  \cp CPP_EEOPTIONS.h_mpi CPP_EEOPTIONS.h  \cp packages.conf_hr packages.conf
 \rm mdsio_readfield.F  
99    
100  cd ..  cd ../build
101  mkdir build  \rm *
102  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  
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  \cp ../inp_thsice/* .  \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, 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.8  
changed lines
  Added in v.1.32
  ViewVC Help
Powered by ViewVC 1.1.22