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

Annotation of /MITgcm/verification/exp4/input/data

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


Revision 1.5 - (hide annotations) (download)
Wed May 5 18:32:39 1999 UTC (25 years ago) by adcroft
Branch: MAIN
CVS Tags: checkpoint21
Changes since 1.4: +4 -3 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= 0.696834, 0.497738, 0.298643, 0.0995477, -0.0995477, -0.298643, -0.497738, -0.696834,
8     sRef= 8*35.,
9     viscAz=1.E-3,
10     viscAh=1.E3,
11     no_slip_sides=.FALSE.,
12     no_slip_bottom=.FALSE.,
13     viscA4=0.E12,
14     diffKhT=1.E3,
15     diffKzT=1.E-5,
16     diffKhS=1.E3,
17     diffKzS=1.E-5,
18     GMkBackground=0.d3,
19     f0=1.e-4,
20     beta=0.E-11,
21     tAlpha=2.E-4,
22     sBeta =0.E-4,
23     gravity=9.81,
24     gBaro=9.81,
25     rigidLid=.FALSE.,
26     implicitFreeSurface=.TRUE.,
27     eosType='LINEAR',
28     hFacMin=0.2,
29     openBoundaries=.TRUE.,
30 adcroft 1.3 nonHydrostatic=.FALSE.,
31 adcroft 1.1 readBinaryPrec=64,
32 adcroft 1.5 globalFiles=.TRUE.,
33 adcroft 1.1 &
34    
35     # Elliptic solver parameters
36     &PARM02
37     cg2dMaxIters=1000,
38     cg2dTargetResidual=1.E-13,
39 adcroft 1.2 cg3dMaxIters=400,
40     cg3dTargetResidual=1.E-13,
41 adcroft 1.1 &
42    
43     # Time stepping parameters
44     &PARM03
45 adcroft 1.2 niter0=0,
46 adcroft 1.5 nTimeSteps=10,
47 adcroft 1.1 deltaT=600.0,
48     abEps=0.1,
49 adcroft 1.5 pChkptFreq=3000.0,
50 adcroft 1.1 chkptFreq=0.0,
51 adcroft 1.5 dumpFreq=6000.0,
52 adcroft 1.1 &
53    
54     # Gridding parameters
55     &PARM04
56     usingCartesianGrid=.TRUE.,
57     usingSphericalPolarGrid=.FALSE.,
58     delX=80*5.e3,
59     delY=42*5.e3,
60     delZ= 8*562.5,
61     &
62    
63     # Input datasets
64     &PARM05
65     bathyFile='topog.bump',
66     &
67    
68     # Open-boundaries
69     &PARM06
70     OB_Jnorth=80*-1,
71     OB_Jsouth=80*1,
72     OB_Ieast=42*-1,
73     OB_Iwest=42*1,
74     &

  ViewVC Help
Powered by ViewVC 1.1.22