Skip to content
Snippets Groups Projects
Commit 39ee37e5 authored by Peters, Wouter's avatar Peters, Wouter
Browse files

other filename to be read

parent 2d5f153c
No related branches found
No related tags found
No related merge requests found
...@@ -29,12 +29,12 @@ def write_mixing_ratios(dacycle): ...@@ -29,12 +29,12 @@ def write_mixing_ratios(dacycle):
Write Sample information to NetCDF files. These files are organized by site and Write Sample information to NetCDF files. These files are organized by site and
have an unlimited time axis to which data is appended each cycle. have an unlimited time axis to which data is appended each cycle.
The needed information is obtained from the sampleinfo.nc files and the original input data files from ObsPack. The needed information is obtained from the sample_auxiliary.nc files and the original input data files from ObsPack.
The steps are: The steps are:
(1) Create a directory to hold timeseries output files (1) Create a directory to hold timeseries output files
(2) Read the sampleinfo.nc file for this cycle and get a list of original files they were obtained from (2) Read the sample_auxiliary.nc file for this cycle and get a list of original files they were obtained from
(3) For each file, copy the original data file from ObsPack (if not yet present) (3) For each file, copy the original data file from ObsPack (if not yet present)
(4) Open the copied file, find the index of each observation, fill in the simulated data (4) Open the copied file, find the index of each observation, fill in the simulated data
...@@ -56,7 +56,7 @@ def write_mixing_ratios(dacycle): ...@@ -56,7 +56,7 @@ def write_mixing_ratios(dacycle):
# Step (1): Get the posterior sample output data file for this cycle # Step (1): Get the posterior sample output data file for this cycle
infile = os.path.join(dacycle['dir.output'], 'sampleinfo_%s.nc' % dacycle['time.sample.stamp']) infile = os.path.join(dacycle['dir.output'], 'sample_auxiliary_%s.nc' % dacycle['time.sample.stamp'])
ncf_in = io.ct_read(infile, 'read') ncf_in = io.ct_read(infile, 'read')
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment