/[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.18 by jmc, Mon Mar 10 01:35:17 2014 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  
31    The experiment uses pkg/gmredi, pkg/kpp, pkg/seaice, and pkg/exf.
32    The test is a 1-cpu, 10-hour integration.   Both the atmospheric
33    state and the open-water surface fluxes are provided by pkg/exf.
34    
35    More pkg/seaice test experiments, configured for low and
36    high-resolution global cube-sphere domains are described
37    in MITgcm_contrib/high_res_cube/README_ice.
38    
39  Using testscript to test sea-ice code  Using testscript to test sea-ice code
40  =====================================  =====================================
41    
42  Running the testscript experiment:  Running the testscript experiment:
43    cd verification    cd MITgcm/verification
44    ./testscript -force lab_sea    ./testreport -t lab_sea [-of my_platform_optionFile]
45    
46  Note that fairly large differences in accuracy occur across different  Note that fairly large differences in accuracy occur across different
47  platforms.  For example, testscript comparisons between g77 (Linux)  platforms.  For example, testscript comparisons between g77 (Linux)
48  and f77 (SGI) generated output gives:  and f77 (SGI) generated output gives:
49    
50     T           S           U           V                  T           S           U           V
51    C D M    c        m  s        m  s        m  s        m  s  C D M    c        m  s        m  s        m  s        m  s
52    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  .
53    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
54    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  .
55    
56    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
57    
58    Instructions for generating and running a 1-CPU experiment
59  Instructions for generating 1-CPU and 2-CPU executables  ==========================================================
60  =======================================================  
61      cd MITgcm/verification/lab_sea
62  Generating 1-CPU executable:  
63    cd ../verification/lab_sea/input  Configure and compile the code:
64    ln -sf ../code/SIZE.h .    cd build
65    ln -sf ../code/CPP_OPTIONS.h .    ../../../tools/genmake2 -mods ../code [-of my_platform_optionFile]
   ln -sf ../code/CPP_EEOPTIONS.h .  
   ../../../tools/genmake -makefile  
     ==> on alhena use:  
     ../../../tools/genmake -platform=o2k_noopt -makefile  
   make clean  
   make depend  
   make  
   mv mitgcmuv mitgcmuv_1  
     
 Generating 2-CPU executable:  
   cd ../../../verification/lab_sea/input  
   ln -sf ../code/SIZE_2x1.h SIZE.h  
   ln -sf ../code/CPP_OPTIONS.h .  
   ln -sf ../code/CPP_EEOPTIONS_MPI.h CPP_EEOPTIONS.h  
   ../../../tools/genmake -mpi -makefile  
     ==> on alhena for comparison purposes use:  
       ../../../tools/genmake -mpi -platform=o2k_noopt -makefile  
     ==> on alhena for fast execution use:  
       ../../../tools/genmake -mpi -platform=o2k -makefile  
   make clean  
66    make depend    make depend
67    make    make
68    mv mitgcmuv mitgcmuv_2x1    cd ..
   
   
 Instructions for running Experiment 1  
 =====================================  
   
 This is a 1-cpu, 10-hour integration used to make sure that all the  
 files are available and that the model compiles and integrates.  It  
 is the default experiment of "verification/testscript lab_sea".  
   
 To run Experiment 1:  
   cd ../../../verification/lab_sea/input  
   ln -sf data.10hours data  
   ln -sf data.seaice.adi data.seaice  
   mitgcmuv_1 >&! output.txt  
   mkdir exp1  
   mv *tave.0000000010.data exp1  
   
 Use matlab script lookat_exp1.m to compare the output  
 of exp1 with that from release1_patch5 sea-ice code:  
   cd ../../../verification/lab_sea/matlab  
   matlab  
   lookat_exp1  
   
   
 Instructions for running Experiment 2  
 =====================================  
   
 This is a 1-cpu test of the LSR solver.  The solution is  
 compared to that of experiment 1, which used the ADI solver.  
   
 To run Experiment 2:  
   cd ../../../verification/lab_sea/input  
   ln -sf data.10hours data  
   ln -sf data.seaice.lsr data.seaice  
   mitgcmuv_1 >&! output.txt  
   mkdir exp2  
   mv *tave.0000000010.data exp2  
   
 Use matlab script lookat_exp2.m to compare  
 the output of exp2 to that of exp1:  
   cd ../../../verification/lab_sea/matlab  
   matlab  
   lookat_exp2  
69    
70    To run:
71      cd run
72      ln -s ../input/* .
73      ln -s ../build/mitgcmuv .
74      ./mitgcmuv > output.txt
75      cd ..
76    
77  Instructions for running Experiment 3  There is comparison output in the directory:
78  =====================================    results/output.txt
   
 This is a test of periodic boundary conditions for LSR  
 and ADI solvers.  The domain has a flat bottom and  
 is periodic both in the x and the y directions.  
 All forcing files are null or constant (u10m = v10m = 5 m/s).  
   
 To run Experiment 3:  
   cd ../../../verification/lab_sea/input  
   ln -sf data.1hour data  
   ln -sf data.seaice.testadi data.seaice  
   mitgcmuv_1 >&! output.txt  
   mkdir exp3a  
   mv *tave.0000000001.data exp3a  
   ln -sf data.seaice.testlsr data.seaice  
   mitgcmuv_1 >&! output.txt  
   mkdir exp3b  
   mv *tave.0000000001.data exp3b  
79    
80  Use matlab script lookat_exp3.m to compare  Use matlab script lookat_ice.m to compare the output
81  the output of exp2 to that of exp1:   with that from checkpoint51f sea-ice code:
82    cd ../../../verification/lab_sea/matlab    cd ../../../verification/lab_sea/matlab
83    matlab    matlab
84    lookat_exp3    lookat_ice
   
 If the dynamic ice solvers are correct, they should converge to a  
 constant solution.  The tests above show that both the ADI and the LSR  
 solver are unable to handle periodic domains and therefore that they  
 cannot be "correctly" parallelized.  
   
   
 Instructions for running Experiment 4  
 =====================================  
85    
86  This is a 2-cpu, 10-hour integration used to test  Instructions for generating and running a 2-CPU experiment
87  tile edges for sea-ice dynamic solvers.  ==========================================================
88    
89  To run Experiment 4:    cd MITgcm/verification/lab_sea
   cd ../../../verification/lab_sea/input  
   ln -sf data.10hours data  
   ln -sf data.seaice.adi data.seaice  
   mpirun -np 2 mitgcmuv_2x1  
   mkdir exp4  
   mv *tave.0000000010.data exp4  
90    
91  Use the matlab script lookat_exp4.m to compare  Configure and compile the code:
92  the 2-cpu output to that of exp1:    cd build
93    cd ../../../verification/lab_sea/matlab    ../../../tools/genmake2 -mpi -mods ../code [-of my_platform_optionFile]
94    matlab    ln -s ../code/SIZE.h_mpi SIZE.h
95    lookat_exp4    make depend
96      make
97  The inaccuracy at the tile boundary can be decreased by increasing    cd ..
 NPSEUDO in data.seaice at the expense of computation time.  
 Also as sea-ice ages (thickens) the difference at the tiles becomes  
 increasingly small.  For solver pkg/seaice/adi.F a value of NPSEUDO=10  
 appears adequate for forward integrations but cannot be used for sea-ice  
 adjoint model.  Work is underway to fix this problem.  
   
   
 Instructions for running Experiment 5  
 =====================================  
   
 This is a 2-cpu, 2-year, test integration.  It illustrates  
 convention used by seaice_get_forcing.F for multi-year forcing.  
   
 To run Experiment 5:  
   cd ../../../verification/lab_sea/input  
   ln -sf      data.2years             data  
   ln -sf  data.seaice.adi      data.seaice  
   ln -sf  evap.labsea1979  evap.labsea1980    
   ln -sf   flo.labsea1979   flo.labsea1980    
   ln -sf   fsh.labsea1979   fsh.labsea1980    
   ln -sf prate.labsea1979 prate.labsea1980  
   ln -sf    qa.labsea1979    qa.labsea1980      
   ln -sf  tair.labsea1979  tair.labsea1980    
   ln -sf    u10m.labsea79    u10m.labsea80      
   ln -sf    v10m.labsea79    v10m.labsea80      
   mpirun -np 2 mitgcmuv_2x1  
   
 Use the matlab script lookat_exp5.m to compare  
 the 2-cpu output to SMMR-SSM/I data:  
   cd ../../../verification/lab_sea/matlab  
   matlab  
   lookat_exp5  
   
 Disclaimer:  
 The comparison here is just for fun, not really supposed  
 to look anything like the data.  
 Otherwise it would put a lot of people out of business :-)  
   
 Instructions for running Experiment 6  
 =====================================  
   
 This is a 1-cpu test of sea-ice thermodynamics (no dynamics).  
 The solution is compared to that of experiment 1.  
   
 To run Experiment 6:  
   cd ../../../verification/lab_sea/input  
   ln -sf data.10hours data  
   ln -sf data.seaice.nodynamics data.seaice  
   mitgcmuv_1 >&! output.txt  
   mkdir exp6  
   mv *tave.0000000010.data exp6  
   
 Use the matlab script lookat_exp6.m to compare  
 the output of exp6 to that of exp1:  
   cd ../verification/lab_sea/matlab  
   matlab  
   lookat_exp6  
   
 Instructions for running Experiment 7  
 =====================================  
98    
99  This is a 2-cpu, 10-hour integration used to test  To run:
100  tile edges for sea-ice thermodynamics (no dynamics).    cd run
101  2-CPU executable from experiment 4 is required.    ln -s ../input/* .
102      mpirun -np 2 ../build/mitgcmuv
103  To run Experiment 7:    cd ..
104    cd ../../../verification/lab_sea/input  
105    ln -sf data.10hours data  Instructions for testing useExfYearlyFields (Note: might not be up-to-date)
106    ln -sf data.seaice.nodynamics data.seaice  ===========================================
107    mpirun -np 2 mitgcmuv_2x1  
108    mkdir exp7    cd MITgcm/verification/lab_sea/build
109    mv *tave.0000000010.data exp7    \rm *
110      \cp ../code/* .
111  Use the matlab script lookat_exp7.m to compare    sed 's/undef EXF_V/define EXF_V/' \
112  the output of exp7 to that of exp6:        ../../../pkg/exf/EXF_OPTIONS.h > EXF_OPTIONS.h
113    cd ../verification/lab_sea/matlab    ../../../tools/genmake2
114    matlab    make depend
115    lookat_exp7    make -j
116      cd ../run
117      \rm *
118      \cp ../input/* .
119      \mv v10m.labsea1979 v10m.labsea_1979
120      \mv u10m.labsea1979 u10m.labsea_1979
121      \mv flo.labsea1979 flo.labsea_1979
122      \mv prate.labsea1979 prate.labsea_1979
123      \mv tair.labsea1979 tair.labsea_1979
124      \mv fsh.labsea1979 fsh.labsea_1979
125      \mv qa.labsea1979 qa.labsea_1979
126      \mv SSS_monthly.labsea1979 SSS_monthly.labsea_1979
127      \ln -sf v10m.labsea_1979 v10m.labsea_1978
128      \ln -sf u10m.labsea_1979 u10m.labsea_1978
129      \ln -sf flo.labsea_1979 flo.labsea_1978
130      \ln -sf prate.labsea_1979 prate.labsea_1978
131      \ln -sf tair.labsea_1979 tair.labsea_1978
132      \ln -sf fsh.labsea_1979 fsh.labsea_1978
133      \ln -sf qa.labsea_1979 qa.labsea_1978
134      \ln -sf SSS_monthly.labsea_1979 SSS_monthly.labsea_1978
135      \ln -sf v10m.labsea_1979 v10m.labsea_1980
136      \ln -sf u10m.labsea_1979 u10m.labsea_1980
137      \ln -sf flo.labsea_1979 flo.labsea_1980
138      \ln -sf prate.labsea_1979 prate.labsea_1980
139      \ln -sf tair.labsea_1979 tair.labsea_1980
140      \ln -sf fsh.labsea_1979 fsh.labsea_1980
141      \ln -sf qa.labsea_1979 qa.labsea_1980
142      \ln -sf SSS_monthly.labsea_1979 SSS_monthly.labsea_1980
143      \mv data.exf_YearlyFields data.exf
144      \mv data_YearlyFields data
145      sed 's/tics = .TRUE./tics = .FALSE./' ../input/data.pkg > data.pkg
146      ../build/mitgcmuv >& output.txt &
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22