[Mitgcm-support] curious

mitgcm-support at dev.mitgcm.org mitgcm-support at dev.mitgcm.org
Wed Jul 9 15:12:41 EDT 2003


running exp5, with the new version of cg2d.F & solve_for_pressure.F
(i.e. with write(0,..  CG2D iters, err =' in solve_for_pressure.F)
give different "Numbers" (differences of the order of truncation error) :

new version, 2nd iter :
 cg2d: Sum(rhs) =          -2.16493489801906E-15 
presious version, 2nd iter :
 cg2d: Sum(rhs) =          -3.49720252756924E-15
And also:
Binary files U.0000000010.001.001.data and ../res_sav/U.0000000010.001.001.data differ

But if I activate the 2 commented "write(0," in the new version of 
cg2d.F (-> as before) (lines 227,228  and lines 390,391)
then I get the same "Numbers" as before. 

Curious, no ? 

Jean-Michel



More information about the MITgcm-support mailing list