[Tgcmgroup] tiegcm1 advec bug
Ben Foster
Ben Foster <foster@hao.ucar.edu>
Tue, 28 Oct 2003 15:04:42 -0700 (MST)
Hi tgcmgroup:
A bug has been found in tiegcm1 advec.F in which 2nd order finite
differencing in latitude is missing (it is doing only 4th order in
latitude). The corrected version is /home/tgcm/tiegcm1/modsrc/advec.F.
Compare this file with /home/tgcm/tiegcm1/advec.F to see the fix.
Calls to advec are made by dt, duv, and for no and n4s composition.
Differences in the results are significant. Please use the correct
version of advec.F for new runs.
This affects tiegcm only. Time-gcm (tgcm24) does not have this bug.
--Ben
-----------------------------------------------------------------------
Ben Foster High Altitude Observatory (HAO)
foster@ucar.edu phone: 303-497-1595 fax: 303-497-1589
Nat. Center for Atmos. Res. P.O. Box 3000 Boulder CO 80307 USA
-----------------------------------------------------------------------