/[MITgcm]/MITgcm_contrib/high_res_cube/README.cs32
ViewVC logotype

Diff of /MITgcm_contrib/high_res_cube/README.cs32

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.2 by dimitri, Sat Apr 7 23:21:41 2007 UTC revision 1.11 by dimitri, Wed Mar 21 23:18:55 2012 UTC
# Line 5  a preliminary sanity check for the cs510 Line 5  a preliminary sanity check for the cs510
5  ##########################################################################  ##########################################################################
6  # getting the code from anonymous CVS server  # getting the code from anonymous CVS server
7    
8    bash or sh shell:  cvs co -P MITgcm_contrib/high_res_cube/README.cs32
9      $ export CVSROOT=':pserver:cvsanon@mitgcm.org:/u/gcmpack'  cvs co -P MITgcm_contrib/high_res_cube/code-mods
10      $ cvs login  cvs co -P MITgcm_contrib/high_res_cube/input
11        ( enter the CVS password: "cvsanon" )  cvs co -P MITgcm_contrib/high_res_cube/results
12    cvs co -P MITgcm_code
13    tcsh or csh shell:  cvs co -P MITgcm/verification/global_ocean.cs32x15/input
14      $ setenv CVSROOT ':pserver:cvsanon@mitgcm.org:/u/gcmpack'  cvs co -P MITgcm/verification/tutorial_held_suarez_cs/input
     $ cvs login  
       ( enter the CVS password: "cvsanon" )  
   
 cvs co MITgcm_contrib/high_res_cube/README.cs32  
 cvs co MITgcm_contrib/high_res_cube/code-mods  
 cvs co MITgcm_contrib/high_res_cube/input  
 cvs co MITgcm_contrib/high_res_cube/results  
 cvs co MITgcm_code  
 cvs co MITgcm/verification/global_ocean.cs32x15  
 cvs co MITgcm/verification/tutorial_held_suarez_cs/input  
15    
16  ##########################################################################  ##########################################################################
17  # 32*32*6*15, 1-cpu cube sphere  # 32*32*6*15, 1-cpu cube sphere
# Line 31  mkdir build run Line 21  mkdir build run
21  cd build  cd build
22  \rm *  \rm *
23  \cp ../../MITgcm_contrib/high_res_cube/code-mods/* .  \cp ../../MITgcm_contrib/high_res_cube/code-mods/* .
 \cp ../utils/exch2/code-mods/s12t_16x32/* .  
24  \rm CPP_EEOPTIONS.h  \rm CPP_EEOPTIONS.h
25  ../tools/genmake2  ../tools/genmake2
26  make depend  make depend
27  make -j  make -j 16
28  cd ../run  cd ../run
29  \rm *  \rm *
30  \cp ../verification/tutorial_held_suarez_cs/input/grid* .  \cp ../verification/tutorial_held_suarez_cs/input/grid* .
31  \cp ../verification/global_ocean.cs32x15/input/lev_*_cs_15k.bin .  \cp ../verification/global_ocean.cs32x15/input/lev_*_cs_15k.bin .
32  \cp ../verification/global_ocean.cs32x15/input/bathy_Hmin50.bin .  \cp ../verification/global_ocean.cs32x15/input/bathy_Hmin50.bin .
33  \cp ../../MITgcm_contrib/high_res_cube/input/* .  \cp ../../MITgcm_contrib/high_res_cube/input/* .
34    \cp data.pkg_lr data.pkg
35    \cp data.diagnostics_lr data.diagnostics
36  ../build/mitgcmuv >& output.txt  ../build/mitgcmuv >& output.txt
37    
38  comparison output is in:  emacs output.txt ../../MITgcm_contrib/high_res_cube/results/output.txt
 ../../MITgcm_contrib/high_res_cube/results/output.txt  
39    
40  to use matlab to look at the output  to use matlab to look at the output:
41  matlab  matlab
   
42  for fld={'SIuice','SIvice','SIheff','SIarea','oceQsw', ...  for fld={'SIuice','SIvice','SIheff','SIarea','oceQsw', ...
43           'surForcT','oceTAUX','oceTAUY','surForcS'}           'surForcT','oceTAUX','oceTAUY','surForcS'}
44      tmp=permute(readbin([fld{1} '.0000000072.data'],[32,6,32]),[1 3 2]);      tmp=permute(readbin([fld{1} '.0000000072.data'],[32,6,32]),[1 3 2]);
# Line 68  subplot(224), mypcolor(vice(:,:,6)'); co Line 57  subplot(224), mypcolor(vice(:,:,6)'); co
57    
58    
59  ##########################################################################  ##########################################################################
60  # 32*32*6*15, 17-tile, 1-cpu cube sphere experiment  # 32*32*6*15, 176-tile, 1-cpu cube sphere experiment
61  # needs fixing; does not work as of checkpoint58x_post  # example that has 16 blank tiles
 # see discussion on MITgcm-devel  
62    
63  cd ../build  cd ../build
64  \rm *  \rm *
65  \cp ../../MITgcm_contrib/high_res_cube/code-mods/* .  \cp ../../MITgcm_contrib/high_res_cube/code-mods/* .
66  \cp ../verification/global_ocean.cs32x15/code_alt/code.176t_8x4/* .  \cp ../utils/exch2/code-mods/s176t_8x4/SIZE.h .
67    ex - SIZE.h >> /dev/null <<EOF
68    /OLx =
69    c
70         &           OLx =   5,
71    .
72    /OLy =
73    c
74         &           OLy =   5,
75    .
76    w
77    q
78    EOF
79  \rm CPP_EEOPTIONS.h  \rm CPP_EEOPTIONS.h
80  ../tools/genmake2  ../tools/genmake2
81  make depend  make depend
82  make -j  make -j
83    cd ..
84    mv run run_orig
85    mkdir run
86    cd run
87    \cp ../verification/tutorial_held_suarez_cs/input/grid* .
88    \cp ../verification/global_ocean.cs32x15/input/lev_*_cs_15k.bin .
89    \cp ../verification/global_ocean.cs32x15/input/bathy_Hmin50.bin .
90    \cp ../utils/exch2/code-mods/data.exch2.16_blk data.exch2
91    \cp ../../MITgcm_contrib/high_res_cube/input/* .
92    ../build/mitgcmuv >& output.txt
93    
94    to use matlab to compare to previous solution:
95    matlab
96    for f={'ETAN','oceTAUY','SIhsnow','surForcS','KPPhbl','PHIBOT', ...
97           'SIuice','surForcT','MXLDEPTH','SIarea','SIvice','UVEL_k1', ...
98           'oceQsw','SIheff','SSS','VVEL_k1','oceTAUX','SIhsalt','SST'}
99     t1=readbin(['../run_orig/' f{1} '.0000000072.data'],[32*6 32]);
100     t2=readbin([f{1} '.0000000072.data'],[32*6 32]);
101     clf, subplot(311), mypcolor(t1'); thincolorbar, title(f{1})
102     subplot(312), mypcolor(t2'); thincolorbar
103     subplot(313), mypcolor(t2'-t1'); thincolorbar, pause
104    end
105    for f={'YG','YC','XG','XC','RAZ','RAW','RAS','RAC','hFacW','hFacS', ...
106         'hFacC','Depth','DYG','DYC','DXG','DXC','AngleSN','AngleCS'}
107     t1=readbin(['../run_orig/' f{1} '.data'],[32*6 32]);
108     t2=readbin([f{1} '.data'],[32*6 32]);
109     clf, subplot(311), mypcolor(t1'); thincolorbar, title(f{1})
110     subplot(312), mypcolor(t2'); thincolorbar
111     subplot(313), mypcolor(t2'-t1'); thincolorbar, pause
112    end
113    
114    
115    ##########################################################################
116    # 32*32*6*15, 2-cpu cube sphere on the NAS altices with ifort
117    
118    cd MITgcm
119    mkdir build run
120    cd build
121    \rm *
122    \cp ../../MITgcm_contrib/high_res_cube/code-mods/* .
123    \mv SIZE.h_mpi SIZE.h
124    ../tools/genmake2 -of ../../MITgcm_contrib/high_res_cube/code-mods/linux_ia64_ifort+mpi_altix_nas
125    make depend
126    make -j 8
127  cd ../run  cd ../run
128  \rm *  \rm *
129  \cp ../verification/tutorial_held_suarez_cs/input/grid* .  \cp ../verification/tutorial_held_suarez_cs/input/grid* .
130  \cp ../verification/global_ocean.cs32x15/input/lev_*_cs_15k.bin .  \cp ../verification/global_ocean.cs32x15/input/lev_*_cs_15k.bin .
131  \cp ../verification/global_ocean.cs32x15/input/bathy_Hmin50.bin .  \cp ../verification/global_ocean.cs32x15/input/bathy_Hmin50.bin .
132  \cp ../../MITgcm_contrib/high_res_cube/input/* .  \cp ../../MITgcm_contrib/high_res_cube/input/* .
133  ../build/mitgcmuv >& output.txt  mpirun -np 2 ../build/mitgcmuv
134    
135    
136  ##########################################################################  ##########################################################################
137  # 32*32*6*15, 2-cpu cube sphere on the altices with ifort  # 32*32*6*15, 2-cpu cube sphere on the JPL altices with ifort
138    
139  cd MITgcm  cd MITgcm
140  mkdir build run  mkdir build run
141  cd build  cd build
142  \rm *  \rm *
143  \cp ../../MITgcm_contrib/high_res_cube/code-mods/* .  \cp ../../MITgcm_contrib/high_res_cube/code-mods/* .
144  \cp ../utils/exch2/code-mods/s12t_16x32/* .  \mv SIZE.h_mpi SIZE.h
145  \cp SIZE.h_mpi SIZE.h  ../tools/genmake2 -of ../tools/build_options/linux_ia64_ifort+mpi_altix_jpl
 ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_ifort+mpi_altix_nas  
146  make depend  make depend
147  make -j  make -j
148  cd ../run  cd ../run
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22