Changeset 94b3788 in flexpart.git

Timestamp:

Apr 8, 2021, 11:35:25 AM (3 years ago)

Author:

Espen Sollum <eso@…>

Branches:

dev

Children:

e2132cf, 8f3cc41

Parents:

b1f28c3

Message:

Fix for reading GFS v16 fields

Location:

src

Files:

• gridcheck_gfs.f90 (modified) (4 diffs)
• readwind_gfs.f90 (modified) (2 diffs)

src/gridcheck_gfs.f90

@@ -73 +73 @@
   !HSO  end
   integer :: ix,jy,i,ifn,ifield,j,k,iumax,iwmax,numskip
-  real :: sizesouth,sizenorth,xauxa,pint
+  real :: sizesouth,sizenorth,xauxa,pint,xsec18
   real :: akm_usort(nwzmax)
   real,parameter :: eps=spacing(2.0_4*360.0_4) 
@@ -149 +149 @@
   call grib_check(iret,gribFunction,gribErrorMsg)
 
+  xsec18 = real(isec1(8))
+
   else ! GRIB Edition 2
 
@@ -168 +170 @@
   isec1(7)=-1
   isec1(8)=-1
+  xsec18 = -1.0
   if ((parCat.eq.2).and.(parNum.eq.2).and.(typSurf.eq.100)) then ! U
     isec1(6)=33          ! indicatorOfParameter
     isec1(7)=100         ! indicatorOfTypeOfLevel
-    isec1(8)=valSurf/100 ! level, convert to hPa
+    xsec18=valSurf/100.0 ! level, convert to hPa
   elseif ((parCat.eq.3).and.(parNum.eq.5).and.(typSurf.eq.1)) then ! TOPO
     isec1(6)=7           ! indicatorOfParameter
     isec1(7)=1           ! indicatorOfTypeOfLevel
-    isec1(8)=0
+    xsec18=real(0)
   elseif ((parCat.eq.0).and.(parNum.eq.0).and.(typSurf.eq.1) &
        .and.(discipl.eq.2)) then ! LSM
     isec1(6)=81          ! indicatorOfParameter
     isec1(7)=1           ! indicatorOfTypeOfLevel
-    isec1(8)=0
+    xsec18=real(0)
   endif
 
@@ -301 +304 @@
   if((isec1(6).eq.33).and.(isec1(7).eq.100)) then ! check for U wind
     iumax=iumax+1
-    pres(iumax)=real(isec1(8))*100.0
+    pres(iumax)=xsec18*100.0
   endif
 

src/readwind_gfs.f90

@@ -217 +217 @@
       isec1(7)=102         ! indicatorOfTypeOfLevel
       xsec18=real(0)
-    elseif ((parCat.eq.3).and.(parNum.eq.0).and.(typSurf.eq.1)) then ! SP
+   elseif ((parCat.eq.3).and.(parNum.eq.0).and.(typSurf.eq.1).and. &
+        (discipl.eq.0)) then ! SP
       isec1(6)=1           ! indicatorOfParameter
       isec1(7)=1           ! indicatorOfTypeOfLevel
@@ -602 +603 @@
         if (temp .eq. 0.0) then 
           write (*, *) i, j, k, n
-          temp = 273.0
+!          temp = 273.0
+          stop
         endif
         plev1=akm(k)+bkm(k)*ps(i,j,1,n)