/[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.3 by dimitri, Wed Apr 30 08:23:08 2003 UTC revision 1.18 by jmc, Mon Mar 10 01:35:17 2014 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    
39    Using testscript to test sea-ice code
40    =====================================
41    
42    Running the testscript experiment:
43      cd MITgcm/verification
44      ./testreport -t lab_sea [-of my_platform_optionFile]
45    
46    Note that fairly large differences in accuracy occur across different
47    platforms.  For example, testscript comparisons between g77 (Linux)
48    and f77 (SGI) generated output gives:
49    
50  Compiling the code:                  T           S           U           V
51    cd verification/lab_sea/input  C D M    c        m  s        m  s        m  s        m  s
52    ln -sf ../code/*.h .  n p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  .
53    ../../../tools/genmake -makefile  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  .
55    
56    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    ==========================================================
60    
61      cd MITgcm/verification/lab_sea
62    
63    Configure and compile the code:
64      cd build
65      ../../../tools/genmake2 -mods ../code [-of my_platform_optionFile]
66    make depend    make depend
67    make    make
68      cd ..
69    
70  Running the code:  To run:
71    ./mitgcmuv >&! output.txt    cd run
72      ln -s ../input/* .
73      ln -s ../build/mitgcmuv .
74      ./mitgcmuv > output.txt
75      cd ..
76    
77  There is comparison output in:  There is comparison output in the directory:
78    diff output.txt ../results/output.txt | more    results/output.txt
79    
80  Matlab comparison script:  Use matlab script lookat_ice.m to compare the output
81     with that from checkpoint51f sea-ice code:
82    cd ../../../verification/lab_sea/matlab    cd ../../../verification/lab_sea/matlab
83    matlab    matlab
84    lookat_exp    lookat_ice
85    
86    Instructions for generating and running a 2-CPU experiment
87    ==========================================================
88    
89      cd MITgcm/verification/lab_sea
90    
91    Configure and compile the code:
92      cd build
93      ../../../tools/genmake2 -mpi -mods ../code [-of my_platform_optionFile]
94      ln -s ../code/SIZE.h_mpi SIZE.h
95      make depend
96      make
97      cd ..
98    
99    To run:
100      cd run
101      ln -s ../input/* .
102      mpirun -np 2 ../build/mitgcmuv
103      cd ..
104    
105    Instructions for testing useExfYearlyFields (Note: might not be up-to-date)
106    ===========================================
107    
108      cd MITgcm/verification/lab_sea/build
109      \rm *
110      \cp ../code/* .
111      sed 's/undef EXF_V/define EXF_V/' \
112          ../../../pkg/exf/EXF_OPTIONS.h > EXF_OPTIONS.h
113      ../../../tools/genmake2
114      make depend
115      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.4.2.3  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22