/[MITgcm]/MITgcm/verification/lab_sea/README
ViewVC logotype

Contents of /MITgcm/verification/lab_sea/README

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.18 - (show annotations) (download)
Mon Mar 10 01:35:17 2014 UTC (6 years, 2 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, checkpoint65, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint65o, checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64w, checkpoint64v, HEAD
Changes since 1.17: +30 -20 lines
update the intructions on how to run (in dir "run" instead of in "build"
 or "input" or elsewhere)

1 Example: Labrador Sea Region with Sea-Ice
2 =========================================
3
4 This example sets up a small (20x16x23) Labrador Sea experiment
5 coupled to a dynamic thermodynamic sea-ice model.
6 A brief description of the sea-ice model is in "seaice.ps".
7
8 The domain of integration spans 280E to 320E and 46N to 78N.
9 Horizontal grid spacing is 2 degrees.
10 The 23 vertical levels and the bathymetry file
11 bathyFile = 'bathy.labsea'
12 are obtained from the the 2-degree ECCO configuration.
13
14 Integration is initialized from annual-mean Levitus climatology
15 hydrogThetaFile = 'LevCli_temp.labsea'
16 hydrogSaltFile = 'LevCli_salt.labsea'
17
18 Surface salinity relaxation is to the monthly mean Levitus climatology
19 saltClimFile = 'SSS.labsea'
20
21 Forcing files are a 1979-1999 monthly climatology computed from the
22 NCEP reanalysis (see pkg/seaice/SEAICE_FFIELDS.h for units and signs)
23 uwindFile = 'u10m.labsea79' # 10-m zonal wind
24 vwindFile = 'v10m.labsea79' # 10-m meridional wind
25 atempFile = 'tair.labsea1979' # 2-m air temperature
26 aqhFile = 'qa.labsea1979' # 2-m specific humidity
27 lwdownFile = 'flo.labsea1979' # downward longwave radiation
28 swdownFile = 'fsh.labsea1979' # downward shortwave radiation
29 precipFile = 'prate.labsea1979' # precipitation
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 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 T S U V
51 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 .
53 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
67 make
68 cd ..
69
70 To run:
71 cd run
72 ln -s ../input/* .
73 ln -s ../build/mitgcmuv .
74 ./mitgcmuv > output.txt
75 cd ..
76
77 There is comparison output in the directory:
78 results/output.txt
79
80 Use matlab script lookat_ice.m to compare the output
81 with that from checkpoint51f sea-ice code:
82 cd ../../../verification/lab_sea/matlab
83 matlab
84 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 &

  ViewVC Help
Powered by ViewVC 1.1.22