--- MITgcm_contrib/high_res_cube/README.cs32	2009/05/25 15:25:08	1.7
+++ MITgcm_contrib/high_res_cube/README.cs32	2012/03/21 23:18:55	1.11
@@ -24,20 +24,21 @@
 \rm CPP_EEOPTIONS.h
 ../tools/genmake2
 make depend
-make -j
+make -j 16
 cd ../run
 \rm *
 \cp ../verification/tutorial_held_suarez_cs/input/grid* .
 \cp ../verification/global_ocean.cs32x15/input/lev_*_cs_15k.bin .
 \cp ../verification/global_ocean.cs32x15/input/bathy_Hmin50.bin .
 \cp ../../MITgcm_contrib/high_res_cube/input/* .
+\cp data.pkg_lr data.pkg
+\cp data.diagnostics_lr data.diagnostics
 ../build/mitgcmuv >& output.txt
 
 emacs output.txt ../../MITgcm_contrib/high_res_cube/results/output.txt
 
-to use matlab to look at the output
+to use matlab to look at the output:
 matlab
-
 for fld={'SIuice','SIvice','SIheff','SIarea','oceQsw', ...
          'surForcT','oceTAUX','oceTAUY','surForcS'}
     tmp=permute(readbin([fld{1} '.0000000072.data'],[32,6,32]),[1 3 2]);
@@ -75,13 +76,14 @@
 w
 q
 EOF
-\cp ../verification/global_ocean.cs32x15/code_alt/code.176t_8x4/* .
 \rm CPP_EEOPTIONS.h
 ../tools/genmake2
 make depend
 make -j
-cd ../run
-\rm *
+cd ..
+mv run run_orig
+mkdir run
+cd run
 \cp ../verification/tutorial_held_suarez_cs/input/grid* .
 \cp ../verification/global_ocean.cs32x15/input/lev_*_cs_15k.bin .
 \cp ../verification/global_ocean.cs32x15/input/bathy_Hmin50.bin .
@@ -89,6 +91,26 @@
 \cp ../../MITgcm_contrib/high_res_cube/input/* .
 ../build/mitgcmuv >& output.txt
 
+to use matlab to compare to previous solution:
+matlab
+for f={'ETAN','oceTAUY','SIhsnow','surForcS','KPPhbl','PHIBOT', ...
+       'SIuice','surForcT','MXLDEPTH','SIarea','SIvice','UVEL_k1', ...
+       'oceQsw','SIheff','SSS','VVEL_k1','oceTAUX','SIhsalt','SST'}
+ t1=readbin(['../run_orig/' f{1} '.0000000072.data'],[32*6 32]);
+ t2=readbin([f{1} '.0000000072.data'],[32*6 32]);
+ clf, subplot(311), mypcolor(t1'); thincolorbar, title(f{1})
+ subplot(312), mypcolor(t2'); thincolorbar
+ subplot(313), mypcolor(t2'-t1'); thincolorbar, pause
+end
+for f={'YG','YC','XG','XC','RAZ','RAW','RAS','RAC','hFacW','hFacS', ...
+     'hFacC','Depth','DYG','DYC','DXG','DXC','AngleSN','AngleCS'}
+ t1=readbin(['../run_orig/' f{1} '.data'],[32*6 32]);
+ t2=readbin([f{1} '.data'],[32*6 32]);
+ clf, subplot(311), mypcolor(t1'); thincolorbar, title(f{1})
+ subplot(312), mypcolor(t2'); thincolorbar
+ subplot(313), mypcolor(t2'-t1'); thincolorbar, pause
+end
+
 
 ##########################################################################
 # 32*32*6*15, 2-cpu cube sphere on the NAS altices with ifort
@@ -101,7 +123,7 @@
 \mv SIZE.h_mpi SIZE.h
 ../tools/genmake2 -of ../../MITgcm_contrib/high_res_cube/code-mods/linux_ia64_ifort+mpi_altix_nas
 make depend
-make -j
+make -j 8
 cd ../run
 \rm *
 \cp ../verification/tutorial_held_suarez_cs/input/grid* .