[MITgcm-support] Model/DIC_forcing time mismatch

[Nicholas Rogan]

Hi all,

I’m having problems with some new forcing/input files I have included in the DIC package. I am using a 1D setup and have spun up my 500m water column for ten years to get a repeating cycle. I have then implemented my perturbations but it seems that the timing of the input files is mismatched relative to the model clock. I am starting my experiment runs from a pickup file that was created following my spin up. I have a perturbation that should be implemented on Day 1 of the new experiment. This is done with a forcing file 1x366 long, i.e. a value for each day of the year. My perturbation is in the new file from days 1-40. But this is not seen in the model run until the end of year 1, as if the timing has shifted over the ten year spin up. 

Is there a way I can run the model from a pickup AND reset the model clock to zero so that I can run it as if it’s the first year?

Alternatively, can anyone give me guidance on my error during my ten year spin up? My DIC forcing period is set to 86400 (Daily) and the model time mode is set to the gregorian calendar, so I thought it should be ok. 

Thanks. I’ll include my data(.dic) files below. Thanks!

Nick

# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
&PARM01
tRef=23*20.,
sRef=23*35.,
viscA4=0.,
viscAh=5.E5,
diffKhT=0.E3,
diffKhS=0.E3,
viscAz=1.E-3,
diffKzT=5.E-5,
diffKzS=5.E-5,
rotationPeriod=86400.,
# beta=1.E-11,
# f0=1.e-4,
tAlpha=2.E-4,
sBeta =7.4E-4,
gravity=9.81,
gBaro=9.81,
rhonil=1035.,
rigidLid=.FALSE.,
implicitFreeSurface=.TRUE.,
eosType='JMD95Z',
implicitDiffusion=.TRUE.,
implicitViscosity=.TRUE.,
ivdc_kappa=0.,
tempAdvection   =.TRUE.,
saltAdvection   =.TRUE.,
momAdvection    =.FALSE.,
tempAdvScheme       = 2,
saltAdvScheme       = 2,
tempStepping        = .TRUE.,
saltStepping        = .TRUE.,
multiDimAdvection   = .TRUE.,
allowFreezing=.TRUE.,
useCDscheme=.FALSE.,
useSingleCpuIO = .TRUE.,
&

# Elliptic solver parameters
&PARM02
cg2dMaxIters=1000,
cg2dTargetResidual=1.E-13,
&

# Time stepping parameters
&PARM03
pickupStrictlyMatch = .FALSE.,
# startTime = 0.,
nIter0 = 87600,
# nIter0 = 0,
# nTimeSteps = 20, 1yr =720 (43200s = 12hrs * 360 d/yr)
nTimeSteps =8040,
deltaTmom = 900.,
# tauCD =     321428.,
deltaTtracer= 3600.,
deltaTClock = 3600.,
abEps = 0.1,
cAdjFreq = 0,
abEps = 0.1,
tracForcingOutAB=1,
# These are 1 year averages, followed by monthly 2592000
#    There are also 5 or 1 day : 432000 / 86400
 pChkptFreq = 432000.,
#  chkptFreq = 432000.,
 dumpFreq = 21600.,
#  taveFreq = 43200.,
#  pChkptFreq = 2592000.,
#  chkptFreq = 2592000.,
#  dumpFreq = 2592000.
#  taveFreq = 2592000.,
#  pChkptFreq = 86400.,
#  chkptFreq = 86400.,
#  dumpFreq = 86400.,
#  taveFreq = 86400.,
monitorFreq = 1.,
tauThetaClimRelax = 5184000.0,
tauSaltClimRelax = 7776000.0,
periodicExternalForcing=.TRUE.,
externForcingPeriod = 2592000.,
externForcingCycle = 31104000.,
&

# Gridding parameters
&PARM04
usingCartesianGrid=.TRUE.,
dXspacing=5000.,
dYspacing=5000.,
# usingSphericalPolarGrid=.TRUE.,
# delX=1*2.E0,
# delY=1*2.E0,
delZ   = 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,
# ygOrigin=65.,
# xgOrigin=280.,
# rSphere = 6371.D3,
&

# Input datasets
&PARM05
bathyFile = 'new_bathy_1x1_500m_testpool',
hydrogThetaFile = 'lev_clim_temp.bin',
hydrogSaltFile = 'lev_clim_salt.bin',
zonalWindFile ='tren_taux.bin',
meridWindFile = 'tren_tauy.bin',
thetaClimFile = 'lev_monthly_temp.bin',
saltClimFile = 'lev_monthly_salt.bin',
surfQFile = 'shi_qnet.bin',
EmPmRFile = 'shi_empmr_year.bin',
the_run_name = 'vII_y',
&

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++


&DIC_FORCING
 DIC_iceFile='day_fice.bin',
 DIC_windFile='day_tren_speed.bin',
 DIC_silicaFile='',
 DIC_ironFile='day_mah_somtimind_molfem2s.bin',
 DIC_sedFile='sedsup_fivfivzer.bin',
 DIC_hydFile='hydro_sup_molHem2s.bin',
#  DIC_sedFile='sedsup_sevninzer.bin',
 DIC_perturbFile='mon1perturb.bin',
 DIC_int1=3,
 DIC_forcingPeriod = 86400,
&