From: Wendi Madsen <Wendi.B.Madsen_at_nyahnyahspammersnyahnyah>
Date: Fri, 29 Sep 2006 10:07:35 -0600

Good Morning

NCEP recently changed the format of their ensemble data, and this change
required changes in the way we get that data from them. I started with
a Perl script to get the data and then let the old processes take care
of it from there, however there were missing ensembles in the netcdf
files. In an effort to fix this problem I wrote a python script using
Nio to read the NCEP ensemble grib files and write the data to netcdf
files directly. This seems to work except for the times that I put into
the file. I am using the python datetime module to convert a date,
year, month and day, to hours since 1/1/1 at 00:00:00. When I compare
the python generated files to the files generated using a Fortran
program, with udunits calls for time conversion, the python times are 24
hours less than the Fortran times. Are there any times functions,
possibly using udunits functionality, that are part of Pyngl/nio?? Or
is there some other python module that will correct this problem. I
have kludged it for the time being by adding 1 to the number of days
since 1/1/1 before the multiply by 24, and then decode it in the reverse.

Thanks for your assistance.

Wendi Madsen

Wendi B. Madsen
Systems Analyst
Physical Sciences Division
NOAA/OAR,Earth System Research Laboratory
325 Broadway
Boulder, Co 80305
voice 303 497-6869
fax: 303 497-6020
pyngl-talk mailing list
Received on Fri Sep 29 2006 - 10:07:35 MDT

This archive was generated by hypermail 2.2.0 : Mon Oct 23 2006 - 15:54:12 MDT