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

src/readwind_nests.f90

@@ -79 +79 @@
   ! OPENING OF DATA FILE (GRIB CODE)
   !
-
+    write(*,*)   'Reading: '//path(numpath+2*(l-1)+1) &
+         (1:length(numpath+2*(l-1)+1))//trim(wfnamen(l,indj))
 5   call grib_open_file(ifile,path(numpath+2*(l-1)+1) &
          (1:length(numpath+2*(l-1)+1))//trim(wfnamen(l,indj)),'r')
@@ -365 +366 @@
   endif
 
-  if(levdiff2.eq.0) then
-    iwmax=nlev_ec+1
+!  if(levdiff2.eq.0) then
+    iwmax=iwmax+1
     do i=0,nxn(l)-1
       do j=0,nyn(l)-1
-        wwhn(i,j,nlev_ec+1,l)=0.
+        wwhn(i,j,iwmax,l)=0.
       end do
     end do
-  endif
+!  endif
 
   do i=0,nxn(l)-1