id summary reporter owner description type status priority milestone component version resolution keywords cc 115 Bug in readspecies.f90 pesei pesei "I found a bug in `readspecies.f90`, probably introduced with the modulation of releases as a function of day of week and hour of day, and more with intro of namelists. '''1. Wrong address for premature end of SPECIES''' In the lines {{{ #!python read(unitspecies,'(a10)',end=22) species(id_pos) ! write(*,*) species(id_pos) }}} etc., the place to jump to must be '''`20`''', otherwise `area_hour, point_hour, area_dow, point_dow` are not intialised. '''2. Attempt to read modulation factors''' If we are already at EOF, we can't use `end=` anymore, so use `err=22`. '''3. Assignment to namelist''' (maybe that should go into #94) {{{ pspecies=species(id_pos) pdecay=decay(id_pos) pweta=weta(id_pos) pwetb=wetb(id_pos) pweta_in=weta_in(id_pos) pwetb_in=wetb_in(id_pos) pwetc_in=wetc_in(id_pos) pwetd_in=wetd_in(id_pos) preldiff=reldiff(id_pos) phenry=henry(id_pos) pf0=f0(id_pos) pdensity=density(id_pos) pdquer=dquer(id_pos) pdsigma=dsigma(id_pos) pdryvel=dryvel(id_pos) pweightmolar=weightmolar(id_pos) pohreact=ohreact(id_pos) pspec_ass=spec_ass(id_pos) pkao=kao(id_pos) }}} has to be split and each line to be put below the corresponding read statement. There is more stuff that should be cleaned up even if it does not give bugs. I started to do so in my wet depo work. " Defect accepted major FLEXPART 9.2 FP coding/compilation FLEXPART 9.0.2