/[MITgcm]/MITgcm/verification/lab_sea/README

Diff of /MITgcm/verification/lab_sea/README

Parent Directory | 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.9 by dimitri,
Thu Jan 15 03:38:22 2004 UTC
 Line 20 
 Surface salinity relaxation is to the mo
  Forcing files are a 1979-1999 monthly climatology computed from the
  NCEP reanalysis (see pkg/seaice/SEAICE_FFIELDS.h for units and signs)
-   gairxFile      = 'u10m.labsea79'    # 10-m zonal wind
+   uwindFile      = 'u10m.labsea79'    # 10-m zonal wind
-   gairyFile      = 'v10m.labsea79'    # 10-m meridional wind
+   vwindFile      = 'v10m.labsea79'    # 10-m meridional wind
-   tairFile       = 'tair.labsea1979'  # 2-m air temperature
+   atempFile      = 'tair.labsea1979'  # 2-m air temperature
-   qaFile         = 'qa.labsea1979'    # 2-m specific humidity
+   aqhFile        = 'qa.labsea1979'    # 2-m specific humidity
-   floFile        = 'flo.labsea1979'   # longwave radiation
+   lwdownFile     = 'flo.labsea1979'   # downward longwave radiation
-   fshFile        = 'fsh.labsea1979'   # shortwave radiation
+   swdownFile     = 'fsh.labsea1979'   # downward shortwave radiation
-   rainFile       = 'prate.labsea1979' # precipitation
+   precipFile     = 'prate.labsea1979' # precipitation
    evapFile       = 'evap.labsea1979'  # evaporation
 Line 36 
 Using testscript to test sea-ice code
  Running the testscript experiment:
    cd verification
    ./testscript -force lab_sea
+ The default experiment is Experiment 8, below.
+ It uses pkg/gmredi, pkg/kpp, pkg/seaice, and pkg/exf.
  Note that fairly large differences in accuracy occur across different
  platforms.  For example, testscript comparisons between g77 (Linux)
  and f77 (SGI) generated output gives:
-    T           S           U           V
+                 T           S           U           V
-   C D M    c        m  s        m  s        m  s        m  s
+ C D M    c        m  s        m  s        m  s        m  s
-   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  .
-   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
-   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  .
-   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
  Instructions for generating 1-CPU and 2-CPU executables
  =======================================================
  Generating 1-CPU executable:
-   cd ../verification/lab_sea/input
+   cd verification/lab_sea/input
+   ln -sf ../code/packages.conf .
    ln -sf ../code/SIZE.h .
-   ln -sf ../code/CPP_OPTIONS.h .
    ln -sf ../code/CPP_EEOPTIONS.h .
-   ../../../tools/genmake -makefile
+   rm -f SEAICE_OPTIONS.h
-     ==> on alhena use:
+   ../../../tools/genmake2
-     ../../../tools/genmake -platform=o2k_noopt -makefile
    make clean
    make depend
    make
    mv mitgcmuv mitgcmuv_1
- Generating 2-CPU executable:
+ Generating 1-CPU executable with SEAICE_EXTERNAL_FORCING:
    cd ../../../verification/lab_sea/input
-   ln -sf ../code/SIZE_2x1.h SIZE.h
+   ln -sf ../code/packages.conf .
+   ln -sf ../code/SIZE.h .
+   ln -sf ../code/CPP_EEOPTIONS.h .
    ln -sf ../code/CPP_OPTIONS.h .
-   ln -sf ../code/CPP_EEOPTIONS_MPI.h CPP_EEOPTIONS.h
+   ln -sf ../code/ECCO_CPPOPTIONS.h .
-   ../../../tools/genmake -mpi -makefile
+   ln -sf ../code/SEAICE_OPTIONS.h .
-     ==> on alhena for comparison purposes use:
+   ../../../tools/genmake2
-       ../../../tools/genmake -mpi -platform=o2k_noopt -makefile
+   make clean
-     ==> on alhena for fast execution use:
+   make depend
-       ../../../tools/genmake -mpi -platform=o2k -makefile
+   make
+   mv mitgcmuv mitgcmuv_exf_forcing
+-CPU executable with SEAICE_EXTERNAL_FORCING, without READ_EVAP
+   cd ../../../verification/lab_sea/input
+   ln -sf ../code/packages.conf .
+   ln -sf ../code/SIZE.h .
+   ln -sf ../code/CPP_EEOPTIONS.h .
+   ln -sf ../code/CPP_OPTIONS.h .
+   ln -sf ../code/ECCO_CPPOPTIONS.BULK.h ECCO_CPPOPTIONS.h
+   ln -sf ../code/SEAICE_OPTIONS.h .
+   ../../../tools/genmake2
+   make clean
+   make depend
+   make
+   mv mitgcmuv mitgcmuv_compute_evap
+ Generating 1-CPU executable with SEAICE_EXTERNAL_FLUXES:
+   cd ../../../verification/lab_sea/input
+   ln -sf ../code/packages.conf .
+   ln -sf ../code/SIZE.h .
+   ln -sf ../code/CPP_EEOPTIONS.h .
+   ln -sf ../code/CPP_OPTIONS.h .
+   ln -sf ../code/ECCO_CPPOPTIONS.BULK.h ECCO_CPPOPTIONS.h
+   ln -sf ../code/SEAICE_OPTIONS.EXF_FLUXES.h SEAICE_OPTIONS.h
+   ../../../tools/genmake2
+   make clean
+   make depend
+   make
+   mv mitgcmuv mitgcmuv_exf_fluxes
+ Generating 1-CPU executable with SEAICE_MULTILEVEL:
+   cd ../../../verification/lab_sea/input
+   ln -sf ../code/packages.conf .
+   ln -sf ../code/SIZE.h .
+   ln -sf ../code/CPP_EEOPTIONS.h .
+   ln -sf ../code/SEAICE_OPTIONS.MULTILEVEL.h SEAICE_OPTIONS.h
+   ../../../tools/genmake2
+   make clean
+   make depend
+   make
+   mv mitgcmuv mitgcmuv_multilevel
+ Generating 2-CPU executable:
+   cd ../../../verification/lab_sea/input
+   ln -sf ../code/packages.conf .
+   ln -sf ../code/SIZE.2x1.h SIZE.h
+   ln -sf ../code/CPP_EEOPTIONS.h_mpi CPP_EEOPTIONS.h
+   rm -f SEAICE_OPTIONS.h
+   ../../../tools/genmake2 -mpi
    make clean
    make depend
    make
-Line 85 
 Generating 2-CPU executable:
+Line 137 
 Generating 2-CPU executable:
  Instructions for running Experiment 1
  =====================================
- 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
- files are available and that the model compiles and integrates.  It
+ the files are available and that the sea-ice model compiles and
- is the default experiment of "verification/testscript lab_sea".
+ runs.  The reading of atmospheric forcing files and the computation
+ of open-water bulk fluxes are carried out within pkg/seaice.
  To run Experiment 1:
    cd ../../../verification/lab_sea/input
-   ln -sf data.10hours data
+   ln -sf DATA.10hours data
-   ln -sf data.seaice.adi data.seaice
+   ln -sf DATA.SEAICE.lsr 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:
+ of exp1 with that from checkpoint51f sea-ice code:
    cd ../../../verification/lab_sea/matlab
    matlab
    lookat_exp1
-Line 107 
 of exp1 with that from release1_patch5 s
+Line 160 
 of exp1 with that from release1_patch5 s
  Instructions for running Experiment 2
  =====================================
- 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
- compared to that of experiment 1, which used the ADI solver.
+ compared to that of experiment 1, which used the LSR solver.
  To run Experiment 2:
    cd ../../../verification/lab_sea/input
-   ln -sf data.10hours data
+   ln -sf DATA.10hours data
-   ln -sf data.seaice.lsr data.seaice
+   ln -sf DATA.SEAICE.adi data.seaice
    mitgcmuv_1 >&! output.txt
    mkdir exp2
    mv *tave.0000000010.data exp2
-Line 135 
 All forcing files are null or constant (
+Line 188 
 All forcing files are null or constant (
  To run Experiment 3:
    cd ../../../verification/lab_sea/input
-   ln -sf data.1hour data
+   ln -sf DATA.1hour data
-   ln -sf data.seaice.testadi data.seaice
+   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
+   ln -sf DATA.SEAICE.testlsr data.seaice
    mitgcmuv_1 >&! output.txt
    mkdir exp3b
    mv *tave.0000000001.data exp3b
- Use matlab script lookat_exp3.m to compare
+ Use matlab script lookat_exp3.m to look at
- the output of exp2 to that of exp1:
+ results from exp3a and exp3b:
    cd ../../../verification/lab_sea/matlab
    matlab
    lookat_exp3
- If the dynamic ice solvers are correct, they should converge to a
+ Both solvers should and do converge to solutions
- constant solution.  The tests above show that both the ADI and the LSR
+ that are constant at each latitude.  But note that
- solver are unable to handle periodic domains and therefore that they
+ the two solutions are not exactly the same because
- cannot be "correctly" parallelized.
+ the two solvers represent slightly different
+ approximations of the viscous/plastic equations.
  Instructions for running Experiment 4
-Line 165 
 tile edges for sea-ice dynamic solvers.
+Line 219 
 tile edges for sea-ice dynamic solvers.
  To run Experiment 4:
    cd ../../../verification/lab_sea/input
-   ln -sf data.10hours data
+   ln -sf DATA.10hours data
-   ln -sf data.seaice.adi data.seaice
+   ln -sf DATA.SEAICE.lsr data.seaice
    mpirun -np 2 mitgcmuv_2x1
    mkdir exp4
    mv *tave.0000000010.data exp4
-Line 177 
 the 2-cpu output to that of exp1:
+Line 231 
 the 2-cpu output to that of exp1:
    matlab
    lookat_exp4
- The inaccuracy at the tile boundary can be decreased by increasing
+ The accuracy of the dynamic solvers is increased by
- NPSEUDO in data.seaice at the expense of computation time.
+ reducing parameter LSR_ERROR for the LSR solver and
- Also as sea-ice ages (thickens) the difference at the tiles becomes
+ by increasing parameter NPSEUDO for the ADI solver.
- increasingly small.  For solver pkg/seaice/adi.F a value of NPSEUDO=10
+ The computational cost/accuracy trade-off is quite
- appears adequate for forward integrations but cannot be used for sea-ice
+ high for the ADI solver but relatively small for
- adjoint model.  Work is underway to fix this problem.
+ the LSR solver, which is the default solver.
  Instructions for running Experiment 5
-Line 193 
 convention used by seaice_get_forcing.F
+Line 247 
 convention used by seaice_get_forcing.F
  To run Experiment 5:
    cd ../../../verification/lab_sea/input
-   ln -sf      data.2years             data
+   ln -sf      DATA.2years             data
-   ln -sf  data.seaice.adi      data.seaice
+   ln -sf  DATA.SEAICE.lsr      data.seaice
    ln -sf  evap.labsea1979  evap.labsea1980
    ln -sf   flo.labsea1979   flo.labsea1980
    ln -sf   fsh.labsea1979   fsh.labsea1980
-Line 211 
 the 2-cpu output to SMMR-SSM/I data:
+Line 265 
 the 2-cpu output to SMMR-SSM/I data:
    matlab
    lookat_exp5
- Disclaimer:
+ The comparison here is just for fun.
- The comparison here is just for fun, not really supposed
+ The configuration is too coarse and artificial
- to look anything like the data.
+ to reproduce realistic sea-ice conditions.
- Otherwise it would put a lot of people out of business :-)
  Instructions for running Experiment 6
  =====================================
-Line 224 
 The solution is compared to that of expe
+Line 278 
 The solution is compared to that of expe
  To run Experiment 6:
    cd ../../../verification/lab_sea/input
-   ln -sf data.10hours data
+   ln -sf DATA.10hours data
-   ln -sf data.seaice.nodynamics data.seaice
+   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
+   cd ../../../verification/lab_sea/matlab
    matlab
    lookat_exp6
  Instructions for running Experiment 7
  =====================================
-Line 245 
 tile edges for sea-ice thermodynamics (n
+Line 300 
 tile edges for sea-ice thermodynamics (n
  To run Experiment 7:
    cd ../../../verification/lab_sea/input
-   ln -sf data.10hours data
+   ln -sf DATA.10hours data
-   ln -sf data.seaice.nodynamics data.seaice
+   ln -sf DATA.SEAICE.nodynamics data.seaice
    mpirun -np 2 mitgcmuv_2x1
    mkdir exp7
    mv *tave.0000000010.data exp7
-Line 256 
 the output of exp7 to that of exp6:
+Line 311 
 the output of exp7 to that of exp6:
    cd ../verification/lab_sea/matlab
    matlab
    lookat_exp7
+ Instructions for running Experiment 8
+ =====================================
+ This is the default experiment which is executed by
+ verification/testscript.  It is a 1-cpu, 10-hour integration
+ used to test CPP option SEAICE_EXTERNAL_FORCING.  Atmospheric
+ state is read-in using pkg/exf, rather than pkg/seaice routines.
+ This experiment uses pkg/gmredi, pkg/kpp, pkg/seaice, and pkg/exf.
+ To run Experiment 8:
+   cd ../../../verification/lab_sea/input
+   ln -sf DATA.10hours_exf data
+   ln -sf DATA.SEAICE.exf data.seaice
+   mitgcmuv_exf_forcing >&! output.txt
+   mkdir exp8
+   mv *tave.0000000010.data exp8
+ Use matlab script lookat_exp8.m to compare
+ the output of exp8 to that of exp1:
+   cd ../../../verification/lab_sea/matlab
+   matlab
+   lookat_exp8
+ Instructions for running Experiment 9
+ =====================================
+ This is a 1-cpu, 10-hour integration used to test CPP option
+ SEAICE_EXTERNAL_FLUXES.  Both the atmospheric state and the
+ open-water surface fluxes are provided by pkg/exf.
+ To run Experiment 9:
+   cd ../../../verification/lab_sea/input
+   ln -sf DATA.10hours_exf data
+   ln -sf DATA.SEAICE.exf data.seaice
+   mitgcmuv_exf_fluxes >&! output.txt
+   mkdir exp9
+   mv *tave.0000000010.data exp9
+ Use matlab script lookat_exp9.m to compare
+ the output of exp9 to that of exp8:
+   cd ../../../verification/lab_sea/matlab
+   matlab
+   lookat_exp9
+ Instructions for running Experiment 10
+ ======================================
+ This is a 1-cpu, 10-hour integration similar to exp8, but
+ with bulk formula computation of evaporation fields.
+ To run Experiment 10:
+   cd ../../../verification/lab_sea/input
+   ln -sf DATA.10hours_exf data
+   ln -sf DATA.SEAICE.exf data.seaice
+   mitgcmuv_compute_evap >&! output.txt
+   mkdir exp10
+   mv *tave.0000000010.data exp10
+ Use matlab script lookat_exp10.m to compare
+ the output of exp10 to that of exp8:
+   cd ../../../verification/lab_sea/matlab
+   matlab
+   lookat_exp10
+ Instructions for running Experiment 11
+ ======================================
+ This is a 1-cpu, 10-hour integration similar to exp1, but
+ with multi-level sea-ice thermodynamics.
+ To run Experiment 10:
+   cd ../../../verification/lab_sea/input
+   ln -sf DATA.10hours data
+   ln -sf DATA.SEAICE.lsr data.seaice
+   mitgcmuv_multilevel >&! output.txt
+   mkdir exp11
+   mv *tave.0000000010.data exp11
+ Use matlab script lookat_exp11.m to compare
+ the output of exp11 to that of exp1:
+   cd ../../../verification/lab_sea/matlab
+   matlab
+   lookat_exp11

 Legend:



Removed from v.1.3
 


changed lines


 
Added in v.1.9
 Legend:



Removed from v.1.3
 


changed lines


 
Added in v.1.9
-Removed from v.1.3
+Added in v.1.9

	ViewVC Help
Powered by ViewVC 1.1.22