/[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.76 by dimitri, Tue Nov 7 15:24:04 2006 UTC revision 1.91 by dimitri, Sat Apr 7 21:52:43 2007 UTC
# Line 1  Line 1 
1    !!!! Instructions in this file are obsolete and will eventually
2    !!!! be removed.  This file has been replaced by
3    !!!! MITgcm_contrib/high_res_cube/README.cs510
4    !!!! MITgcm_contrib/high_res_cube/README.cs32
5    
6    =================
7    
8  Instructions for setting up a cube-sphere integration with sea-ice.  Instructions for setting up a cube-sphere integration with sea-ice.
9    
10    Forcing fields and other input files needed for the CS510
11    configuration can be obtained from
12    http://ecco2.org/products/model_setups/cs510.html
13    
14  =================  =================
15    
16  UPDATES  UPDATES
17    
18    February 6, 2006: cube56, same as cube 50 with SEAICEuseFlooding
19    
20    February 1, 2006: back to cube50 default (bottomDragQuadratic = 0.002)
21                      removed UVEL and VVEL from diagnostics
22    
23    January 5, 2006: cube51, like cube50 but with bottomDragQuadratic = 0.0012
24    
25    January 5, 2006: cube50, like cube48 but with snow advection
26    
27    December 30, 2006: updates for cube49, the first
28    "blend" integration.
29    
30    December 11, 2006: preparing for cube47, same as cube43
31    but with Arctic monthly river runoff climatology.
32    
33    December 5, 2006: 216-CPU code for cube46
34    Changed diagnostics to accomodate JMCs latest and some
35    changes to code-mods/CPP_EEOPTIONS.h and W2_OPTIONS.h
36    to speed up code.
37    
38  November 7, 2006: 216-CPU code for cube43  November 7, 2006: 216-CPU code for cube43
39  Bug fixes in growth.F routine and new set of diagnostics  Bug fixes in growth.F routine and new set of diagnostics
40    
# Line 24  for online T/S profile diagnostics all t Line 56  for online T/S profile diagnostics all t
56  see pkg/ecco/ecco_cost_init_fixed.F  see pkg/ecco/ecco_cost_init_fixed.F
57      pkg/ecco/cost_averagesfields.F      pkg/ecco/cost_averagesfields.F
58    
59    ##########################################################################
60    ##########################################################################
61    # Get forcing files, etc.
62    # Example below is for cube46 on /nobackup2a/menemenl
63    
64    cd /nobackup2a/menemenl
65    mkdir cube47
66    cd cube47
67    ln -sf ../run_template .
68    ln -sf ../ncep_rgau .
69    
70    
71  ##########################################################################  ##########################################################################
72  ##########################################################################  ##########################################################################
73  # getting the code from anonymous CVS server  # getting the code from anonymous CVS server
 #  
74    
75    bash or sh shell:    bash or sh shell:
76      $ 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 94  cvs co MITgcm/verification/tutorial_held
94  ##########################################################################  ##########################################################################
95  ##########################################################################  ##########################################################################
96  # 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  
97    
98  cd MITgcm/verification/global_ocean.cs32x15/code  cd MITgcm/verification/global_ocean.cs32x15/code
99  \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 103  cd MITgcm/verification/global_ocean.cs32
103    
104  cd ../build  cd ../build
105  \rm *  \rm *
106  ../../../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
107  make depend  make depend
108  make -j  make -j
109    
# Line 76  cd run Line 116  cd run
116  ln -sf ../../../../run_template/* .  ln -sf ../../../../run_template/* .
117  \rm pickup* eedata.mth data.pkg data.seaice data.exf_clim data.exf data.mnc data pr*  \rm pickup* eedata.mth data.pkg data.seaice data.exf_clim data.exf data.mnc data pr*
118  \rm bathy_Hmin50.bin *192_94* *cs32* *.m lev* tren* job_c* job_altix150 job_altix92_54t  \rm bathy_Hmin50.bin *192_94* *cs32* *.m lev* tren* job_c* job_altix150 job_altix92_54t
119  ln -sf ../../../../run_template/pickup.0000000216.cube38 pickup.0000000216  ln -sf ../../../../run_template/pickup.0000000216.cube49 pickup.0000000216
120  ln -sf  ../../../../run_template/pickup_seaice.0000000216.cube38 pickup_seaice.0000000216  ln -sf  ../../../../run_template/pickup_seaice.0000000216.cube48 pickup_seaice.0000000216
121  \cp ../build/mitgcmuv .  \cp ../build/mitgcmuv .
122  qsub job_altix  qsub job_altix
123    
# Line 89  qsub job_altix Line 129  qsub job_altix
129  cd MITgcm/verification/global_ocean.cs32x15/code  cd MITgcm/verification/global_ocean.cs32x15/code
130  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .  \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
131  \cp ../../../utils/exch2/code-mods/s12t_16x32/* .  \cp ../../../utils/exch2/code-mods/s12t_16x32/* .
132    \rm CPP_EEOPTIONS.h
133  cd ../build  cd ../build
134  \rm *  \rm *
135  ../../../tools/genmake2 -mods=../code  ../../../tools/genmake2 -mods=../code
# Line 109  comparison output is in: Line 150  comparison output is in:
150  to use matlab to look at the output  to use matlab to look at the output
151  matlab  matlab
152    
153  for fld={'VICEtave','UICEtave','HEFFtave','AREAtave','QSWtave', ...  for fld={'SIuice','SIvice','SIheff','SIarea','oceQsw', ...
154           'QNETtave','FVtave','FUtave','EmPmRtave'}           'surForcT','oceTAUX','oceTAUY','surForcS'}
155      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]);
156      cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause      cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause
157  end  end
158    
159  heff=permute(readbin(['HEFFtave.0000000020.data'],[32,6,32]),[1 3 2]);  heff=permute(readbin(['SIheff.0000000072.data'],[32,6,32]),[1 3 2]);
160  mask=0*heff; mask(find(heff>.001))=1;  mask=0*heff; mask(find(heff>.001))=1;
161  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]);
162  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]);
163  clf  clf
164  subplot(221), mypcolor(uice(:,:,3)'); colorbar, title('uice, tile 3')  subplot(221), mypcolor(uice(:,:,3)'); colorbar, title('uice, tile 3')
165  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.76  
changed lines
  Added in v.1.91
  ViewVC Help
Powered by ViewVC 1.1.22