[Mitgcm-support] bug in gad_advection.F ?
mitgcm-support at dev.mitgcm.org
mitgcm-support at dev.mitgcm.org
Wed Jul 9 15:31:01 EDT 2003
In gad_advection there are the following lines:
IF (useCubedSphereExchange) THEN
nipass=3
ELSE
nipass=1
ENDIF
nipass=1
Is the last line (unconditional nipass=1) intended?
Also, there's an issue with the multi-dimensional advection scheme.
It's ipass loop requires storing to a new type of common blocks.
This can be done, but in the absence of an #ifdef
this common block will ALWAYS be created, thus unnecessarily
increasing the code size if multiDimAdvection is not wanted.
Patrick
_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/
Patrick Heimbach FON: +1 / 617 / 253-5259
Massachusetts Institute of Technology FAX: +1 / 617 / 253-4464
EAPS, Room 54-1518
77 Massachusetts Avenue mailto:heimbach at mit.edu
Cambridge MA 02139, U.S.A. http://www.mit.edu/~heimbach/
>>>>>>> http://web.mit.edu/cms/reconstructions <<<<<<<
More information about the MITgcm-support
mailing list