Changeset 94b3788 in flexpart.git for src/readwind_gfs.f90

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

File:

• src/readwind_gfs.f90 (modified) (2 diffs)

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)