[92fab65] | 1 | ! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt |
---|
| 2 | ! SPDX-License-Identifier: GPL-3.0-or-later |
---|
[332fbbd] | 3 | |
---|
[e200b7a] | 4 | subroutine readcommand |
---|
| 5 | |
---|
| 6 | !***************************************************************************** |
---|
| 7 | ! * |
---|
| 8 | ! This routine reads the user specifications for the current model run. * |
---|
| 9 | ! * |
---|
| 10 | ! Author: A. Stohl * |
---|
| 11 | ! * |
---|
| 12 | ! 18 May 1996 * |
---|
[b4d29ce] | 13 | ! HSO, 1 July 2014 * |
---|
| 14 | ! Added optional namelist input * |
---|
[e200b7a] | 15 | ! * |
---|
| 16 | !***************************************************************************** |
---|
| 17 | ! * |
---|
| 18 | ! Variables: * |
---|
| 19 | ! bdate beginning date as Julian date * |
---|
| 20 | ! ctl factor by which time step must be smaller than * |
---|
| 21 | ! Lagrangian time scale * |
---|
| 22 | ! ibdate,ibtime beginnning date and time (YYYYMMDD, HHMISS) * |
---|
| 23 | ! ideltas [s] modelling period * |
---|
| 24 | ! iedate,ietime ending date and time (YYYYMMDD, HHMISS) * |
---|
| 25 | ! ifine reduction factor for vertical wind time step * |
---|
| 26 | ! outputforeachrel for forward runs it is possible either to create * |
---|
| 27 | ! one outputfield or several for each releasepoint * |
---|
| 28 | ! iflux switch to turn on (1)/off (0) flux calculations * |
---|
| 29 | ! iout 1 for conc. (residence time for backward runs) output,* |
---|
| 30 | ! 2 for mixing ratio output, 3 both, 4 for plume * |
---|
| 31 | ! trajectory output, 5 = options 1 and 4 * |
---|
| 32 | ! ipin 1 continue simulation with dumped particle data, 0 no * |
---|
| 33 | ! ipout 0 no particle dump, 1 every output time, 3 only at end* |
---|
[0a94e13] | 34 | ! ipoutfac increase particle dump interval by factor (default 1) * |
---|
[e200b7a] | 35 | ! itsplit [s] time constant for particle splitting * |
---|
| 36 | ! loutaver [s] concentration output is an average over loutaver * |
---|
| 37 | ! seconds * |
---|
| 38 | ! loutsample [s] average is computed from samples taken every [s] * |
---|
| 39 | ! seconds * |
---|
| 40 | ! loutstep [s] time interval of concentration output * |
---|
| 41 | ! lsynctime [s] synchronisation time interval for all particles * |
---|
| 42 | ! lagespectra switch to turn on (1)/off (0) calculation of age * |
---|
| 43 | ! spectra * |
---|
| 44 | ! lconvection value of either 0 and 1 indicating mixing by * |
---|
| 45 | ! convection * |
---|
| 46 | ! = 0 .. no convection * |
---|
| 47 | ! + 1 .. parameterisation of mixing by subgrid-scale * |
---|
| 48 | ! convection = on * |
---|
| 49 | ! lsubgrid switch to turn on (1)/off (0) subgrid topography * |
---|
| 50 | ! parameterization * |
---|
| 51 | ! method method used to compute the particle pseudovelocities * |
---|
| 52 | ! mdomainfill 1 use domain-filling option, 0 not, 2 use strat. O3 * |
---|
| 53 | ! * |
---|
| 54 | ! Constants: * |
---|
| 55 | ! unitcommand unit connected to file COMMAND * |
---|
| 56 | ! * |
---|
| 57 | !***************************************************************************** |
---|
| 58 | |
---|
| 59 | use par_mod |
---|
| 60 | use com_mod |
---|
| 61 | |
---|
| 62 | implicit none |
---|
| 63 | |
---|
| 64 | real(kind=dp) :: juldate |
---|
| 65 | character(len=50) :: line |
---|
| 66 | logical :: old |
---|
[f13406c] | 67 | integer :: readerror |
---|
| 68 | |
---|
| 69 | namelist /command/ & |
---|
[4fbe7a5] | 70 | ldirect, & |
---|
| 71 | ibdate,ibtime, & |
---|
| 72 | iedate,ietime, & |
---|
| 73 | loutstep, & |
---|
| 74 | loutaver, & |
---|
| 75 | loutsample, & |
---|
| 76 | itsplit, & |
---|
| 77 | lsynctime, & |
---|
| 78 | ctl, & |
---|
| 79 | ifine, & |
---|
| 80 | iout, & |
---|
| 81 | ipout, & |
---|
[0a94e13] | 82 | ipoutfac, & |
---|
[4fbe7a5] | 83 | lsubgrid, & |
---|
| 84 | lconvection, & |
---|
| 85 | lagespectra, & |
---|
| 86 | ipin, & |
---|
| 87 | ioutputforeachrelease, & |
---|
| 88 | iflux, & |
---|
| 89 | mdomainfill, & |
---|
| 90 | ind_source, & |
---|
| 91 | ind_receptor, & |
---|
| 92 | mquasilag, & |
---|
| 93 | nested_output, & |
---|
| 94 | linit_cond, & |
---|
[8a65cb0] | 95 | lnetcdfout, & |
---|
| 96 | surf_only, & |
---|
[78e62dc] | 97 | cblflag, & |
---|
[2eefa58] | 98 | linversionout, & |
---|
[25b4532] | 99 | ohfields_path, & |
---|
| 100 | d_trop, & |
---|
| 101 | d_strat |
---|
[f13406c] | 102 | |
---|
| 103 | ! Presetting namelist command |
---|
[b4d29ce] | 104 | ldirect=0 |
---|
[f13406c] | 105 | ibdate=20000101 |
---|
| 106 | ibtime=0 |
---|
| 107 | iedate=20000102 |
---|
| 108 | ietime=0 |
---|
| 109 | loutstep=10800 |
---|
| 110 | loutaver=10800 |
---|
| 111 | loutsample=900 |
---|
| 112 | itsplit=999999999 |
---|
| 113 | lsynctime=900 |
---|
| 114 | ctl=-5.0 |
---|
| 115 | ifine=4 |
---|
| 116 | iout=3 |
---|
| 117 | ipout=0 |
---|
[0a94e13] | 118 | ipoutfac=1 |
---|
[f13406c] | 119 | lsubgrid=1 |
---|
| 120 | lconvection=1 |
---|
| 121 | lagespectra=0 |
---|
| 122 | ipin=1 |
---|
| 123 | ioutputforeachrelease=1 |
---|
| 124 | iflux=1 |
---|
| 125 | mdomainfill=0 |
---|
| 126 | ind_source=1 |
---|
| 127 | ind_receptor=1 |
---|
| 128 | mquasilag=0 |
---|
| 129 | nested_output=0 |
---|
| 130 | linit_cond=0 |
---|
[8a65cb0] | 131 | lnetcdfout=0 |
---|
[f13406c] | 132 | surf_only=0 |
---|
[c04b739] | 133 | cblflag=0 ! if using old-style COMMAND file, set to 1 here to use mc cbl routine |
---|
[2eefa58] | 134 | linversionout=0 |
---|
[78e62dc] | 135 | ohfields_path="../../flexin/" |
---|
[e200b7a] | 136 | |
---|
[b5127f9] | 137 | !Af set release-switch |
---|
| 138 | WETBKDEP=.false. |
---|
| 139 | DRYBKDEP=.false. |
---|
| 140 | |
---|
[e200b7a] | 141 | ! Open the command file and read user options |
---|
[f13406c] | 142 | ! Namelist input first: try to read as namelist file |
---|
| 143 | !************************************************************************** |
---|
[b4d29ce] | 144 | open(unitcommand,file=path(1)(1:length(1))//'COMMAND',status='old',form='formatted',err=999) |
---|
[e200b7a] | 145 | |
---|
[b4d29ce] | 146 | ! try namelist input (default) |
---|
| 147 | read(unitcommand,command,iostat=readerror) |
---|
[f13406c] | 148 | close(unitcommand) |
---|
[e200b7a] | 149 | |
---|
[b4d29ce] | 150 | ! distinguish namelist from fixed text input |
---|
| 151 | if ((readerror.ne.0).or.(ldirect.eq.0)) then ! parse as text file format |
---|
| 152 | |
---|
| 153 | open(unitcommand,file=path(1)(1:length(1))//'COMMAND',status='old', err=999) |
---|
[f13406c] | 154 | |
---|
| 155 | ! Check the format of the COMMAND file (either in free format, |
---|
| 156 | ! or using formatted mask) |
---|
| 157 | ! Use of formatted mask is assumed if line 10 contains the word 'DIRECTION' |
---|
| 158 | !************************************************************************** |
---|
| 159 | |
---|
| 160 | call skplin(9,unitcommand) |
---|
| 161 | read (unitcommand,901) line |
---|
| 162 | 901 format (a) |
---|
| 163 | if (index(line,'LDIRECT') .eq. 0) then |
---|
| 164 | old = .false. |
---|
[4c64400] | 165 | if (lroot) write(*,*) 'COMMAND in old short format, & |
---|
| 166 | &please update to namelist format' |
---|
[f13406c] | 167 | else |
---|
| 168 | old = .true. |
---|
[4c64400] | 169 | if (lroot) write(*,*) 'COMMAND in old long format, & |
---|
| 170 | &please update to namelist format' |
---|
[f13406c] | 171 | endif |
---|
| 172 | rewind(unitcommand) |
---|
| 173 | |
---|
[4fbe7a5] | 174 | |
---|
[f13406c] | 175 | ! Read parameters |
---|
| 176 | !**************** |
---|
| 177 | |
---|
| 178 | call skplin(7,unitcommand) |
---|
| 179 | if (old) call skplin(1,unitcommand) |
---|
| 180 | read(unitcommand,*) ldirect |
---|
| 181 | if (old) call skplin(3,unitcommand) |
---|
| 182 | read(unitcommand,*) ibdate,ibtime |
---|
| 183 | if (old) call skplin(3,unitcommand) |
---|
| 184 | read(unitcommand,*) iedate,ietime |
---|
| 185 | if (old) call skplin(3,unitcommand) |
---|
| 186 | read(unitcommand,*) loutstep |
---|
| 187 | if (old) call skplin(3,unitcommand) |
---|
| 188 | read(unitcommand,*) loutaver |
---|
| 189 | if (old) call skplin(3,unitcommand) |
---|
| 190 | read(unitcommand,*) loutsample |
---|
| 191 | if (old) call skplin(3,unitcommand) |
---|
| 192 | read(unitcommand,*) itsplit |
---|
| 193 | if (old) call skplin(3,unitcommand) |
---|
| 194 | read(unitcommand,*) lsynctime |
---|
| 195 | if (old) call skplin(3,unitcommand) |
---|
| 196 | read(unitcommand,*) ctl |
---|
| 197 | if (old) call skplin(3,unitcommand) |
---|
| 198 | read(unitcommand,*) ifine |
---|
| 199 | if (old) call skplin(3,unitcommand) |
---|
| 200 | read(unitcommand,*) iout |
---|
| 201 | if (old) call skplin(3,unitcommand) |
---|
| 202 | read(unitcommand,*) ipout |
---|
| 203 | if (old) call skplin(3,unitcommand) |
---|
| 204 | read(unitcommand,*) lsubgrid |
---|
| 205 | if (old) call skplin(3,unitcommand) |
---|
| 206 | read(unitcommand,*) lconvection |
---|
| 207 | if (old) call skplin(3,unitcommand) |
---|
| 208 | read(unitcommand,*) lagespectra |
---|
| 209 | if (old) call skplin(3,unitcommand) |
---|
| 210 | read(unitcommand,*) ipin |
---|
| 211 | if (old) call skplin(3,unitcommand) |
---|
| 212 | read(unitcommand,*) ioutputforeachrelease |
---|
| 213 | if (old) call skplin(3,unitcommand) |
---|
| 214 | read(unitcommand,*) iflux |
---|
| 215 | if (old) call skplin(3,unitcommand) |
---|
| 216 | read(unitcommand,*) mdomainfill |
---|
| 217 | if (old) call skplin(3,unitcommand) |
---|
| 218 | read(unitcommand,*) ind_source |
---|
| 219 | if (old) call skplin(3,unitcommand) |
---|
| 220 | read(unitcommand,*) ind_receptor |
---|
| 221 | if (old) call skplin(3,unitcommand) |
---|
| 222 | read(unitcommand,*) mquasilag |
---|
| 223 | if (old) call skplin(3,unitcommand) |
---|
| 224 | read(unitcommand,*) nested_output |
---|
| 225 | if (old) call skplin(3,unitcommand) |
---|
| 226 | read(unitcommand,*) linit_cond |
---|
[4fbe7a5] | 227 | if (old) call skplin(3,unitcommand) |
---|
| 228 | read(unitcommand,*) surf_only |
---|
[c04b739] | 229 | ! Removed for backwards compatibility. |
---|
| 230 | ! if (old) call skplin(3,unitcommand) !added by mc |
---|
| 231 | ! read(unitcommand,*) cblflag !added by mc |
---|
[78e62dc] | 232 | |
---|
[f13406c] | 233 | close(unitcommand) |
---|
| 234 | |
---|
| 235 | endif ! input format |
---|
| 236 | |
---|
| 237 | ! write command file in namelist format to output directory if requested |
---|
[8a65cb0] | 238 | if (nmlout.and.lroot) then |
---|
[f13406c] | 239 | open(unitcommand,file=path(2)(1:length(2))//'COMMAND.namelist',err=1000) |
---|
| 240 | write(unitcommand,nml=command) |
---|
| 241 | close(unitcommand) |
---|
[e200b7a] | 242 | endif |
---|
| 243 | |
---|
| 244 | ifine=max(ifine,1) |
---|
| 245 | |
---|
| 246 | ! Determine how Markov chain is formulated (for w or for w/sigw) |
---|
| 247 | !*************************************************************** |
---|
[8a65cb0] | 248 | if (cblflag.eq.1) then !---- added by mc to properly set parameters for CBL simulations |
---|
[e200b7a] | 249 | turbswitch=.true. |
---|
[8a65cb0] | 250 | if (lsynctime>maxtl) lsynctime=maxtl !maxtl defined in com_mod.f90 |
---|
| 251 | if (ctl.lt.5) then |
---|
| 252 | print *,'WARNING: CBL flag active the ratio of TLu/dt has been set to 5' |
---|
| 253 | ctl=5. |
---|
| 254 | end if |
---|
| 255 | if (ifine*ctl.lt.50) then |
---|
| 256 | ifine=int(50./ctl)+1 |
---|
| 257 | |
---|
| 258 | print *,'WARNING: CBL flag active the ratio of TLW/dt was < 50, ifine has been re-set to',ifine |
---|
| 259 | !pause |
---|
| 260 | endif |
---|
| 261 | print *,'WARNING: CBL flag active the ratio of TLW/dt is ',ctl*ifine |
---|
| 262 | print *,'WARNING: CBL flag active lsynctime is ',lsynctime |
---|
| 263 | else !added by mc |
---|
| 264 | if (ctl.ge.0.1) then |
---|
| 265 | turbswitch=.true. |
---|
| 266 | else |
---|
| 267 | turbswitch=.false. |
---|
| 268 | ifine=1 |
---|
| 269 | endif |
---|
| 270 | endif !added by mc |
---|
[e200b7a] | 271 | fine=1./real(ifine) |
---|
| 272 | ctl=1./ctl |
---|
| 273 | |
---|
| 274 | ! Set the switches required for the various options for input/output units |
---|
| 275 | !************************************************************************* |
---|
| 276 | !AF Set the switches IND_REL and IND_SAMP for the release and sampling |
---|
| 277 | !Af switches for the releasefile: |
---|
| 278 | !Af IND_REL = 1 : xmass * rho |
---|
| 279 | !Af IND_REL = 0 : xmass * 1 |
---|
| 280 | |
---|
| 281 | !Af switches for the conccalcfile: |
---|
| 282 | !AF IND_SAMP = 0 : xmass * 1 |
---|
| 283 | !Af IND_SAMP = -1 : xmass / rho |
---|
| 284 | |
---|
| 285 | !AF IND_SOURCE switches between different units for concentrations at the source |
---|
| 286 | !Af NOTE that in backward simulations the release of computational particles |
---|
| 287 | !Af takes place at the "receptor" and the sampling of p[articles at the "source". |
---|
| 288 | !Af 1 = mass units |
---|
| 289 | !Af 2 = mass mixing ratio units |
---|
| 290 | !Af IND_RECEPTOR switches between different units for concentrations at the receptor |
---|
| 291 | !Af 1 = mass units |
---|
| 292 | !Af 2 = mass mixing ratio units |
---|
[5deb48c] | 293 | ! 3 = wet deposition in outputfield |
---|
| 294 | ! 4 = dry deposition in outputfield |
---|
[e200b7a] | 295 | |
---|
| 296 | if ( ldirect .eq. 1 ) then ! FWD-Run |
---|
| 297 | !Af set release-switch |
---|
| 298 | if (ind_source .eq. 1 ) then !mass |
---|
| 299 | ind_rel = 0 |
---|
| 300 | else ! mass mix |
---|
| 301 | ind_rel = 1 |
---|
| 302 | endif |
---|
| 303 | !Af set sampling switch |
---|
| 304 | if (ind_receptor .eq. 1) then !mass |
---|
| 305 | ind_samp = 0 |
---|
| 306 | else ! mass mix |
---|
| 307 | ind_samp = -1 |
---|
| 308 | endif |
---|
| 309 | elseif (ldirect .eq. -1 ) then !BWD-Run |
---|
| 310 | !Af set sampling switch |
---|
| 311 | if (ind_source .eq. 1 ) then !mass |
---|
| 312 | ind_samp = -1 |
---|
| 313 | else ! mass mix |
---|
| 314 | ind_samp = 0 |
---|
| 315 | endif |
---|
[54cbd6c] | 316 | select case (ind_receptor) |
---|
| 317 | case (1) ! 1 .. concentration at receptor |
---|
[e200b7a] | 318 | ind_rel = 1 |
---|
[54cbd6c] | 319 | case (2) ! 2 .. mixing ratio at receptor |
---|
[e200b7a] | 320 | ind_rel = 0 |
---|
[54cbd6c] | 321 | case (3) ! 3 .. wet deposition in outputfield |
---|
| 322 | ind_rel = 3 |
---|
[20963b1] | 323 | if (lroot) then |
---|
| 324 | write(*,*) ' #### FLEXPART WET DEPOSITION BACKWARD MODE #### ' |
---|
| 325 | write(*,*) ' #### Releaseheight is forced to 0 - 20km #### ' |
---|
| 326 | write(*,*) ' #### Release is performed above ground lev #### ' |
---|
| 327 | end if |
---|
[54cbd6c] | 328 | WETBKDEP=.true. |
---|
| 329 | allocate(xscav_frac1(maxpart,maxspec)) |
---|
| 330 | case (4) ! 4 .. dry deposition in outputfield |
---|
| 331 | ind_rel = 4 |
---|
[20963b1] | 332 | if (lroot) then |
---|
| 333 | write(*,*) ' #### FLEXPART DRY DEPOSITION BACKWARD MODE #### ' |
---|
| 334 | write(*,*) ' #### Releaseheight is forced to 0 - 2*href #### ' |
---|
| 335 | write(*,*) ' #### Release is performed above ground lev #### ' |
---|
| 336 | end if |
---|
[54cbd6c] | 337 | DRYBKDEP=.true. |
---|
| 338 | allocate(xscav_frac1(maxpart,maxspec)) |
---|
| 339 | end select |
---|
[e200b7a] | 340 | endif |
---|
| 341 | |
---|
| 342 | !************************************************************* |
---|
| 343 | ! Check whether valid options have been chosen in file COMMAND |
---|
| 344 | !************************************************************* |
---|
| 345 | |
---|
| 346 | ! Check options for initial condition output: Switch off for forward runs |
---|
| 347 | !************************************************************************ |
---|
| 348 | |
---|
| 349 | if (ldirect.eq.1) linit_cond=0 |
---|
| 350 | if ((linit_cond.lt.0).or.(linit_cond.gt.2)) then |
---|
| 351 | write(*,*) ' #### FLEXPART MODEL ERROR! INVALID OPTION #### ' |
---|
| 352 | write(*,*) ' #### FOR LINIT_COND IN FILE "COMMAND". #### ' |
---|
| 353 | stop |
---|
| 354 | endif |
---|
| 355 | |
---|
| 356 | ! Check input dates |
---|
| 357 | !****************** |
---|
| 358 | |
---|
| 359 | if (iedate.lt.ibdate) then |
---|
| 360 | write(*,*) ' #### FLEXPART MODEL ERROR! BEGINNING DATE #### ' |
---|
| 361 | write(*,*) ' #### IS LARGER THAN ENDING DATE. CHANGE #### ' |
---|
| 362 | write(*,*) ' #### EITHER POINT 2 OR POINT 3 IN FILE #### ' |
---|
| 363 | write(*,*) ' #### "COMMAND". #### ' |
---|
| 364 | stop |
---|
| 365 | else if (iedate.eq.ibdate) then |
---|
| 366 | if (ietime.lt.ibtime) then |
---|
| 367 | write(*,*) ' #### FLEXPART MODEL ERROR! BEGINNING TIME #### ' |
---|
| 368 | write(*,*) ' #### IS LARGER THAN ENDING TIME. CHANGE #### ' |
---|
| 369 | write(*,*) ' #### EITHER POINT 2 OR POINT 3 IN FILE #### ' |
---|
| 370 | write(*,*) ' #### "COMMAND". #### ' |
---|
| 371 | stop |
---|
| 372 | endif |
---|
| 373 | endif |
---|
| 374 | |
---|
| 375 | |
---|
| 376 | ! Determine kind of dispersion method |
---|
| 377 | !************************************ |
---|
| 378 | |
---|
| 379 | if (ctl.gt.0.) then |
---|
| 380 | method=1 |
---|
| 381 | mintime=minstep |
---|
| 382 | else |
---|
| 383 | method=0 |
---|
| 384 | mintime=lsynctime |
---|
| 385 | endif |
---|
| 386 | |
---|
[a9cf4b1] | 387 | ! Check for netcdf output switch |
---|
| 388 | !******************************* |
---|
[8a65cb0] | 389 | if (iout.ge.8) then |
---|
| 390 | lnetcdfout = 1 |
---|
| 391 | iout = iout - 8 |
---|
[a9cf4b1] | 392 | #ifndef USE_NCF |
---|
| 393 | write(*,*) 'ERROR: netcdf output not activated during compile time but used in COMMAND file!' |
---|
| 394 | write(*,*) 'Please recompile with netcdf library (`make [...] ncf=yes`) or use standard output format.' |
---|
| 395 | stop |
---|
| 396 | #endif |
---|
[8a65cb0] | 397 | endif |
---|
| 398 | |
---|
[e200b7a] | 399 | ! Check whether a valid option for gridded model output has been chosen |
---|
| 400 | !********************************************************************** |
---|
| 401 | |
---|
[2eefa58] | 402 | if ((iout.lt.1).or.(iout.gt.5)) then |
---|
[e200b7a] | 403 | write(*,*) ' #### FLEXPART MODEL ERROR! FILE COMMAND: #### ' |
---|
[54cbd6c] | 404 | write(*,*) ' #### IOUT MUST BE 1, 2, 3, 4 OR 5 FOR #### ' |
---|
[8a65cb0] | 405 | write(*,*) ' #### STANDARD FLEXPART OUTPUT OR 9 - 13 #### ' |
---|
| 406 | write(*,*) ' #### FOR NETCDF OUTPUT #### ' |
---|
[e200b7a] | 407 | stop |
---|
| 408 | endif |
---|
| 409 | |
---|
| 410 | !AF check consistency between units and volume mixing ratio |
---|
| 411 | if ( ((iout.eq.2).or.(iout.eq.3)).and. & |
---|
| 412 | (ind_source.gt.1 .or.ind_receptor.gt.1) ) then |
---|
| 413 | write(*,*) ' #### FLEXPART MODEL ERROR! FILE COMMAND: #### ' |
---|
| 414 | write(*,*) ' #### VOLUME MIXING RATIO ONLY SUPPORTED #### ' |
---|
| 415 | write(*,*) ' #### FOR MASS UNITS (at the moment) #### ' |
---|
| 416 | stop |
---|
| 417 | endif |
---|
| 418 | |
---|
| 419 | |
---|
| 420 | ! For quasilag output for each release is forbidden |
---|
| 421 | !***************************************************************************** |
---|
| 422 | |
---|
| 423 | if ((ioutputforeachrelease.eq.1).and.(mquasilag.eq.1)) then |
---|
| 424 | write(*,*) '#### FLEXPART MODEL ERROR! FILE COMMAND: ####' |
---|
| 425 | write(*,*) '#### OUTPUTFOREACHRELEASE AND QUASILAGRANGIAN####' |
---|
| 426 | write(*,*) '#### MODE IS NOT POSSIBLE ! ####' |
---|
| 427 | stop |
---|
| 428 | endif |
---|
| 429 | |
---|
| 430 | |
---|
| 431 | ! For quasilag backward is forbidden |
---|
| 432 | !***************************************************************************** |
---|
| 433 | |
---|
| 434 | if ((ldirect.lt.0).and.(mquasilag.eq.1)) then |
---|
| 435 | write(*,*) '#### FLEXPART MODEL ERROR! FILE COMMAND: ####' |
---|
| 436 | write(*,*) '#### FOR BACKWARD RUNS, QUASILAGRANGIAN MODE ####' |
---|
| 437 | write(*,*) '#### IS NOT POSSIBLE ! ####' |
---|
| 438 | stop |
---|
| 439 | endif |
---|
| 440 | |
---|
| 441 | |
---|
| 442 | ! For backward runs one releasefield for all releases makes no sense, |
---|
| 443 | ! For quasilag and domainfill ioutputforechrelease is forbidden |
---|
| 444 | !***************************************************************************** |
---|
| 445 | |
---|
| 446 | if ((ldirect.lt.0).and.(ioutputforeachrelease.eq.0)) then |
---|
| 447 | write(*,*) '#### FLEXPART MODEL ERROR! FILE COMMAND: ####' |
---|
| 448 | write(*,*) '#### FOR BACKWARD RUNS, IOUTPUTFOREACHRLEASE ####' |
---|
| 449 | write(*,*) '#### MUST BE SET TO ONE! ####' |
---|
| 450 | stop |
---|
| 451 | endif |
---|
| 452 | |
---|
| 453 | |
---|
| 454 | ! For backward runs one releasefield for all releases makes no sense, |
---|
| 455 | ! and is "forbidden" |
---|
| 456 | !***************************************************************************** |
---|
| 457 | |
---|
| 458 | if ((mdomainfill.eq.1).and.(ioutputforeachrelease.eq.1)) then |
---|
| 459 | write(*,*) '#### FLEXPART MODEL ERROR! FILE COMMAND: ####' |
---|
| 460 | write(*,*) '#### FOR DOMAIN FILLING RUNS OUTPUT FOR ####' |
---|
| 461 | write(*,*) '#### EACH RELEASE IS FORBIDDEN ! ####' |
---|
| 462 | stop |
---|
| 463 | endif |
---|
| 464 | |
---|
[2eefa58] | 465 | ! Inversion output format only for backward runs |
---|
| 466 | !***************************************************************************** |
---|
| 467 | |
---|
| 468 | if ((linversionout.eq.1).and.(ldirect.eq.1)) then |
---|
| 469 | write(*,*) '#### FLEXPART MODEL ERROR! FILE COMMAND: ####' |
---|
| 470 | write(*,*) '#### INVERSION OUTPUT FORMAT ONLY FOR ####' |
---|
| 471 | write(*,*) '#### BACKWARD RUNS ####' |
---|
| 472 | stop |
---|
| 473 | endif |
---|
| 474 | |
---|
[e200b7a] | 475 | |
---|
| 476 | ! For domain-filling trajectories, a plume centroid trajectory makes no sense, |
---|
| 477 | ! For backward runs, only residence time output (iout=1) or plume trajectories (iout=4), |
---|
| 478 | ! or both (iout=5) makes sense; other output options are "forbidden" |
---|
| 479 | !***************************************************************************** |
---|
| 480 | |
---|
| 481 | if (ldirect.lt.0) then |
---|
| 482 | if ((iout.eq.2).or.(iout.eq.3)) then |
---|
| 483 | write(*,*) '#### FLEXPART MODEL ERROR! FILE COMMAND: ####' |
---|
| 484 | write(*,*) '#### FOR BACKWARD RUNS, IOUT MUST BE 1,4,OR 5####' |
---|
| 485 | stop |
---|
| 486 | endif |
---|
| 487 | endif |
---|
| 488 | |
---|
| 489 | |
---|
| 490 | ! For domain-filling trajectories, a plume centroid trajectory makes no sense, |
---|
| 491 | ! and is "forbidden" |
---|
| 492 | !***************************************************************************** |
---|
| 493 | |
---|
| 494 | if (mdomainfill.ge.1) then |
---|
| 495 | if ((iout.eq.4).or.(iout.eq.5)) then |
---|
| 496 | write(*,*) '#### FLEXPART MODEL ERROR! FILE COMMAND: ####' |
---|
| 497 | write(*,*) '#### FOR DOMAIN-FILLING TRAJECTORY OPTION, ####' |
---|
| 498 | write(*,*) '#### IOUT MUST NOT BE SET TO 4 OR 5. ####' |
---|
| 499 | stop |
---|
| 500 | endif |
---|
| 501 | endif |
---|
| 502 | |
---|
| 503 | |
---|
| 504 | |
---|
| 505 | ! Check whether a valid options for particle dump has been chosen |
---|
| 506 | !**************************************************************** |
---|
| 507 | |
---|
[0a94e13] | 508 | if ((ipout.ne.0).and.(ipout.ne.1).and.(ipout.ne.2).and.(ipout.ne.3)) then |
---|
[e200b7a] | 509 | write(*,*) ' #### FLEXPART MODEL ERROR! FILE COMMAND: #### ' |
---|
[0a94e13] | 510 | write(*,*) ' #### IPOUT MUST BE 0, 1, 2 OR 3! #### ' |
---|
[e200b7a] | 511 | stop |
---|
| 512 | endif |
---|
| 513 | |
---|
[f13406c] | 514 | if(lsubgrid.ne.1.and.verbosity.eq.0) then |
---|
[e200b7a] | 515 | write(*,*) ' ---------------- ' |
---|
| 516 | write(*,*) ' INFORMATION: SUBGRIDSCALE TERRAIN EFFECT IS' |
---|
| 517 | write(*,*) ' NOT PARAMETERIZED DURING THIS SIMULATION. ' |
---|
| 518 | write(*,*) ' ---------------- ' |
---|
| 519 | endif |
---|
| 520 | |
---|
| 521 | |
---|
| 522 | ! Check whether convection scheme is either turned on or off |
---|
| 523 | !*********************************************************** |
---|
| 524 | |
---|
| 525 | if ((lconvection.ne.0).and.(lconvection.ne.1)) then |
---|
| 526 | write(*,*) ' #### FLEXPART MODEL ERROR! FILE COMMAND: #### ' |
---|
| 527 | write(*,*) ' #### LCONVECTION MUST BE SET TO EITHER 1 OR 0#### ' |
---|
| 528 | stop |
---|
| 529 | endif |
---|
| 530 | |
---|
| 531 | |
---|
| 532 | ! Check whether synchronisation interval is sufficiently short |
---|
| 533 | !************************************************************* |
---|
| 534 | |
---|
| 535 | if (lsynctime.gt.(idiffnorm/2)) then |
---|
| 536 | write(*,*) ' #### FLEXPART MODEL ERROR! SYNCHRONISATION #### ' |
---|
| 537 | write(*,*) ' #### TIME IS TOO LONG. MAKE IT SHORTER. #### ' |
---|
| 538 | write(*,*) ' #### MINIMUM HAS TO BE: ', idiffnorm/2 |
---|
| 539 | stop |
---|
| 540 | endif |
---|
| 541 | |
---|
| 542 | |
---|
| 543 | ! Check consistency of the intervals, sampling periods, etc., for model output |
---|
| 544 | !***************************************************************************** |
---|
| 545 | |
---|
| 546 | if (loutaver.eq.0) then |
---|
| 547 | write(*,*) ' #### FLEXPART MODEL ERROR! TIME AVERAGE OF #### ' |
---|
| 548 | write(*,*) ' #### CONCENTRATION FIELD OUTPUT MUST NOT BE #### ' |
---|
| 549 | write(*,*) ' #### ZERO. #### ' |
---|
| 550 | write(*,*) ' #### CHANGE INPUT IN FILE COMMAND. #### ' |
---|
| 551 | stop |
---|
| 552 | endif |
---|
| 553 | |
---|
| 554 | if (loutaver.gt.loutstep) then |
---|
| 555 | write(*,*) ' #### FLEXPART MODEL ERROR! TIME AVERAGE OF #### ' |
---|
| 556 | write(*,*) ' #### CONCENTRATION FIELD OUTPUT MUST NOT BE #### ' |
---|
| 557 | write(*,*) ' #### GREATER THAN INTERVAL OF OUTPUT. #### ' |
---|
| 558 | write(*,*) ' #### CHANGE INPUT IN FILE COMMAND. #### ' |
---|
| 559 | stop |
---|
| 560 | endif |
---|
| 561 | |
---|
| 562 | if (loutsample.gt.loutaver) then |
---|
| 563 | write(*,*) ' #### FLEXPART MODEL ERROR! SAMPLING TIME OF #### ' |
---|
| 564 | write(*,*) ' #### CONCENTRATION FIELD OUTPUT MUST NOT BE #### ' |
---|
| 565 | write(*,*) ' #### GREATER THAN TIME AVERAGE OF OUTPUT. #### ' |
---|
| 566 | write(*,*) ' #### CHANGE INPUT IN FILE COMMAND. #### ' |
---|
| 567 | stop |
---|
| 568 | endif |
---|
| 569 | |
---|
| 570 | if (mod(loutaver,lsynctime).ne.0) then |
---|
| 571 | write(*,*) ' #### FLEXPART MODEL ERROR! AVERAGING TIME OF #### ' |
---|
| 572 | write(*,*) ' #### CONCENTRATION FIELD MUST BE A MULTIPLE #### ' |
---|
| 573 | write(*,*) ' #### OF THE SYNCHRONISATION INTERVAL #### ' |
---|
| 574 | stop |
---|
| 575 | endif |
---|
| 576 | |
---|
| 577 | if ((loutaver/lsynctime).lt.2) then |
---|
| 578 | write(*,*) ' #### FLEXPART MODEL ERROR! AVERAGING TIME OF #### ' |
---|
| 579 | write(*,*) ' #### CONCENTRATION FIELD MUST BE AT LEAST #### ' |
---|
| 580 | write(*,*) ' #### TWICE THE SYNCHRONISATION INTERVAL #### ' |
---|
| 581 | stop |
---|
| 582 | endif |
---|
| 583 | |
---|
| 584 | if (mod(loutstep,lsynctime).ne.0) then |
---|
| 585 | write(*,*) ' #### FLEXPART MODEL ERROR! INTERVAL BETWEEN #### ' |
---|
| 586 | write(*,*) ' #### CONCENTRATION FIELDS MUST BE A MULTIPLE #### ' |
---|
| 587 | write(*,*) ' #### OF THE SYNCHRONISATION INTERVAL #### ' |
---|
| 588 | stop |
---|
| 589 | endif |
---|
| 590 | |
---|
| 591 | if ((loutstep/lsynctime).lt.2) then |
---|
| 592 | write(*,*) ' #### FLEXPART MODEL ERROR! INTERVAL BETWEEN #### ' |
---|
| 593 | write(*,*) ' #### CONCENTRATION FIELDS MUST BE AT LEAST #### ' |
---|
| 594 | write(*,*) ' #### TWICE THE SYNCHRONISATION INTERVAL #### ' |
---|
| 595 | stop |
---|
| 596 | endif |
---|
| 597 | |
---|
| 598 | if (mod(loutsample,lsynctime).ne.0) then |
---|
| 599 | write(*,*) ' #### FLEXPART MODEL ERROR! SAMPLING TIME OF #### ' |
---|
| 600 | write(*,*) ' #### CONCENTRATION FIELD MUST BE A MULTIPLE #### ' |
---|
| 601 | write(*,*) ' #### OF THE SYNCHRONISATION INTERVAL #### ' |
---|
| 602 | stop |
---|
| 603 | endif |
---|
| 604 | |
---|
| 605 | if (itsplit.lt.loutaver) then |
---|
| 606 | write(*,*) ' #### FLEXPART MODEL ERROR! SPLITTING TIME FOR#### ' |
---|
| 607 | write(*,*) ' #### PARTICLES IS TOO SHORT. PLEASE INCREASE #### ' |
---|
| 608 | write(*,*) ' #### SPLITTING TIME CONSTANT. #### ' |
---|
| 609 | stop |
---|
| 610 | endif |
---|
| 611 | |
---|
| 612 | if ((mquasilag.eq.1).and.(iout.ge.4)) then |
---|
| 613 | write(*,*) ' #### FLEXPART MODEL ERROR! CONFLICTING #### ' |
---|
| 614 | write(*,*) ' #### OPTIONS: IF MQUASILAG=1, PLUME #### ' |
---|
| 615 | write(*,*) ' #### TRAJECTORY OUTPUT IS IMPOSSIBLE. #### ' |
---|
| 616 | stop |
---|
| 617 | endif |
---|
| 618 | |
---|
| 619 | ! Compute modeling time in seconds and beginning date in Julian date |
---|
| 620 | !******************************************************************* |
---|
| 621 | |
---|
| 622 | outstep=real(abs(loutstep)) |
---|
| 623 | if (ldirect.eq.1) then |
---|
| 624 | bdate=juldate(ibdate,ibtime) |
---|
| 625 | edate=juldate(iedate,ietime) |
---|
| 626 | ideltas=nint((edate-bdate)*86400.) |
---|
| 627 | else if (ldirect.eq.-1) then |
---|
| 628 | loutaver=-1*loutaver |
---|
| 629 | loutstep=-1*loutstep |
---|
| 630 | loutsample=-1*loutsample |
---|
| 631 | lsynctime=-1*lsynctime |
---|
| 632 | bdate=juldate(iedate,ietime) |
---|
| 633 | edate=juldate(ibdate,ibtime) |
---|
| 634 | ideltas=nint((edate-bdate)*86400.) |
---|
| 635 | else |
---|
| 636 | write(*,*) ' #### FLEXPART MODEL ERROR! DIRECTION IN #### ' |
---|
| 637 | write(*,*) ' #### FILE "COMMAND" MUST BE EITHER -1 OR 1. #### ' |
---|
| 638 | stop |
---|
| 639 | endif |
---|
| 640 | |
---|
| 641 | return |
---|
| 642 | |
---|
| 643 | 999 write(*,*) ' #### FLEXPART MODEL ERROR! FILE "COMMAND" #### ' |
---|
| 644 | write(*,*) ' #### CANNOT BE OPENED IN THE DIRECTORY #### ' |
---|
| 645 | write(*,'(a)') path(1)(1:length(1)) |
---|
| 646 | stop |
---|
| 647 | |
---|
[f13406c] | 648 | 1000 write(*,*) ' #### FLEXPART MODEL ERROR! FILE "COMMAND" #### ' |
---|
[b4d29ce] | 649 | write(*,*) ' #### CANNOT BE OPENED IN THE DIRECTORY #### ' |
---|
[01e6052] | 650 | write(*,'(a)') path(2)(1:length(2)) |
---|
[b4d29ce] | 651 | stop |
---|
[e200b7a] | 652 | end subroutine readcommand |
---|