Reading a restart file takes ~50s on my system, which is the main bottleneck of the FTLE computation. (Compared to this, reading the file in Numpy takes about 2 seconds).
The relevant code:
|
diskfield=OPEN(field_name) |
|
WITH diskfield: |
|
LOOP FOR iV=0 TO 2 AND iz=0 TO nzd-1 |
|
ARRAY(0..ny, 0..nxd-1) OF REAL buf = fieldimage(iV,iz,*,*); |
|
DO V(ix,iz,iy,iV) = buf(iy,ix) FOR iy=0 TO ny AND ix=0 TO nxd-1 |
|
REPEAT |
|
CLOSE diskfield |
Reading a restart file takes ~50s on my system, which is the main bottleneck of the FTLE computation. (Compared to this, reading the file in Numpy takes about 2 seconds).
The relevant code:
activeBarriers/aftle/xCompactInterface.cpl
Lines 30 to 36 in f76bc1c