/[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.4 by dimitri, Sun Sep 28 13:44:33 2008 UTC revision 1.9 by dimitri, Thu Jun 18 18:29:51 2009 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/* .  
 ex - SIZE.h >> /dev/null <<EOF  
 /OLx =  
 c  
      &           OLx =   5,  
 .  
 /OLy =  
 c  
      &           OLy =   5,  
 .  
 w  
 q  
 EOF  
24  \rm CPP_EEOPTIONS.h  \rm CPP_EEOPTIONS.h
25  ../tools/genmake2  ../tools/genmake2
26  make depend  make depend
# Line 58  cd ../run Line 35  cd ../run
35    
36  emacs output.txt ../../MITgcm_contrib/high_res_cube/results/output.txt  emacs output.txt ../../MITgcm_contrib/high_res_cube/results/output.txt
37    
38  to use matlab to look at the output  to use matlab to look at the output:
39  matlab  matlab
   
40  for fld={'SIuice','SIvice','SIheff','SIarea','oceQsw', ...  for fld={'SIuice','SIvice','SIheff','SIarea','oceQsw', ...
41           'surForcT','oceTAUX','oceTAUY','surForcS'}           'surForcT','oceTAUX','oceTAUY','surForcS'}
42      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 79  subplot(224), mypcolor(vice(:,:,6)'); co Line 55  subplot(224), mypcolor(vice(:,:,6)'); co
55    
56    
57  ##########################################################################  ##########################################################################
58  # 32*32*6*15, 17-tile, 1-cpu cube sphere experiment  # 32*32*6*15, 176-tile, 1-cpu cube sphere experiment
59  # needs fixing; does not work as of checkpoint58x_post  # example that has 16 blank tiles
 # see discussion on MITgcm-devel  
60    
61  cd ../build  cd ../build
62  \rm *  \rm *
63  \cp ../../MITgcm_contrib/high_res_cube/code-mods/* .  \cp ../../MITgcm_contrib/high_res_cube/code-mods/* .
64  \cp ../verification/global_ocean.cs32x15/code_alt/code.176t_8x4/* .  \cp ../utils/exch2/code-mods/s176t_8x4/SIZE.h .
65    ex - SIZE.h >> /dev/null <<EOF
66    /OLx =
67    c
68         &           OLx =   5,
69    .
70    /OLy =
71    c
72         &           OLy =   5,
73    .
74    w
75    q
76    EOF
77  \rm CPP_EEOPTIONS.h  \rm CPP_EEOPTIONS.h
78  ../tools/genmake2  ../tools/genmake2
79  make depend  make depend
80  make -j  make -j
81  cd ../run  cd ..
82  \rm *  mv run run_orig
83    mkdir run
84    cd run
85  \cp ../verification/tutorial_held_suarez_cs/input/grid* .  \cp ../verification/tutorial_held_suarez_cs/input/grid* .
86  \cp ../verification/global_ocean.cs32x15/input/lev_*_cs_15k.bin .  \cp ../verification/global_ocean.cs32x15/input/lev_*_cs_15k.bin .
87  \cp ../verification/global_ocean.cs32x15/input/bathy_Hmin50.bin .  \cp ../verification/global_ocean.cs32x15/input/bathy_Hmin50.bin .
88    \cp ../utils/exch2/code-mods/data.exch2.16_blk data.exch2
89  \cp ../../MITgcm_contrib/high_res_cube/input/* .  \cp ../../MITgcm_contrib/high_res_cube/input/* .
90  ../build/mitgcmuv >& output.txt  ../build/mitgcmuv >& output.txt
91    
92    to use matlab to compare to previous solution:
93    matlab
94    for f={'ETAN','oceTAUY','SIhsnow','surForcS','KPPhbl','PHIBOT', ...
95           'SIuice','surForcT','MXLDEPTH','SIarea','SIvice','UVEL_k1', ...
96           'oceQsw','SIheff','SSS','VVEL_k1','oceTAUX','SIhsalt','SST'}
97     t1=readbin(['../run_orig/' f{1} '.0000000072.data'],[32*6 32]);
98     t2=readbin([f{1} '.0000000072.data'],[32*6 32]);
99     clf, subplot(311), mypcolor(t1'); thincolorbar, title(f{1})
100     subplot(312), mypcolor(t2'); thincolorbar
101     subplot(313), mypcolor(t2'-t1'); thincolorbar, pause
102    end
103    for f={'YG','YC','XG','XC','RAZ','RAW','RAS','RAC','hFacW','hFacS', ...
104         'hFacC','Depth','DYG','DYC','DXG','DXC','AngleSN','AngleCS'}
105     t1=readbin(['../run_orig/' f{1} '.data'],[32*6 32]);
106     t2=readbin([f{1} '.data'],[32*6 32]);
107     clf, subplot(311), mypcolor(t1'); thincolorbar, title(f{1})
108     subplot(312), mypcolor(t2'); thincolorbar
109     subplot(313), mypcolor(t2'-t1'); thincolorbar, pause
110    end
111    
112    
113  ##########################################################################  ##########################################################################
114  # 32*32*6*15, 2-cpu cube sphere on the NAS altices with ifort  # 32*32*6*15, 2-cpu cube sphere on the NAS altices with ifort
# Line 108  mkdir build run Line 118  mkdir build run
118  cd build  cd build
119  \rm *  \rm *
120  \cp ../../MITgcm_contrib/high_res_cube/code-mods/* .  \cp ../../MITgcm_contrib/high_res_cube/code-mods/* .
121  \cp ../utils/exch2/code-mods/s12t_16x32/* .  \mv SIZE.h_mpi SIZE.h
 \cp SIZE.h_mpi SIZE.h  
122  ../tools/genmake2 -of ../../MITgcm_contrib/high_res_cube/code-mods/linux_ia64_ifort+mpi_altix_nas  ../tools/genmake2 -of ../../MITgcm_contrib/high_res_cube/code-mods/linux_ia64_ifort+mpi_altix_nas
123  make depend  make depend
124  make -j  make -j 8
125  cd ../run  cd ../run
126  \rm *  \rm *
127  \cp ../verification/tutorial_held_suarez_cs/input/grid* .  \cp ../verification/tutorial_held_suarez_cs/input/grid* .
# Line 130  mkdir build run Line 139  mkdir build run
139  cd build  cd build
140  \rm *  \rm *
141  \cp ../../MITgcm_contrib/high_res_cube/code-mods/* .  \cp ../../MITgcm_contrib/high_res_cube/code-mods/* .
142  \cp ../utils/exch2/code-mods/s12t_16x32/* .  \mv SIZE.h_mpi SIZE.h
 \cp SIZE.h_mpi SIZE.h  
143  ../tools/genmake2 -of ../tools/build_options/linux_ia64_ifort+mpi_altix_jpl  ../tools/genmake2 -of ../tools/build_options/linux_ia64_ifort+mpi_altix_jpl
144  make depend  make depend
145  make -j  make -j
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22