/[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.6 by dimitri,
Wed Apr 30 07:04:08 2003 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/SIZE.h .
-   ln -sf ../code/CPP_OPTIONS.h .
    ln -sf ../code/CPP_EEOPTIONS.h .
+   ln -sf ../code/CPP_OPTIONS.NOEXF.h CPP_OPTIONS.h
+   rm -f SEAICE_OPTIONS.h
    ../../../tools/genmake -makefile
      ==> on alhena use:
      ../../../tools/genmake -platform=o2k_noopt -makefile
-Line 65 
 Generating 1-CPU executable:
+Line 68 
 Generating 1-CPU executable:
    make depend
    make
    mv mitgcmuv mitgcmuv_1
+ Generating 1-CPU executable with SEAICE_EXTERNAL_FORCING:
+   cd ../../../verification/lab_sea/input
+   ln -sf ../code/SIZE.h .
+   ln -sf ../code/CPP_EEOPTIONS.h .
+   ln -sf ../code/CPP_OPTIONS.h .
+   ln -sf ../code/ECCO_CPPOPTIONS.h .
+   ln -sf ../code/SEAICE_OPTIONS.h .
+   ../../../tools/genmake -makefile
+     ==> on alhena use:
+     ../../../tools/genmake -platform=o2k_noopt -makefile
+   make clean
+   make depend
+   make
+   mv mitgcmuv mitgcmuv_exf_forcing
+-CPU executable with SEAICE_EXTERNAL_FORCING, without READ_EVAP
+   cd ../../../verification/lab_sea/input
+   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/genmake -makefile
+     ==> on alhena use:
+     ../../../tools/genmake -platform=o2k_noopt -makefile
+   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/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/genmake -makefile
+     ==> on alhena use:
+     ../../../tools/genmake -platform=o2k_noopt -makefile
+   make clean
+   make depend
+   make
+   mv mitgcmuv mitgcmuv_exf_fluxes
  Generating 2-CPU executable:
    cd ../../../verification/lab_sea/input
-   ln -sf ../code/SIZE_2x1.h SIZE.h
+   ln -sf ../code/SIZE.2x1.h SIZE.h
-   ln -sf ../code/CPP_OPTIONS.h .
+   ln -sf ../code/CPP_EEOPTIONS.MPI.h CPP_EEOPTIONS.h
-   ln -sf ../code/CPP_EEOPTIONS_MPI.h CPP_EEOPTIONS.h
+   ln -sf ../code/CPP_OPTIONS.NOEXF.h CPP_OPTIONS.h
+   rm -f SEAICE_OPTIONS.h
    ../../../tools/genmake -mpi -makefile
      ==> on alhena for comparison purposes use:
        ../../../tools/genmake -mpi -platform=o2k_noopt -makefile
-Line 85 
 Generating 2-CPU executable:
+Line 134 
 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.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:
+ of exp1 with that from release1_p10 sea-ice code:
    cd ../../../verification/lab_sea/matlab
    matlab
    lookat_exp1
-Line 112 
 compared to that of experiment 1, which
+Line 162 
 compared to that of experiment 1, which
  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.lsr data.seaice
    mitgcmuv_1 >&! output.txt
    mkdir exp2
    mv *tave.0000000010.data exp2
-Line 135 
 All forcing files are null or constant (
+Line 185 
 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
- constant solution.  The tests above show that both the ADI and the LSR
+ solution that is constant at each latitude.
- solver are unable to handle periodic domains and therefore that they
- cannot be "correctly" parallelized.
  Instructions for running Experiment 4
-Line 165 
 tile edges for sea-ice dynamic solvers.
+Line 213 
 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.adi data.seaice
    mpirun -np 2 mitgcmuv_2x1
    mkdir exp4
    mv *tave.0000000010.data exp4
-Line 182 
 NPSEUDO in data.seaice at the expense of
+Line 230 
 NPSEUDO in data.seaice at the expense of
  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.
+ adjoint model.  With NPSEUDO=300 or greater, the tile disappears, but the
+ computational cost is very large (11 times that of teh forward model).
  Instructions for running Experiment 5
-Line 193 
 convention used by seaice_get_forcing.F
+Line 242 
 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.adi      data.seaice
    ln -sf  evap.labsea1979  evap.labsea1980
    ln -sf   flo.labsea1979   flo.labsea1980
    ln -sf   fsh.labsea1979   fsh.labsea1980
-Line 224 
 The solution is compared to that of expe
+Line 273 
 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
-Line 245 
 tile edges for sea-ice thermodynamics (n
+Line 294 
 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 305 
 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.adi_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.adi_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.adi_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

 Legend:



Removed from v.1.3
 


changed lines


 
Added in v.1.6
 Legend:



Removed from v.1.3
 


changed lines


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

	ViewVC Help
Powered by ViewVC 1.1.22