Scripts to create emissions files in the right format for PREP PGD
The following scripts are designed to convert the original emission files taken from the GEIA/ECCAD international program (http://eccad.sedoo.fr) to the latitude-longitude-value ascii format needed by the PREP_PGD .
Emissions files are downloaded from the ECCAD web site in the netcdf format.
The script emiss_conv_nc_to_asc.ncl converts netcdf format to lat-lon-val with final emissions in molec/cm2/s.
The script needs to be adapted to the user requirements:
- The molecular mass has to be changed according to the considered molecule.
The following script is an example of how to change the molecular mass for the conversion BC_conversion_nc_to_asc_O3R.ncl (source: S. Strada)
- Emissions in the PREP_PGD program should be in
- ppp m/s
- All the emissions files have to be in the same units (i.e. you cannot mix emissions in molec/cm2/s and in ppp m/s in the same PRE_PGD1.nam)
When running the model, the routine build_emisstabn.F90 in SURFEX will automatically make the conversion of the above mentioned units to molec/m2/s according to the keyword provided in the chemistry file under the EMISUNIT keyword: MOL (if emissions are provided in microMol/m2/day), MIX (if emissions are provided in ppp m/s), CON (emissions are provided in molec/cm2/s).
EMISUNIT emission fluxes in microMole/m2/day MOL
EMISUNIT emission fluxes in ppp m/s MIX
EMISUNIT emission fluxes in molec/cm2/s CON