gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

PYTHON: Add refinement 'level' to be read as dataset.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2019-02-18 11:27:47 -03:00
parent 81b94a8009
commit 0ceca17a95
1 changed files with 6 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
python/amun.py
View File

@ -236,6 +236,7 @@ def amun_dataset(fname, vname, shrink = 1, progress = False):
# add derived variables if possible # add derived variables if possible
# #
variables.append('level')
if 'velx' in variables and 'vely' in variables and 'velz' in variables: if 'velx' in variables and 'vely' in variables and 'velz' in variables:
variables.append('velo') variables.append('velo')
variables.append('divv') variables.append('divv')
@ -311,7 +312,11 @@ def amun_dataset(fname, vname, shrink = 1, progress = False):
dz = amun_coordinate(lname, 'dz') dz = amun_coordinate(lname, 'dz')
f = h5.File(lname, 'r') f = h5.File(lname, 'r')
g = f['variables'] g = f['variables']
if vname == 'velo': if vname == 'level':
dataset = np.zeros(g[variables[0]].shape)
for l in range(dblocks):
dataset[:,:,:,l] = levels[l]
elif vname == 'velo':
dataset = np.sqrt(g['velx'][:,:,:,:]**2 \ dataset = np.sqrt(g['velx'][:,:,:,:]**2 \
+ g['vely'][:,:,:,:]**2 \ + g['vely'][:,:,:,:]**2 \
+ g['velz'][:,:,:,:]**2) + g['velz'][:,:,:,:]**2)