cube13: 1992-2003, same as cube12 but all albedos increased by 0.1
==================================================================
# SEAICE_dryIceAlb = 0.75 _d 0
# SEAICE_wetIceAlb = 0.66 _d 0
# SEAICE_drySnowAlb = 0.84 _d 0
# SEAICE_wetSnowAlb = 0.7 _d 0
# SEAICE_waterAlbedo = 0.15
SEAICE_dryIceAlb = 0.85
SEAICE_wetIceAlb = 0.76
SEAICE_drySnowAlb = 0.94
SEAICE_wetSnowAlb = 0.8
SEAICE_waterAlbedo = 0.25
Note that cube13 crashed in 2001, so starting in 2001 switched to:
viscA4GridMin=3.e-4, viscA4GridMax=4.e-3, viscC4Leith=0.2.
Driven by NCEP reanalysis, 6-hourly, atmospheric state.
Uses bulk formulae and dynamic/thermodynamic sea-ice model.
Includes sea-surface salinity relaxation.
Output files are in
lou.nas.nasa.gov:/u/bcheng/c6ncp10
Example matlab scripts for reading and visualizing data are in
lou.nas.nasa.gov:/ecco/menemenl/cube12/matlab/lookat.m.
All output fields are computed at the tracer grid points (XC,YC),
with following exceptions.
uVeltave, FUtave, uFluxtave, UUtave: Southwest C-grid U position (XG,YC)
vVeltave, FVtave, vFluxtave, VVtave: Southwest C-grid V position (XC,YG)
UICEtave, VICEtave : South-West B-grid U,V point (XG,YG)
Available output files are 3-day averages:
------------------------------------------
ETAtave.*.data : sea-surface height (m)
uVeltave.*.data : uVel (m/s)
vVeltave.*.data : vVel (m/s)
wVeltave.*.data : wVel (m/s)
Ttave.*.data : theta (deg C)
Stave.*.data : salinity (PSU)
AREAtave.*.data : sea-ice area
HEFFtave.*.data : sea-ice thickness
UICEtave.*.data : sea-ice uVel
VICEtave.*.data : sea-ice vVel
EmPmRtave.*.data : following are surface flux diagnostics,
QNETtave.*.data : some redundant, some not, will describe
QSWtave.*.data : in detail later ...
FUtave.*.data :
FVtave.*.data :
sFluxtave.*.data :
tFluxtave.*.data :
uFluxtave.*.data :
vFluxtave.*.data :
===============================================================
== input data file
===============================================================
# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
&PARM01
tRef = 3*23.,3*22.,21.,2*20.,19.,2*18.,17.,2*16.,15.,14.,13.,
12.,11.,2*9.,8.,7.,2*6.,2*5.,3*4.,3*3.,4*2.,12*1.,
sRef = 50*34.5,
viscAr= 1.77e-5,
no_slip_sides=.FALSE.,
no_slip_bottom=.TRUE.,
diffKhT=0.,
diffK4T=0.,
diffKrT=1.51e-5,
diffKhS=0.,
diffK4S=0.,
diffKrS=1.51e-5,
tAlpha=2.E-4,
sBeta =7.4E-4,
gravity=9.81,
eosType='JMD95Z',
readBinaryPrec=32,
writeBinaryPrec=32,
useSingleCPUio=.TRUE.
globalFiles=.TRUE.
hFacMinDr=50.,
hFacMin=0.3,
hFacInf=0.1,
hFacSup=5.,
# z* - begin
select_rStar=2,
nonlinFreeSurf=4,
# * - end
implicitDiffusion=.true.,
implicitViscosity=.true.,
allowFreezing=.false.,
viscAh=0.,
# viscAhMax=unset
viscAhGrid=0.,
viscC2Leith=0.,
viscA4=0.,
# viscA4Max=0.,
viscA4Grid=0.,
viscA4GridMin=3.e-4,
viscA4GridMax=4.e-3,
viscC4Leith=0.2,
bottomDragQuadratic = 0.002,
tempAdvScheme=33,
saltAdvScheme=33,
StaggerTimeStep=.true.,
multiDimAdvection=.true.,
vectorInvariantMomentum=.TRUE.,
rigidLid=.FALSE.,
implicitFreeSurface=.TRUE.,
exactConserv=.TRUE.,
debuglevel=0
&
# Elliptic solver parameters
&PARM02
cg2dMaxIters=100,
cg2dTargetResidual=1.E-5,
&
# Time stepping parameters
&PARM03
niter0=315360,
endTime=379728000,
abEps=0.1,
deltaT=1200,
cAdjFreq=0.,
pChkptFreq=379728000,
chkptFreq=2592000
tavefreq=259200
tauSaltClimRelax = 3842380.,
monitorFreq=86400
forcing_In_AB = .FALSE.
&
# Gridding parameters
&PARM04
usingCartesianGrid=.FALSE.,
usingSphericalPolarGrid=.FALSE.,
usingCurvilinearGrid=.TRUE.,
delR = 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.01,
10.03, 10.11, 10.32, 10.80, 11.76, 13.42, 16.04 , 19.82, 24.85,
31.10, 38.42, 46.50, 55.00, 63.50, 71.58, 78.90, 85.15, 90.18,
93.96, 96.58, 98.25, 99.25,100.01,101.33,104.56,111.33,122.83,
139.09,158.94,180.83,203.55,226.50,249.50,272.50,295.50,318.50,
341.50,364.50,387.50,410.50,433.50,456.50,
&
# Input datasets
&PARM05
bathyFile ='ETOPO2_510x6x510_ver5e_filled.bin',
hydrogThetaFile='LEVT01_JAN_510x6x510x50.bin',
hydrogSaltFile ='LEVS01_JAN_510x6x510x50.bin',
&