[MITgcm-support] time step problem for cs510?
Jean-Michel Campin
jmc at ocean.mit.edu
Wed Aug 24 22:36:39 EDT 2011
Hi Zhaomin,
May be a little bit of 'context' type information might help
(type of set-up, type of problem you want to addess, ...).
And also some numerical details about your set-up could
also be useful (in general, the usefulness of the response you
could expect from this list is largely function of the
how precise the problem is described).
Cheers,
Jean-Michel
On Wed, Aug 24, 2011 at 03:19:04PM +0100, Wang, Zhaomin wrote:
> Hi,
>
> We have one run probably stopped by numerical instability. We used cube sphere 510 configuration, and initially used deltaT=1200. When the run was stopped, we reduced deltaT to 600. When the one was stopped again, we used deltaT=300. However, the run was also stopped even for this very small time step. We did not try to reduce deltaT further, because the instability may not be caused by a too large time step. We got the following error message:
>
> ***********************
> WARNING: r*FacW < hFacInf at 8 pts : bi,bj,Thid,Iter= 1 1 1 1694086
> e.g. at i,j= 24 24 ; rStarFac,H,eta = NaN 4.281800E+02 NaN -1.058899E+00
> STOP in CALC_R_STAR : too SMALL rStarFacW !
> WARNING: r*FacS < hFacInf at 8 pts : bi,bj,Thid,Iter= 1 1 1 1694086
> e.g. at i,j= 23 25 ; rStarFac,H,eta = NaN 4.281800E+02 NaN -8.253981E-01
> STOP in CALC_R_STAR : too SMALL rStarFacS !
> WARNING: r*FacC < hFacInf at 5 pts : bi,bj,Thid,Iter= 1 1 1 1694086
> e.g. at i,j= 23 24 ; rStarFac,H,eta = NaN 4.281800E+02 NaN
> STOP in CALC_R_STAR : too SMALL rStarFacC !
> ********************
>
> Thanks.
>
> Zhaomin
>
> British Antarctic Survey
> High Cross, Madingley Road
> Cambridge, CB3 0ET
> UK
> Tel: +44(0)1223221488
> Fax: +44(0)1223221279
>
> --
> This message (and any attachments) is for the recipient only. NERC
> is subject to the Freedom of Information Act 2000 and the contents
> of this email and any reply you make may be disclosed by NERC unless
> it is exempt from release under the Act. Any material supplied to
> NERC may be stored in an electronic records management system.
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
More information about the MITgcm-support
mailing list