Newer
Older
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!! Settings for S5p/SO2 processing.
!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!-----------------------------------------------------------------------
! user specific settings:
!-----------------------------------------------------------------------
! include user specfic settings:
#include cso-user-settings.rc
!======================================================================
!===
!===
!======================================================================
!-----------------------------------------------------------------------
! inquire SciHub portal
!-----------------------------------------------------------------------
! Obtain names of all available S5p files.
! Stored as csv with processing type, processor version, filenames, etc.
! inquire full time range:
cso.s5p.so2.inquire-table-scihub.timerange.start : ${my.full-timerange.start}
cso.s5p.so2.inquire-table-scihub.timerange.end : ${my.full-timerange.end}
!
! server url ;
! provide login/password in ~/.netrc:
!
! machine s5phub.copernicus.eu login s5pguest password s5pguest
!
cso.s5p.so2.inquire-table-scihub.url : https://s5phub.copernicus.eu/dhus
!!~ empty for no limitation:
!cso.s5p.so2.inquire-table-scihub.area :
!~ domain specified as: west,south:east,north
cso.s5p.so2.inquire-table-scihub.area : ${my.region.west},${my.region.south}:${my.region.east},${my.region.north}
! search query, obtained from interactive download:
!
! platformname : Sentinel-5
! producttype : L2__SO2___ | L2__CO____ (always 10 characters!)
! processinglevel : L2
!
cso.s5p.so2.inquire-table-scihub.query : platformname:Sentinel-5 AND \
producttype:L2__SO2___ AND \
processinglevel:L2
! output table, date of today:
cso.s5p.so2.inquire-table-scihub.output.file : ${my.work}/Copernicus-inquire/Copernicus_S5p_SO2_scihub__%Y-%m-%d.csv
!-----------------------------------------------------------
! overview plot of versions
!-----------------------------------------------------------
! renew existing plots?
cso.s5p.so2.inquire-plot.renew : True
! listing files:
cso.s5p.so2.inquire-plot.file : ${cso.s5p.so2.inquire-table-scihub.output.file}
!!~ specify dates ("yyyy-mm-dd") to use historic tables,
!! default is table of today:
!cso.s5p.so2.inquire-plot.filedate : 2023-08-07
! annote:
cso.s5p.so2.inquire-plot.title : S5p/SO2 %Y-%m-%d
! output figure, date of today:
cso.s5p.so2.inquire-plot.output.file : ${my.work}/Copernicus-inquire/Copernicus_S5p_SO2_scihub__%Y-%m-%d.png
!======================================================================
!===
!===
!======================================================================
! renew existing files (True|False) ?
cso.s5p.so2.convert.renew : True
!cso.s5p.so2.convert.renew : False
! time range:
cso.s5p.so2.convert.timerange.start : ${my.timerange.start}
cso.s5p.so2.convert.timerange.end : ${my.timerange.end}
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! input files
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
! listing of available source files,
! created by 'inquire-scihub' job:
cso.s5p.so2.convert.inquire.file : ${my.work}/Copernicus-inquire/Copernicus_S5p_SO2_scihub__%Y-%m-%d.csv
!!~ historic inquire ...
!cso.s5p.so2.convert.inquire.filedate : ${cso.s5p.so2.inquire-plot.filedate}
! selection keyword:
my.s5p.so2.selection : C03
! Provide ';' seperated list of to decide if a particular orbit file should be processed.
! If more than one file is available for a particular orbit (from "OFFL" and "RPRO" processing),
! the file with the first match will be used.
! The expressions should include templates '%{header}' for the column values.
! Example to select files from collection '03', preferably from processing 'RPRO' but otherwise from 'OFFL':
! (%{collection} == '03') and (%{processing} == 'RPRO') ; \
! (%{collection} == '03') and (%{processing} == 'OFFL')
!
#if "${my.s5p.so2.selection}" == "C03"
cso.s5p.so2.convert.selection : (%{collection} == '03') and (%{processing} == 'RPRO') ; \
(%{collection} == '03') and (%{processing} == 'OFFL')
#else
#error unsupported my.s5p.so2.selection "${my.s5p.so2.selection}"
#endif
! input directory;
! files are searched here or downloaded to if not present yet;
! supported templates:
! %{processing}
cso.s5p.so2.convert.input.dir : ${my.work}/Copernicus/S5P/%{processing}/SO2/%Y/%m
! remove downloaded input files after convert?
cso.s5p.so2.convert.downloads.cleanup : False
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
! selection names:
cso.s5p.so2.convert.filters : lons lats valid quality sza ground_pixel cloud_fraction
! filter settings:
cso.s5p.so2.convert.filter.lons.var : PRODUCT/longitude
cso.s5p.so2.convert.filter.lons.type : minmax
cso.s5p.so2.convert.filter.lons.minmax : ${my.region.west} ${my.region.east}
cso.s5p.so2.convert.filter.lons.units : degrees_east
! filter settings:
cso.s5p.so2.convert.filter.lats.var : PRODUCT/latitude
cso.s5p.so2.convert.filter.lats.type : minmax
cso.s5p.so2.convert.filter.lats.minmax : ${my.region.south} ${my.region.north}
cso.s5p.so2.convert.filter.lats.units : degrees_north
! filter ground pixel sides, more uncertain observations 50-400, here a few more:
!Theys, N., Hedelt, P., De Smedt, I. et al. Global monitoring of volcanic SO2
! degassing with unprecedented resolution from TROPOMI onboard Sentinel-5 Precursor. Sci Rep 9, 2643 (2019).
! https://doi.org/10.1038/s41598-019-39279-y
cso.s5p.so2.convert.filter.ground_pixel.var : PRODUCT/ground_pixel
cso.s5p.so2.convert.filter.ground_pixel.type : minmax
cso.s5p.so2.convert.filter.ground_pixel.minmax : 25 425
cso.s5p.so2.convert.filter.ground_pixel.units : 1
! 2 pixels on the west side of the track are "no-data",
! also "qa_value" is not filled here;
! therefore need to remove these:
cso.s5p.so2.convert.filter.valid.var : PRODUCT/sulfurdioxide_total_vertical_column
cso.s5p.so2.convert.filter.valid.type : valid
! Comment in "PRODUCT/qa_value" variable:
! "A continuous quality descriptor,
! varying between 0 (no data) and 1 (full quality data).
! Recommend to ignore data with qa_value < 0.5"
cso.s5p.so2.convert.filter.quality.var : PRODUCT/qa_value
cso.s5p.so2.convert.filter.quality.type : min
cso.s5p.so2.convert.filter.quality.min : 0.5
cso.s5p.so2.convert.filter.quality.units : 1
! no cloudy pixels:
cso.s5p.so2.convert.filter.cloud_fraction.var : PRODUCT/SUPPORT_DATA/INPUT_DATA/cloud_fraction_crb
cso.s5p.so2.convert.filter.cloud_fraction.type : max
cso.s5p.so2.convert.filter.cloud_fraction.max : 0.2
cso.s5p.so2.convert.filter.cloud_fraction.units : 1
! detection flag for SO2:
! 0 = no detection,
! 1 = detection
! 2 = clear detection close to known volcano
! 3 = clear detection close to known anthropogenic source
! 4 = detection at high SZA
cso.s5p.so2.convert.filter.detection.var : PRODUCT/SUPPORT_DATA/DETAILED_RESULTS/sulfurdioxide_detection_flag
cso.s5p.so2.convert.filter.detection.units : 1
cso.s5p.so2.convert.filter.detection.type : minmax
cso.s5p.so2.convert.filter.detection.minmax : 1 3
! this removes the pixels at very high solar-zenit-angle:
cso.s5p.so2.convert.filter.sza.var : PRODUCT/SUPPORT_DATA/GEOLOCATIONS/solar_zenith_angle
cso.s5p.so2.convert.filter.sza.type : max
cso.s5p.so2.convert.filter.sza.max : 70
cso.s5p.so2.convert.filter.sza.units : degree
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! output files
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! output directory and filename:
! - times are taken from mid of selection, rounded to hours
! - use '%{processing}' for the processing name
! - use '%{orbit}' for orbit number
cso.s5p.so2.convert.output.filename : ${my.work}/CSO-data/${my.region}/S5p/SO2/${my.s5p.so2.selection}/%Y/%m/S5p_SO2_%{orbit}.nc
! pack variables on output:
cso.s5p.so2.convert.output.packed : True
! zlib compression level, 0 for no compression:
cso.s5p.so2.convert.output.complevel : 1
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
! global attributes:
cso.s5p.so2.convert.output.attrs : format Conventions \
author institution email
!
cso.s5p.so2.convert.output.attr.format : ${my.cso.format}
cso.s5p.so2.convert.output.attr.Conventions : ${my.cso.convention}
cso.s5p.so2.convert.output.attr.author : ${my.attr.author}
cso.s5p.so2.convert.output.attr.institution : ${my.attr.institution}
cso.s5p.so2.convert.output.attr.email : ${my.attr.email}
! no need to swap layes:
cso.s5p.so2.convert.swap_layers : False
! ~ variables
! which fields to be put out ?
cso.s5p.so2.convert.output.vars : longitude longitude_bounds \
latitude latitude_bounds \
track_longitude track_longitude_bounds \
track_latitude track_latitude_bounds \
time \
pressure kernel qa_value \
vcd vcd_errvar \
detection_flag cloud_fraction solar_zenith_angle ground_pixel
!
! Describe per variable:
! * .dims : dimensions list:
! pixel : selected pixels
! corner : number of footprint bounds (probably 4)
! layer : number of layers in atmospheric profile (layers in kernel)
! layeri : number of layer interfaces in atmospheric profile (layer+1)
! retr : number of layers in retrieval product (1 for columns) ;
! for error covariance use (retr,retr0) to avoid repeated dimensions
! track_scan : original scan index in 2D track
! track_pixel : original ground pixel in 2D track
! * .specal : keyword to select special processing
! * None : no special processing (default)
! * track_longitude : longitudes at centers of original 2D track
! * track_latitude : latitudes at centers of original 2D track
! * track_longitude_bounds : longitude bounds at centers of original 2D track
! * track_latitude_bounds : latitude bounds at centers of original 2D track
! * .units : target units if different from original
! * .oper : special postprocessing, currently supported:
! * square : fill variable with squared valued (used to form variance from standard deviation)
! In case no special processing is needed:
! * .from : original variable (group path and variable name)
!
! center longitudes; remove bounds attribute, no coordinate ...
cso.s5p.so2.convert.output.var.longitude.dims : pixel
cso.s5p.so2.convert.output.var.longitude.from : PRODUCT/longitude
cso.s5p.so2.convert.output.var.longitude.attrs : { 'bounds' : None }
! center latitudes; remove bounds attribute, no coordinate ...
cso.s5p.so2.convert.output.var.latitude.dims : pixel
cso.s5p.so2.convert.output.var.latitude.from : PRODUCT/latitude
cso.s5p.so2.convert.output.var.latitude.attrs : { 'bounds' : None }
! corner longitudes; no units in file:
cso.s5p.so2.convert.output.var.longitude_bounds.dims : pixel corner
cso.s5p.so2.convert.output.var.longitude_bounds.from : PRODUCT/SUPPORT_DATA/GEOLOCATIONS/longitude_bounds
cso.s5p.so2.convert.output.var.longitude_bounds.units : degrees_east
! corner latitudes, no units in file:
cso.s5p.so2.convert.output.var.latitude_bounds.dims : pixel corner
cso.s5p.so2.convert.output.var.latitude_bounds.from : PRODUCT/SUPPORT_DATA/GEOLOCATIONS/latitude_bounds
cso.s5p.so2.convert.output.var.latitude_bounds.units : degrees_north
! original track:
!~ center lon; remove bounds attribute, no coordinate ...
cso.s5p.so2.convert.output.var.track_longitude.dims : track_scan track_pixel
cso.s5p.so2.convert.output.var.track_longitude.special : track_longitude
cso.s5p.so2.convert.output.var.track_longitude.from : PRODUCT/longitude
cso.s5p.so2.convert.output.var.track_longitude.attrs : { 'bounds' : None }
!~ center lat; remove bounds attribute, no coordinate ...
cso.s5p.so2.convert.output.var.track_latitude.dims : track_scan track_pixel
cso.s5p.so2.convert.output.var.track_latitude.special : track_latitude
cso.s5p.so2.convert.output.var.track_latitude.from : PRODUCT/latitude
cso.s5p.so2.convert.output.var.track_latitude.attrs : { 'bounds' : None }
!~ corner lons
cso.s5p.so2.convert.output.var.track_longitude_bounds.dims : track_scan track_pixel corner
cso.s5p.so2.convert.output.var.track_longitude_bounds.special : track_longitude_bounds
cso.s5p.so2.convert.output.var.track_longitude_bounds.from : PRODUCT/SUPPORT_DATA/GEOLOCATIONS/longitude_bounds
cso.s5p.so2.convert.output.var.track_longitude_bounds.units : degrees_east
!~ corner lats
cso.s5p.so2.convert.output.var.track_latitude_bounds.dims : track_scan track_pixel corner
cso.s5p.so2.convert.output.var.track_latitude_bounds.special : track_latitude_bounds
cso.s5p.so2.convert.output.var.track_latitude_bounds.from : PRODUCT/SUPPORT_DATA/GEOLOCATIONS/latitude_bounds
cso.s5p.so2.convert.output.var.track_latitude_bounds.units : degrees_north
!~ solar zenith angle
cso.s5p.so2.convert.output.var.solar_zenith_angle.dims : pixel
cso.s5p.so2.convert.output.var.solar_zenith_angle.from : PRODUCT/SUPPORT_DATA/GEOLOCATIONS/solar_zenith_angle
cso.s5p.so2.convert.output.var.solar_zenith_angle.units : degree
! time value per scan line
cso.s5p.so2.convert.output.var.time.dims : pixel
cso.s5p.so2.convert.output.var.time.special : time-delta
cso.s5p.so2.convert.output.var.time.tref : PRODUCT/time
cso.s5p.so2.convert.output.var.time.dt : PRODUCT/delta_time
! ground position across track
cso.s5p.so2.convert.output.var.ground_pixel.dims : track_pixel
cso.s5p.so2.convert.output.var.ground_pixel.special : ground_pixel
cso.s5p.so2.convert.output.var.ground_pixel.from : PRODUCT/ground_pixel
! vertical column density:
cso.s5p.so2.convert.output.var.vcd.dims : pixel retr
cso.s5p.so2.convert.output.var.vcd.from : PRODUCT/sulfurdioxide_total_vertical_column
! error variance in vertical column density (after application of kernel),
! fill with single element 'covariance matrix', from square of standard error:
! use dims with different names to avoid that cf-checker complains:
cso.s5p.so2.convert.output.var.vcd_errvar.dims : pixel retr retr0
cso.s5p.so2.convert.output.var.vcd_errvar.special : square
cso.s5p.so2.convert.output.var.vcd_errvar.from : PRODUCT/sulfurdioxide_total_vertical_column_precision
!~ skip standard name, modifier "standard_error" is not valid anymore:
cso.s5p.so2.convert.output.var.vcd_errvar.attrs : { 'standard_name' : None }
! Convert from hybride coefficient bounds in (2,nlev) aray to 3D half level pressure:
cso.s5p.so2.convert.output.var.pressure.dims : pixel layeri
cso.s5p.so2.convert.output.var.pressure.special : hymid_to_pressure
cso.s5p.so2.convert.output.var.pressure.sp : PRODUCT/SUPPORT_DATA/INPUT_DATA/surface_pressure
cso.s5p.so2.convert.output.var.pressure.hyam : PRODUCT/SUPPORT_DATA/INPUT_DATA/tm5_constant_a
cso.s5p.so2.convert.output.var.pressure.hybm : PRODUCT/SUPPORT_DATA/INPUT_DATA/tm5_constant_b
cso.s5p.so2.convert.output.var.pressure.units : Pa
cso.s5p.so2.convert.output.var.kernel.dims : pixel layer retr
cso.s5p.so2.convert.output.var.kernel.from : PRODUCT/SUPPORT_DATA/DETAILED_RESULTS/averaging_kernel
! cloud property:
cso.s5p.so2.convert.output.var.cloud_fraction.dims : pixel
cso.s5p.so2.convert.output.var.cloud_fraction.from : PRODUCT/SUPPORT_DATA/INPUT_DATA/cloud_fraction_crb
cso.s5p.so2.convert.output.var.cloud_fraction.attrs : { 'coordinates' : None, 'source' : None }
!! cloud property:
!cso.s5p.so2.convert.output.var.cloud_radiance_fraction.dims : pixel
!cso.s5p.so2.convert.output.var.cloud_radiance_fraction.from : PRODUCT/SUPPORT_DATA/DETAILED_RESULTS/cloud_radiance_fraction_nitrogendioxide_window
!cso.s5p.so2.convert.output.var.cloud_radiance_fraction.attrs : { 'coordinates' : None, 'ancillary_variables' : None }
! quality flag:
cso.s5p.so2.convert.output.var.qa_value.dims : pixel
cso.s5p.so2.convert.output.var.qa_value.from : PRODUCT/qa_value
!~ skip some attributes, cf-checker complains ...
cso.s5p.so2.convert.output.var.qa_value.attrs : { 'valid_min' : None, 'valid_max' : None }
! detection flag:
cso.s5p.so2.convert.output.var.detection_flag.dims : pixel
cso.s5p.so2.convert.output.var.detection_flag.from : PRODUCT/SUPPORT_DATA/DETAILED_RESULTS/sulfurdioxide_detection_flag
cso.s5p.so2.convert.output.var.detection_flag.attrs : { 'coordinates' : None }
cso.s5p.so2.convert.output.var.detection_flag.dtype : i1
!======================================================================
!===
!=== listing
!===
!======================================================================
! csv file that will hold records per file with:
! - timerange of pixels in file
! - orbit number
cso.s5p.so2.listing.file : ${my.work}/CSO-data/${my.region}/S5p/SO2/${my.s5p.so2.selection}__listing.csv
! renew table if file already exists?
cso.s5p.so2.listing.renew : True
! time range:
cso.s5p.so2.listing.timerange.start : ${my.timerange.start}
cso.s5p.so2.listing.timerange.end : ${my.timerange.end}
! filename filters relative to listing file that should be scanned for orbit files;
! names could include time templates ;
! if same orbit is found in multiple directories, the first found is used;
! remove existing table for safety to ensure that this is done correctly ...
cso.s5p.so2.listing.patterns : ${my.s5p.so2.selection}/%Y/%m/S5p_*.nc
! extra columns to be added, read from global attributes:
cso.s5p.so2.listing.xcolumns : orbit
!======================================================================
!===
!=== catalogue
!===
!======================================================================
! listing file with filenames/timerange.
cso.s5p.so2.catalogue.input.listing : ${cso.s5p.so2.listing.file}
! time range:
cso.s5p.so2.catalogue.timerange.start : ${my.timerange.start}
cso.s5p.so2.catalogue.timerange.end : ${my.timerange.end}
! target filenames; templates:
! - time values
! - %{orbit} : from listing
! - %{varname} for variable
cso.s5p.so2.catalogue.output.file : ${my.work}/CSO-data-catalogue/${my.region}/S5p/SO2/${my.s5p.so2.selection}/%Y/%m/%d/S5p_SO2_%{orbit}__%{varname}.png
! map domain (west east south north)
cso.s5p.so2.catalogue.domain : ${my.region.west} ${my.region.east} ${my.region.south} ${my.region.north}
! figure size (inches), default is (8,6):
cso.s5p.so2.catalogue.figsize : ${my.region.figsize}
! renew existing files?
!cso.s5p.so2.catalogue.renew : False
cso.s5p.so2.catalogue.renew : True
! variables to be plotted:
cso.s5p.so2.catalogue.vars : vcd vcd_errvar qa_value \
cloud_fraction detection_flag solar_zenith_angle
!! color for no-data values in track, default '0.8' (gray):
!cso.s5p.so2.catalogue.color_nan : white
!! extra keyword arguments for map:
!cso.s5p.so2.catalogue.map : resolution='h'
! convert units:
cso.s5p.so2.catalogue.var.vcd.units : umol/m2
! style:
cso.s5p.so2.catalogue.var.vcd.vmin : 0.0
cso.s5p.so2.catalogue.var.vcd.vmax : 100.0
! show error as std.dev, convert to vcd units:
cso.s5p.so2.catalogue.var.vcd_errvar.units : umol/m2
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
! style:
cso.s5p.so2.catalogue.var.vcd_errvar.vmax : 100.0
! style:
cso.s5p.so2.catalogue.var.qa_value.vmin : 0.5
cso.s5p.so2.catalogue.var.qa_value.vmax : 1.0
cso.s5p.so2.catalogue.var.qa_value.colors : ['red','yellow','green']
! style:
cso.s5p.so2.catalogue.var.solar_zenith_angle.vmin : 0
cso.s5p.so2.catalogue.var.solar_zenith_angle.vmax : 90
! style:
cso.s5p.so2.catalogue.var.detection_flag.vmin : 0
cso.s5p.so2.catalogue.var.detection_flag.vmax : 4
cso.s5p.so2.catalogue.var.detection_flag.colors : ['purple','red','yellow','green','blue']
! style:
cso.s5p.so2.catalogue.var.cloud_fraction.vmax : 1.0
cso.s5p.so2.catalogue.var.cloud_fraction.colors : ['blue','cyan','white']
!! style:
!cso.s5p.so2.catalogue.var.cloud_radiance_fraction.vmax : 1.0
!cso.s5p.so2.catalogue.var.cloud_radiance_fraction.colors : ['blue','cyan','white']
!-----------------------------------------------------------
! catalogue index
!-----------------------------------------------------------
! target location:
cso.s5p.so2.catalogue-index.outdir : ${my.work}/CSO-data-catalogue/${my.region}/S5p/SO2/${my.s5p.so2.selection}
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
! title:
cso.s5p.so2.catalogue-index.header : CSO catalogue
! show info on created page?
cso.s5p.so2.catalogue-index.info : True
! create new page for each value?
cso.s5p.so2.catalogue-index.newpage : True
! content type:
cso.s5p.so2.catalogue-index.type : list
! define row values:
cso.s5p.so2.catalogue-index.name : date
cso.s5p.so2.catalogue-index.values : TimeSeries( ${my.timerange.start}, ${my.timerange.end}, '1 day', '%Y%m%d' )
! create new page for each value:
cso.s5p.so2.catalogue-index.date.newpage : True
! content type:
cso.s5p.so2.catalogue-index.date.type : table-row
! define row values:
cso.s5p.so2.catalogue-index.date.name : orbit
cso.s5p.so2.catalogue-index.date.values : CsvFile( '%{date[0:4]}/%{date[4:6]}/%{date[6:8]}/orbits.csv' )
! content type:
cso.s5p.so2.catalogue-index.date.orbit.type : table-col
! define column values:
cso.s5p.so2.catalogue-index.date.orbit.name : var
cso.s5p.so2.catalogue-index.date.orbit.values : ${cso.s5p.so2.catalogue.vars}
! content type:
cso.s5p.so2.catalogue-index.date.orbit.var.type : img
! define image:
cso.s5p.so2.catalogue-index.date.orbit.var.img : %{date[0:4]}/%{date[4:6]}/%{date[6:8]}/S5p_SO2_%{orbit}__%{var}.png
cso.s5p.so2.catalogue-index.date.orbit.var.kwargs : height=300
!======================================================================
!===
!===
!======================================================================
!-----------------------------------------------------------
!-----------------------------------------------------------
! time range:
cso.s5p.so2.gridded.timerange.start : ${my.timerange.start}
cso.s5p.so2.gridded.timerange.end : ${my.timerange.end}
! create one gridded file per hour:
cso.s5p.so2.gridded.timerange.step : hour
! renew existing files?
cso.s5p.so2.gridded.renew : True
! target directory, incl. subdir for resolution and filters:
my.gridded.dir : ${my.work}/CSO-gridded/${my.region}__r01x01/S5p/SO2/${my.s5p.so2.selection}__${my.s5p.so2.gridded-selection}
!~
! grid definition:
!~ same as pixel selection on conversion:
cso.s5p.so2.gridded.grid.west : ${my.region.west}
cso.s5p.so2.gridded.grid.east : ${my.region.east}
cso.s5p.so2.gridded.grid.south : ${my.region.south}
cso.s5p.so2.gridded.grid.north : ${my.region.north}
! resolution:
cso.s5p.so2.gridded.grid.dlon : 0.1
cso.s5p.so2.gridded.grid.dlat : 0.1
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
! level of recursive splitting of footprint into triangles,
! and assignment of centroids to grid cells;
! for 4-corner footprints, number of centroids is:
! 1 (levels=0), 4 (1), 8 (2), 16 (3), 64 (5), 256 (7)
cso.s5p.so2.gridded.mapping.levels : 5
!~
! keywords for source files;
! the first one should have the footprints;
! here the only source are the converted files:
cso.s5p.so2.gridded.sources : data
! input files for each source type:
!~ here: specify listing file,
! this is only supported for a single source
cso.s5p.so2.gridded.source.data.listing : ${cso.s5p.so2.listing.file}
!!~ alternative: filename patterns with time templates,
!! but here the source files do not have time stamps
!cso.s5p.so2.gridded.source.data.filenames : ${my.work}/CSO-data/S5p/RPRO/SO2/${my.region}/%Y/%m/S5p_*.nc
!~
! filter description:
my.s5p.so2.gridded-selection : qa08
! switch:
#if "${my.s5p.so2.gridded-selection}" == "qa08"
! keywords for filters:
cso.s5p.so2.gridded.filters : quality
! minimum quality value required:
cso.s5p.so2.gridded.filter.quality.var : qa_value
cso.s5p.so2.gridded.filter.quality.type : min
cso.s5p.so2.gridded.filter.quality.min : 0.8
cso.s5p.so2.gridded.filter.quality.units : 1
#else
#error unsupported my.s5p.so2.gridded-selection "${my.s5p.so2.gridded-selection}"
#endif
!~
! target file, might contain templates:
! %Y,%m,etc : time values
cso.s5p.so2.gridded.output.file : ${my.gridded.dir}/%Y/%m/S5p_SO2_%Y%m%d_%H%M_gridded.nc
! pack floats as short values?
cso.s5p.so2.gridded.output.packed : True
! zlib compression level (default 1, 0 for no compression):
cso.s5p.so2.gridded.output.complevel : 1
! data variables to be created:
cso.s5p.so2.gridded.output.vars : yr
! input variables:
! data:vcd : variable "vcd" from source "data"
cso.s5p.so2.gridded.output.yr.source : data:vcd
!-----------------------------------------------------------
! catalogue of gridded fields
!-----------------------------------------------------------
! time range:
cso.s5p.so2.gridded-catalogue.timerange.start : ${my.timerange.start}
cso.s5p.so2.gridded-catalogue.timerange.end : ${my.timerange.end}
! create one gridded file per hour:
cso.s5p.so2.gridded-catalogue.timerange.step : hour
! renew existing files?
!cso.s5p.so2.gridded-catalogue.renew : True
cso.s5p.so2.gridded-catalogue.renew : False
! target directory for catalogue:
cso.s5p.so2.gridded-catalogue.output.dir : ${my.gridded.dir}/catalogue
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
cso.s5p.so2.gridded-catalogue.input.file : ${my.gridded.dir}/%Y/%m/S5p_SO2_%Y%m%d_%H%M_gridded.nc
!!~ idem for daily average:
!cso.s5p.so2.gridded-catalogue.input.file : ${my.gridded.dir}/S5p_SO2_%Y%m%d_aver_gridded.nc
! target files, time tempates and variable name are replaced:
cso.s5p.so2.gridded-catalogue.output.file : ${my.gridded.dir}/catalogue/%Y/%m/%d/S5p_SO2_%Y%m%d_%H%M_gridded_%{var}.png
!!~ idem for daily average:
!cso.s5p.so2.gridded-catalogue.output.file : ${my.gridded.dir}/catalogue/%Y/%m/%d/S5p_SO2_%Y%m%d_aver_gridded_%{var}.png
! figure size (inches), default is (8,6):
cso.s5p.so2.gridded-catalogue.figsize : (6,6)
! variables to be plotted:
cso.s5p.so2.gridded-catalogue.vars : yr
! variable:
cso.s5p.so2.gridded-catalogue.var.yr.source : yr
! convert units:
cso.s5p.so2.gridded-catalogue.var.yr.units : umol/m2
! style:
cso.s5p.so2.gridded-catalogue.var.yr.vmin : 0.0
cso.s5p.so2.gridded-catalogue.var.yr.vmax : 100.0
!-----------------------------------------------------------
! gridded catalogue index
!-----------------------------------------------------------
! target location:
cso.s5p.so2.gridded-catalogue-index.outdir : ${cso.s5p.so2.gridded-catalogue.output.dir}
! title:
cso.s5p.so2.gridded-catalogue-index.header : CSO catalogue
! show info on created page?
cso.s5p.so2.gridded-catalogue-index.info : True
! create new page for each value?
cso.s5p.so2.gridded-catalogue-index.newpage : True
! content type:
cso.s5p.so2.gridded-catalogue-index.type : list
! define row values:
cso.s5p.so2.gridded-catalogue-index.name : date
cso.s5p.so2.gridded-catalogue-index.values : TimeSeries( ${my.timerange.start}, ${my.timerange.end}, '1 day', '%Y%m%d' )
! create new page for each value:
cso.s5p.so2.gridded-catalogue-index.date.newpage : True
! content type:
cso.s5p.so2.gridded-catalogue-index.date.type : table-row
! define row values:
cso.s5p.so2.gridded-catalogue-index.date.name : time
cso.s5p.so2.gridded-catalogue-index.date.values : Range( 0, 23, 1, '%2.2i00' )
! content type:
cso.s5p.so2.gridded-catalogue-index.date.time.type : table-col
! define column values:
cso.s5p.so2.gridded-catalogue-index.date.time.name : var
cso.s5p.so2.gridded-catalogue-index.date.time.values : ${cso.s5p.so2.gridded-catalogue.vars}
! content type:
cso.s5p.so2.gridded-catalogue-index.date.time.var.type : img
! define image:
cso.s5p.so2.gridded-catalogue-index.date.time.var.img : %{date[0:4]}/%{date[4:6]}/%{date[6:8]}/S5p_SO2_%{date}_%{time}_gridded_%{var}.png
cso.s5p.so2.gridded-catalogue-index.date.time.var.kwargs : height=300
!======================================================================
!===
!=== simulation catalogue
!===
!======================================================================
!-----------------------------------------------------------
! simulation catalogue
!-----------------------------------------------------------
! time range:
cso.s5p.so2.sim-catalogue.timerange.start : ${my.timerange.start}
cso.s5p.so2.sim-catalogue.timerange.end : ${my.timerange.end}
cso.s5p.so2.sim-catalogue.timerange.step : hour
! input files:
cso.s5p.so2.sim-catalogue.input.data.file : ${my.work}/CSO-oper/CSO_output_%Y%m%d_%H%M_data.nc
cso.s5p.so2.sim-catalogue.input.state.file : ${my.work}/CSO-oper/CSO_output_%Y%m%d_%H%M_state.nc
! target files, time tempates and variable name are replaced:
cso.s5p.so2.sim-catalogue.output.file : ${my.work}/CSO-data-sim-catalogue/S5p/SO2/${my.region}/%Y/%m/%d/S5p_SO2_%Y%m%d_%H%M_%{var}.png
! map domain used for simulations (west east south north):
cso.s5p.so2.sim-catalogue.domain : ${my.region.west} ${my.region.east} ${my.region.south} ${my.region.north}
! figure size (inches), default is (8,6):
cso.s5p.so2.sim-catalogue.figsize : ${my.region.figsize}
cso.s5p.so2.sim-catalogue.renew : False
!cso.s5p.so2.sim-catalogue.renew : True
! variables to be plotted:
cso.s5p.so2.sim-catalogue.vars : yr ys
! variable:
cso.s5p.so2.sim-catalogue.var.yr.source : data:yr
! convert units:
cso.s5p.so2.sim-catalogue.var.yr.units : umol/m2
! style:
cso.s5p.so2.sim-catalogue.var.yr.vmin : 0.0
cso.s5p.so2.sim-catalogue.var.yr.vmax : 100.0
! variable:
cso.s5p.so2.sim-catalogue.var.ys.source : state:ys
! convert units:
cso.s5p.so2.sim-catalogue.var.ys.units : umol/m2
! style:
cso.s5p.so2.sim-catalogue.var.ys.vmin : 0.0
cso.s5p.so2.sim-catalogue.var.ys.vmax : 100.0
!-----------------------------------------------------------
! sim catalogue index
!-----------------------------------------------------------
! target location:
cso.s5p.so2.sim-catalogue-index.outdir : ${my.work}/CSO-sim-catalogue
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
! title:
cso.s5p.so2.sim-catalogue-index.header : CSO catalogue
! show info on created page?
cso.s5p.so2.sim-catalogue-index.info : True
! create new page for each value?
cso.s5p.so2.sim-catalogue-index.newpage : True
! content type:
cso.s5p.so2.sim-catalogue-index.type : list
! define row values:
cso.s5p.so2.sim-catalogue-index.name : date
cso.s5p.so2.sim-catalogue-index.values : TimeSeries( ${my.timerange.start}, ${my.timerange.end}, '1 day', '%Y%m%d' )
! create new page for each value:
cso.s5p.so2.sim-catalogue-index.date.newpage : True
! content type:
cso.s5p.so2.sim-catalogue-index.date.type : table-row
! define row values:
cso.s5p.so2.sim-catalogue-index.date.name : time
cso.s5p.so2.sim-catalogue-index.date.values : Range( 0, 23, 1, '%2.2i00' )
! content type:
cso.s5p.so2.sim-catalogue-index.date.time.type : table-col
! define column values:
cso.s5p.so2.sim-catalogue-index.date.time.name : var
cso.s5p.so2.sim-catalogue-index.date.time.values : ${cso.s5p.so2.sim-catalogue.vars}
! content type:
cso.s5p.so2.sim-catalogue-index.date.time.var.type : img
! define image:
cso.s5p.so2.sim-catalogue-index.date.time.var.img : %{date[0:4]}/%{date[4:6]}/%{date[6:8]}/S5p_SO2_%{date}_%{time}_%{var}.png
cso.s5p.so2.sim-catalogue-index.date.time.var.kwargs : height=300
!======================================================================
!===
!===
!======================================================================
!-----------------------------------------------------------
! gridded orbits
!-----------------------------------------------------------
! time range:
cso.s5p.so2.sim-gridded.timerange.start : ${my.timerange.start}
cso.s5p.so2.sim-gridded.timerange.end : ${my.timerange.end}
! keyword used for subdir for resolution and filters:
my.gridded.subdir : ${my.region}__r01x01__qa08
!~
! testing ..
my.gridded.region : ${my.region}
! grid definition:
!~ same as pixel selection on conversion:
cso.s5p.so2.sim-gridded.grid.west : ${my.region.west}
cso.s5p.so2.sim-gridded.grid.east : ${my.region.east}
cso.s5p.so2.sim-gridded.grid.south : ${my.region.south}
cso.s5p.so2.sim-gridded.grid.north : ${my.region.north}
cso.s5p.so2.sim-gridded.grid.dlon : 0.1
cso.s5p.so2.sim-gridded.grid.dlat : 0.1
!~
! keywords for source files;
! the first one should have the footprints:
cso.s5p.so2.sim-gridded.sources : data state
! input files for each source type:
cso.s5p.so2.sim-gridded.source.data : ${my.work}/CSO-oper/CSO_output_%Y%m%d_%H%M_data.nc
cso.s5p.so2.sim-gridded.source.state : ${my.work}/CSO-oper/CSO_output_%Y%m%d_%H%M_state.nc
!~
! keywords for filters:
cso.s5p.so2.sim-gridded.filters : quality detection
! minimum quality value required:
cso.s5p.so2.sim-gridded.filter.quality.var : qa_value
cso.s5p.so2.sim-gridded.filter.quality.type : min
cso.s5p.so2.sim-gridded.filter.quality.min : 0.75
cso.s5p.so2.sim-gridded.filter.quality.units : 1
! detection flag:
cso.s5p.so2.gridded.filter.detection.var : detection_flag
cso.s5p.so2.gridded.filter.detection.type : minmax
cso.s5p.so2.gridded.filter.detection.min : 1 3
cso.s5p.so2.gridded.filter.detection.units : 1
!~
! target file, might contain templates:
! %Y,%m,etc : time values
! %{basename} : basename (without extension) of first source file
cso.s5p.so2.sim-gridded.output.file : ${my.work}/CSO-gridded/${my.gridded.subdir}/%Y/%m/CSO_output_%Y%m%d_%H%M_gridded.nc
! data variables to be created:
! input variables:
! data:yr : from data file
! state:ys : from state file
cso.s5p.so2.sim-gridded.output.yr.source : data:yr
cso.s5p.so2.sim-gridded.output.ys.source : state:ys
!-----------------------------------------------------------
! catalogue of gridded fields
!-----------------------------------------------------------
! time range:
cso.s5p.so2.sim-gridded-catalogue.timerange.start : ${my.timerange.start}
cso.s5p.so2.sim-gridded-catalogue.timerange.end : ${my.timerange.end}
cso.s5p.so2.sim-gridded-catalogue.renew : True
!cso.s5p.so2.sim-gridded-catalogue.renew : False
! target directory for catalogue:
cso.s5p.so2.sim-gridded-catalogue.output.dir : ${my.work}/CSO-gridded/${my.gridded.subdir}/catalogue
cso.s5p.so2.sim-gridded-catalogue.input.file : ../%Y/%m/CSO_output_%Y%m%d_%H%M_gridded.nc
!cso.s5p.so2.sim-gridded-catalogue.input.file : ../CSO_output_%Y%m%d_aver_gridded.nc
! target files, time tempates and variable name are replaced:
cso.s5p.so2.sim-gridded-catalogue.output.file : %Y/%m/%d/S5p_RPRO_SO2_%Y%m%d_%H%M_gridded_%{var}.png
!cso.s5p.so2.sim-gridded-catalogue.output.file : %Y/%m/%d/S5p_RPRO_SO2_%Y%m%d_aver_gridded_%{var}.png
! figure size (inches), default is (8,6):
cso.s5p.so2.sim-gridded-catalogue.figsize : (6,6)
cso.s5p.so2.sim-gridded-catalogue.vars : yr ys
cso.s5p.so2.sim-gridded-catalogue.var.yr.source : yr
cso.s5p.so2.sim-gridded-catalogue.var.yr.units : umol/m2
cso.s5p.so2.sim-gridded-catalogue.var.yr.vmin : 0.0
cso.s5p.so2.sim-gridded-catalogue.var.yr.vmax : 100.0
cso.s5p.so2.sim-gridded-catalogue.var.ys.source : ys
cso.s5p.so2.sim-gridded-catalogue.var.ys.units : umol/m2
cso.s5p.so2.sim-gridded-catalogue.var.ys.vmin : 0.0
cso.s5p.so2.sim-gridded-catalogue.var.ys.vmax : 100.0
!-----------------------------------------------------------
! gridded catalogue index
!-----------------------------------------------------------
! target location:
cso.s5p.so2.sim-gridded-catalogue-index.outdir : ${cso.s5p.so2.sim-gridded-catalogue.output.dir}
cso.s5p.so2.sim-gridded-catalogue-index.header : CSO catalogue
cso.s5p.so2.sim-gridded-catalogue-index.info : True
! create new page for each value?
cso.s5p.so2.sim-gridded-catalogue-index.newpage : True
cso.s5p.so2.sim-gridded-catalogue-index.type : list
cso.s5p.so2.sim-gridded-catalogue-index.name : date
cso.s5p.so2.sim-gridded-catalogue-index.values : TimeSeries( ${my.timerange.start}, ${my.timerange.end}, '1 day', '%Y%m%d' )
! create new page for each value:
cso.s5p.so2.sim-gridded-catalogue-index.date.newpage : True
cso.s5p.so2.sim-gridded-catalogue-index.date.type : table-row
cso.s5p.so2.sim-gridded-catalogue-index.date.name : time
cso.s5p.so2.sim-gridded-catalogue-index.date.values : Range( 0, 23, 1, '%2.2i00' )
cso.s5p.so2.sim-gridded-catalogue-index.date.time.type : table-col
cso.s5p.so2.sim-gridded-catalogue-index.date.time.name : var
cso.s5p.so2.sim-gridded-catalogue-index.date.time.values : ${cso.s5p.so2.sim-gridded-catalogue.vars}
cso.s5p.so2.sim-gridded-catalogue-index.date.time.var.type : img
cso.s5p.so2.sim-gridded-catalogue-index.date.time.var.img : %{date[0:4]}/%{date[4:6]}/%{date[6:8]}/S5p_RPRO_SO2_%{date}_%{time}_gridded_%{var}.png
cso.s5p.so2.sim-gridded-catalogue-index.date.time.var.kwargs : height=300
!======================================================================
!===
!=== end
!===
!======================================================================