Re: [Pflexible-users] pflexible example on introduction page

[Janusz M. <jan...@gm...>]

Hello Everyone,

I have been trying to get pflexible to run for almost a day and a half now
according to the example set out in:

http://pflexible.readthedocs.org/en/latest/getting_started.html

unsuccessfully unfortuantely. Could someone please help me with this...

This is what I did..

wget http://niflheim.nilu.no/~burkhart/sharing/pflexible_testdata.tgz (1 GB)

To start off I noticed that the above test data is not available... so I
downloaded

wget http://niflheim.nilu.no/~burkhart/sharing/pflexpart_testdata.tgz (1.7 GB)

I tried the following instructions with the Pure python version and the
FortFlex module:


import pflexible as pf

H = pf.Header('/home_local/jmichalik/example_data')

FD = pf.read_grid(H,time_ret=0,nspec_ret=0)


D=H.fill_backward(nspec=(0,1))

T = pf.read_trajectories(H)

for s,k in H.D:
   data = H.D[(s,k)]
   TC = pf.plot_totalcolumn(H,data,region='POLARCAT',FIGURE=TC)
   TC = pf.plot_trajectory(H,T,k,FIGURE=TC)

With Pure Python I get the following error after the for loop:

Traceback (most recent call last):


  File "<input>", line 1, in <module>
  File "/usr/local/lib/python2.7/dist-packages/pflexible-0.9.0-py2.7-linux-x86_64.egg/pflexible.py",
line 4667, in __ge
tattr__
    return self[attr]


KeyError: 'D'

With the FortFlex module I get the following error after

getting grid for:  ['20070202150000']
Assumed V8 Flexpart
720 180 3 0 0 0 1 88
20070202150000
Traceback (most recent call last):
  File "<input>", line 1, in <module>
  File "/usr/local/lib/python2.7/dist-packages/pflexible-0.9.0-py2.7-linux-x86_64.egg/pflexible.py",
line 1861, in read
_grid
    return readgridV8(H, **kwargs)
  File "/usr/local/lib/python2.7/dist-packages/pflexible-0.9.0-py2.7-linux-x86_64.egg/pflexible.py",
line 2066, in read
gridV8
    OPS.scaledepo, OPS.scaleconc, H.decayconstant, npspec_int)
TypeError: Required argument 'scaleconc' (pos 11) not found

I did try scaleconc=1 as per issue 14 in bitbucket but I still get the same
error.

Any help is appreciated, I need to be able to postprocess binary Flexpart
data...

Thank you in advance.

Best,
Janusz