Description: Fixes for gfortran-10 transition
 Promote some arrays from REAL(4) to REAL(8) to agree with interfaces
Author: Alastair McKinstry <mckinstry@debian.org>
Last-Updated: 2020-04-23
Forwarded: no

Index: flexpart-9.02/readwind_gfs.f90
===================================================================
--- flexpart-9.02.orig/readwind_gfs.f90
+++ flexpart-9.02/readwind_gfs.f90
@@ -95,7 +95,7 @@ subroutine readwind(indj,n,uuh,vvh,wwh)
   real(kind=4) :: xaux,yaux,xaux0,yaux0
   real(kind=8) :: xauxin,yauxin
   real,parameter :: eps=1.e-4
-  real(kind=4) :: ewss(0:nxmax-1,0:nymax-1),nsss(0:nxmax-1,0:nymax-1)
+  real(kind=8) :: ewss(0:nxmax-1,0:nymax-1),nsss(0:nxmax-1,0:nymax-1)
   real :: plev1,hlev1,ff10m,fflev1
 
   logical :: hflswitch,strswitch
Index: flexpart-9.02/readwind.f90
===================================================================
--- flexpart-9.02.orig/readwind.f90
+++ flexpart-9.02/readwind.f90
@@ -74,9 +74,9 @@ subroutine readwind(indj,n,uuh,vvh,wwh)
   integer :: gotGrid
   !HSO  end
 
-  real(kind=4) :: uuh(0:nxmax-1,0:nymax-1,nuvzmax)
-  real(kind=4) :: vvh(0:nxmax-1,0:nymax-1,nuvzmax)
-  real(kind=4) :: wwh(0:nxmax-1,0:nymax-1,nwzmax)
+  real(kind=8) :: uuh(0:nxmax-1,0:nymax-1,nuvzmax)
+  real(kind=8) :: vvh(0:nxmax-1,0:nymax-1,nuvzmax)
+  real(kind=8) :: wwh(0:nxmax-1,0:nymax-1,nwzmax)
   integer :: indj,i,j,k,n,levdiff2,ifield,iumax,iwmax
 
   ! VARIABLES AND ARRAYS NEEDED FOR GRIB DECODING
@@ -92,7 +92,7 @@ subroutine readwind(indj,n,uuh,vvh,wwh)
   real(kind=4) :: xaux,yaux,xaux0,yaux0
   real(kind=8) :: xauxin,yauxin
   real,parameter :: eps=1.e-4
-  real(kind=4) :: nsss(0:nxmax-1,0:nymax-1),ewss(0:nxmax-1,0:nymax-1)
+  real(kind=8) :: nsss(0:nxmax-1,0:nymax-1),ewss(0:nxmax-1,0:nymax-1)
   real :: plev1,pmean,tv,fu,hlev1,ff10m,fflev1
 
   logical :: hflswitch,strswitch
