/[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.3 by dimitri, Thu Dec 5 08:43:03 2002 UTC revision 1.5 by dimitri, Tue Feb 18 05:33:56 2003 UTC
# Line 20  Surface salinity relaxation is to the mo Line 20  Surface salinity relaxation is to the mo
20    
21  Forcing files are a 1979-1999 monthly climatology computed from the  Forcing files are a 1979-1999 monthly climatology computed from the
22  NCEP reanalysis (see pkg/seaice/SEAICE_FFIELDS.h for units and signs)  NCEP reanalysis (see pkg/seaice/SEAICE_FFIELDS.h for units and signs)
23    gairxFile      = 'u10m.labsea79'    # 10-m zonal wind    uwindFile      = 'u10m.labsea79'    # 10-m zonal wind
24    gairyFile      = 'v10m.labsea79'    # 10-m meridional wind    vwindFile      = 'v10m.labsea79'    # 10-m meridional wind
25    tairFile       = 'tair.labsea1979'  # 2-m air temperature    atempFile      = 'tair.labsea1979'  # 2-m air temperature
26    qaFile         = 'qa.labsea1979'    # 2-m specific humidity    aqhFile        = 'qa.labsea1979'    # 2-m specific humidity
27    floFile        = 'flo.labsea1979'   # longwave radiation    lwdownFile     = 'flo.labsea1979'   # downward longwave radiation
28    fshFile        = 'fsh.labsea1979'   # shortwave radiation    swdownFile     = 'fsh.labsea1979'   # downward shortwave radiation
29    rainFile       = 'prate.labsea1979' # precipitation    precipFile     = 'prate.labsea1979' # precipitation
30    evapFile       = 'evap.labsea1979'  # evaporation    evapFile       = 'evap.labsea1979'  # evaporation
31    
32    
# Line 36  Using testscript to test sea-ice code Line 36  Using testscript to test sea-ice code
36  Running the testscript experiment:  Running the testscript experiment:
37    cd verification    cd verification
38    ./testscript -force lab_sea    ./testscript -force lab_sea
39    The default experiment is Experiment 8, below.
40    It uses pkg/gmredi, pkg/kpp, pkg/seaice, and pkg/exf.
41    
42  Note that fairly large differences in accuracy occur across different  Note that fairly large differences in accuracy occur across different
43  platforms.  For example, testscript comparisons between g77 (Linux)  platforms.  For example, testscript comparisons between g77 (Linux)
44  and f77 (SGI) generated output gives:  and f77 (SGI) generated output gives:
45    
46     T           S           U           V                  T           S           U           V
47    C D M    c        m  s        m  s        m  s        m  s  C D M    c        m  s        m  s        m  s        m  s
48    n p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  .  n p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  .
49    f n k u  2  i  a  a  d  i  a  a  d  i  a  a  d  i  a  a  d  f n k u  2  i  a  a  d  i  a  a  d  i  a  a  d  i  a  a  d
50    g d e n  d  n  x  n  .  n  x  n  .  n  x  n  .  n  x  n  .  g d e n  d  n  x  n  .  n  x  n  .  n  x  n  .  n  x  n  .
51    
52    Y Y Y Y  5  5  7  7  7  8 10  9  6  6  6  6  7  5  7  5  7 FAIL  lab_sea  Y Y Y Y  8 10  9 11 10  9 11 13 10  9  8  8  9  8  9  8  9 FAIL  lab_sea
53    
54    
55  Instructions for generating 1-CPU and 2-CPU executables  Instructions for generating 1-CPU and 2-CPU executables
56  =======================================================  =======================================================
57    
58  Generating 1-CPU executable:  Generating 1-CPU executable:
59    cd ../verification/lab_sea/input    cd verification/lab_sea/input
60    ln -sf ../code/SIZE.h .    ln -sf ../code/SIZE.h .
   ln -sf ../code/CPP_OPTIONS.h .  
61    ln -sf ../code/CPP_EEOPTIONS.h .    ln -sf ../code/CPP_EEOPTIONS.h .
62      ln -sf ../code/CPP_OPTIONS.NOEXF.h CPP_OPTIONS.h
63      rm -f SEAICE_OPTIONS.h
64    ../../../tools/genmake -makefile    ../../../tools/genmake -makefile
65      ==> on alhena use:      ==> on alhena use:
66      ../../../tools/genmake -platform=o2k_noopt -makefile      ../../../tools/genmake -platform=o2k_noopt -makefile
# Line 65  Generating 1-CPU executable: Line 68  Generating 1-CPU executable:
68    make depend    make depend
69    make    make
70    mv mitgcmuv mitgcmuv_1    mv mitgcmuv mitgcmuv_1
71    
72    Generating 1-CPU executable with SEAICE_EXTERNAL_FORCING:
73      cd ../../../verification/lab_sea/input
74      ln -sf ../code/SIZE.h .
75      ln -sf ../code/CPP_EEOPTIONS.h .
76      ln -sf ../code/CPP_OPTIONS.h .
77      ln -sf ../code/ECCO_CPPOPTIONS.h .
78      ln -sf ../code/SEAICE_OPTIONS.h .
79      ../../../tools/genmake -makefile
80        ==> on alhena use:
81        ../../../tools/genmake -platform=o2k_noopt -makefile
82      make clean
83      make depend
84      make
85      mv mitgcmuv mitgcmuv_exf_forcing
86    
87    1-CPU executable with SEAICE_EXTERNAL_FORCING, without READ_EVAP
88      cd ../../../verification/lab_sea/input
89      ln -sf ../code/SIZE.h .
90      ln -sf ../code/CPP_EEOPTIONS.h .
91      ln -sf ../code/CPP_OPTIONS.h .
92      ln -sf ../code/ECCO_CPPOPTIONS.BULK.h ECCO_CPPOPTIONS.h
93      ln -sf ../code/SEAICE_OPTIONS.h .
94      ../../../tools/genmake -makefile
95        ==> on alhena use:
96        ../../../tools/genmake -platform=o2k_noopt -makefile
97      make clean
98      make depend
99      make
100      mv mitgcmuv mitgcmuv_compute_evap
101    
102    Generating 1-CPU executable with SEAICE_EXTERNAL_FLUXES:
103      cd ../../../verification/lab_sea/input
104      ln -sf ../code/SIZE.h .
105      ln -sf ../code/CPP_EEOPTIONS.h .
106      ln -sf ../code/CPP_OPTIONS.h .
107      ln -sf ../code/ECCO_CPPOPTIONS.BULK.h ECCO_CPPOPTIONS.h
108      ln -sf ../code/SEAICE_OPTIONS.EXF_FLUXES.h SEAICE_OPTIONS.h
109      ../../../tools/genmake -makefile
110        ==> on alhena use:
111        ../../../tools/genmake -platform=o2k_noopt -makefile
112      make clean
113      make depend
114      make
115      mv mitgcmuv mitgcmuv_exf_fluxes
116        
117  Generating 2-CPU executable:  Generating 2-CPU executable:
118    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
119    ln -sf ../code/SIZE_2x1.h SIZE.h    ln -sf ../code/SIZE.2x1.h SIZE.h
120    ln -sf ../code/CPP_OPTIONS.h .    ln -sf ../code/CPP_EEOPTIONS.MPI.h CPP_EEOPTIONS.h
121    ln -sf ../code/CPP_EEOPTIONS_MPI.h CPP_EEOPTIONS.h    ln -sf ../code/CPP_OPTIONS.NOEXF.h CPP_OPTIONS.h
122      rm -f SEAICE_OPTIONS.h
123    ../../../tools/genmake -mpi -makefile    ../../../tools/genmake -mpi -makefile
124      ==> on alhena for comparison purposes use:      ==> on alhena for comparison purposes use:
125        ../../../tools/genmake -mpi -platform=o2k_noopt -makefile        ../../../tools/genmake -mpi -platform=o2k_noopt -makefile
# Line 85  Generating 2-CPU executable: Line 134  Generating 2-CPU executable:
134  Instructions for running Experiment 1  Instructions for running Experiment 1
135  =====================================  =====================================
136    
137  This is a 1-cpu, 10-hour integration used to make sure that all the  This is a 1-cpu, 10-hour integration used to make sure that all
138  files are available and that the model compiles and integrates.  It  the files are available and that the sea-ice model compiles and
139  is the default experiment of "verification/testscript lab_sea".  runs.  The reading of atmospheric forcing files and the computation
140    of open-water bulk fluxes are carried out within pkg/seaice.
141    
142  To run Experiment 1:  To run Experiment 1:
143    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
144    ln -sf data.10hours data    ln -sf DATA.10hours data
145    ln -sf data.seaice.adi data.seaice    ln -sf DATA.SEAICE.adi data.seaice
146    mitgcmuv_1 >&! output.txt    mitgcmuv_1 >&! output.txt
147    mkdir exp1    mkdir exp1
148    mv *tave.0000000010.data exp1    mv *tave.0000000010.data exp1
149    
150  Use matlab script lookat_exp1.m to compare the output  Use matlab script lookat_exp1.m to compare the output
151  of exp1 with that from release1_patch5 sea-ice code:  of exp1 with that from release1_p10 sea-ice code:
152    cd ../../../verification/lab_sea/matlab    cd ../../../verification/lab_sea/matlab
153    matlab    matlab
154    lookat_exp1    lookat_exp1
# Line 112  compared to that of experiment 1, which Line 162  compared to that of experiment 1, which
162    
163  To run Experiment 2:  To run Experiment 2:
164    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
165    ln -sf data.10hours data    ln -sf DATA.10hours data
166    ln -sf data.seaice.lsr data.seaice    ln -sf DATA.SEAICE.lsr data.seaice
167    mitgcmuv_1 >&! output.txt    mitgcmuv_1 >&! output.txt
168    mkdir exp2    mkdir exp2
169    mv *tave.0000000010.data exp2    mv *tave.0000000010.data exp2
# Line 135  All forcing files are null or constant ( Line 185  All forcing files are null or constant (
185    
186  To run Experiment 3:  To run Experiment 3:
187    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
188    ln -sf data.1hour data    ln -sf DATA.1hour data
189    ln -sf data.seaice.testadi data.seaice    ln -sf DATA.SEAICE.testadi data.seaice
190    mitgcmuv_1 >&! output.txt    mitgcmuv_1 >&! output.txt
191    mkdir exp3a    mkdir exp3a
192    mv *tave.0000000001.data exp3a    mv *tave.0000000001.data exp3a
193    ln -sf data.seaice.testlsr data.seaice    ln -sf DATA.SEAICE.testlsr data.seaice
194    mitgcmuv_1 >&! output.txt    mitgcmuv_1 >&! output.txt
195    mkdir exp3b    mkdir exp3b
196    mv *tave.0000000001.data exp3b    mv *tave.0000000001.data exp3b
# Line 165  tile edges for sea-ice dynamic solvers. Line 215  tile edges for sea-ice dynamic solvers.
215    
216  To run Experiment 4:  To run Experiment 4:
217    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
218    ln -sf data.10hours data    ln -sf DATA.10hours data
219    ln -sf data.seaice.adi data.seaice    ln -sf DATA.SEAICE.adi data.seaice
220    mpirun -np 2 mitgcmuv_2x1    mpirun -np 2 mitgcmuv_2x1
221    mkdir exp4    mkdir exp4
222    mv *tave.0000000010.data exp4    mv *tave.0000000010.data exp4
# Line 193  convention used by seaice_get_forcing.F Line 243  convention used by seaice_get_forcing.F
243    
244  To run Experiment 5:  To run Experiment 5:
245    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
246    ln -sf      data.2years             data    ln -sf      DATA.2years             data
247    ln -sf  data.seaice.adi      data.seaice    ln -sf  DATA.SEAICE.adi      data.seaice
248    ln -sf  evap.labsea1979  evap.labsea1980      ln -sf  evap.labsea1979  evap.labsea1980  
249    ln -sf   flo.labsea1979   flo.labsea1980      ln -sf   flo.labsea1979   flo.labsea1980  
250    ln -sf   fsh.labsea1979   fsh.labsea1980      ln -sf   fsh.labsea1979   fsh.labsea1980  
# Line 224  The solution is compared to that of expe Line 274  The solution is compared to that of expe
274    
275  To run Experiment 6:  To run Experiment 6:
276    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
277    ln -sf data.10hours data    ln -sf DATA.10hours data
278    ln -sf data.seaice.nodynamics data.seaice    ln -sf DATA.SEAICE.nodynamics data.seaice
279    mitgcmuv_1 >&! output.txt    mitgcmuv_1 >&! output.txt
280    mkdir exp6    mkdir exp6
281    mv *tave.0000000010.data exp6    mv *tave.0000000010.data exp6
# Line 245  tile edges for sea-ice thermodynamics (n Line 295  tile edges for sea-ice thermodynamics (n
295    
296  To run Experiment 7:  To run Experiment 7:
297    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
298    ln -sf data.10hours data    ln -sf DATA.10hours data
299    ln -sf data.seaice.nodynamics data.seaice    ln -sf DATA.SEAICE.nodynamics data.seaice
300    mpirun -np 2 mitgcmuv_2x1    mpirun -np 2 mitgcmuv_2x1
301    mkdir exp7    mkdir exp7
302    mv *tave.0000000010.data exp7    mv *tave.0000000010.data exp7
# Line 256  the output of exp7 to that of exp6: Line 306  the output of exp7 to that of exp6:
306    cd ../verification/lab_sea/matlab    cd ../verification/lab_sea/matlab
307    matlab    matlab
308    lookat_exp7    lookat_exp7
309    
310    
311    Instructions for running Experiment 8
312    =====================================
313    
314    This is the default experiment which is executed by
315    verification/testscript.  It is a 1-cpu, 10-hour integration
316    used to test CPP option SEAICE_EXTERNAL_FORCING.  Atmospheric
317    state is read-in using pkg/exf, rather than pkg/seaice routines.
318    This experiment uses pkg/gmredi, pkg/kpp, pkg/seaice, and pkg/exf.
319    
320    To run Experiment 8:
321      cd ../../../verification/lab_sea/input
322      ln -sf DATA.10hours_exf data
323      ln -sf DATA.SEAICE.adi_exf data.seaice
324      mitgcmuv_exf_forcing >&! output.txt
325      mkdir exp8
326      mv *tave.0000000010.data exp8
327    
328    Use matlab script lookat_exp8.m to compare
329    the output of exp8 to that of exp1:
330      cd ../../../verification/lab_sea/matlab
331      matlab
332      lookat_exp8
333    
334    
335    Instructions for running Experiment 9
336    =====================================
337    
338    This is a 1-cpu, 10-hour integration used to test CPP option
339    SEAICE_EXTERNAL_FLUXES.  Both the atmospheric state and the
340    open-water surface fluxes are provided by pkg/exf.
341    
342    To run Experiment 9:
343      cd ../../../verification/lab_sea/input
344      ln -sf DATA.10hours_exf data
345      ln -sf DATA.SEAICE.adi_exf data.seaice
346      mitgcmuv_exf_fluxes >&! output.txt
347      mkdir exp9
348      mv *tave.0000000010.data exp9
349    
350    Use matlab script lookat_exp9.m to compare
351    the output of exp9 to that of exp8:
352      cd ../../../verification/lab_sea/matlab
353      matlab
354      lookat_exp9
355    
356    
357    Instructions for running Experiment 10
358    ======================================
359    
360    This is a 1-cpu, 10-hour integration similar to exp8, but
361    with bulk formula computation of evaporation fields.
362    
363    To run Experiment 10:
364      cd ../../../verification/lab_sea/input
365      ln -sf DATA.10hours_exf data
366      ln -sf DATA.SEAICE.adi_exf data.seaice
367      mitgcmuv_compute_evap >&! output.txt
368      mkdir exp10
369      mv *tave.0000000010.data exp10
370    
371    Use matlab script lookat_exp10.m to compare
372    the output of exp10 to that of exp8:
373      cd ../../../verification/lab_sea/matlab
374      matlab
375      lookat_exp10
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22