[MITgcm-support] Freshwater OBCS problem
spanishcastlemagik at gmail.com
Sat Dec 22 16:24:13 EST 2012
Yes, it dies actually during the first time step. I do get negative salinities. I have tried changing my OBCS (river) S=1 psu and S=10 psu with the same results. It looks like my advection scheme is creating a undershoot. I have tried advection schemes 7,33,and 77 with the same results.
cg2d starts increasing positively from the start and then becomes NaN's. cg3d looks fine.
Yes, hopefully someone with experience with river or freshwater input will help out here as I am stuck on this.
> Does it die after one time step or one save step? You can set your MonitorFreq to equal your time step to get diagnostics each step. I often set dumpFreq equal to my time step as well to save each step if I have something blow up early in the run.
> Any idea where your NaN's are? An uneducated guess would be that your salinities go slightly negative and the model doesn't accept negative salinity for some reason and NaN's those salinities. Have you tried setting S=1 psu (or 10 psu) just to see if it can run for a while? I don't think the advection scheme is up to keeping a zero salinity front strictly greater than zero.
> Good luck - hopefully someone else with experience putting a river into their model speaks up.
> Cheers, Jody
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the MITgcm-support