The W3FB04 routine converts the coordinates of a location on
earth from the natural coordinate system of
latitude/longitude to the grid (i,j) coordinate
system overlaid on a polar stereographic map
projection true at 60 degrees n or s latitude.
USAGE: CALL W3FB04 (ALAT, ALONG, XMESHL, ORIENT, XI, XJ)
ALAT LATITUDE IN DEGREES (<0 IF SH)
ALONG WEST LONGITUDE IN DEGREES
XMESHL MESH LENGTH OF GRID IN KM AT 60 DEG LAT(<0 IF SH)
(190.5 LFM GRID, 381.0 NH PE GRID,-381.0 SH PE GRID)
ORIENT ORIENTATION WEST LONGITUDE OF THE GRID
(105.0 LFM GRID, 80.0 NH PE GRID, 260.0 SH PE GRID)
XI I OF THE POINT RELATIVE TO NORTH OR SOUTH POLE
XJ J OF THE POINT RELATIVE TO NORTH OR SOUTH POLE
All parameters in the calling statement must be real.
The range of allowable latitudes is from a pole to 30
degrees into the opposite hemisphere. The grid used in
this subroutine has its origin (i=0,j=0) at the pole in
either hemisphere, so if the user's grid has its
origin at a point other than the pole, a translation is
needed to get i and j. The gridlines of i=constant are
parallel to a longitude designated by the user.
The earth's radius is taken to be 6371.2 km.
This code is not vectorized. You will have put
it line to vectorize it.
||Library contains Fortran 90 decoder/encoder
routines for GRIB edition 1. (Fortran90)
Date posted: 2/22/2007