Changeset dba4221 in flexpart.git for src/readspecies.f90

Timestamp:

Nov 9, 2022, 4:50:39 PM (18 months ago)

Author:

Pirmin Kaufmann <pirmin.kaufmann@…>

Branches:

bugfixes+enhancements

Children:

8ad70c7

Parents:

49e63b9

Message:

Bugfixes:

• options/SPECIES/SPECIES_009: corrected wrong number format, replaced comma by decimal point
• options/SPECIES/SPECIES_028: corrected wrong number format, moved sign of exponent to after the E
• options/SPECIES/specoverview.f90: added namelist parameters that appear in SPECIES files but were missing here
• src/FLEXPART.f90: replaced compiler-specific command line argument routines by standard Fortran intrinsic routines
• src/FLEXPART_MPI.f90: ditto
• src/gridcheck_ecmwf.f90: corrected handling of vertical levels when input files do not contain uppermost layers
• src/gridcheck_nests.f90: ditto
• src/readwind_ecmwf.f90: corrected handling of vertical levels when input files do not contain uppermost layers
• readwind_ecmwf_mpi.f90: ditto

Code enhancements:

• options/OUTGRID: added comments describing contents
• options/SPECIES/SPECIES_*: aligned comments
• options/SPECIES/specoverview.f90: removed commented lines, rectified lines indenting
• src/FLEXPART.f90: rectified lines indenting, updated date in version string
• src/FLEXPART_MPI.f90: ditto, and realigned code with src/FLEXPART.f90
• src/gridcheck_*.f90: added code to write out name of file before it is opened (helps a lot when an input file causes troubles)
• src/par_mod.f90: added comment explaining relevance of nuvzmax for GRIB input
• src/readreleases.f90: write out warning if too few particles are used to randomize release
• src/readspecies.f90: write out name of SPECIES file before it is read
• src/readwind_*.f90: write out name of input file before opening it
• src/writeheader_txt.f90: removed wrong comment

File:

• src/readspecies.f90 (modified) (1 diff)

src/readspecies.f90

@@ -104 +104 @@
   specnum(pos_spec)=id_spec
   write(aspecnumb,'(i3.3)') specnum(pos_spec)
+  write(*,*) 'Reading: '// &
+                        path(1)(1:length(1))//'SPECIES/SPECIES_'//aspecnumb
   open(unitspecies,file=path(1)(1:length(1))//'SPECIES/SPECIES_'//aspecnumb,status='old',form='formatted',err=998)
 !write(*,*) 'reading SPECIES',specnum(pos_spec)