/[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.5 by dimitri, Tue Feb 18 05:33:56 2003 UTC revision 1.8 by dimitri, Thu Aug 7 02:31:29 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 checkpoint51f 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 195  To run Experiment 3: Line 209  To run Experiment 3:
209    mkdir exp3b    mkdir exp3b
210    mv *tave.0000000001.data exp3b    mv *tave.0000000001.data exp3b
211    
212  Use matlab script lookat_exp3.m to compare  Use matlab script lookat_exp3.m to look at
213  the output of exp2 to that of exp1:  results from exp3a and exp3b:
214    cd ../../../verification/lab_sea/matlab    cd ../../../verification/lab_sea/matlab
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  constant solution.  The tests above show that both the ADI and the LSR  that are constant at each latitude.  But note that
220  solver are unable to handle periodic domains and therefore that they  the two solutions are not exactly the same because
221  cannot be "correctly" parallelized.  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 216  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 227  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.  Work is underway to fix this problem.  the LSR solver, which is the default solver.
251    
252    
253  Instructions for running Experiment 5  Instructions for running Experiment 5
# Line 244  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 261  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 282  To run Experiment 6: Line 297  To run Experiment 6:
297    
298  Use the matlab script lookat_exp6.m to compare  Use the matlab script lookat_exp6.m to compare
299  the output of exp6 to that of exp1:  the output of exp6 to that of exp1:
300    cd ../verification/lab_sea/matlab    cd ../../../verification/lab_sea/matlab
301    matlab    matlab
302    lookat_exp6    lookat_exp6
303    
304    
305  Instructions for running Experiment 7  Instructions for running Experiment 7
306  =====================================  =====================================
307    
# Line 320  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 342  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 363  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 373  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.5  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22