/[MITgcm]/MITgcm/verification/lab_sea/input_ad/data
ViewVC logotype

Annotation of /MITgcm/verification/lab_sea/input_ad/data

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

Revision 1.1 - (hide annotations) (download)
Fri Dec 15 18:11:59 2006 UTC (17 years, 4 months ago) by heimbach
Branch: MAIN 
Initial check-in of lab_sea adjoint
(very similar to natl_box_adjoint/code_seaice)
* used C-grid and EVP solver
* used KPP, GM/Redi
* uses inAdExact=.TRUE.
* gradient checks fail for both ADM and TLM
  due to sensitive snow melt code in seaice_grwoth (look for "JZ")
  f.d. gradients actually seem more problematic than ADM, TLM gradients;
  issue is related to forward mismatch between g77 an ifort
1 heimbach 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     # f0 - Reference coriolis parameter,
16     # south edge of f on beta plane (1/s)
17     # beta - df/dy (s^-1.m^-1)
18     # tAlpha - Linear EOS thermal expansion coefficient (1/oC)
19     # sBeta - Linear EOS haline contraction coefficient (1/ppt)
20     # gravity - Acceleration due to gravity (m/s^2)
21     # gBaro - Accel. due to gravity used in barotropic equation (m/s^2)
22     # rigidLid - Set to true to use rigid lid
23     # implicitFreeSurface - Set to true to use implicit free surface
24     # eosType - Flag for linear or polynomial equation of state
25     # momAdvection - On/Off flag for momentum self transport
26     # momViscosity - On/Off flag for momentum mixing
27     #
28     &PARM01
29     tRef= 24.0 , 23.0 , 22.0 , 21.0 , 20.0 ,
30     19.0 , 18.0 , 17.0 , 16.0 , 15.0 ,
31     14.0 , 13.0 , 12.0 , 11.0 , 10.0 ,
32     9.0 , 8.0 , 7.0 , 6.0, 5.0 ,
33     4.0 , 3.0 , 2.0 ,
34     sRef= 34.65, 34.75, 34.82, 34.87, 34.90,
35     34.90, 34.86, 34.78, 34.69, 34.60,
36     34.58, 34.62, 34.68, 34.72, 34.73,
37     34.74, 34.73, 34.73, 34.72, 34.72,
38     34.71, 34.70, 34.69,
39     no_slip_sides=.FALSE.,
40     no_slip_bottom=.TRUE.,
41     viscAz=1.93e-5,
42     viscAh=5.E4,
43     diffKhT=0.0,
44     diffKzT=1.46e-5,
45     diffKhS=0.0,
46     diffKzS=1.46e-5,
47     rotationPeriod=86400.,
48     beta=1.E-11,
49     tAlpha=2.E-4,
50     sBeta =7.4E-4,
51     rigidLid=.FALSE.,
52     implicitFreeSurface=.TRUE.,
53     eosType='JMD95Z',
54     readBinaryPrec=32,
55     writeBinaryPrec=32,
56     saltStepping=.TRUE.,
57     tempStepping=.TRUE.,
58     momStepping=.TRUE.,
59     implicitDiffusion=.TRUE.,
60     implicitViscosity=.TRUE.,
61     #globalFiles=.TRUE.,
62     #- not safe to use globalFiles in multi-processors runs; set instead useSingleCpuIO
63     useSingleCpuIO=.FALSE.,
64     allowFreezing=.FALSE.,
65     HeatCapacity_Cp = 3986.D0
66     gravity = 9.8156,
67     gBaro = 9.8156,
68     rhoNil = 1027.D0
69     useCDscheme=.TRUE.,
70     useNHMTerms=.TRUE.,
71     inAdExact=.TRUE.,
72     &
73    
74     # Elliptic solver parameters
75     #
76     # cg2dMaxIters - Maximum number of 2d solver iterations
77     # cg2dTargetResidual - Solver target residual
78     #
79     &PARM02
80     cg2dMaxIters=500,
81     cg2dTargetResidual=1.D-8
82     &
83    
84     # Time stepping parameters
85     #
86     # startTime - Integration starting time (s)
87     # endTime - Integration ending time (s)
88     # tauCD - CD scheme coupling timescale (s)
89     # deltaTMom - Timestep for momemtum equations (s)
90     # deltaTtracer - Tracer timestep (s)
91     # deltaTClock - Timestep used as model "clock" (s)
92     # abEps - Adams-Bashforth stabilising factor
93     # pChkPtFreq - Frequency of permanent check pointing (s)
94     # chkPtFreq - Frequency of rolling check pointing (s)
95     # dumpFreq - Frequency at which model state is stored (s)
96     # tauThetaClimRelax - Relaxation to climatology time scale (s)
97     # tauSaltClimRelax - Relaxation to climatology time scale (s)
98     #
99     &PARM03
100     tauCD=172800.,
101     startTime=0.0,
102     nTimeSteps=10,
103     deltaTmom=3600.0,
104     deltaTtracer=3600.0,
105     deltaTClock =3600.0,
106     cAdjFreq=0.,
107     abEps=0.1,
108     pChkptFreq=36000.,
109     chkptFreq= 0.,
110     dumpFreq = 0.,
111     tavefreq = 36000.,
112     tauSaltClimRelax = 4142330.,
113     tauThetaClimRelax = 0.,
114     monitorFreq=1.,
115     adjMonitorFreq=1.,
116     #forcing_In_AB = .FALSE.
117     tracForcingOutAB=1,
118     &
119    
120     # Gridding parameters
121     #
122     # usingSphericalPolarGrid - On/Off flag for spherical polar coordinates
123     # usingCartesianGrid - On/Off flag for selecting cartesian coordinates
124     # delX - Zonal grid spacing (degrees)
125     # delY - Meridional grid spacing (degrees)
126     # delZ - Vertical grid spacing (m)
127     # phiMin - Southern boundary latitude (degrees)
128     #
129     &PARM04
130     usingCartesianGrid=.FALSE.,
131     usingSphericalPolarGrid=.TRUE.,
132     delX=20*2.E0,
133     delY=16*2.E0,
134     delZ= 10., 10., 15., 20., 20., 25., 35., 50., 75.,
135     100., 150., 200., 275., 350., 415., 450.,
136     500., 500., 500., 500., 500., 500., 500.
137     phimin=46.,
138     thetamin=280.,
139     rSphere = 6371.D3
140     &
141    
142     # Input datasets
143     #
144     # bathyFile - File containing bathymetry
145     # hydrogThetaFile - File containing initial potential temperature data
146     # hydrogSaltFile - File containing initial salinity data
147     # zonalWindFile - File containing zonal wind data
148     # meridWindFile - File containing meridional wind data
149     # thetaClimFile - File containing theta climatology used for relaxation
150     # saltClimFile - File containing salt climatology used for relaxation
151     #
152     &PARM05
153     bathyFile = 'bathy.labsea1979',
154     hydrogThetaFile = 'LevCli_temp.labsea1979',
155     hydrogSaltFile = 'LevCli_salt.labsea1979',
156     &
  ViewVC Help
Powered by ViewVC 1.1.22