--- MITgcm_contrib/high_res_cube/README.cs32 2007/04/07 21:56:47 1.1 +++ MITgcm_contrib/high_res_cube/README.cs32 2009/06/18 18:29:51 1.9 @@ -5,51 +5,38 @@ ########################################################################## # getting the code from anonymous CVS server - bash or sh shell: - $ export CVSROOT=':pserver:cvsanon@mitgcm.org:/u/gcmpack' - $ cvs login - ( enter the CVS password: "cvsanon" ) - - tcsh or csh shell: - $ setenv CVSROOT ':pserver:cvsanon@mitgcm.org:/u/gcmpack' - $ 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 +cvs co -P MITgcm_contrib/high_res_cube/README.cs32 +cvs co -P MITgcm_contrib/high_res_cube/code-mods +cvs co -P MITgcm_contrib/high_res_cube/input +cvs co -P MITgcm_contrib/high_res_cube/results +cvs co -P MITgcm_code +cvs co -P MITgcm/verification/global_ocean.cs32x15/input +cvs co -P MITgcm/verification/tutorial_held_suarez_cs/input ########################################################################## # 32*32*6*15, 1-cpu cube sphere -cd MITgcm/verification/global_ocean.cs32x15/code -\cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* . -\cp ../../../utils/exch2/code-mods/s12t_16x32/* . -\rm CPP_EEOPTIONS.h -cd ../build +cd MITgcm +mkdir build run +cd build \rm * -../../../tools/genmake2 -mods=../code +\cp ../../MITgcm_contrib/high_res_cube/code-mods/* . +\rm CPP_EEOPTIONS.h +../tools/genmake2 make depend -make -cd .. -mkdir run -cd run +make -j +cd ../run \rm * -\cp ../../tutorial_held_suarez_cs/input/grid* . -\cp ../input/* . -\cp ../../../../MITgcm_contrib/high_res_cube/input/* . +\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/* . ../build/mitgcmuv >& output.txt -comparison output is in: -../../../../MITgcm_contrib/high_res_cube/results/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]); @@ -66,45 +53,100 @@ subplot(223), mypcolor(vice(:,:,3)'); colorbar, title('vice, tile 3') subplot(224), mypcolor(vice(:,:,6)'); colorbar, title('vice, tile 6') + ########################################################################## -# 32*32*6*15, 17-tile, 1-cpu cube sphere experiment +# 32*32*6*15, 176-tile, 1-cpu cube sphere experiment +# example that has 16 blank tiles -cd MITgcm/verification/global_ocean.cs32x15/code -\cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* . -\cp ../code_alt/code.176t_8x4/* . cd ../build \rm * -../../../tools/genmake2 -mods=../code +\cp ../../MITgcm_contrib/high_res_cube/code-mods/* . +\cp ../utils/exch2/code-mods/s176t_8x4/SIZE.h . +ex - SIZE.h >> /dev/null <& output.txt -########################################################################## -# 32*32*6*15, 2-cpu cube sphere on the altices with ifort +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 -cd MITgcm/verification/global_ocean.cs32x15/code -\cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* . -\cp ../../../utils/exch2/code-mods/s12t_16x32/* . -\cp SIZE.h_mpi SIZE.h -cd ../build +########################################################################## +# 32*32*6*15, 2-cpu cube sphere on the NAS altices with ifort + +cd MITgcm +mkdir build run +cd build +\rm * +\cp ../../MITgcm_contrib/high_res_cube/code-mods/* . +\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 8 +cd ../run \rm * -../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_ifort+mpi_altix_nas +\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/* . +mpirun -np 2 ../build/mitgcmuv + + +########################################################################## +# 32*32*6*15, 2-cpu cube sphere on the JPL altices with ifort + +cd MITgcm +mkdir build run +cd build +\rm * +\cp ../../MITgcm_contrib/high_res_cube/code-mods/* . +\mv SIZE.h_mpi SIZE.h +../tools/genmake2 -of ../tools/build_options/linux_ia64_ifort+mpi_altix_jpl make depend make -j - -cd .. -mkdir run -cd run +cd ../run \rm * -\cp ../../tutorial_held_suarez_cs/input/grid* . -\cp ../input/* . -\cp ../../../../MITgcm_contrib/high_res_cube/input/* . +\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/* . mpirun -np 2 ../build/mitgcmuv