/[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.4.2.2.2.1 by dimitri, Wed Mar 12 15:07:22 2003 UTC revision 1.16 by dimitri, Tue Jan 29 05:57:10 2008 UTC
# Line 5  This example sets up a small (20x16x23) Line 5  This example sets up a small (20x16x23)
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 "seaice.ps".  A brief description of the sea-ice model is in "seaice.ps".
7    
 This experiment 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.  
   
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.
10  The 23 vertical levels and the bathymetry file  The 23 vertical levels and the bathymetry file
# Line 32  NCEP reanalysis (see pkg/seaice/SEAICE_F Line 27  NCEP reanalysis (see pkg/seaice/SEAICE_F
27    lwdownFile     = 'flo.labsea1979'   # downward longwave radiation    lwdownFile     = 'flo.labsea1979'   # downward longwave radiation
28    swdownFile     = 'fsh.labsea1979'   # downward shortwave radiation    swdownFile     = 'fsh.labsea1979'   # downward shortwave radiation
29    precipFile     = 'prate.labsea1979' # precipitation    precipFile     = 'prate.labsea1979' # precipitation
   evapFile       = 'evap.labsea1979'  # evaporation  
30    
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  Instructions for running default Experiment  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    
 Compiling the code:  
   cd verification/lab_sea/input  
   ln -sf ../code/*.h .  
   ../../../tools/genmake -makefile  
   make depend  
   make  
39    
40  Running the code:  Using testscript to test sea-ice code
41    ./mitgcmuv >&! output.txt  =====================================
42    
43    Running the testscript experiment:
44      cd MITgcm/verification
45      ./testreport -t lab_sea
46    
47    Note that fairly large differences in accuracy occur across different
48    platforms.  For example, testscript comparisons between g77 (Linux)
49    and f77 (SGI) generated output gives:
50    
51  There is comparison output in:                  T           S           U           V
52    diff output.txt ../results/output.txt | more  C D M    c        m  s        m  s        m  s        m  s
53    n p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  .
54    f n k u  2  i  a  a  d  i  a  a  d  i  a  a  d  i  a  a  d
55    g d e n  d  n  x  n  .  n  x  n  .  n  x  n  .  n  x  n  .
56    
57    Y Y Y Y  8 10  9 11 10  9 11 13 10  9  8  8  9  8  9  8  9 FAIL  lab_sea
58    
59    
60    Instructions for generating and running a 1-CPU experiment
61    ==========================================================
62    
63      cd MITgcm/verification/lab_sea
64      mkdir build
65      cd build
66      cp ../code/*.h ../code/packages.conf .
67      ../../../tools/genmake2
68      make depend
69      make
70      cd ../input
71      ../build/mitgcmuv > output.txt
72    
73  Matlab comparison script:  Use matlab script lookat_exp1.m to compare the output
74    of exp1 with that from checkpoint51f sea-ice code:
75    cd ../../../verification/lab_sea/matlab    cd ../../../verification/lab_sea/matlab
76    matlab    matlab
77    lookat_exp    lookat_ice
78    
79    
80    Instructions for generating and running a 2-CPU experiment on kalpana
81    =====================================================================
82    
83      cd MITgcm/verification/lab_sea
84      mkdir build
85      cd build
86      \cp ../code/*.h ../code/packages.conf .
87      \cp ../code/SIZE.h_mpi SIZE.h
88      ../../../tools/genmake2 -of ../../../tools/build_options/linux_ia64_efc+mpi_altix
89      make depend
90      make
91      cd ../input
92      mpirun -np 2 ../build/mitgcmuv
93    
94    
95    Instructions for testing useExfYearlyFields
96    ===========================================
97    
98      cd MITgcm/verification/lab_sea/build
99      \rm *
100      \cp ../code/* .
101      sed 's/undef EXF_V/define EXF_V/' \
102          ../../../pkg/exf/EXF_OPTIONS.h > EXF_OPTIONS.h
103      ../../../tools/genmake2
104      make depend
105      make -j
106      cd ../run
107      \rm *
108      \cp ../input/* .
109      \mv v10m.labsea1979 v10m.labsea_1979
110      \mv u10m.labsea1979 u10m.labsea_1979
111      \mv flo.labsea1979 flo.labsea_1979
112      \mv prate.labsea1979 prate.labsea_1979
113      \mv tair.labsea1979 tair.labsea_1979
114      \mv fsh.labsea1979 fsh.labsea_1979
115      \mv qa.labsea1979 qa.labsea_1979
116      \mv SSS_monthly.labsea1979 SSS_monthly.labsea_1979
117      \ln -sf v10m.labsea_1979 v10m.labsea_1978
118      \ln -sf u10m.labsea_1979 u10m.labsea_1978
119      \ln -sf flo.labsea_1979 flo.labsea_1978
120      \ln -sf prate.labsea_1979 prate.labsea_1978
121      \ln -sf tair.labsea_1979 tair.labsea_1978
122      \ln -sf fsh.labsea_1979 fsh.labsea_1978
123      \ln -sf qa.labsea_1979 qa.labsea_1978
124      \ln -sf SSS_monthly.labsea_1979 SSS_monthly.labsea_1978
125      \mv data.exf_YearlyFields data.exf
126      \mv data_YearlyFields data
127      sed 's/tics = .TRUE./tics = .FALSE./' ../input/data.pkg > data.pkg
128      ../build/mitgcmuv >& output.txt &

Legend:
Removed from v.1.4.2.2.2.1  
changed lines
  Added in v.1.16

  ViewVC Help
Powered by ViewVC 1.1.22