W3FB00

The W3FB00 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 the polar
           stereographic map projection true at 60 N.

USAGE:  CALL W3FB00 (ALAT, ALONG, XMESHL, XI, XJ)

   Input variables:
      ALAT     LATITUDE IN DEG.  (-20.0(S. HEMIS) ) ALAT ) 90.0)
      ALONG    WEST LONGITUDE IN DEGREES
      XMESHL   MESH LENGTH OF GRID IN KILOMETERS AT 60N

   Output variables:
      XI       I OF THE POINT RELATIVE TO NORTH POLE
      XJ       J OF THE POINT RELATIVE TO NORTH POLE

REMARKS:
        The grid used in this subroutine has its origin (i=0,j=0)
        at the north pole, so if the user's grid has its origin at a
        point other than the north pole, a translation is required to
        get i and j. The subroutine grid is oriented so that longitude
        80W is parallel to the gridlines of i=constant. the radius of
        the earth is taken to be 6371.2 km.
        All parameters in the call statement must be real. 
        This code will not vectorize on a cray. You will have put
        it line to vectorize it.

 


  
W3lib.tar
  
Library contains Fortran 90 decoder/encoder
   routines for GRIB edition 1. (Fortran90) 

   Date posted: 2/22/2007