[MITgcm-support] Large runoff causes negative salinity
Mattia Almansi
mattia.almansi at jhu.edu
Sun Jun 19 14:56:25 EDT 2016
Hi Ryan,
I'm currently using the advection scheme 33.
Attached are my data and data.exf files.
Thanks,
Mattia
- Here is the contents of the data file attached:
#
# ******************************
# Model parameters
# Continuous equation parameters
# ******************************
&PARM01
viscAz =1.E-5,
diffK4T=0.0,
diffKzT=1.E-5,
diffK4S=0.0,
diffKzS=1.E-5,
viscAhGridMax=0.25,
viscC2Leith=0.0,
viscAhGrid=0.0,
useAnisotropicViscAGridMax=.TRUE.,
viscA4GridMax=0.25,
viscC4Leith=1.0,
viscA4Grid=0.0,
no_slip_sides=.TRUE.,
no_slip_bottom=.TRUE.,
bottomDragQuadratic = 0.001,
rigidLid=.FALSE.,
implicitFreeSurface=.TRUE.,
#--- z* - begin
nonlinFreeSurf = 4,
exactConserv = .TRUE.,
select_rStar = 2,
hFacInf = 0.01,
hFacSup = 2.5,
#--- z* - end
#--- - debugMode in eedata
debugLevel = 3,
eosType='JMD95Z',
readBinaryPrec=64,
writeBinaryPrec=64,
useSingleCpuIo=.TRUE.,
globalFiles=.TRUE.,
hFacMin=0.1,
hFacMinDz=5.,
allowFreezing=.FALSE.,
HeatCapacity_Cp = 3986.D0,
gravity = 9.8156,
gBaro = 9.8156,
rhoNil = 1027.D0,
implicitDiffusion=.TRUE.,
implicitViscosity=.TRUE.,
vectorInvariantMomentum=.TRUE.,
staggerTimeStep=.TRUE.,
tempAdvScheme=33,
saltAdvScheme=33,
MultiDimAdvection=.TRUE.,
&
# **************************
# Elliptic solver parameters
# **************************
&PARM02
cg2dMaxIters=1000,
cg2dTargetResidual=1.E-13,
&
# ************************
# Time stepping parameters
# ************************
&PARM03
nIter0=0,
endTime=259200.0,
deltaT=30.0,
tauCD=0.,
cAdjFreq=0.,
abEps=0.1,
forcing_In_AB=.FALSE.,
chkptFreq=604800.,
pChkptFreq=864000.0,
dumpFreq=21600.,
diagFreq=0.,
monitorFreq=3600.,
&
# *******************
# Gridding parameters
# *******************
&PARM04
usingCartesianGrid=.FALSE.,
usingSphericalPolarGrid=.TRUE.,
delZ= 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12,
13, 14, 203*15,
ygOrigin=+56.79340,
xgOrigin=-46.95949,
delY= [see attached file]
delX= [see attached file]
rSphere = 6371.D3,
&
# **********
# Data Files
# **********
&PARM05
bathyFile='Bat_irm2kmGlacAuto.bin',
hydrogThetaFile='Temp_Initial_irm2kmGlacAuto.bin',
hydrogSaltFile='Salt_Initial_irm2kmGlacAuto.bin',
uVelInitFile='Uvel_Initial_irm2kmGlacAuto.bin',
vVelInitFile='Vvel_Initial_irm2kmGlacAuto.bin',
pSurfInitFile='Eta_Initial_irm2kmGlacAuto.bin',
# set to false to avoid checking for theta=0 and salt=0 when we start with
sea ice temp and salt
checkIniTemp=.TRUE.,
checkIniSalt=.TRUE.,
&
# Open boundaries
&PARM06
&
# KPP Mixing scheme
&PARM07
&
- Here is the contents of the data.exf file attached:
#
# *********************
# External Forcing Data
# *********************
&EXF_NML_01
#
useExfCheckRange = .TRUE.,
exf_iprec = 64,
exf_yftype = 'RL',
&
#
# *********************
&EXF_NML_02
# format
# startdate1=YYYYMMDD
# startdate2=HHMMSS
# *********************
#
atempstartdate1 = 20061231,
atempstartdate2 = 000000,
atempperiod = 10800.0,
#
aqhstartdate1 = 20061231,
aqhstartdate2 = 000000,
aqhperiod = 10800.0,
#
evapstartdate1 = 20061231,
evapstartdate2 = 013000,
evapperiod = 10800.0,
#
precipstartdate1 = 20061231,
precipstartdate2 = 013000,
precipperiod = 10800.0,
#
uwindstartdate1 = 20061231,
uwindstartdate2 = 000000,
uwindperiod = 10800.0,
#
vwindstartdate1 = 20061231,
vwindstartdate2 = 000000,
vwindperiod = 10800.0,
#
# swfluxstartdate1 = 20061231,
# swfluxstartdate2 = 013000,
# swfluxperiod = 10800.0,
#
# lwfluxstartdate1 = 20061231,
# lwfluxstartdate2 = 013000,
# lwfluxperiod = 10800.0,
#
swdownstartdate1 = 20061231,
swdownstartdate2 = 013000,
swdownperiod = 10800.0,
#
lwdownstartdate1 = 20061231,
lwdownstartdate2 = 013000,
lwdownperiod = 10800.0,
#
climsststartdate1 = 20061231,
climsststartdate2 = 120000,
climsstperiod = 86400.0,
# restoring timescale for sea sur. temperature
climsstTauRelax = 86400.0,
#
runoffstartdate1 = 20061231,
runoffstartdate2 = 000000,
runoffperiod = 86400.0,
#
atempfile = 'atemp_EraInterim3h_irm2kmGlacAuto_20061231_20070112.bin',
aqhfile = 'sphum_EraInterim3h_irm2kmGlacAuto_20061231_20070112.bin',
uwindfile = 'uwind_EraInterim3h_irm2kmGlacAuto_20061231_20070112.bin',
vwindfile = 'vwind_EraInterim3h_irm2kmGlacAuto_20061231_20070112.bin',
evapfile = 'evapr_EraInterim3h_irm2kmGlacAuto_20061231_20070112.bin',
precipfile = 'preci_EraInterim3h_irm2kmGlacAuto_20061231_20070112.bin',
# lwfluxfile = 'nlwfl_EraInterim3h_irm2kmGlacAuto_20061231_20070112.bin',
# swfluxfile = 'nswfl_EraInterim3h_irm2kmGlacAuto_20061231_20070112.bin',
lwdownfile = 'dlwfl_EraInterim3h_irm2kmGlacAuto_20061231_20070112.bin',
swdownfile = 'dswfl_EraInterim3h_irm2kmGlacAuto_20061231_20070112.bin',
climsstfile = 'OstiaDailySst_irm2kmGlacAuto_20061231_20070112.bin',
runoffFile = 'runoff_Glacier_irm2kmGlacAuto_20061231_20070112.bin',
runoftempFile = 'runoftemp_Glacier_irm2kmGlacAuto_20061231_20070112.bin',
&
#
# *********************
&EXF_NML_03
&
#
# *********************
&EXF_NML_04
&
#
# *********************
&EXF_NML_OBCS
obcsNstartdate1 = 20061231,
obcsNstartdate2 = 000000,
obcsNperiod = 86400.0,
#
obcsSstartdate1 = 20061231,
obcsSstartdate2 = 000000,
obcsSperiod = 86400.0,
#
obcsEstartdate1 = 20061231,
obcsEstartdate2 = 000000,
obcsEperiod = 86400.0,
#
obcsWstartdate1 = 20061231,
obcsWstartdate2 = 000000,
obcsWperiod = 86400.0,
#
##### SEA ICE OBCS
#
siobNstartdate1 = 20061215,
siobNstartdate2 = 000000,
siobNperiod = 2635200.0,
#
siobSstartdate1 = 20061215,
siobSstartdate2 = 000000,
siobSperiod = 2635200.0,
#
siobEstartdate1 = 20061215,
siobEstartdate2 = 000000,
siobEperiod = 2635200.0,
#
siobWstartdate1 = 20061215,
siobWstartdate2 = 000000,
siobWperiod = 2635200.0,
&
2016-06-19 13:18 GMT-04:00 Ryan Abernathey <ryan.abernathey at gmail.com>:
> Hi Mattia,
>
> This could be related to your advection scheme. Please provide your full
> data and data.exf files in order to help people diagnose the problem.
>
> -Ryan
>
> On Sun, Jun 19, 2016 at 1:02 PM, Mattia Almansi <mattia.almansi at jhu.edu>
> wrote:
>
>> Dear all,
>>
>> I'm setting up a configuration which includes part of Greenland's East
>> coast.
>> I'm having some issues with the runoff: the salinity output field is
>> negative where the runoff is large. There are some large runoff values
>> (e.g. 1.e-4 m/s into the fjords) because I want to include both glacier
>> runoff and solid ice discharge.
>> Is there any way to avoid this problem?
>>
>> Thanks in advance,
>> Mattia
>>
>> --
>> Mattia Almansi,
>> Ph.D. Student,
>> Johns Hopkins University,
>> Department of Earth and Planetary Sciences,
>> 315 Olin Hall - 3400 N. Charles Street, Baltimore, MD 21218
>>
>> e-mail: mattia.almansi at jhu.edu
>>
>> _______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>>
>>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>
>
--
Mattia Almansi,
Ph.D. Student,
Johns Hopkins University,
Department of Earth and Planetary Sciences,
315 Olin Hall - 3400 N. Charles Street, Baltimore, MD 21218
e-mail: mattia.almansi at jhu.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20160619/641d35c7/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: data.exf
Type: application/octet-stream
Size: 3247 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20160619/641d35c7/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: data
Type: application/octet-stream
Size: 9010 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20160619/641d35c7/attachment-0003.obj>
More information about the MITgcm-support
mailing list