/[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.69 by dimitri, Tue Oct 10 23:27:04 2006 UTC revision 1.87 by dimitri, Sat Jan 6 07:09:50 2007 UTC
# Line 1  Line 1 
1  Instructions for setting up a cube-sphere integration with sea-ice.  Instructions for setting up a cube-sphere integration with sea-ice.
2    
3    Forcing fields and other input files needed for the CS510
4    configuration can be obtained from
5    http://ecco2.org/products/model_setups/cs510.html
6    
7    =================
8    UPDATES
9    
10    January 5, 2006: cube51, like cube50 but with bottomDragQuadratic = 0.0012
11    
12    January 5, 2006: cube50, like cube48 but with snow advection
13    
14    December 30, 2006: updates for cube49, the first
15    "blend" integration.
16    
17    December 11, 2006: preparing for cube47, same as cube43
18    but with Arctic monthly river runoff climatology.
19    
20    December 5, 2006: 216-CPU code for cube46
21    Changed diagnostics to accomodate JMCs latest and some
22    changes to code-mods/CPP_EEOPTIONS.h and W2_OPTIONS.h
23    to speed up code.
24    
25    November 7, 2006: 216-CPU code for cube43
26    Bug fixes in growth.F routine and new set of diagnostics
27    
28    October 15, 2006: 216-CPU code for cube40
29    
30  September 28, 2006: 216-CPU code for cube38  September 28, 2006: 216-CPU code for cube38
31    
32  September 4, 2006  September 4, 2006
# Line 7  Modified 32*32*6*15 experiments to use g Line 34  Modified 32*32*6*15 experiments to use g
34   verification/tutorial_held_suarez_cs/input instead of   verification/tutorial_held_suarez_cs/input instead of
35   tile* files from verification/global_ocean.cs32x15/input   tile* files from verification/global_ocean.cs32x15/input
36    
37  August 8, 2006  August 8, 2006: 216-CPU code and input used for cube37
 216-CPU code and input used for integrating cube37  
38  using "June 10, 2006" MITgcm code.  using "June 10, 2006" MITgcm code.
39    
40  for online T/S profile diagnostics all that is needed is  for online T/S profile diagnostics all that is needed is
# Line 17  for online T/S profile diagnostics all t Line 43  for online T/S profile diagnostics all t
43  see pkg/ecco/ecco_cost_init_fixed.F  see pkg/ecco/ecco_cost_init_fixed.F
44      pkg/ecco/cost_averagesfields.F      pkg/ecco/cost_averagesfields.F
45    
46    ##########################################################################
47    ##########################################################################
48    # Get forcing files, etc.
49    # Example below is for cube46 on /nobackup2a/menemenl
50    
51    cd /nobackup2a/menemenl
52    mkdir cube47
53    cd cube47
54    ln -sf ../run_template .
55    ln -sf ../ncep_rgau .
56    
57    
58  ##########################################################################  ##########################################################################
59  ##########################################################################  ##########################################################################
60  # getting the code from anonymous CVS server  # getting the code from anonymous CVS server
 #  
61    
62    bash or sh shell:    bash or sh shell:
63      $ export CVSROOT=':pserver:cvsanon@mitgcm.org:/u/gcmpack'      $ export CVSROOT=':pserver:cvsanon@mitgcm.org:/u/gcmpack'
# Line 44  cvs co MITgcm/verification/tutorial_held Line 80  cvs co MITgcm/verification/tutorial_held
80    
81  ##########################################################################  ##########################################################################
82  ##########################################################################  ##########################################################################
83    # 510*510*6*50, 216-cpu cube sphere on altix with ifort and pkg/diagnostics
84    
85    cd MITgcm/verification/global_ocean.cs32x15/code
86    \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
87    \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* .
88    \cp SIZE.h_216 SIZE.h
89    \cp packages.conf_hr packages.conf
90    
91    cd ../build
92    \rm *
93    ../../../tools/genmake2 -mods=../code -of ../code/linux_ia64_ifort+mpi_altix_nas
94    make depend
95    make -j
96    
97    cd ..
98    mkdir run
99    cd run
100    \rm *
101    \cp ../input/* .
102    \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
103    ln -sf ../../../../run_template/* .
104    \rm pickup* eedata.mth data.pkg data.seaice data.exf_clim data.exf data.mnc data pr*
105    \rm bathy_Hmin50.bin *192_94* *cs32* *.m lev* tren* job_c* job_altix150 job_altix92_54t
106    ln -sf ../../../../run_template/pickup.0000000216.cube49 pickup.0000000216
107    ln -sf  ../../../../run_template/pickup_seaice.0000000216.cube48 pickup_seaice.0000000216
108    \cp ../build/mitgcmuv .
109    qsub job_altix
110    
111    
112    ##########################################################################
113    ##########################################################################
114  # 32*32*6*15, 1-cpu cube sphere  # 32*32*6*15, 1-cpu cube sphere
115    
116  cd MITgcm/verification/global_ocean.cs32x15/code  cd MITgcm/verification/global_ocean.cs32x15/code
# Line 69  comparison output is in: Line 136  comparison output is in:
136  to use matlab to look at the output  to use matlab to look at the output
137  matlab  matlab
138    
139  for fld={'VICEtave','UICEtave','HEFFtave','AREAtave','QSWtave', ...  for fld={'SIuice','SIvice','SIheff','SIarea','oceQsw', ...
140           'QNETtave','FVtave','FUtave','EmPmRtave'}           'surForcT','oceTAUX','oceTAUY','surForcS'}
141      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]);
142      cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause      cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause
143  end  end
144    
145  heff=permute(readbin(['HEFFtave.0000000020.data'],[32,6,32]),[1 3 2]);  heff=permute(readbin(['SIheff.0000000072.data'],[32,6,32]),[1 3 2]);
146  mask=0*heff; mask(find(heff>.001))=1;  mask=0*heff; mask(find(heff>.001))=1;
147  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]);
148  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]);
149  clf  clf
150  subplot(221), mypcolor(uice(:,:,3)'); colorbar, title('uice, tile 3')  subplot(221), mypcolor(uice(:,:,3)'); colorbar, title('uice, tile 3')
151  subplot(222), mypcolor(uice(:,:,6)'); colorbar, title('uice, tile 6')  subplot(222), mypcolor(uice(:,:,6)'); colorbar, title('uice, tile 6')
# Line 87  subplot(224), mypcolor(vice(:,:,6)'); co Line 154  subplot(224), mypcolor(vice(:,:,6)'); co
154    
155    
156  ##########################################################################  ##########################################################################
 ##########################################################################  
 # 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  
   
 cd MITgcm/verification/global_ocean.cs32x15/code  
 \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .  
 \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* .  
 \cp SIZE.h_216 SIZE.h  
 \cp packages.conf_hr packages.conf  
   
 cd ../build  
 \rm *  
 ../../../tools/genmake2 -mods=../code -of  
 ../../../tools/build_options/linux_ia\  
 64_ifort+mpi_altix_nas  
 make depend  
 make -j  
   
 cd ..  
 mkdir run  
 cd run  
 \rm *  
 \cp ../input/* .  
 \cp ../../../../MITgcm_contrib/high_res_cube/input/* .  
 ln -sf ../../../../run_template/* .  
 \rm pickup* data.pkg data.seaice data.exf_clim data.exf data bathy_Hmin50.bin  
 \rm *192_94* *cs32* *.m lev* tren* job_c* job_altix150 job_altix92_54t  
 ln -sf ../../../../run_template/pickup.0000000216.cube38 pickup.0000000216  
 ln -sf  ../../../../run_template/pickup_seaice.0000000216.cube38 pickup_seaice.0000000216  
 \cp ../build/mitgcmuv .  
 qsub job_altix  
   
   
 ##########################################################################  
157  ##########################################################################  ##########################################################################
158  # 32*32*6*15, 17-tile, 1-cpu cube sphere experiment  # 32*32*6*15, 17-tile, 1-cpu cube sphere experiment
159    
Legend:
Removed from v.1.69  
changed lines
  Added in v.1.87
  ViewVC Help
Powered by ViewVC 1.1.22