/[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.1.2.1 by heimbach, Mon Nov 11 22:03:47 2002 UTC revision 1.1.2.2 by dimitri, Thu Dec 5 08:06:18 2002 UTC
# Line 3  Example: Labrador Sea Region with Sea-Ic Line 3  Example: Labrador Sea Region with Sea-Ic
3    
4  This example sets up a small (20x16x23) Labrador Sea experiment  This example sets up a small (20x16x23) Labrador Sea experiment
5  coupled to a dynamic thermodynamic sea-ice model.  coupled to a dynamic thermodynamic sea-ice model.
6  A brief description of the sea-ice model is in "doc/seaice.ps".  A brief description of the sea-ice model is in "seaice.ps".
7    
8  The domain of integration spans 280E to 320E and 46N to 78N.  The domain of integration spans 280E to 320E and 46N to 78N.
9  Horizontal grid spacing is 2 degrees.  Horizontal grid spacing is 2 degrees.
# Line 30  NCEP reanalysis (see pkg/seaice/SEAICE_F Line 30  NCEP reanalysis (see pkg/seaice/SEAICE_F
30    evapFile       = 'evap.labsea1979'  # evaporation    evapFile       = 'evap.labsea1979'  # evaporation
31    
32    
33  Instructions for running Experiment 1  Using testscript to test sea-ice code
34  =====================================  =====================================
35    
36  This is a 1-cpu, 10-hour integration used to make sure that all the  Running the testscript experiment:
37  files are available and that the model compiles and integrates.    cd verification
38      ./testscript -force lab_sea
39  To configure and compile the code:  
40    cd MITgcmUV  Note that fairly large differences in accuracy occur across different
41    mkdir bin exe  platforms.  For example, testscript comparisons between g77 (Linux)
42    cd bin  and f77 (SGI) generated output gives:
43    ln -sf ../verification/lab_sea/code/SIZE.h .  
44    ln -sf ../verification/lab_sea/code/CPP_OPTIONS.h .     T           S           U           V
45    ln -sf ../verification/lab_sea/code/CPP_EEOPTIONS.h .    C D M    c        m  s        m  s        m  s        m  s
46    ../tools/genmake -makefile    n p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  .
47      f n k u  2  i  a  a  d  i  a  a  d  i  a  a  d  i  a  a  d
48      g d e n  d  n  x  n  .  n  x  n  .  n  x  n  .  n  x  n  .
49    
50      Y Y Y Y  5  5  7  7  7  8 10  9  6  6  6  6  7  5  7  5  7 FAIL  lab_sea
51    
52    
53    Instructions for generating 1-CPU and 2-CPU executables
54    =======================================================
55    
56    Generating 1-CPU executable:
57      cd ../verification/lab_sea/input
58      ln -sf ../code/SIZE.h .
59      ln -sf ../code/CPP_OPTIONS.h .
60      ln -sf ../code/CPP_EEOPTIONS.h .
61      ../../../tools/genmake -makefile
62        ==> on alhena use:
63        ../../../tools/genmake -platform=o2k_noopt -makefile
64      make clean
65      make depend
66      make
67      mv mitgcmuv mitgcmuv_1
68      
69    Generating 2-CPU executable:
70      cd ../../../verification/lab_sea/input
71      ln -sf ../code/SIZE_2x1.h SIZE.h
72      ln -sf ../code/CPP_OPTIONS.h .
73      ln -sf ../code/CPP_EEOPTIONS_MPI.h CPP_EEOPTIONS.h
74      ../../../tools/genmake -mpi -makefile
75        ==> on alhena for comparison purposes use:
76          ../../../tools/genmake -mpi -platform=o2k_noopt -makefile
77        ==> on alhena for fast execution use:
78          ../../../tools/genmake -mpi -platform=o2k -makefile
79      make clean
80    make depend    make depend
81    make    make
82      mv mitgcmuv mitgcmuv_2x1
83    
84    
85    Instructions for running Experiment 1
86    =====================================
87    
88  To run:  This is a 1-cpu, 10-hour integration used to make sure that all the
89    cd ../exe  files are available and that the model compiles and integrates.  It
90    cp ../verification/lab_sea/input/* .  is the default experiment of "verification/testscript lab_sea".
91    
92    To run Experiment 1:
93      cd ../../../verification/lab_sea/input
94    ln -sf data.10hours data    ln -sf data.10hours data
95    ln -sf data.seaice.adi data.seaice    ln -sf data.seaice.adi data.seaice
96    mv mitgcmuv mitgcmuv1    mitgcmuv_1 >&! output.txt
   mitgcmuv1 >&! output.txt  
97    mkdir exp1    mkdir exp1
98    mv *tave.0000000010.data exp1    mv *tave.0000000010.data exp1
99    
100  There is comparison output in:  Use matlab script lookat_exp1.m to compare the output
101    diff output.txt ../verification/lab_sea/results/output.txt | more  of exp1 with that from release1_patch5 sea-ice code:
102      cd ../../../verification/lab_sea/matlab
 Use the matlab script lookat_exp1.m to compare the output  
 of exp1 with that from release1_beta1 sea-ice code:  
   cd ../verification/lab_sea/matlab  
103    matlab    matlab
104    lookat_exp1    lookat_exp1
105    
106    
107  Instructions for running Experiment 2  Instructions for running Experiment 2
108  =====================================  =====================================
109    
110  This is a 1-cpu test of the LSR solver.  The solution is  This is a 1-cpu test of the LSR solver.  The solution is
111  compared to that of experiment 1, which used the ADI solver.  compared to that of experiment 1, which used the ADI solver.
 1-CPU executable from experiment 1 is required.  
112    
113  To run:  To run Experiment 2:
114    cd ../exe    cd ../../../verification/lab_sea/input
   cp ../verification/lab_sea/input/* .  
115    ln -sf data.10hours data    ln -sf data.10hours data
116    ln -sf data.seaice.lsr data.seaice    ln -sf data.seaice.lsr data.seaice
117    mitgcmuv1 >&! output.txt    mitgcmuv_1 >&! output.txt
118    mkdir exp2    mkdir exp2
119    mv *tave.0000000010.data exp2    mv *tave.0000000010.data exp2
120    
121  Use the matlab script lookat_exp2.m to compare  Use matlab script lookat_exp2.m to compare
122  the output of exp2 to that of exp1:  the output of exp2 to that of exp1:
123    cd ../verification/lab_sea/matlab    cd ../../../verification/lab_sea/matlab
124    matlab    matlab
125    lookat_exp2    lookat_exp2
126    
127    
128  Instructions for running Experiment 3  Instructions for running Experiment 3
129  =====================================  =====================================
130    
131  This is a test of periodic boundary conditions for LSR  This is a test of periodic boundary conditions for LSR
132  and ADI solvers.  The domain has a flat bottom and  and ADI solvers.  The domain has a flat bottom and
133  is periodic both in the x and the y directions.  is periodic both in the x and the y directions.
134  All forcing files are null or constant  All forcing files are null or constant (u10m = v10m = 5 m/s).
 (u10m = v10m = 5 m/s).  
   
 1-CPU executable from experiment 1 is required.  
135    
136  To run:  To run Experiment 3:
137    cd ../exe    cd ../../../verification/lab_sea/input
   cp ../verification/lab_sea/input/* .  
138    ln -sf data.1hour data    ln -sf data.1hour data
139    ln -sf data.seaice.testadi data.seaice    ln -sf data.seaice.testadi data.seaice
140    mitgcmuv1 >&! output.txt    mitgcmuv_1 >&! output.txt
141    mkdir exp3a    mkdir exp3a
142    mv *tave.0000000001.data exp3a    mv *tave.0000000001.data exp3a
143    ln -sf data.seaice.testlsr data.seaice    ln -sf data.seaice.testlsr data.seaice
144    mitgcmuv1 >&! output.txt    mitgcmuv_1 >&! output.txt
145    mkdir exp3b    mkdir exp3b
146    mv *tave.0000000001.data exp3b    mv *tave.0000000001.data exp3b
147    
148  Use the matlab script lookat_exp3.m to compare  Use matlab script lookat_exp3.m to compare
149  the output of exp2 to that of exp1:  the output of exp2 to that of exp1:
150    cd ../verification/lab_sea/matlab    cd ../../../verification/lab_sea/matlab
151    matlab    matlab
152    lookat_exp3    lookat_exp3
153    
# Line 123  constant solution.  The tests above show Line 156  constant solution.  The tests above show
156  solver are unable to handle periodic domains and therefore that they  solver are unable to handle periodic domains and therefore that they
157  cannot be "correctly" parallelized.  cannot be "correctly" parallelized.
158    
159    
160  Instructions for running Experiment 4  Instructions for running Experiment 4
161  =====================================  =====================================
162    
163  This is a 2-cpu, 10-hour integration used to test  This is a 2-cpu, 10-hour integration used to test
164  tile edges for sea-ice dynamic solvers.  tile edges for sea-ice dynamic solvers.
165    
166  To configure and compile the code:  To run Experiment 4:
167    cd ../bin    cd ../../../verification/lab_sea/input
   rm -rf *.f  
   rm -rf *.o  
   ln -sf ../verification/lab_sea/code/SIZE_2x1.h SIZE.h  
   ln -sf ../verification/lab_sea/code/CPP_OPTIONS.h .  
   ln -sf ../verification/lab_sea/code/CPP_EEOPTIONS_MPI.h CPP_EEOPTIONS.h  
   ../tools/genmake -mpi -makefile  
      ( on alhena use: ../tools/genmake -mpi -platform=o2k -makefile )  
   make depend  
   make  
   
 To run:  
   cd ../exe  
   cp ../verification/lab_sea/input/* .  
168    ln -sf data.10hours data    ln -sf data.10hours data
169    ln -sf data.seaice.adi data.seaice    ln -sf data.seaice.adi data.seaice
   mv mitgcmuv mitgcmuv_2x1  
170    mpirun -np 2 mitgcmuv_2x1    mpirun -np 2 mitgcmuv_2x1
171    mkdir exp4    mkdir exp4
172    mv *tave.0000000010.data exp4    mv *tave.0000000010.data exp4
173    
174  Use the matlab script lookat_exp4.m to compare  Use the matlab script lookat_exp4.m to compare
175  the 2-cpu output to that of exp1:  the 2-cpu output to that of exp1:
176    cd ../verification/lab_sea/matlab    cd ../../../verification/lab_sea/matlab
177    matlab    matlab
178    lookat_exp4    lookat_exp4
179    
# Line 169  Instructions for running Experiment 5 Line 189  Instructions for running Experiment 5
189  =====================================  =====================================
190    
191  This is a 2-cpu, 2-year, test integration.  It illustrates  This is a 2-cpu, 2-year, test integration.  It illustrates
192  convention used for forcing data files by seaice_get_forcing.F.  convention used by seaice_get_forcing.F for multi-year forcing.
193    
194  To run:  To run Experiment 5:
195    cd ../exe    cd ../../../verification/lab_sea/input
196    cp ../verification/lab_sea/input/* .    ln -sf      data.2years             data
197    cp eedata_mpi eedata    ln -sf  data.seaice.adi      data.seaice
198    ln -sf data.2years data    ln -sf  evap.labsea1979  evap.labsea1980  
199    ln -sf data.seaice.adi data.seaice    ln -sf   flo.labsea1979   flo.labsea1980  
200    ln -sf evap.labsea1979  evap.labsea1980      ln -sf   fsh.labsea1979   fsh.labsea1980  
   ln -sf flo.labsea1979   flo.labsea1980    
   ln -sf fsh.labsea1979   fsh.labsea1980    
201    ln -sf prate.labsea1979 prate.labsea1980    ln -sf prate.labsea1979 prate.labsea1980
202    ln -sf qa.labsea1979    qa.labsea1980        ln -sf    qa.labsea1979    qa.labsea1980    
203    ln -sf tair.labsea1979  tair.labsea1980      ln -sf  tair.labsea1979  tair.labsea1980  
204    ln -sf u10m.labsea79    u10m.labsea80        ln -sf    u10m.labsea79    u10m.labsea80    
205    ln -sf v10m.labsea79    v10m.labsea80        ln -sf    v10m.labsea79    v10m.labsea80    
206    mpirun -np 2 mitgcmuv_2x1    mpirun -np 2 mitgcmuv_2x1
207    
208  Use the matlab script lookat_exp3.m to compare  Use the matlab script lookat_exp5.m to compare
209  the 2-cpu output to SMMR-SSM/I data:  the 2-cpu output to SMMR-SSM/I data:
210    cd ../verification/lab_sea/matlab    cd ../../../verification/lab_sea/matlab
211    matlab    matlab
212    lookat_exp5    lookat_exp5
213    
# Line 203  Instructions for running Experiment 6 Line 221  Instructions for running Experiment 6
221    
222  This is a 1-cpu test of sea-ice thermodynamics (no dynamics).  This is a 1-cpu test of sea-ice thermodynamics (no dynamics).
223  The solution is compared to that of experiment 1.  The solution is compared to that of experiment 1.
 1-CPU executable from experiment 1 is required.  
224    
225  To run:  To run Experiment 6:
226    cd ../exe    cd ../../../verification/lab_sea/input
   cp ../verification/lab_sea/input/* .  
227    ln -sf data.10hours data    ln -sf data.10hours data
228    ln -sf data.seaice.nodynamics data.seaice    ln -sf data.seaice.nodynamics data.seaice
229    mitgcmuv1 >&! output.txt    mitgcmuv_1 >&! output.txt
230    mkdir exp6    mkdir exp6
231    mv *tave.0000000010.data exp6    mv *tave.0000000010.data exp6
232    
233  Use the matlab script lookat_exp6.m to compare  Use the matlab script lookat_exp6.m to compare
234  the output of exp2 to that of exp1:  the output of exp6 to that of exp1:
235    cd ../verification/lab_sea/matlab    cd ../verification/lab_sea/matlab
236    matlab    matlab
237    lookat_exp6    lookat_exp6
238    
239    Instructions for running Experiment 7
240    =====================================
241    
242    This is a 2-cpu, 10-hour integration used to test
243    tile edges for sea-ice thermodynamics (no dynamics).
244    2-CPU executable from experiment 4 is required.
245    
246    To run Experiment 7:
247      cd ../../../verification/lab_sea/input
248      ln -sf data.10hours data
249      ln -sf data.seaice.nodynamics data.seaice
250      mpirun -np 2 mitgcmuv_2x1
251      mkdir exp7
252      mv *tave.0000000010.data exp7
253    
254    Use the matlab script lookat_exp7.m to compare
255    the output of exp7 to that of exp6:
256      cd ../verification/lab_sea/matlab
257      matlab
258      lookat_exp7

Legend:
Removed from v.1.1.2.1  
changed lines
  Added in v.1.1.2.2

  ViewVC Help
Powered by ViewVC 1.1.22