readcommand Subroutine

public subroutine readcommand()

                                                                        *
 This routine reads the user specifications for the current model run.  *
                                                                        *
 Author: A. Stohl                                                       *
                                                                        *
 18 May 1996                                                            *
 HSO, 1 July 2014                                                       *
 Added optional namelist input                                          *
                                                                        * 
                                                                        *
 January 2024 Rona Thompson                                             *
 Added new variables for LCM                                            *
                                                                        *

                                                                        *

Variables: * bdate beginning date as Julian date * ctl factor by which time step must be smaller than * Lagrangian time scale * ibdate,ibtime beginnning date and time (YYYYMMDD, HHMISS) * ideltas [s] modelling period * iedate,ietime ending date and time (YYYYMMDD, HHMISS) * ifine reduction factor for vertical wind time step * outputforeachrel for forward runs it is possible either to create * one outputfield or several for each releasepoint * iflux switch to turn on (1)/off (0) flux calculations * iout 1 for conc. (residence time for backward runs) output, 2 for mixing ratio output, 3 both, 4 for plume * trajectory output, 5 = options 1 and 4 * ipin 1 continue simulation with restart.bin file, * 2 continue simulaion with dumped particle data, 0 no * 3 use self-defined initial conditions in netcdf * 4 initial run using option 3, restart from restart.bin ipout 0 no particle dump, 1 every output time, 3 only at end* ipoutfac increase particle dump interval by factor (default 1) * loutaver [s] concentration output is an average over loutaver * seconds * loutsample [s] average is computed from samples taken every [s] * seconds * loutstep [s] time interval of concentration output * lrecoutstep [s] time interval of receptor output * lrecoutaver [s] receptor output is an average of lrecoutaver seconds * lrecoutsample [s] average is computed from samples taken every [s] * lsynctime [s] synchronisation time interval for all particles * lagespectra switch to turn on (1)/off (0) calculation of age * spectra * lconvection value of either 0 and 1 indicating mixing by * convection * = 0 .. no convection * + 1 .. parameterisation of mixing by subgrid-scale * convection = on * lsubgrid switch to turn on (1)/off (0) subgrid topography * parameterization * method method used to compute the particle pseudovelocities * mdomainfill 1 use domain-filling option, 0 not, 2 use strat. O3 * * Constants: * unitcommand unit connected to file COMMAND * *

AF Set the switches IND_REL and IND_SAMP for the release and sampling Af switches for the releasefile: Af IND_REL = 1 : xmass * rho Af IND_REL = 0 : xmass * 1

MDOMAINFILL = 1 | LLCMOUTPUT = true IND_SOURCE = 1 | IND_SAMP = 0 IND_RECEPTOR = 1 | calculates mass ratio mixing ratio IOUT = 2 | as ratio species_mass to airtracer_mass

Check whether valid options have been chosen in file COMMAND

if (((iout.eq.4).or.(iout.eq.5)).and.((ipin.eq.3).or.(ipin.eq.4))) then write(,) ' #### FLEXPART MODEL ERROR! FILE COMMAND: #### ' write(,) ' #### IOUT CANNOT BE 4 or 5 (plume) WHEN #### ' write(,) ' #### READING FROM part_ic.nc (ipin=4/5) #### ' error stop endif

Arguments