First, Delia, if you're not a member of this group, would you please join
so I don't have to remember to keep cc'ing you?! :) I'll send you listserv
info later.
Before the bug report - I'm not clear on how development is being performed
for this software. Would you prefer that issues be reported via this list,
or the issues portion of the bitbucket repository (I note that nobody but
me has used this for the pflexible software yet :)). Or, would you even
prefer that we make our own changes to the source and submit them to
bitbucket version control for consideration?
OK, here's the minor bug.....
In trying to adapt some of Delia's "nesting" plots to pflexible data
extraction, I have found a discrepancy in the documentation for the read
header routines. They say to specify that you wanted to use nested output
by inserting the "nest=" flag in function call, but, in fact, it seems like
it should be "nested=". Here is what I wrote up on the bitbucket site:
---------------------------
The documentation for the various read header methods says to use "nest="
to specify nesting or no nesting. However, the correct use should be
"nested="
The following code demonstrates that, if using "nest=1", the nested key
returned with the header dictionary is 0, but if using "nested=1" it is
correct:
#!/usr/bin/env python
FLEXPARTOutputDir =
'/home/morton/Research/CTBTO/FLEXPART_Benchmarking/NestExperiments/OutputNestOnlyWithDryWetConv/output'
# From NILU pflexible project
# Users will need to add the location of pflexible to $PYTHONPATH
import pflexible as pf
# Read the general header information common to all files
H = pf.Header(FLEXPARTOutputDir, nested=1)
isNested = H.nested # 0 = not nested, 1 = nested
print 'isNested: ' + str(isNested)
H = pf.Header(FLEXPARTOutputDir, nest=1)
isNested = H.nested # 0 = not nested, 1 = nested
print 'isNested: ' + str(isNested)
The output of the above segment (which I forgot to post into the bitbucket
issues) is:
Header read:
/home/morton/Research/CTBTO/FLEXPART_Benchmarking/NestExperiments/OutputNestOnlyWithDryWetConv/output/header_nest
isNested: 1
Header read:
/home/morton/Research/CTBTO/FLEXPART_Benchmarking/NestExperiments/OutputNestOnlyWithDryWetConv/output/header
isNested: 0
Thanks,
Don Morton
--
Voice: +1 907 450 8679
Arctic Region Supercomputing Center
http://weather.arsc.edu/
http://people.arsc.edu/~morton/ <http://www.arsc.edu/%7Emorton/>
|