/[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.4 by dimitri, Sat Dec 28 10:11:11 2002 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    lwfluxFile     = 'flo.labsea1979'   # longwave radiation
28    fshFile        = 'fsh.labsea1979'   # shortwave radiation    swfluxFile     = 'fsh.labsea1979'   # 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    Generating 1-CPU executable with SEAICE_EXTERNAL_FLUXES:
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.EXF_FLUXES.h 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_exf_fluxes
101        
102  Generating 2-CPU executable:  Generating 2-CPU executable:
103    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
104    ln -sf ../code/SIZE_2x1.h SIZE.h    ln -sf ../code/SIZE.2x1.h SIZE.h
105    ln -sf ../code/CPP_OPTIONS.h .    ln -sf ../code/CPP_EEOPTIONS.MPI.h CPP_EEOPTIONS.h
106    ln -sf ../code/CPP_EEOPTIONS_MPI.h CPP_EEOPTIONS.h    ln -sf ../code/CPP_OPTIONS.NOEXF.h CPP_OPTIONS.h
107      rm -f SEAICE_OPTIONS.h
108    ../../../tools/genmake -mpi -makefile    ../../../tools/genmake -mpi -makefile
109      ==> on alhena for comparison purposes use:      ==> on alhena for comparison purposes use:
110        ../../../tools/genmake -mpi -platform=o2k_noopt -makefile        ../../../tools/genmake -mpi -platform=o2k_noopt -makefile
# Line 85  Generating 2-CPU executable: Line 119  Generating 2-CPU executable:
119  Instructions for running Experiment 1  Instructions for running Experiment 1
120  =====================================  =====================================
121    
122  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
123  files are available and that the model compiles and integrates.  It  the files are available and that the sea-ice model compiles and
124  is the default experiment of "verification/testscript lab_sea".  runs.  The reading of atmospheric forcing files and the computation
125    of open-water bulk fluxes are carried out within pkg/seaice.
126    
127  To run Experiment 1:  To run Experiment 1:
128    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
129    ln -sf data.10hours data    ln -sf DATA.10hours data
130    ln -sf data.seaice.adi data.seaice    ln -sf DATA.SEAICE.adi data.seaice
131    mitgcmuv_1 >&! output.txt    mitgcmuv_1 >&! output.txt
132    mkdir exp1    mkdir exp1
133    mv *tave.0000000010.data exp1    mv *tave.0000000010.data exp1
134    
135  Use matlab script lookat_exp1.m to compare the output  Use matlab script lookat_exp1.m to compare the output
136  of exp1 with that from release1_patch5 sea-ice code:  of exp1 with that from release1_p10 sea-ice code:
137    cd ../../../verification/lab_sea/matlab    cd ../../../verification/lab_sea/matlab
138    matlab    matlab
139    lookat_exp1    lookat_exp1
# Line 112  compared to that of experiment 1, which Line 147  compared to that of experiment 1, which
147    
148  To run Experiment 2:  To run Experiment 2:
149    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
150    ln -sf data.10hours data    ln -sf DATA.10hours data
151    ln -sf data.seaice.lsr data.seaice    ln -sf DATA.SEAICE.lsr data.seaice
152    mitgcmuv_1 >&! output.txt    mitgcmuv_1 >&! output.txt
153    mkdir exp2    mkdir exp2
154    mv *tave.0000000010.data exp2    mv *tave.0000000010.data exp2
# Line 135  All forcing files are null or constant ( Line 170  All forcing files are null or constant (
170    
171  To run Experiment 3:  To run Experiment 3:
172    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
173    ln -sf data.1hour data    ln -sf DATA.1hour data
174    ln -sf data.seaice.testadi data.seaice    ln -sf DATA.SEAICE.testadi data.seaice
175    mitgcmuv_1 >&! output.txt    mitgcmuv_1 >&! output.txt
176    mkdir exp3a    mkdir exp3a
177    mv *tave.0000000001.data exp3a    mv *tave.0000000001.data exp3a
178    ln -sf data.seaice.testlsr data.seaice    ln -sf DATA.SEAICE.testlsr data.seaice
179    mitgcmuv_1 >&! output.txt    mitgcmuv_1 >&! output.txt
180    mkdir exp3b    mkdir exp3b
181    mv *tave.0000000001.data exp3b    mv *tave.0000000001.data exp3b
# Line 165  tile edges for sea-ice dynamic solvers. Line 200  tile edges for sea-ice dynamic solvers.
200    
201  To run Experiment 4:  To run Experiment 4:
202    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
203    ln -sf data.10hours data    ln -sf DATA.10hours data
204    ln -sf data.seaice.adi data.seaice    ln -sf DATA.SEAICE.adi data.seaice
205    mpirun -np 2 mitgcmuv_2x1    mpirun -np 2 mitgcmuv_2x1
206    mkdir exp4    mkdir exp4
207    mv *tave.0000000010.data exp4    mv *tave.0000000010.data exp4
# Line 193  convention used by seaice_get_forcing.F Line 228  convention used by seaice_get_forcing.F
228    
229  To run Experiment 5:  To run Experiment 5:
230    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
231    ln -sf      data.2years             data    ln -sf      DATA.2years             data
232    ln -sf  data.seaice.adi      data.seaice    ln -sf  DATA.SEAICE.adi      data.seaice
233    ln -sf  evap.labsea1979  evap.labsea1980      ln -sf  evap.labsea1979  evap.labsea1980  
234    ln -sf   flo.labsea1979   flo.labsea1980      ln -sf   flo.labsea1979   flo.labsea1980  
235    ln -sf   fsh.labsea1979   fsh.labsea1980      ln -sf   fsh.labsea1979   fsh.labsea1980  
# Line 224  The solution is compared to that of expe Line 259  The solution is compared to that of expe
259    
260  To run Experiment 6:  To run Experiment 6:
261    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
262    ln -sf data.10hours data    ln -sf DATA.10hours data
263    ln -sf data.seaice.nodynamics data.seaice    ln -sf DATA.SEAICE.nodynamics data.seaice
264    mitgcmuv_1 >&! output.txt    mitgcmuv_1 >&! output.txt
265    mkdir exp6    mkdir exp6
266    mv *tave.0000000010.data exp6    mv *tave.0000000010.data exp6
# Line 245  tile edges for sea-ice thermodynamics (n Line 280  tile edges for sea-ice thermodynamics (n
280    
281  To run Experiment 7:  To run Experiment 7:
282    cd ../../../verification/lab_sea/input    cd ../../../verification/lab_sea/input
283    ln -sf data.10hours data    ln -sf DATA.10hours data
284    ln -sf data.seaice.nodynamics data.seaice    ln -sf DATA.SEAICE.nodynamics data.seaice
285    mpirun -np 2 mitgcmuv_2x1    mpirun -np 2 mitgcmuv_2x1
286    mkdir exp7    mkdir exp7
287    mv *tave.0000000010.data exp7    mv *tave.0000000010.data exp7
# Line 256  the output of exp7 to that of exp6: Line 291  the output of exp7 to that of exp6:
291    cd ../verification/lab_sea/matlab    cd ../verification/lab_sea/matlab
292    matlab    matlab
293    lookat_exp7    lookat_exp7
294    
295    
296    Instructions for running Experiment 8
297    =====================================
298    
299    This is the default experiment which is executed by
300    verification/testscript.  It is a 1-cpu, 10-hour integration
301    used to test CPP option SEAICE_EXTERNAL_FORCING.  Atmospheric
302    state is read-in using pkg/exf, rather than pkg/seaice routines.
303    This experiment uses pkg/gmredi, pkg/kpp, pkg/seaice, and pkg/exf.
304    
305    To run Experiment 8:
306      cd ../../../verification/lab_sea/input
307      ln -sf DATA.10hours_exf data
308      ln -sf DATA.SEAICE.adi_exf data.seaice
309      mitgcmuv_exf_forcing >&! output.txt
310      mkdir exp8
311      mv *tave.0000000010.data exp8
312    
313    Use matlab script lookat_exp8.m to compare
314    the output of exp8 to that of exp1:
315      cd ../../../verification/lab_sea/matlab
316      matlab
317      lookat_exp8
318    
319    
320    Instructions for running Experiment 9
321    =====================================
322    
323    This is a 1-cpu, 10-hour integration used to test CPP option
324    SEAICE_EXTERNAL_FLUXES.  Both the atmospheric state and the
325    open-water surface fluxes are provided by pkg/exf.
326    
327    To run Experiment 9:
328      cd ../../../verification/lab_sea/input
329      ln -sf DATA.10hours_exf data
330      ln -sf DATA.SEAICE.adi_exf data.seaice
331      mitgcmuv_exf_fluxes >&! output.txt
332      mkdir exp9
333      mv *tave.0000000010.data exp9
334    
335    Use matlab script lookat_exp9.m to compare
336    the output of exp9 to that of exp8:
337      cd ../../../verification/lab_sea/matlab
338      matlab
339      lookat_exp9
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22