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

Annotation of /MITgcm/verification/seaice_obcs/input/data

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


Revision 1.9 - (hide annotations) (download)
Thu May 12 22:58:27 2011 UTC (13 years ago) by jmc
Branch: MAIN
CVS Tags: checkpoint63h, checkpoint63i, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint63, checkpoint62z, checkpoint62y, checkpoint62x
Changes since 1.8: +4 -3 lines
turn on "useEnergyConservingCoriolis" (better discretisation)

1 dimitri 1.1 # ====================
2     # | Model parameters |
3     # ====================
4     #
5     # Continuous equation parameters
6     #
7     # tRef - Reference vertical potential temperature (deg C)
8     # sRef - Reference vertical salinity (PSU)
9     # viscAh - Horizontal eddy viscosity coefficient (m^2/s)
10     # viscAz - Vertical eddy viscosity coefficient (m^2/s)
11     # diffKhT - Horizontal temperature diffusivity (m^2/s)
12     # diffKzT - Vertical temperature diffusivity (m^2/s)
13     # diffKhS - Horizontal salt diffusivity (m^2/s)
14     # diffKzS - Vertical salt diffusivity (m^2/s)
15     # gravity - Acceleration due to gravity (m/s^2)
16     # rigidLid - Set to true to use rigid lid
17     # implicitFreeSurface - Set to true to use implicit free surface
18     # eosType - Flag for linear or polynomial equation of state
19     # momAdvection - On/Off flag for momentum self transport
20     # momViscosity - On/Off flag for momentum mixing
21     #
22     &PARM01
23     tRef= 24.0 , 23.0 , 22.0 , 21.0 , 20.0 ,
24     19.0 , 18.0 , 17.0 , 16.0 , 15.0 ,
25     14.0 , 13.0 , 12.0 , 11.0 , 10.0 ,
26     9.0 , 8.0 , 7.0 , 6.0, 5.0 ,
27     4.0 , 3.0 , 2.0 ,
28     sRef= 34.65, 34.75, 34.82, 34.87, 34.90,
29     34.90, 34.86, 34.78, 34.69, 34.60,
30     34.58, 34.62, 34.68, 34.72, 34.73,
31     34.74, 34.73, 34.73, 34.72, 34.72,
32     34.71, 34.70, 34.69,
33     no_slip_sides=.FALSE.,
34     no_slip_bottom=.TRUE.,
35     viscAz=1.93e-5,
36     viscAh=5.E4,
37     diffKhT=0.0,
38     diffKzT=1.46e-5,
39     diffKhS=0.0,
40     diffKzS=1.46e-5,
41 jmc 1.8 #- put small value (<< stab.limit ~ 2.e15) only to test biharmonic-diffusivity:
42     diffK4T=1.E11,
43     diffK4S=1.E11,
44     #-
45 dimitri 1.1 rigidLid=.FALSE.,
46     implicitFreeSurface=.TRUE.,
47     eosType='JMD95Z',
48     saltStepping=.TRUE.,
49     tempStepping=.TRUE.,
50     momStepping=.TRUE.,
51 jmc 1.9 useEnergyConservingCoriolis=.TRUE.,
52 dimitri 1.1 implicitDiffusion=.TRUE.,
53     implicitViscosity=.TRUE.,
54 jmc 1.5 staggerTimeStep=.TRUE.,
55     tempAdvScheme =7,
56     saltAdvScheme =7,
57     convertFW2Salt=-1,
58 dimitri 1.1 allowFreezing=.FALSE.,
59     HeatCapacity_Cp = 3986.D0,
60     gravity = 9.8156,
61     rhoNil = 1027.D0,
62 ifenty 1.7 rhoConstFresh = 999.8,
63 jmc 1.9 readBinaryPrec=32,
64     writeBinaryPrec=32,
65     useSingleCpuIO=.TRUE.,
66 dimitri 1.1 &
67    
68     # Elliptic solver parameters
69     #
70     # cg2dMaxIters - Maximum number of 2d solver iterations
71     # cg2dTargetResidual - Solver target residual
72     #
73     &PARM02
74     cg2dMaxIters=500,
75     cg2dTargetResidual=1.E-12,
76     &
77    
78     # Time stepping parameters
79     #
80     # startTime - Integration starting time (s)
81     # endTime - Integration ending time (s)
82     # deltaTMom - Timestep for momemtum equations (s)
83     # deltaTtracer - Tracer timestep (s)
84     # deltaTClock - Timestep used as model "clock" (s)
85     # abEps - Adams-Bashforth stabilising factor
86     # pChkPtFreq - Frequency of permanent check pointing (s)
87     # chkPtFreq - Frequency of rolling check pointing (s)
88     # dumpFreq - Frequency at which model state is stored (s)
89     # tauThetaClimRelax - Relaxation to climatology time scale (s)
90     # tauSaltClimRelax - Relaxation to climatology time scale (s)
91     #
92     &PARM03
93     startTime=0.0,
94 jmc 1.8 #- with current input files, cannot run longer than 10 hours
95 dimitri 1.2 endTime=21600.,
96 dimitri 1.1 deltaTmom=3600.0,
97     deltaTtracer=3600.0,
98     deltaTClock =3600.0,
99     cAdjFreq=0.,
100     abEps=0.1,
101 jmc 1.5 forcing_In_AB = .FALSE.,
102 dimitri 1.1 pChkptFreq=3600000.,
103 jmc 1.5 dumpFreq = 86400.,
104 dimitri 1.1 monitorFreq=1.,
105     &
106    
107     # Gridding parameters
108     #
109     # usingSphericalPolarGrid - On/Off flag for spherical polar coordinates
110     # delX - Zonal grid spacing (degrees)
111     # delY - Meridional grid spacing (degrees)
112     # delZ - Vertical grid spacing (m)
113 jmc 1.4 # ygOrigin - Southern boundary latitude (degrees)
114 dimitri 1.1 #
115     &PARM04
116     usingSphericalPolarGrid=.TRUE.,
117     delX=10*2.E0,
118     delY= 8*2.E0,
119     delZ= 10., 10., 15., 20., 20., 25., 35., 50., 75.,
120     100., 150., 200., 275., 350., 415., 450.,
121     500., 500., 500., 500., 500., 500., 500.,
122 jmc 1.4 ygOrigin=46.,
123     xgOrigin=280.,
124 dimitri 1.1 rSphere = 6371.D3,
125     &
126    
127     # Input datasets
128     #
129     # bathyFile - File containing bathymetry
130     # hydrogThetaFile - File containing initial potential temperature data
131     # hydrogSaltFile - File containing initial salinity data
132     # zonalWindFile - File containing zonal wind data
133     # meridWindFile - File containing meridional wind data
134     # thetaClimFile - File containing theta climatology used for relaxation
135     # saltClimFile - File containing salt climatology used for relaxation
136     #
137     &PARM05
138     bathyFile = 'bathy.seaice_obcs',
139     hydrogThetaFile = 'T.seaice_obcs',
140     hydrogSaltFile = 'S.seaice_obcs',
141     uVelInitFile = 'U.seaice_obcs',
142     vVelInitFile = 'V.seaice_obcs',
143     &

  ViewVC Help
Powered by ViewVC 1.1.22