/[MITgcm]/MITgcm/verification/exp2/input/data
ViewVC logotype

Annotation of /MITgcm/verification/exp2/input/data

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


Revision 1.2 - (hide annotations) (download)
Wed May 5 18:32:36 1999 UTC (25 years ago) by adcroft
Branch: MAIN
CVS Tags: checkpoint21
Changes since 1.1: +5 -5 lines
Implemented new I/O package (mdsio.F). This package does parallel
I/O in much the same way as dfile.F used to except it uses "direct
access" rather than (f77) unformatted sequential access.

Problems with dfile.F package included:
  o unnecessary memory use (each process had two global sized buffers)
  o inability to read the files it had written without post-processing
  o "tiled" files were tiled by process/thread rather than actual tiles
  o created huge numbers of files with no alternatives

Features of the mdsio.F package:
  o direct-access binary writes
  o no excessive memory use
  o ability to read/write from multiple record files
  o "tiled" files are based on "WRAPPER" tiles so that the number
    and content of files is independent of the number of threads
    and/or processes
  o option to create single "global" files rather than "tiled" files
  o ability to read both "global" and "tiled" files
    [Caveat: the tiling of files must match the model tiles]
  o checkpoints now use a single file per model section
    ie.  one file for the hydrostatic model core, one file
    for the non-hydrostatic extensions and one file for the C-D
    extensions
  o the mid-level I/O routines now is broken into more source files
    read_write_fld.F supplies basic I/O routines with the same interface
                     as the original I/O package
    read_write_rec.F supplies I/O routines which allow multiple records
    write_state.F    writes the model state
    checkpoint.F     supplies the read/write checkpoint routines

All the example input data has had to be modified to be direct-access.
Otherwise only routines that used I/O have been affected and not
all of those have been due to the continuity of arguments in
the read_write_fld.F routines.

What needs to be done?  We have to create a suite of conversion
utilities for users with old-style data. Also supply the option
for using old-style I/O, not just for die-hards but for reading
data too extensive to be converted. And more...

1 adcroft 1.1 # ====================
2     # | Model parameters |
3     # ====================
4     #
5     # Continuous equation parameters
6     &PARM01
7     tRef= 16.0 , 15.2 , 14.5 , 13.9 , 13.3 ,
8     12.4 , 11.3 , 9.9 , 8.4 , 6.7 ,
9     5.2 , 3.8 , 2.9 , 2.3 , 1.8 ,
10     1.5 , 1.1 , 0.8 , 0.66, 0.63,
11     sRef= 34.65, 34.75, 34.82, 34.87, 34.90,
12     34.90, 34.86, 34.78, 34.69, 34.60,
13     34.58, 34.62, 34.68, 34.72, 34.73,
14     34.74, 34.73, 34.73, 34.72, 34.72,
15     viscAz=1.E-3,
16     viscAh=5.E5,
17     no_slip_sides=.FALSE.,
18     no_slip_bottom=.TRUE.,
19     diffKhT=1.E3,
20     diffKzT=3.E-5,
21     diffKhS=1.E3,
22     diffKzS=3.E-5,
23     beta=1.E-11,
24     f0=1.e-4,
25     tAlpha=2.E-4,
26     sBeta =7.4E-4,
27     gravity=9.81,
28     rigidLid=.FALSE.,
29     implicitFreeSurface=.TRUE.,
30     eosType='LINEAR',
31     GMkBackground=0000.d0,
32     readBinaryPrec=32,
33     &
34    
35     # Elliptic solver parameters
36     &PARM02
37     cg2dMaxIters=1000,
38     cg2dTargetResidual=1.E-13,
39     &
40    
41     # Time stepping parameters
42     &PARM03
43     startTime=0,
44     #endTime=3110400000,
45     endTime=2808000.,
46     deltaTmom=2400.0,
47     tauCD=321428.,
48     deltaTtracer=108000.0,
49     deltaTClock =108000.0,
50     abEps=0.1,
51     pChkptFreq=0.0,
52     chkptFreq=0.0,
53     dumpFreq=2592000.0,
54     tauThetaClimRelax=2592000.0,
55     tauSaltClimRelax=2592000.0,
56     &
57    
58     # Gridding parameters
59     &PARM04
60     usingCartesianGrid=.FALSE.,
61     usingSphericalPolarGrid=.TRUE.,
62     delZ= 5.000000e+01, 5.000000e+01, 5.500000e+01, 6.000000e+01,
63     6.500000e+01, 7.000000e+01, 8.000000e+01, 9.500000e+01,
64     1.200000e+02, 1.550000e+02, 2.000000e+02, 2.600000e+02,
65     3.200000e+02, 4.000000e+02, 4.800000e+02, 5.700000e+02,
66     6.550000e+02, 7.250000e+02, 7.750000e+02, 8.150000e+02,
67     phiMin=-80.,
68     delY=40*4,
69     delX=90*4,
70     &
71    
72     # Input datasets
73     &PARM05
74 adcroft 1.2 hydrogThetaFile='theta.bin',
75     hydrogSaltFile='salt.bin',
76     bathyFile='topog.bin',
77     zonalWindFile='windx.bin',
78     meridWindFile='windy.bin',
79 adcroft 1.1 thetaClimFile='SST.bin'
80     saltClimFile='SSS.bin'
81     &

  ViewVC Help
Powered by ViewVC 1.1.22