Changeset 21170 for issm/trunk-jpl/test/Par/SquareShelf2.py

Timestamp:

08/25/16 21:53:46 (9 years ago)

Author:

agscott1

Message:

BUG: Fixed archread incorrectly reading the order of a matrix, and fixed Par python files to correctly read arch files

File:

• issm/trunk-jpl/test/Par/SquareShelf2.py (modified) (3 diffs)

issm/trunk-jpl/test/Par/SquareShelf2.py ¶

@@ -1 +1 @@
 import os.path
 import inspect
-import netCDF4
+from arch import *
 import numpy
 from verbose import verbose
@@ -21 +21 @@
 
 #Initial velocity and pressure
-iVelF = netCDF4.Dataset('../Data/SquareShelf.nc','r')
-
-x=numpy.reshape(iVelF.variables['x'][:],(-1))
-y=numpy.reshape(iVelF.variables['y'][:],(-1))
-vx=iVelF.variables['vx'][:]
-vy=iVelF.variables['vy'][:]
-index=iVelF.variables['index'][:].astype(int)
-
+x         = numpy.array(archread('../Data/SquareShelf.arch','x'))
+y         = numpy.array(archread('../Data/SquareShelf.arch','y'))
+vx        = numpy.array(archread('../Data/SquareShelf.arch','vx'));
+vy        = numpy.array(archread('../Data/SquareShelf.arch','vy'));
+index     = archread('../Data/SquareShelf.arch','index').astype(int);
 #dbg - begin
 # #print 'vars in SquareShelf.nc:'
@@ -34 +31 @@
 # #     print v
 #dbg - end 
-
-iVelF.close()
 
 [md.initialization.vx]=InterpFromMeshToMesh2d(index,x,y,vx,md.mesh.x,md.mesh.y)