Re: Converting v[3,N] to vx[N], vy[N], vz[N] without additional memory [message #19885] |
Tue, 02 May 2000 00:00 |
colinr
Messages: 30 Registered: July 1999
|
Member |
|
|
On Mon, 01 May 2000 11:21:52 -0700,
Ricardo Fonseca <zamb@physics.ucla.edu> wrote:
> Hi everyone
>
> I need to transform an array v[3,Nx,Ny,Nz] to three arrays vx[Nx,Ny,Nz],
> vy[Nx,Ny,Nz], vz[Nx,Ny,Nz] for use with the interpolate function. The simple
> way to do it is
>
> vx = reform(v[0,*,*,*])
>
> and so on, and then freeing the memory used by v. But the problem is v is a
> very large (>180 Mb) array so the additional memory needed to do this
> becomes a problem. Is there another way of doing this, that doesn't require
> additional memory?
>
> Thanks in advance, Ricardo
A dirty way would be to write out each of v(0,*,*,*) to three separate
files, free the memory from v, and then read them in separately as vx,
vy and vz.
--
Colin Rosenthal
Astrophysics Institute
University of Oslo
|
|
|