Follow instructions for testing edge problems and other tiling issues with 32x32x6 cube sphere configuration. Baseline 6-tile integration =========================== cd MITgcm/verification/global_ocean.cs32x15/build rm * ../../../tools/genmake2 -mods=../code make depend rm SIZE.h cp ../../../../MITgcm_contrib/edge_problem/SIZE.h . make cd .. mkdir run6 cd run6 rm * cp ../input/* . cp ../../../../MITgcm_contrib/edge_problem/data.base data cp ../../../../MITgcm_contrib/edge_problem/tile* . cp ../build/mitgcmuv . mitgcmuv > output.txt Baseline 12-tile integration ============================ cd ../build rm * ../../../tools/genmake2 -mods="../code_alt/code.12t_32x16 ../code" make depend rm SIZE.h cp ../../../../MITgcm_contrib/edge_problem/SIZE.h_12 SIZE.h make cd .. mkdir run12 cd run12 rm * cp ../input/* . cp ../../../../MITgcm_contrib/edge_problem/data.base data cp ../../../../MITgcm_contrib/edge_problem/tile* . cp ../build/mitgcmuv . mitgcmuv > output.txt cd .. Use matlab to compare Eta for 6-tile and 12-tile integrations ============================================================= matlab -nojvm -nosplash f1='run6/Eta.'; f2='run12/Eta.'; tm=72003; ix=.5:191.5; iy=.5:31.5; tmp1=readbin([f1 myint2str(tm,10) '.data'],[32*6,32]); tmp2=readbin([f2 myint2str(tm,10) '.data'],[32*6,32]); clf, subplot(211) image(ix,iy,tmp1','CDataMapping','scaled'); colorbar, grid set(gca,'ydir','normal','xtick',0:32:192,'ytick',0:16:32) title('Eta, time step 3, 6-tile integration') subplot(212) image(ix,iy,tmp2'-tmp1','CDataMapping','scaled'); colorbar, grid set(gca,'ydir','normal','xtick',0:32:192,'ytick',0:16:32) title('Difference between 12-tile and 6-tile integrations') quit Now repeat with visca4grid ========================== cd run6 cp ../../../../MITgcm_contrib/edge_problem/data.visca4grid data mitgcmuv > output.txt cd ../run12 cp ../../../../MITgcm_contrib/edge_problem/data.visca4grid data mitgcmuv > output.txt cd .. Now repeat with advscheme=33 ============================ cd run6 cp ../../../../MITgcm_contrib/edge_problem/data.tempadv data mitgcmuv > output.txt cd ../run12 cp ../../../../MITgcm_contrib/edge_problem/data.tempadv data mitgcmuv > output.txt cd ..