| 283 | | |
| 284 | | |
| | 283 | === GRIB2FLEXPART === |
| | 284 | |
| | 285 | GRIB2FLEXPART is an add-on utility that, ''currently only for non-depositing species'', allows the user to pre-process each of the meteorological files and have it available for multiple runs. The utility makes use of FLEXPART source code (getfields.f90 and related subroutines) and generates a set of intermediate files in netCDF4 that can be optionally used to run FLEXPART. The utility consists of the following: |
| | 286 | |
| | 287 | * ''grib2flexpart'' directory: includes the main routines that a) will handle the calling of the FLEXPART routines to convert calculate and interpolate the meteorological data needed to drive FLEXPART, b) the routines to produce the intermediate output in netCDF4 format and c) a test_preprocmet_mod routine that includes testing routines. |
| | 288 | * ''preproc_gridcheck.f90'' and ''preproc_getfields.f90'' (in the usual src directory of the FLEXPART distribution): these two routines substitute the usual ''getfields.f90'' and ''readwind.f90'' ones and directly read in the intermediate files generated by the user in a former step |
| | 289 | * An adapted makefile that includes the dependencies for the compilation of the preprocessing tool and those to use the pre-processed fields. |
| | 290 | |
| | 291 | The utility is thought to process one grib file at a time. The user should generate a batch script for generating multiple files or a sequence of files. The tool is called by: |
| | 292 | |
| | 293 | {{{ |
| | 294 | grib2flexpart <path to input gribfile> <path to output NetCDF4 file> |
| | 295 | }}} |
| | 296 | |
| | 297 | Once the intermediate files are generated, the user will need to prepare the usual AVAILABLE in the usual format and simply adapt the pathnames file to point to where the netCDF4 files are stored and include in COMMAND the option “1” for the preproc_met field. In the namelist format, preprocmet is set to 0 by default. |
| | 298 | |
| | 299 | Limitations: |
| | 300 | |
| | 301 | * In FLEXPART, the internal vertical coordinate system is initialised from the first meteorological file read in (that is, from its surface pressure and the orography field at the lower left corner of the met domain). If this is not treated in a consistent way, and the first file read in is different every time, then the vertical coordinate system is defined differently with every read in and, when passed through the testing environment, differences will appear. We have opted for looking for a trade-off solution and assume the potential l inconsistency in the treatment of the vertical coordinates. In this option we guarantee that only one single day of data needs to be processed (for backward and forward) nevertheless, although the number of internal vertical levels is kept the same, each time a file is generated a new system is initialised. For this case, each meteorological file (or in case of nests, mother and nest at the same time) may be processed via command line. |
| | 302 | * Depositing species cannot be handled in this approach mainly because of the deposition parameters calculated in calcpar without actually reading information about the species to be simulated. |
| | 303 | |
| | 304 | === CTBTO specific software === |
| | 305 | |
| | 306 | As explained in Ticket on CTBTO specific code in Fp10 . For Version 9.3 the implementation of such code was made using some compile IFDEF directives. However, for version 10 and, after consultation with the Flexpart developers, a better solution has been found: a new logical variable “lctbto” has been added in ''par_mod.f90'' file. The variable is by default set to False and only those who want to use the CTBTO specific code will have to set it to True and recompile. The specific code includes the following: |
| | 307 | |
| | 308 | * Add the possibility to generate receptors output in backward runs |
| | 309 | * Add the possibility to obtain concentration as backwards unit (when using IND_SOURCE and IND_RECEPTOR set to one) - this includes the 10**12 factor as well. |
| | 310 | * Allow for aggregated output in backward mode (allowing to have IFOR = 0 for bwd runs) |
| | 311 | * Small modification in the welcome string indicating you are using the CTBTO specific version |
| | 312 | |
| | 313 | == TODOs == |
| | 314 | Dynamical allocation of many of the meteorological variables has been tested and prototyped. This has not yet been fully implemented and remains as a TODO. |
| | 315 | |
