Add capacity of detecting scaning mode north2south vs south2north in …

…gaussian grid from Dusan. (#52)
NOAA-EMC · Jan 17, 2020 · 927708e · 927708e
1 parent 5f081ef
commit 927708e
Showing 1 changed file with 5 additions and 0 deletions.
diff --git a/sorc/ncep_post.fd/grib2_module.f b/sorc/ncep_post.fd/grib2_module.f
@@ -1394,6 +1394,11 @@ subroutine getgds(ldfgrd,len3,ifield3len,igds,ifield3)
        ifield3(16) = lonlast 
        ifield3(17) = NINT(360./(IM)*1000000.)
        ifield3(18) = NINT(JM/2.0)
+       if( latstart < latlast ) then
+        ifield3(19) = 64      !for SN scan
+       else
+        ifield3(19) = 0       !for NS scan
+       endif
 !
 !** Latlon grid
       ELSE IF(MAPTYPE == 0 ) THEN