W3FI75
The W3FI75 routine packs a grib field and forms octets (1-11)
of the binary data section (BDS).
USAGE: CALL W3FI75 (IBITL,ITYPE,ITOSS,FLD,IFLD,IBMAP,IBDSFL,
NPTS,BDS11,IPFLD,PFLD,LEN,LENBDS,IBERR,PDS,IGDS)
Input argument list:
IBITL - 0, COMPUTER COMPUTES PACKING LENGTH FROM POWER
OF 2 THAT BEST FITS THE DATA.
8, 12, ETC. COMPUTER RESCALES DATA TO FIT INTO
SET NUMBER OF BITS.
ITYPE - 0 = IF INPUT DATA IS FLOATING POINT (FLD)
1 = IF INPUT DATA IS INTEGER (IFLD)
ITOSS - 0 = NO BIT MAP IS INCLUDED (DON'T TOSS DATA)
1 = TOSS NULL DATA ACCORDING TO IBMAP
FLD - REAL ARRAY OF DATA TO BE PACKED IF ITYPE=0
IFLD - INTEGER ARRAY TO BE PACKED IF ITYPE=1
IBMAP - BIT MAP SUPPLIED FROM USER
IBDSFL - INTEGER ARRAY CONTAINING TABLE 11 FLAG INFO
BDS OCTET 4:
(1) 0 = GRID POINT DATA
1 = SPHERICAL HARMONIC COEFFICIENTS
(2) 0 = SIMPLE PACKING
1 = SECOND ORDER PACKING
(3) 0 = ORIGINAL DATA WERE FLOATING POINT VALUES
1 = ORIGINAL DATA WERE INTEGER VALUES
(4) 0 = NO ADDITIONAL FLAGS AT OCTET 14
1 = OCTET 14 CONTAINS FLAG BITS 5-12
(5) 0 = RESERVED - ALWAYS SET TO 0
(6) 0 = SINGLE DATUM AT EACH GRID POINT
1 = MATRIX OF VALUES AT EACH GRID POINT
(7) 0 = NO SECONDARY BIT MAPS
1 = SECONDARY BIT MAPS PRESENT
(8) 0 = SECOND ORDER VALUES HAVE CONSTANT WIDTH
1 = SECOND ORDER VALUES HAVE DIFFERENT WIDTHS
NPTS - NUMBER OF GRIDPOINTS IN ARRAY TO BE PACKED
IGDS - ARRAY OF GDS INFORMATION
Output argument list:
BDS11 - FIRST 11 OCTETS OF BDS
PFLD - PACKED GRIB FIELD
LEN - LENGTH OF PFLD
LENBDS - LENGTH OF BDS
IBERR - 1, ERROR CONVERTING IEEE F.P. NUMBER TO IBM370 F.P.
REMARKS:
Subprogram can be called from a multiprocessing environment.
W3lib.tar |
Library contains Fortran 90 decoder/encoder
routines for GRIB edition 1. (Fortran90)
Date posted: 2/22/2007 |