[MITgcm-support] the floating-point exception issue

[han]

Dear all,


I am a new user of MITgcm.Recently, I've been  trying to establish an ocean simulation model for a certain area (similar to tutorial_global_oce_latlon case). 
After I used the initial bathymetry(bathyFile), temperature(hydrogThetaFile), salinity(hydrogSaltFile),velocity(uVelInitFile,vVelInitFile), wind stress(zonalWindFile,meridWindFile) files generated by myself and modified some parameters, I found that the model could run normally, but there would be floating-point exceptions, as shown in the following code


Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG
STOP NORMAL END
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG
STOP NORMAL END
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG
STOP NORMAL END
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_OVERFLOW_FLAG
STOP NORMAL END




In the STDOUT file, it can be seen that the calculation has diverged and NAN values appear, like this


(PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector =    F
cg2d: Sum(rhs),rhsMax =   1.24232198821943E-02  4.08393649685423E+00
...
cg2d: Sum(rhs),rhsMax =                    NaN                   NaN
cg2d: Sum(rhs),rhsMax =                    NaN                   NaN
cg2d: Sum(rhs),rhsMax =                    NaN                   NaN
cg2d: Sum(rhs),rhsMax =                    NaN                   NaN
cg2d: Sum(rhs),rhsMax =                    NaN                   NaN
...




I am looking to you for some guidance and advice.

Thanks very much in advance!



Best regards,
Han

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20250423/dadcb990/attachment.html>