Diff of /MITgcm_contrib/high_res_cube/README.cs32

Parent Directory | Revision Log | Revision Graph | Patch

--- MITgcm_contrib/high_res_cube/README.cs32	2009/05/25 15:25:08	1.7
+++ MITgcm_contrib/high_res_cube/README.cs32	2009/05/26 00:25:12	1.8
@@ -35,9 +35,8 @@
 
 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 +74,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 +89,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