Convert 25 word array to GRIB pds. More...

Functions/Subroutines
subroutine	w3fi68 (id, pds)
	Converts an array of 25, or 27 integer words into a grib product definition section (pds) of 28 bytes , or 30 bytes.

Detailed Description

Convert 25 word array to GRIB pds.

Definition in file w3fi68.f.

Function/Subroutine Documentation

Converts an array of 25, or 27 integer words into a grib product definition section (pds) of 28 bytes , or 30 bytes.

if pds bytes > 30, they are set to zero.

Program history log:

Ralph Jones 1991-05-08
Ralph Jones 1992-09-25 Change to 25 words of input, level can be in two words. (10,11)
Ralph Jones 1993-01-08 Change for time range indicator if 10, store time p1 in pds bytes 19-20.
Ralph Jones 1993-01-26 Correction for fixed height above ground level
Ralph Jones 1993-03-29 Add save statement
Bill Cavanaugh 1993-06-24 Modified program to allow for generation of pds greater than 28 bytes (the desired pds size is in id(1).
Farley 1993-09-30 Change to allow for subcenter id; put id(24) into pds(26).
Ralph Jones 1993-10-12 Changes for on388 rev. oct 9,1993, new levels 125, 200, 201.
Ralph Jones 1994-02-23 Take out sbytes, replace with do loop
Ralph Jones 1994-04-14 Changes for on388 rev. mar 24,1994, new levels 115,116.
Ralph Jones 1994-12-04 Change to add id words 26, 27 for pds bytes 29 and 30.
Ralph Jones 1995-09-07 Change for new level 117, 119.
Mark Iredell 1995-10-31 REmoved saves and prints
Ebisuzaki 1998-06-30 Linux port
Stephen Gilbert 2001-06-05 Changed fortran intrinsic function OR() to f90 standard intrinsic IOR().
Mark Iredell 2003-02-25 Recognize level type 126
D. C. Stokes 2005-05-06 Recognize level types 235, 237, 238

Parameters

[in]	ID	25,27 word integer array.
[out]	PDS	28 30 or greater character pds for edition 1.

Note

Layout of 'id' array:

ID(1) = Number of bytes in product definition section (pds)
ID(2) = Parameter table version number
ID(3) = Identification of originating center
ID(4) = Model identification (allocated by originating center)
ID(5) = Grid identification
ID(6) = 0 if no gds section, 1 if gds section is included
ID(7) = 0 if no bms section, 1 if bms section is included
ID(8) = Indicator of parameter and units (table 2)
ID(9) = Indicator of type of level (table 3)
ID(10) = Value 1 of level (0 for 1-100,102,103,105,107 109,111,113,115,117,119,125,126,160,200,201,235,237,238 level is in id word 11)
ID(11) = Value 2 of level
ID(12) = Year of century
ID(13) = Month of year
ID(14) = Day of month
ID(15) = Hour of day
ID(16) = Minute of hour (in most cases set to 0)
ID(17) = Fcst time unit
ID(18) = P1 period of time
ID(19) = P2 period of time
ID(20) = Time range indicator
ID(21) = Number included in average
ID(22) = Number missing from averages
ID(23) = Century (20, change to 21 on jan. 1, 2001)
ID(24) = Subcenter identification
ID(25) = Scaling power of 10

ID(26) = Flag byte, 8 on/off flags

BIT NUMBER	VALUE	ID(26)	DEFINITION
1	0	0	FULL FCST FIELD
	1	128	FCST ERROR FIELD
2	0	0	ORIGINAL FCST FIELD
	1	64	BIAS CORRECTED FCST FIELD
3	0	0	ORIGINAL RESOLUTION RETAINED
	1	32	SMOOTHED FIELD

ID(26) can be the sum of bits 1, 2, 3. bits 4-8 not used, set to zero if ID(1) is 28, you do not need ID(26) and ID(27).

Definition at line 84 of file w3fi68.f.