[MITgcm-support] Re: single-layer configuration problem update
Patrick Heimbach
heimbach at MIT.EDU
Tue Aug 30 14:01:45 EDT 2005
Hi David,
no, if you're using exf,
the parameters mentioned below are ignored,
and data.exf, data.exf_clim get used.
The parameter you're looking for replacing externForcingCycle
is repeatPeriod
-p.
> -----Original Message-----
> From: mitgcm-support-bounces at mitgcm.org
> [mailto:mitgcm-support-bounces at mitgcm.org] On Behalf Of David Wang
> Sent: Tuesday, August 30, 2005 1:46 PM
> To: mitgcm-support at mitgcm.org
> Subject: Re: [MITgcm-support] Re: single-layer
> configuration problem update
>
>
> Hey all,
>
> I try to catch up this thread and throw out my question.
>
> If I use exf package to specify my wind stress forcing (monthly
> climatology) and restoring SST/SSS (also monthly
> climatology), do the
> following parameters in data namelist still count? I've not
> found the
> counterpart of externForcingCycle (supposed to be 1 year in
> my case) in
> data.exf and data.exf_clim namelists.
>
> periodicExternalForcing
> externForcingPeriod
> externForcingCycle
> tauThetaClimRelax
> tauSaltClimRelax
>
> Thanks,
> - D.W.
>
> Patrick Heimbach wrote:
>
> >Sergey,
> >
> >several things:
> >1. I thought you were using the SINGLE_LAYER_MODE flag(?)
> >2. The std. output has NAN's in cg2d after a few timesteps,
> > but you're saying the model doesn't explode.
> > Can you clarify?
> >3. The alternative, simplere forcing would be via
> > S/R external_fields_load.F
> > and not using exf (i.e. remove 'exf' from packages.conf)
> > An example is in verification/global_ocean.90x40x15/code/
> >
> >-Patrick
> >
> >
> >
> > > -----Original Message-----
> > > From: mitgcm-support-bounces at mitgcm.org
> > > [mailto:mitgcm-support-bounces at mitgcm.org] On Behalf Of
> > > Sergey Vinogradov
> > > Sent: Tuesday, August 30, 2005 10:50 AM
> > > To: mitgcm-support at mitgcm.org
> > > Cc: Rui Ponte
> > > Subject: Re: [MITgcm-support] Re: single-layer
> > > configuration problem update
> > >
> > >
> > > Martin:
> > > I don't know of any other way than EXF to introduce the
> > > variable forcing
> > > into the mitgcm. The resulting field is not exploding
> > > though, it stays
> > > equal to i.c.
> > > Anyway, SINGLE_LAYER_MODE is probably what I was
> looking for in a
> > > context of my experiments.
> > > Thank you,
> > > Sergey
> > >
> > > Martin Losch wrote:
> > >
> > > > Sergey,
> > > > that's really puzzling. Have you tried the same
> thing without the
> > > > EXF-package? (I have had unexpected behavior with this
> > > package before).
> > > >
> > > > However, if --- as you suspect --- no forcing fields
> > > introduced by EXF
> > > > is not "felt" by the model, then it wouldn't explode, but
> > > would just
> > > > keep quiet without moving at all. So I think, that on the
> > > contrary,
> > > > the model feels too much of the forcing; or
> something completely
> > > > different is going on.
> > > >
> > > > I can't remember your configuration, and obviously, I am
> > > too stupid to
> > > > search the MITgcm-support email list effectively so I
> > > didn't find your
> > > > previous emails, but were you using the CPP-flag
> > > SINGLE_LAYER_MODE
> > > > (kicks out most of thermodynamics)? However, even if you
> > > were using
> > > > it, I don't see where it should interfere with the
> EXF package.
> > > >
> > > > Martin
> > > >
> > > > On Aug 30, 2005, at 3:49 PM, Sergey Vinogradov wrote:
> > > >
> > > >> Hi all,
> > > >>
> > > >> I just want to update you on my quest in solving the
> > > problem with
> > > >> single-layer run under dynamic forcing (and maybe get
> > > some useful
> > > >> feedback from you).
> > > >>
> > > >> The issue was erroneous computation of ocean dynamics
> > > (sea-level,
> > > >> bottom pressure) in a single layer configuration with a
> > > use of EXF
> > > >> package interface for the time-varying atm.
> pressure loading.
> > > >>
> > > >> Appearently, it is not a problem of just pressure
> > > loading. It is far
> > > >> more general problem, as I tried similar experiment with
> > > wind stress.
> > > >> It looks like any forcing that is introduced through EXF
> > > in one-layer
> > > >> run, is not "felt" by the model.
> > > >>
> > > >> Below is the exerpt from the model output with one
> > > layer under wind
> > > >> stress at the beginning of forward integration where you
> > > can see the
> > > >> error accumulating in cg solver with each step. Similar
> > > run with 23
> > > >> layers computes reasonable fields.
> > > >>
> > > >> I'd appreciate any suggestions/comments/tips very much!
> > > >> Regards,
> > > >> Sergey
> > > >>
> > > >> (PID.TID 0000.0001) //
> > > >> =======================================================
> > > >> (PID.TID 0000.0001) // End MONITOR dynamic field statistics
> > > >> (PID.TID 0000.0001) //
> > > >> =======================================================
> > > >> (PID.TID 0000.0001) MDSREADFIELD: opening global file:
> > > >> /anvil/scratch/day07/svv/forcing/NCEP_R6_ustr_1999
> > > >> (PID.TID 0000.0001) MDSREADFIELD: opening global file:
> > > >> /anvil/scratch/day07/svv/forcing/NCEP_R6_ustr_1999
> > > >> cg2d: Sum(rhs),rhsMax = 5.83105612406154E-15
> > > 2.52789606828280E+01
> > > >> cg2d: Sum(rhs),rhsMax = 4.66381008031730E-09
> > > 2.03846096766287E+01
> > > >> cg2d: Sum(rhs),rhsMax = -7.72491057666330E-11
> > > 1.63318316358843E+02
> > > >> cg2d: Sum(rhs),rhsMax = 7.82384650122037E-10
> > > 3.94671184216402E+02
> > > >> cg2d: Sum(rhs),rhsMax = 1.84482164570718E-10
> > > 1.39786242903642E+03
> > > >> cg2d: Sum(rhs),rhsMax = -3.86795120889058E-10
> > > 5.51987500042437E+03
> > > >> cg2d: Sum(rhs),rhsMax = 5.60930842031203E-10
> > > 2.96468380272229E+04
> > > >> cg2d: Sum(rhs),rhsMax = -4.60258786649532E-11
> > > 1.03144921478204E+06
> > > >> cg2d: Sum(rhs),rhsMax = 1.93486835808804E-12
> > > 3.16417950241384E+08
> > > >> cg2d: Sum(rhs),rhsMax = -2.47236726757771E-12
> > > 2.35353664301718E+13
> > > >> cg2d: Sum(rhs),rhsMax = -3.50420466366502E-16
> > > 2.44341352738558E+23
> > > >> cg2d: Sum(rhs),rhsMax = -2.42051059455915E-16
> > > 1.68257839407874E+43
> > > >> cg2d: Sum(rhs),rhsMax = -2.02970152146811E-16
> > > 1.42078931888867E+83
> > > >> cg2d: Sum(rhs),rhsMax = -4.29685153894816E-16
> > > 9.27740314802081+162
> > > >> cg2d: Sum(rhs),rhsMax = NAN
> > > INF
> > > >> cg2d: Sum(rhs),rhsMax = NAN
> > > 0.00000000000000E+00
> > > >> cg2d: Sum(rhs),rhsMax = NAN
> > > 0.00000000000000E+00
> > > >> cg2d: Sum(rhs),rhsMax = NAN
> > > 0.00000000000000E+00
> > > >> (PID.TID 0000.0001) MDSREADFIELD: opening global file:
> > > >> /anvil/scratch/day07/svv/forcing/NCEP_R6_ustr_1999
> > > >> cg2d: Sum(rhs),rhsMax = NAN
> > > 0.00000000000000E+00
> > > >> cg2d: Sum(rhs),rhsMax = NAN
> > > 0.00000000000000E+00
> > > >>
> > > >
> > > >
> > > > _______________________________________________
> > > > 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
> > >
> >
> >_______________________________________________
> >MITgcm-support mailing list
> >MITgcm-support at mitgcm.org
> >http://mitgcm.org/mailman/listinfo/mitgcm-support
> >
> >
>
>
> --
> Daiwei (David) Wang Graduate Student
> Lamont-Doherty Earth Observatory of Columbia University
> 103B Oceanography
> 61 Route 9W, Palisades, NY 10964-8000
> Voice: (845) 365-8694
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>
More information about the MITgcm-support
mailing list