User Tools

Site Tools


codes:calibration

Calibration Functions

This page provides a brief overview of the functions avaliable in Eubrewnet to determine Brewer configuration parameters.


Langley

  • Function: langley
  • Description: Determine the ETC and the ETC filter corrections applying the Langley-plot method to each half day.
  • Input: a complete list of input arguments is provided next, but we strongly recommed to use the web interface at https://rbcce-test.aemet.es/eubrewnet/caldata/langley to access this function
    1. format: see Format
    2. data: input data used, available options are o3l15 (default if the option is not passed in the url), o3l15v2, o3l1, sco3 for the O3L1.5, O3L1.5V2, O3L1 and sun scan measurements, respectively. For the sco3 option, the step option with the value of the micrometer step must be also passed in the url. Individual observation (not summaries) are used in all cases. Note you can override the default configuration used to process the data in Eubrewnet using the urlconfig method (eg, to pass new ETC filter correction methods add &urlconfig&fc=1,2,3,4,5 to the url)
    3. calc: calculation to perform, available options are etc (Langley-plot with all the data, default option), filters (the Langley-plot fit is done using dummy variables for filters 3, 4, 5, and thus allows to determine the ETC filter correction of those filters with respect to the ETC obtained for the set measurements with filters 0, 1, and 2), and filters_etc (this option performs two fits: first, a filters run to determine the ETC filter corrections, and then an etc calculation using these ETC filter corrections). For etc runs, two different fits are performed: ms9 vs airmass (brewer method) and ms9/airmass vs 1/airmass (dobson method)
    4. minairmass: airmass minimum in etc calculations (defaults to 1.70 if the option is not passed in the url)
    5. maxairmass: airmass maximum in etc calculations (defaults to 3.75 if the option is not passed in the url)
    6. minobs: minimum number of observations in the half-day to perform the Langley in etc calculations (default value: 25)
    7. filters_minairmass: airmass minimum in filters calculations (defaults to 1 if the option is not passed in the url)
    8. filters_maxairmass: airmass maximum in filters calculations (defaults to 5 if the option is not passed in the url)
    9. filters_minobs: minimum number of observations in the half-day to perform the Langley in filters calculations (default value: 15)
    10. maxstdo3: don't perform the langley if the maximum of the standard deviation of the half-day ozone is above this value (default: 2.5)
    11. minrsquared: don't accept the langley results if the r2 of the fits are below this value (default: 0.999)
  • Output:
    1. date: date and time of each half-day langley, morning langleys correspond to 00 time, and afternoon ones, to 12
    2. etc: calculated ETC
    3. o3xabs: calculated ozone multiplied by the ozone absorption coefficient
    4. rsquared: r2 of the langley fit
    5. nobs: number of observations used in the fit
    6. stdo3: standard deviation of the half-day ozone
    7. meano3: mean ozone of the half-day
    8. filter_corr_3,filter_corr_4,filter_corr_5: ETC filter corrections, either calculated in filters run or supplied passing &urlconfig&fc=1,2,3,4,5 in the url
    9. info: a brief description of the calculation performed

ETC filter corrections from ozone differences

  • Name: fcor_o3dif
  • Description: this function returns the ETC filter corrections (with respect to the mean of filters 1 and 2) as calculated from the ozone differences of near-simultaneous measurements taken with different filters. Mathematically, for two measurements taken with different filters at almost the same airmass μ and with a ΔO3L1 difference in the (uncorrected) ozone, Δfiltercorr=ΔO3L1∗μ∗abscoeff∗10
  • Input:
    1. date: see Date in common inputs. If not specified, the calculation will be carried out for all the date ranges between the configurations available in Eubrewnet
    2. enddate: see Enddate in common inputs. If not specified, the current date is used.
    3. format: see Format
    4. maxminutes: maximum difference in minutes between consecutive measurements, default value: 15
    5. airmass: maximum airmass between consecutive measurements, default value: 0.1
  • Output:
    1. config_f1_cor, config_f2_cor, config_f3_cor, config_f4_cor , config_f5_cor: ETC filter corrections set in the configuration at Eubrewnet
    2. calc_f1_cor, calc_f2_cor, calc_f3_cor, calc_f4_cor, calc_f5_cor: median of all the ETC filter corrections calculated in the period, referred to the mean of filters 1 and 2
    3. calc_f1_sd, calc_f2_sd, calc_f3_sd, calc_f4_sd, calc_f5_sd: standard deviation of the ETC filter corrections in the period
    4. calc_f01_nobs, calc_f12_nobs, calc_f23_nobs, calc_f34_nobs, calc_f45_nobs: number of near simultaneous observations with the indicated filters

ETC filter corrections from ozone differences -- raw data

  • Name: fcor_o3dif_raw
  • Description: this function returns the raw data used in the fcor_o3dif function. That is, it returns Δfiltercorr=ΔO3L1∗μ∗abscoeff∗10 for each pair of measurements, without referring the result to filters 1 and 2 nor calculating any median over a date range.
  • Input:
    1. date: see Date in common inputs. If not specified, the calculation will be carried out for all the date ranges between the configurations available in Eubrewnet
    2. enddate: see Enddate in common inputs. If not specified, the current date is used.
    3. format: see Format
    4. maxminutes: maximum difference in minutes between consecutive measurements, default value: 15
    5. airmass: maximum airmass between consecutive measurements, default value: 0.1
  • Output:
    1. previous_gmt: date and time of the first measurement in the pair
    2. gmt: date and time of the second measurement
    3. previous_airmass: airmass of the first measurement
    4. airmass: airmass of the second measurement
    5. previous_filt: filter position for the first measurement
    6. filt: filter position for the second measurement
    7. previous_o3_l1: level 1 ozone for the first measurement
    8. o3_l1: level 1 ozone for the second measurement
    9. configdate: date of the operative configuration
    10. o3o3rate: ozone absorption coefficient
    11. fcor: filter correction
    12. fchange: a single string with both filter positions

ETC filter correction from FIOAVG data

  • Name: fcor_fioavg
  • Description: this function returns the ETC filter corrections (with respect to the mean of filters 1 and 2) as calculated from the results of the Brewer filter tests with the standard lamp, as found in the FIOAVG file.
  • Input:
    1. date: see Date in common inputs. If not specified, the calculation will be carried out for all the date ranges between the configurations available in Eubrewnet
    2. enddate: see Enddate in common inputs. If not specified, the current date is used.
    3. format: see Format
  • Output:
    1. config_f1_cor, config_f2_cor, config_f3_cor, config_f4_cor , config_f5_cor: ETC filter corrections set in the configuration at Eubrewnet
    2. calc_f1_cor, calc_f2_cor, calc_f3_cor, calc_f4_cor, calc_f5_cor: median of all the ETC filter corrections calculated in the period, referred to the mean of filters 1 and 2
    3. calc_f1_sd, calc_f2_sd, calc_f3_sd, calc_f4_sd, calc_f5_sd: standard deviation of the ETC filter corrections in the period
    4. calc_nobs: number of measurements in the FIOAVG file

codes/calibration.txt · Last modified: 2023/01/27 11:18 by jlsolano