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

Contents of /MITgcm/verification/seaice_obcs/input/data

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


Revision 1.8 - (show annotations) (download)
Fri Apr 29 20:41:43 2011 UTC (13 years ago) by jmc
Branch: MAIN
Changes since 1.7: +5 -0 lines
only to test biharmonic-diffusivity+OBCS: add small diffK4T,S=1.e11

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 #- put small value (<< stab.limit ~ 2.e15) only to test biharmonic-diffusivity:
42 diffK4T=1.E11,
43 diffK4S=1.E11,
44 #-
45 rigidLid=.FALSE.,
46 implicitFreeSurface=.TRUE.,
47 eosType='JMD95Z',
48 readBinaryPrec=32,
49 writeBinaryPrec=32,
50 saltStepping=.TRUE.,
51 tempStepping=.TRUE.,
52 momStepping=.TRUE.,
53 implicitDiffusion=.TRUE.,
54 implicitViscosity=.TRUE.,
55 staggerTimeStep=.TRUE.,
56 tempAdvScheme =7,
57 saltAdvScheme =7,
58 convertFW2Salt=-1,
59 useSingleCpuIO=.TRUE.,
60 allowFreezing=.FALSE.,
61 HeatCapacity_Cp = 3986.D0,
62 gravity = 9.8156,
63 rhoNil = 1027.D0,
64 rhoConstFresh = 999.8,
65 &
66
67 # Elliptic solver parameters
68 #
69 # cg2dMaxIters - Maximum number of 2d solver iterations
70 # cg2dTargetResidual - Solver target residual
71 #
72 &PARM02
73 cg2dMaxIters=500,
74 cg2dTargetResidual=1.E-12,
75 &
76
77 # Time stepping parameters
78 #
79 # startTime - Integration starting time (s)
80 # endTime - Integration ending time (s)
81 # deltaTMom - Timestep for momemtum equations (s)
82 # deltaTtracer - Tracer timestep (s)
83 # deltaTClock - Timestep used as model "clock" (s)
84 # abEps - Adams-Bashforth stabilising factor
85 # pChkPtFreq - Frequency of permanent check pointing (s)
86 # chkPtFreq - Frequency of rolling check pointing (s)
87 # dumpFreq - Frequency at which model state is stored (s)
88 # tauThetaClimRelax - Relaxation to climatology time scale (s)
89 # tauSaltClimRelax - Relaxation to climatology time scale (s)
90 #
91 &PARM03
92 startTime=0.0,
93 #- with current input files, cannot run longer than 10 hours
94 endTime=21600.,
95 deltaTmom=3600.0,
96 deltaTtracer=3600.0,
97 deltaTClock =3600.0,
98 cAdjFreq=0.,
99 abEps=0.1,
100 forcing_In_AB = .FALSE.,
101 pChkptFreq=3600000.,
102 dumpFreq = 86400.,
103 monitorFreq=1.,
104 &
105
106 # Gridding parameters
107 #
108 # usingSphericalPolarGrid - On/Off flag for spherical polar coordinates
109 # delX - Zonal grid spacing (degrees)
110 # delY - Meridional grid spacing (degrees)
111 # delZ - Vertical grid spacing (m)
112 # ygOrigin - Southern boundary latitude (degrees)
113 #
114 &PARM04
115 usingSphericalPolarGrid=.TRUE.,
116 delX=10*2.E0,
117 delY= 8*2.E0,
118 delZ= 10., 10., 15., 20., 20., 25., 35., 50., 75.,
119 100., 150., 200., 275., 350., 415., 450.,
120 500., 500., 500., 500., 500., 500., 500.,
121 ygOrigin=46.,
122 xgOrigin=280.,
123 rSphere = 6371.D3,
124 &
125
126 # Input datasets
127 #
128 # bathyFile - File containing bathymetry
129 # hydrogThetaFile - File containing initial potential temperature data
130 # hydrogSaltFile - File containing initial salinity data
131 # zonalWindFile - File containing zonal wind data
132 # meridWindFile - File containing meridional wind data
133 # thetaClimFile - File containing theta climatology used for relaxation
134 # saltClimFile - File containing salt climatology used for relaxation
135 #
136 &PARM05
137 bathyFile = 'bathy.seaice_obcs',
138 hydrogThetaFile = 'T.seaice_obcs',
139 hydrogSaltFile = 'S.seaice_obcs',
140 uVelInitFile = 'U.seaice_obcs',
141 vVelInitFile = 'V.seaice_obcs',
142 &

  ViewVC Help
Powered by ViewVC 1.1.22