[ARPSSUPPORT] ext2arps using grib2 issue
Evans.Randy
evans.randy at ensco.com
Fri Sep 4 10:37:34 CDT 2009
Hi Kevin,
Thanks for the response. We are using half-degree GFS data and I have dumped out the g2print output for the first record (see below) and while I don't see anything that indicates it is grid #4, I am pretty sure that it is. We are running arps version 5.2.11.
Thanks for your help
Randy Evans
ENSCO, Inc.
ungrib - grib edition num 2
reading from grib file =
/ensco/archive/weatherdata/gfs-halfdegree/2008/2008091500.grb2
GRIB2 SECTION 0 - INDICATOR SECTION:
Discipline : 0
GRIB Edition Number : 2
GRIB length : 132582
GRIB2 SECTION 1 - IDENTIFICATION SECTION:
Length of Section : 13
Originating Center ID : 7
Subcenter ID : 0
GRIB Master Table Version : 2
GRIB Local Table Version : 1
Significance of Reference Time : 1
Year : 2008
Month : 9
Day : 15
Hour : 0
Minute : 0
Second : 0
Production Status of data : 0
Type of processed data : 1
GRIB2 SECTION 2 - LOCAL SECTION:
Length of Section 2 : 0
GRIB2 SECTION 3 - GRID DEFINITION SECTION:
Source of grid definition : 0
Number of grid points : 259920
Number of octets for addnl points : 0
Interpretation list : 0
Grid Definition Template Number : 0
Lat/Lon or Cylindrical Equidistant Grid
Shape of the Earth : 6
Scale factor of spher. Earth : 0
Scaled value of spher. Earth : 0
Scale factor of major axis : 0
Scaled value of major axis : 0
Scale factor of minor axis : 0
Scaled value of minor axis : 0
Ni - points along a parallel : 720
Nj - points along a meridian : 361
Basic angle of initial domain : 0
Subdivisions of basic angle : 0
La1 : 90000000
Lo1 : 0
Resolution and Component : 00110000
La2 : -90000000
Lo2 : 359500000
Di - i-dir increment : 500000
Dj - j-dir increment : 500000
Scanning mode : 0
GRIB2 SECTION 4 - PRODUCT DEFINITION SECTION:
Product Definition Template Number : 0
Template Entry 1 : 3
Template Entry 2 : 5
Template Entry 3 : 2
Template Entry 4 : 0
Template Entry 5 : 81
Template Entry 6 : 0
Template Entry 7 : 0
Template Entry 8 : 1
Template Entry 9 : 0
Template Entry 10 : 100
Template Entry 11 : 0
Template Entry 12 : 1000
Template Entry 13 : 255
Template Entry 14 : 0
Template Entry 15 : 0
Product Abbreviated Name : HGT
GRIB2 SECTION 5 - DATA REPRESENTATION SECTION:
Data Representation Template Number : 40
Template Entry 1 : 1243990156
Template Entry 2 : 3
Template Entry 3 : 2
Template Entry 4 : 16
Template Entry 5 : 0
Template Entry 6 : 0
Template Entry 7 : 255
GRIB2 SECTION 6 - BIT-MAP SECTION:
Bit-map indicator : 255
GRIB2 SECTION 7 - DATA SECTION:
Minimum Data Value : 27161.94922
Maximum Data Value : 31355.06836
-----Original Message-----
From: Kevin W. Thomas [mailto:kwthomas at wizard.caps.ou.edu]
Sent: Friday, September 04, 2009 10:23 AM
To: arpssupport at caps.ou.edu; Evans.Randy
Subject: Re: [ARPSSUPPORT] ext2arps using grib2 issue
>We are trying to run ext2arps using a gfs grib2 file as input and we are getting the following core dump. When we read the grib2 file using wgrib2 the file appears to be a good file. When we built ext2arps we built with all the required grib2 libraries (libjasper, libpng, libz). Any suggestions?
>
>&extdfile
> dir_extd = '/ensco/home/randy/',
> extdname = 'gfs',
> extdopt = 20,
> extdfmt = 1,
> nextdfil = 1,
> extdtime(01) = '2008-09-15.00:00:00+000:00:00',
> iorder = 2,
> intropt = 1,
> nsmooth = 1,
> hydradj = 0,
> wndadj = 2,
> obropt = 12,
> obrzero = 12000.0000,
> ext_lbc = 1,
> ext_vbc = 1,
> exttrnopt = 0,
> extntmrg = 7,
> extsfcopt = 0,
> grdbasopt = 0,
> i2dfmt = 0,
> outheader = 'data2d/nam.2008030300',
> &END
>
> ARPS domain NW: ( 67.43,-142.34) NE: ( 45.14, -24.68)
> SW: ( 12.65,-122.74) SE: ( 2.38, -65.17).
>
>
> The base state is formed from a mean sounding created
> with a 50.00 meter interval.
>
> nx_ext = 720, ny_ext = 361, nz_ext = 26, nzsoil_ext = 2
>
> Calling getncepgfs_grb2 looking for 000:00:00 hour forecast
> initialized at 2008-09-15.00:00:00 UTC
> Julian filename: 082590000
> Which is 1537056000 abs seconds or
> 0 seconds from the ARPS initial time.
> rdnmcgrib2: Opening GRIB /ensco/home/randy/gfs.2008091500f00
> Fortran I/O unit 39 picked from the free list.
>Segmentation fault (core dumped)
Hi...
There are two types of GFS data. One is GFS grid 3, with 1 deg x 1 deg
resolution. The other is GFS grid 4, with 0.5 deg x 0.5 deg resolution.
You specified "extdopt" of 20, so you should be using the grid 4 data.
Please confirm that is what you are using. If you are using the correct
data, please let us know what version of ARPS you are using, so we can look
into the problem.
Although we don't use GFS in our daily forecasts, we bring in both the grid 3
and grid 4 data.
Kevin W. Thomas
Center for Analysis and Prediction of Storms
University of Oklahoma
Norman, Oklahoma
Email: kwthomas at ou.edu
The information contained in this email message is intended only for the use of the individual(s) to whom it is addressed and may contain information that is privileged and sensitive. If you are not the intended recipient, or otherwise have received this communication in error, please notify the sender immediately by email at the above referenced address and note that any further dissemination, distribution or copying of this communication is strictly prohibited.
The U.S. Export Control Laws regulate the export and re-export of technology originating in the United States. This includes the electronic transmission of information and software to foreign countries and to certain foreign nationals. Recipient agrees to abide by these laws and their regulations -- including the U.S. Department of Commerce Export Administration Regulations and the U.S. Department of State International Traffic in Arms Regulations -- and not to transfer, by electronic transmission or otherwise, any content derived from this email to either a foreign national or a foreign destination in violation of such laws.
More information about the ARPSSUPPORT
mailing list