/[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.75 by dimitri, Fri Nov 3 23:17:42 2006 UTC revision 1.80 by dimitri, Wed Dec 6 01:07:44 2006 UTC
# Line 3  Instructions for setting up a cube-spher Line 3  Instructions for setting up a cube-spher
3  =================  =================
4  UPDATES  UPDATES
5    
6  November 3, 2006: 216-CPU code for cube43  December 5, 2006: 216-CPU code for cube46
7  Martin's growth.F routine and Patrick's pkg/seaice diagnostics  Changed diagnostics to accomodate JMCs latest and some
8    changes to code-mods/CPP_EEOPTIONS.h and W2_OPTIONS.h
9    to speed up code.
10    
11    November 7, 2006: 216-CPU code for cube43
12    Bug fixes in growth.F routine and new set of diagnostics
13    
14  October 15, 2006: 216-CPU code for cube40  October 15, 2006: 216-CPU code for cube40
15    
# Line 24  for online T/S profile diagnostics all t Line 29  for online T/S profile diagnostics all t
29  see pkg/ecco/ecco_cost_init_fixed.F  see pkg/ecco/ecco_cost_init_fixed.F
30      pkg/ecco/cost_averagesfields.F      pkg/ecco/cost_averagesfields.F
31    
32    ##########################################################################
33    ##########################################################################
34    # Get forcing files, etc.
35    # Example below is for cube46 on /nobackup2a/menemenl
36    
37    cd /nobackup2a/menemenl
38    mkdir cube46
39    cd cube46
40    ln -sf ../run_template .
41    ln -sf ../ncep_rgau .
42    
43    
44  ##########################################################################  ##########################################################################
45  ##########################################################################  ##########################################################################
46  # getting the code from anonymous CVS server  # getting the code from anonymous CVS server
 #  
47    
48    bash or sh shell:    bash or sh shell:
49      $ export CVSROOT=':pserver:cvsanon@mitgcm.org:/u/gcmpack'      $ export CVSROOT=':pserver:cvsanon@mitgcm.org:/u/gcmpack'
# Line 52  cvs co MITgcm/verification/tutorial_held Line 67  cvs co MITgcm/verification/tutorial_held
67  ##########################################################################  ##########################################################################
68  ##########################################################################  ##########################################################################
69  # 510*510*6*50, 216-cpu cube sphere on altix with ifort and pkg/diagnostics  # 510*510*6*50, 216-cpu cube sphere on altix with ifort and pkg/diagnostics
 # requires input files in sudirectories run_template and ncep_rgau  
 # these can be obained from lou:/u/menemenl/cube  
70    
71  cd MITgcm/verification/global_ocean.cs32x15/code  cd MITgcm/verification/global_ocean.cs32x15/code
72  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
# Line 63  cd MITgcm/verification/global_ocean.cs32 Line 76  cd MITgcm/verification/global_ocean.cs32
76    
77  cd ../build  cd ../build
78  \rm *  \rm *
79  ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_ifort+mpi_altix_nas  ../../../tools/genmake2 -mods=../code -of ../code/linux_ia64_ifort+mpi_altix_nas
80  make depend  make depend
81  make -j  make -j
82    
# Line 109  comparison output is in: Line 122  comparison output is in:
122  to use matlab to look at the output  to use matlab to look at the output
123  matlab  matlab
124    
125  for fld={'VICEtave','UICEtave','HEFFtave','AREAtave','QSWtave', ...  for fld={'SIuice','SIvice','SIheff','SIarea','oceQsw', ...
126           'QNETtave','FVtave','FUtave','EmPmRtave'}           'surForcT','oceTAUX','oceTAUY','surForcS'}
127      tmp=permute(readbin([fld{1} '.0000000020.data'],[32,6,32]),[1 3 2]);      tmp=permute(readbin([fld{1} '.0000000072.data'],[32,6,32]),[1 3 2]);
128      cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause      cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause
129  end  end
130    
131  heff=permute(readbin(['HEFFtave.0000000020.data'],[32,6,32]),[1 3 2]);  heff=permute(readbin(['SIheff.0000000072.data'],[32,6,32]),[1 3 2]);
132  mask=0*heff; mask(find(heff>.001))=1;  mask=0*heff; mask(find(heff>.001))=1;
133  uice=mask.*permute(readbin(['UICEtave.0000000020.data'],[32,6,32]),[1 3 2]);  uice=mask.*permute(readbin(['SIuice.0000000072.data'],[32,6,32]),[1 3 2]);
134  vice=mask.*permute(readbin(['VICEtave.0000000020.data'],[32,6,32]),[1 3 2]);  vice=mask.*permute(readbin(['SIvice.0000000072.data'],[32,6,32]),[1 3 2]);
135  clf  clf
136  subplot(221), mypcolor(uice(:,:,3)'); colorbar, title('uice, tile 3')  subplot(221), mypcolor(uice(:,:,3)'); colorbar, title('uice, tile 3')
137  subplot(222), mypcolor(uice(:,:,6)'); colorbar, title('uice, tile 6')  subplot(222), mypcolor(uice(:,:,6)'); colorbar, title('uice, tile 6')
Legend:
Removed from v.1.75  
changed lines
  Added in v.1.80
  ViewVC Help
Powered by ViewVC 1.1.22