/[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.7 by dimitri, Mon May 25 15:25:08 2009 UTC revision 1.8 by dimitri, Tue May 26 00:25:12 2009 UTC
# Line 35  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 75  c Line 74  c
74  w  w
75  q  q
76  EOF  EOF
 \cp ../verification/global_ocean.cs32x15/code_alt/code.176t_8x4/* .  
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 .
# Line 89  cd ../run Line 89  cd ../run
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
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22