Changeset 8a65cb0 in flexpart.git for src/readwind_gfs.f90

Timestamp:

Mar 2, 2015, 3:11:55 PM (9 years ago)

Author:

Espen Sollum ATMOS <espen@…>

Branches:

master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie

Children:

1d207bb

Parents:

60403cd

Message:

Added code, makefile for dev branch

File:

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

src/readwind_gfs.f90

@@ -143 +143 @@
   !first see if we read GRIB1 or GRIB2
   call grib_get_int(igrib,'editionNumber',gribVer,iret)
-  call grib_check(iret,gribFunction,gribErrorMsg)
+!  call grib_check(iret,gribFunction,gribErrorMsg)
 
   if (gribVer.eq.1) then ! GRIB Edition 1
@@ -149 +149 @@
   !read the grib1 identifiers
   call grib_get_int(igrib,'indicatorOfParameter',isec1(6),iret)
-  call grib_check(iret,gribFunction,gribErrorMsg)
+!  call grib_check(iret,gribFunction,gribErrorMsg)
   call grib_get_int(igrib,'indicatorOfTypeOfLevel',isec1(7),iret)
-  call grib_check(iret,gribFunction,gribErrorMsg)
+!  call grib_check(iret,gribFunction,gribErrorMsg)
   call grib_get_int(igrib,'level',isec1(8),iret)
-  call grib_check(iret,gribFunction,gribErrorMsg)
+!  call grib_check(iret,gribFunction,gribErrorMsg)
 
   else ! GRIB Edition 2
@@ -159 +159 @@
   !read the grib2 identifiers
   call grib_get_int(igrib,'discipline',discipl,iret)
-  call grib_check(iret,gribFunction,gribErrorMsg)
+!  call grib_check(iret,gribFunction,gribErrorMsg)
   call grib_get_int(igrib,'parameterCategory',parCat,iret)
-  call grib_check(iret,gribFunction,gribErrorMsg)
+!  call grib_check(iret,gribFunction,gribErrorMsg)
   call grib_get_int(igrib,'parameterNumber',parNum,iret)
-  call grib_check(iret,gribFunction,gribErrorMsg)
+!  call grib_check(iret,gribFunction,gribErrorMsg)
   call grib_get_int(igrib,'typeOfFirstFixedSurface',typSurf,iret)
-  call grib_check(iret,gribFunction,gribErrorMsg)
+!  call grib_check(iret,gribFunction,gribErrorMsg)
   call grib_get_int(igrib,'scaledValueOfFirstFixedSurface', &
        valSurf,iret)
-  call grib_check(iret,gribFunction,gribErrorMsg)
+!  call grib_check(iret,gribFunction,gribErrorMsg)
 
   !convert to grib1 identifiers
@@ -262 +262 @@
   !  get the size and data of the values array
     call grib_get_real4_array(igrib,'values',zsec4,iret)
-    call grib_check(iret,gribFunction,gribErrorMsg)
+!    call grib_check(iret,gribFunction,gribErrorMsg)
   endif
 
@@ -271 +271 @@
   call grib_get_int(igrib,'numberOfPointsAlongAParallel', &
        isec2(2),iret)
-  call grib_check(iret,gribFunction,gribErrorMsg)
+!  call grib_check(iret,gribFunction,gribErrorMsg)
   call grib_get_int(igrib,'numberOfPointsAlongAMeridian', &
        isec2(3),iret)
-  call grib_check(iret,gribFunction,gribErrorMsg)
+!  call grib_check(iret,gribFunction,gribErrorMsg)
   call grib_get_real8(igrib,'longitudeOfFirstGridPointInDegrees', &
        xauxin,iret)
-  call grib_check(iret,gribFunction,gribErrorMsg)
+!  call grib_check(iret,gribFunction,gribErrorMsg)
   call grib_get_real8(igrib,'latitudeOfLastGridPointInDegrees', &
        yauxin,iret)
-  call grib_check(iret,gribFunction,gribErrorMsg)
+!  call grib_check(iret,gribFunction,gribErrorMsg)
   xaux=xauxin+real(nxshift)*dx
   yaux=yauxin