[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
-----------------------------------------------------------------------