|
|
Blitz Support : |
From: John Bray (home_at_[hidden])
Date: 2005-02-16 15:59:38
Thanks for all the advice on defining null arrays and then resizing
them. That worked fine, and I've now replicated my MPP environment with
bounds swapping routines across processors.
I would like to have Fortran-like array ordering though, but it seems
that resize cannot take an ordering argument, so
GeneralArrayStorage<3> storage;
storage.ordering() = firstDim, secondDim, thirdDim;
a.resize(3,4,5,storage);
doesn't work, and I can't define the storage order of the array in
Array<float,3> a;
The reason I want to do this is I have a weather forecasting domain,
where x,y,z correspond to longitude,latitude, vertical level, and most
of the calculations are in the lat-long plane, so I'd like the longitude
values contiguous, not the vertical levels.
John