Changeset dba4221 in flexpart.git for src/readwind_ecmwf.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/readwind_ecmwf.f90 (modified) (2 diffs)

src/readwind_ecmwf.f90

@@ -103 +103 @@
 ! OPENING OF DATA FILE (GRIB CODE)
 !
+  write(*,*)   'Reading: '//path(3)(1:length(3)) &
+       //trim(wfname(indj))
 5 call grib_open_file(ifile,path(3)(1:length(3)) &
        //trim(wfname(indj)),'r',iret)
@@ -386 +388 @@
   endif
 
-  if(levdiff2.eq.0) then
-    iwmax=nlev_ec+1
+!  if(levdiff2.eq.0) then
+    iwmax=iwmax+1
     do i=0,nxmin1
       do j=0,nymin1
-        wwh(i,j,nlev_ec+1)=0.
+        wwh(i,j,iwmax)=0.
       end do
     end do
-  endif
+!  endif
 
 ! For global fields, assign the leftmost data column also to the rightmost