[MITgcm-support] cg3d problem
Nicolas Grisouard
nicolas.grisouard at grenoble-inp.fr
Fri Nov 14 06:38:30 EST 2008
Dear all,
I am using a configuration, in which the code is slowly converging and
am wondering if there is something to do about it. I tried to run a test
with the -ftrace option activated, the resulting file is attached.
To sum up, the cg3d routine is using 75% of the time, which is maybe
usual, I don't know. This is my first question : I am using a 2D
configuration, in this case the size of my grid is nx*ny*nz =
1000*1*223. Why is there a 3d process involved here? Is it the fact that
the boundary conditions north and south are periodic? (OBCs are on)
Now, is there a way to avoid such a long convergence of the 3D solver?
Actually, in the data file, I set these lines:
# Elliptic solver parameters
&PARM02
cg2dMaxIters=10000,
cg2dTargetResidual=1.E-13,
cg3dMaxIters=10000,
cg3dTargetResidual=1.E-13,
&
And in spite of these high values of MaxIters, the cg3d routine does not
converge, as can be read in the STDOUT file :
(PID.TID 0000.0001) cg2d_init_res =
1.40478312751614E+00
(PID.TID 0000.0001) cg2d_iters = 1778
(PID.TID 0000.0001) cg2d_res =
9.75881835078007E-14
cg3d: Sum(rhs),rhsMax = -1.12229402786101E-05 1.57274284756183E-12
(PID.TID 0000.0001) cg3d_init_res =
4.88713802210862E+00
(PID.TID 0000.0001) cg3d_iters = 10000
(PID.TID 0000.0001) cg3d_res =
2.38471462841502E-05
I hope there is something to do about it, instead the cost of these
simulations will get prohibitive.
Thanks for the help,
Nicolas
--
Nicolas GRISOUARD
PhD Student - ERES research group
Laboratoire des Ecoulements Geophysiques et Industriels
BP 53
38041 Grenoble cedex 9 France
tel : +33 (0)476 825 037 - fax : +33 (0)476 827 022
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: mar_tst02.o16716
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20081114/92c856f2/attachment.el>
More information about the MITgcm-support
mailing list