User Tools

Site Tools


tools:actools:cmaq_reduced

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
tools:actools:cmaq_reduced [2017/02/01 09:39]
ctena [Content]
tools:actools:cmaq_reduced [2017/02/14 15:41] (current)
ctena [8h_max]
Line 10: Line 10:
   * Reduced specifications   * Reduced specifications
  
-=== Code ===+==== Code ====
 //cmaq_reduce.py// is fully in Python 2.7.x language and uses some libraries: //cmaq_reduce.py// is fully in Python 2.7.x language and uses some libraries:
   * ConfigParser -> To read the configuration file.   * ConfigParser -> To read the configuration file.
Line 22: Line 22:
 To run the script is needed two more files to specify What to do (Configuration file) and How to do it (Reduced specifications). To run the script is needed two more files to specify What to do (Configuration file) and How to do it (Reduced specifications).
  
-=== Configuration file ===+==== Configuration file ====
 //cmaq_reduce.conf// is the configuration file needed to run the script. You must give the configuration file as an argument <code>./cmaq_reduced.py cmaq_reduced.conf</code> //cmaq_reduce.conf// is the configuration file needed to run the script. You must give the configuration file as an argument <code>./cmaq_reduced.py cmaq_reduced.conf</code>
  
Line 60: Line 60:
   * output_directory -> Path to store the outputs. In this path will be generated the conventions folders to store all the outputs.   * output_directory -> Path to store the outputs. In this path will be generated the conventions folders to store all the outputs.
  
 +
 +==== Reduced specifications ====
 +The General section of the configuration file have the csv_path parameter that have to point to this file.
 +
 +As you are supposing these specifications are under a CSV file like this:
 +<code>
 +var,var_cmaq,filename,factor,unit,timeint_stat
 +sconco3,O3,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,1,ppm,1h 8h_max day_mean day_max
 +sconcco,CO,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,1,ppm,1h 8h_max day_mean day_max
 +sconcno2,NO2,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,1,ppm,1h 8h_max day_mean day_max
 +sconcso2,SO2,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,1,ppm,1h 8h_max day_mean day_max
 +sconcno,NO,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,1,ppm,1h
 +sconchono,HONO,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,1,ppm,1h
 +sconcnh3,NH3,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,1,ppm,1h
 +sconchno3,HNO3,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,1,ppm,1h
 +pm2p5,ASO4J ASO4I ANH4J ANH4I ANO3J ANO3I AORGPAJ AORGPAI AALKJ AXYL1J AXYL2J AXYL3J ATOL1J ATOL2J ATOL3J ABNZ1J ABNZ2J ABNZ3J AOLGAJ ATRP1J ATRP2J AISO1J AISO2J AISO3J ASQTJ AOLGBJ AORGCJ AECJ AECI A25J ANAJ ACLJ ACLI,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,0.00000001,kgm-3,1h day_mean day_max
 +pm2p5so4,ASO4I ASO4J,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,0.00000001,kgm-3,1h
 +pm2p5ppm,A25J AECJ AECI AORGPAJ AORGPAI AECI AECJ,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,0.00000001,kgm-3,1h
 +pm2p5oc,AORGPAJ AORGPAI AALKJ AXYL1J AXYL2J AXYL3J ATOL1J ATOL2J ATOL3J ABNZ1J ABNZ2J ABNZ3J AOLGAJ ATRP1J ATRP2J AISO1J AISO2J AISO3J ASQTJ AOLGBJ AORGCJ,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,0.00000001,kgm-3,1h
 +pm2p5no3,ANO3I ANO3J,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,0.00000001,kgm-3,1h
 +pm2p5nh4,ANH4I ANH4J,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,0.00000001,kgm-3,1h
 +pm2p5na,ANAJ,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,0.00000001,kgm-3,1h
 +pm2p5cl,ACLJ ACLI,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,0.00000001,kgm-3,1h
 +pm2p5bc,AECI AECJ,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,0.00000001,kgm-3,1h
 +pm2p5du,,,,,
 +pm2p5soa,AALKJ AXYL1J AXYL2J AXYL3J ATOL1J ATOL2J ATOL3J ABNZ1J ABNZ2J ABNZ3J AOLGAJ,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,0.00000001,kgm-3,1h
 +pm10,ASO4J ASO4I ANH4J ANH4I ANO3J ANO3I AORGPAJ AORGPAI AALKJ AXYL1J AXYL2J AXYL3J ATOL1J ATOL2J ATOL3J ABNZ1J ABNZ2J ABNZ3J AOLGAJ ATRP1J ATRP2J AISO1J AISO2J AISO3J ASQTJ AOLGBJ AORGCJ AECJ AECI A25J ANAJ ACLJ ACLI ASO4K ANH4K ANO3K ANAK ACLK ACORS ASOIL,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,0.00000001,kgm-3,1h day_mean day_max
 +pm10so4,ASO4I ASO4J,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,0.00000001,kgm-3,1h
 +pm10ppm,A25J AECJ AECI AORGPAJ AORGPAI ACORS ASOIL AECI AECJ,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,0.00000001,kgm-3,1h
 +pm10oc,AALKJ AXYL1J AXYL2J AXYL3J ATOL1J ATOL2J ATOL3J ABNZ1J ABNZ2J ABNZ3J AOLGAJ ATRP1J ATRP2J AISO1J AISO2J AISO3J ASQTJ AOLGBJ AORGCJ AORGPAJ AORGPAI,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,0.00000001,kgm-3,1h
 +pm10no3,ANO3I ANO3J ANO3K,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,0.00000001,kgm-3,1h
 +pm10nh4,ANH4I ANH4J ANH4K,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,0.00000001,kgm-3,1h
 +pm10na,ANAJ ANAK,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,0.00000001,kgm-3,1h
 +pm10cl,ACLJ ACLI ACLK,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,0.00000001,kgm-3,1h
 +pm10bc,AECI AECJ,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,0.00000001,kgm-3,1h
 +pm10du,,,,,
 +pm10soa,AALKJ AXYL1J AXYL2J AXYL3J ATOL1J ATOL2J ATOL3J ABNZ1J ABNZ2J ABNZ3J AOLGAJ ATRP1J ATRP2J AISO1J AISO2J AISO3J ASQTJ AOLGBJ AORGCJ,CCTM_FORECAST_<domain>.ACONC.FORECAST_<domain>,0.00000001,kgm-3,1h
 +vmro3,O3,CCTM_FORECAST_<domain>.CONC.FORECAST_<domain>,1,mol/mol,1h
 +vmrno2,NO2,CCTM_FORECAST_<domain>.CONC.FORECAST_<domain>,1,mol/mol,1h
 +vmrso2,SO2,CCTM_FORECAST_<domain>.CONC.FORECAST_<domain>,1,mol/mol,1h
 +t2,TEMP2,METCRO2D_<domain>,1,K,1h
 +pblh,PBL,METCRO2D_<domain>,1,m,1h
 +wspd10,WSPD10,METCRO2D_<domain>,1,ms-1,1h
 +sdir10,WDIR10,METCRO2D_<domain>,1,deg,1h
 +psfc,PRSFC,METCRO2D_<domain>,1,Pa,1h
 +ustar,USTAR,METCRO2D_<domain>,1,ms-1,1h
 +t,TA,METCRO3D_<domain>,1,K,1h
 +pint,PRES,METCRO3D_<domain>,1,Pa,1h
 +z,ZH,METCRO3D_<domain>,1,m,1h
 +zf,ZF,METCRO3D_<domain>,1,m,1h
 +</code>
 +
 +In the first line it has the header that specifies the content of the following lines.
 +<code>
 +var,var_cmaq,filename,factor,unit,timeint_stat
 +</code>
 +  * var -> Name of the output variable following the conventions.
 +  * var_cmaq -> Variable or List of variables of the CMAQ file to reduce.
 +  * filename -> CMAQ file to reduce.
 +  * factor -> Conversion factor to store the value in the correct units.
 +  * unit -> Units to save the outputs.
 +  * timeint_stat -> Type (or list of types) of the reduceds to do. the Reduced types will be explained in the REDUCED TYPES SECTION.
 +
 +===== Reduced types =====
 +All the output types are saved following the NetCDF conventions with the lat lon localization and the datetime.
 +==== 1h ====
 +This output type is directly the hourly values.
 +
 +The output type will have the same time steps of the input.
 +==== 8h_max ====
 +That's the most complex reduced type.
 +
 +The result is one value for each cell and each day.
 +This result is the maximum daily 8 hour mean.
 +
 +DIRECTIVE 2008/50/EC OF THE EUROPEAN PARLIAMENT AND OF THE COUNCIL
 +
 +Maximum daily eight-hour mean: the maximum daily eight-hour mean concentration shall be selected by examining eight-hour running averages, calculated from hourly data and updated each hour. Each eight -hour average so calculated shall be assigned to the day on which it ends. i.e. the first calculation period for any one day will be the period from 17:00 on the previous day to 01:00 on that day; the last calculation  period for any one day will be the period from 16:00 to 24:00 on the day.
 +
 +==== day_mean ====
 +The mean of each day
 +==== day_max ====
 +The maximum value of each day
tools/actools/cmaq_reduced.1485941996.txt.gz ยท Last modified: 2017/02/01 09:39 by ctena