/[MITgcm]/MITgcm/verification/lab_sea/README
ViewVC logotype

Diff of /MITgcm/verification/lab_sea/README

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

revision 1.6 by dimitri, Wed Apr 30 07:04:08 2003 UTC revision 1.7 by dimitri, Fri May 23 20:19:16 2003 UTC
# Line 113  Generating 1-CPU executable with SEAICE_ Line 113  Generating 1-CPU executable with SEAICE_
113    make depend    make depend
114    make    make
115    mv mitgcmuv mitgcmuv_exf_fluxes    mv mitgcmuv mitgcmuv_exf_fluxes
116      
117    Generating 1-CPU executable with SEAICE_MULTILEVEL:
118      cd ../../../verification/lab_sea/input
119      ln -sf ../code/SIZE.h .
120      ln -sf ../code/CPP_EEOPTIONS.h .
121      ln -sf ../code/CPP_OPTIONS.NOEXF.h CPP_OPTIONS.h
122      ln -sf ../code/SEAICE_OPTIONS.MULTILEVEL.h SEAICE_OPTIONS.h
123      ../../../tools/genmake -makefile
124        ==> on alhena use:
125        ../../../tools/genmake -platform=o2k_noopt -makefile
126      make clean
127      make depend
128      make
129      mv mitgcmuv mitgcmuv_multilevel
130    
131  Generating 2-CPU executable:  Generating 2-CPU executable:
132    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
133    ln -sf ../code/SIZE.2x1.h SIZE.h    ln -sf ../code/SIZE.2x1.h SIZE.h
# Line 142  of open-water bulk fluxes are carried ou Line 156  of open-water bulk fluxes are carried ou
156  To run Experiment 1:  To run Experiment 1:
157    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
158    ln -sf DATA.10hours data    ln -sf DATA.10hours data
159    ln -sf DATA.SEAICE.adi data.seaice    ln -sf DATA.SEAICE.lsr data.seaice
160    mitgcmuv_1 >&! output.txt    mitgcmuv_1 >&! output.txt
161    mkdir exp1    mkdir exp1
162    mv *tave.0000000010.data exp1    mv *tave.0000000010.data exp1
163    
164  Use matlab script lookat_exp1.m to compare the output  Use matlab script lookat_exp1.m to compare the output
165  of exp1 with that from release1_p10 sea-ice code:  of exp1 with that from release1_p16 sea-ice code:
166    cd ../../../verification/lab_sea/matlab    cd ../../../verification/lab_sea/matlab
167    matlab    matlab
168    lookat_exp1    lookat_exp1
# Line 157  of exp1 with that from release1_p10 sea- Line 171  of exp1 with that from release1_p10 sea-
171  Instructions for running Experiment 2  Instructions for running Experiment 2
172  =====================================  =====================================
173    
174  This is a 1-cpu test of the LSR solver.  The solution is  This is a 1-cpu test of the ADI solver.  The solution is
175  compared to that of experiment 1, which used the ADI solver.  compared to that of experiment 1, which used the LSR solver.
176    
177  To run Experiment 2:  To run Experiment 2:
178    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
179    ln -sf DATA.10hours data    ln -sf DATA.10hours data
180    ln -sf DATA.SEAICE.lsr data.seaice    ln -sf DATA.SEAICE.adi data.seaice
181    mitgcmuv_1 >&! output.txt    mitgcmuv_1 >&! output.txt
182    mkdir exp2    mkdir exp2
183    mv *tave.0000000010.data exp2    mv *tave.0000000010.data exp2
# Line 201  results from exp3a and exp3b: Line 215  results from exp3a and exp3b:
215    matlab    matlab
216    lookat_exp3    lookat_exp3
217    
218  If the dynamic ice solvers are correct, they should converge to a  Both solvers should and do converge to solutions
219  solution that is constant at each latitude.  that are constant at each latitude.  But note that
220    the two solutions are not exactly the same because
221    the two solvers represent slightly different
222    approximations of the viscous/plastic equations.
223    
224    
225  Instructions for running Experiment 4  Instructions for running Experiment 4
# Line 214  tile edges for sea-ice dynamic solvers. Line 231  tile edges for sea-ice dynamic solvers.
231  To run Experiment 4:  To run Experiment 4:
232    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
233    ln -sf DATA.10hours data    ln -sf DATA.10hours data
234    ln -sf DATA.SEAICE.adi data.seaice    ln -sf DATA.SEAICE.lsr data.seaice
235    mpirun -np 2 mitgcmuv_2x1    mpirun -np 2 mitgcmuv_2x1
236    mkdir exp4    mkdir exp4
237    mv *tave.0000000010.data exp4    mv *tave.0000000010.data exp4
# Line 225  the 2-cpu output to that of exp1: Line 242  the 2-cpu output to that of exp1:
242    matlab    matlab
243    lookat_exp4    lookat_exp4
244    
245  The inaccuracy at the tile boundary can be decreased by increasing  The accuracy of the dynamic solvers is increased by
246  NPSEUDO in data.seaice at the expense of computation time.  reducing parameter LSR_ERROR for the LSR solver and
247  Also as sea-ice ages (thickens) the difference at the tiles becomes  by increasing parameter NPSEUDO for the ADI solver.
248  increasingly small.  For solver pkg/seaice/adi.F a value of NPSEUDO=10  The computational cost/accuracy trade-off is quite
249  appears adequate for forward integrations but cannot be used for sea-ice  high for the ADI solver but relatively small for
250  adjoint model.  With NPSEUDO=300 or greater, the tile disappears, but the  the LSR solver, which is the default solver.
 computational cost is very large (11 times that of teh forward model).  
251    
252    
253  Instructions for running Experiment 5  Instructions for running Experiment 5
# Line 243  convention used by seaice_get_forcing.F Line 259  convention used by seaice_get_forcing.F
259  To run Experiment 5:  To run Experiment 5:
260    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
261    ln -sf      DATA.2years             data    ln -sf      DATA.2years             data
262    ln -sf  DATA.SEAICE.adi      data.seaice    ln -sf  DATA.SEAICE.lsr      data.seaice
263    ln -sf  evap.labsea1979  evap.labsea1980      ln -sf  evap.labsea1979  evap.labsea1980  
264    ln -sf   flo.labsea1979   flo.labsea1980      ln -sf   flo.labsea1979   flo.labsea1980  
265    ln -sf   fsh.labsea1979   fsh.labsea1980      ln -sf   fsh.labsea1979   fsh.labsea1980  
# Line 260  the 2-cpu output to SMMR-SSM/I data: Line 276  the 2-cpu output to SMMR-SSM/I data:
276    matlab    matlab
277    lookat_exp5    lookat_exp5
278    
279  Disclaimer:  The comparison here is just for fun.
280  The comparison here is just for fun, not really supposed  The configuration is too coarse and artificial
281  to look anything like the data.  to reproduce realistic sea-ice conditions.
282  Otherwise it would put a lot of people out of business :-)  
283    
284  Instructions for running Experiment 6  Instructions for running Experiment 6
285  =====================================  =====================================
# Line 285  the output of exp6 to that of exp1: Line 301  the output of exp6 to that of exp1:
301    matlab    matlab
302    lookat_exp6    lookat_exp6
303    
304    
305  Instructions for running Experiment 7  Instructions for running Experiment 7
306  =====================================  =====================================
307    
# Line 319  This experiment uses pkg/gmredi, pkg/kpp Line 336  This experiment uses pkg/gmredi, pkg/kpp
336  To run Experiment 8:  To run Experiment 8:
337    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
338    ln -sf DATA.10hours_exf data    ln -sf DATA.10hours_exf data
339    ln -sf DATA.SEAICE.adi_exf data.seaice    ln -sf DATA.SEAICE.exf data.seaice
340    mitgcmuv_exf_forcing >&! output.txt    mitgcmuv_exf_forcing >&! output.txt
341    mkdir exp8    mkdir exp8
342    mv *tave.0000000010.data exp8    mv *tave.0000000010.data exp8
# Line 341  open-water surface fluxes are provided b Line 358  open-water surface fluxes are provided b
358  To run Experiment 9:  To run Experiment 9:
359    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
360    ln -sf DATA.10hours_exf data    ln -sf DATA.10hours_exf data
361    ln -sf DATA.SEAICE.adi_exf data.seaice    ln -sf DATA.SEAICE.exf data.seaice
362    mitgcmuv_exf_fluxes >&! output.txt    mitgcmuv_exf_fluxes >&! output.txt
363    mkdir exp9    mkdir exp9
364    mv *tave.0000000010.data exp9    mv *tave.0000000010.data exp9
# Line 362  with bulk formula computation of evapora Line 379  with bulk formula computation of evapora
379  To run Experiment 10:  To run Experiment 10:
380    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
381    ln -sf DATA.10hours_exf data    ln -sf DATA.10hours_exf data
382    ln -sf DATA.SEAICE.adi_exf data.seaice    ln -sf DATA.SEAICE.exf data.seaice
383    mitgcmuv_compute_evap >&! output.txt    mitgcmuv_compute_evap >&! output.txt
384    mkdir exp10    mkdir exp10
385    mv *tave.0000000010.data exp10    mv *tave.0000000010.data exp10
# Line 372  the output of exp10 to that of exp8: Line 389  the output of exp10 to that of exp8:
389    cd ../../../verification/lab_sea/matlab    cd ../../../verification/lab_sea/matlab
390    matlab    matlab
391    lookat_exp10    lookat_exp10
392    
393    
394    Instructions for running Experiment 11
395    ======================================
396    
397    This is a 1-cpu, 10-hour integration similar to exp1, but
398    with multi-level sea-ice thermodynamics.
399    
400    To run Experiment 10:
401      cd ../../../verification/lab_sea/input
402      ln -sf DATA.10hours data
403      ln -sf DATA.SEAICE.lsr data.seaice
404      mitgcmuv_multilevel >&! output.txt
405      mkdir exp11
406      mv *tave.0000000010.data exp11
407    
408    Use matlab script lookat_exp11.m to compare
409    the output of exp11 to that of exp1:
410      cd ../../../verification/lab_sea/matlab
411      matlab
412      lookat_exp11
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22