/[MITgcm]/MITgcm_contrib/high_res_cube/README_ice
ViewVC logotype

Diff of /MITgcm_contrib/high_res_cube/README_ice

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.36 by dimitri, Sun Mar 28 09:28:21 2004 UTC revision 1.37 by dimitri, Tue May 11 04:03:55 2004 UTC
# Line 24  cvs co -r hrcube5 MITgcm/verification/gl Line 24  cvs co -r hrcube5 MITgcm/verification/gl
24  cd MITgcm/verification/global_ocean.cs32x15/code  cd MITgcm/verification/global_ocean.cs32x15/code
25  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
26  \cp ../../../utils/exch2/code-mods/s12t_16x32/* .  \cp ../../../utils/exch2/code-mods/s12t_16x32/* .
   
27  cd ../build  cd ../build
28    \rm *
29  ../../../tools/genmake2 -mods=../code  ../../../tools/genmake2 -mods=../code
30  make depend  make depend
31  make  make
32    cd ..
33    mkdir run
34    cd run
35    \rm *
36    \cp ../input/* .
37    \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
38    ../build/mitgcmuv >& output.txt
39    
40    comparison output is in:
41    ../../../../MITgcm_contrib/high_res_cube/results/output.txt
42    
43    to use matlab to look at the output
44    matlab
45    
46    for fld={'Eta','VICE','UICE','HEFF','AREA','Qsw','Qnet','FV','FU','EmPmR'}
47        tmp=permute(readbin([fld{1} '.0000000020.data'],[32,6,32]),[1 3 2]);
48        cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause
49    end
50    
51    heff=permute(readbin(['HEFF.0000000020.data'],[32,6,32]),[1 3 2]);
52    mask=0*heff; mask(find(heff>.001))=1;
53    uice=mask.*permute(readbin(['UICE.0000000020.data'],[32,6,32]),[1 3 2]);
54    vice=mask.*permute(readbin(['VICE.0000000020.data'],[32,6,32]),[1 3 2]);
55    clf
56    subplot(221), mypcolor(uice(:,:,3)'), colorbar, title('uice, tile 3')
57    subplot(222), mypcolor(uice(:,:,6)'), colorbar, title('uice, tile 6')
58    subplot(223), mypcolor(vice(:,:,3)'), colorbar, title('vice, tile 3')
59    subplot(224), mypcolor(vice(:,:,6)'), colorbar, title('vice, tile 6')
60    
61    
62    ##########################################################################
63    ##########################################################################
64    # 32*32*6*15, 17-tile, 1-cpu cube sphere experiment
65    
66    cd MITgcm/verification/global_ocean.cs32x15/code
67    \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
68    \cp ../code_alt/code.176t_8x4/* .
69    cd ../build
70    \rm *
71    ../../../tools/genmake2 -mods=../code
72    make depend
73    make
74  cd ..  cd ..
75  mkdir run  mkdir run
76  cd run  cd run
77    \rm *
78  \cp ../input/* .  \cp ../input/* .
79  \cp ../../../../MITgcm_contrib/high_res_cube/input/* .  \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
80  ../build/mitgcmuv >& output.txt  ../build/mitgcmuv >& output.txt
Legend:
Removed from v.1.36  
changed lines
  Added in v.1.37
  ViewVC Help
Powered by ViewVC 1.1.22