W3FT05
The W3FT05 routine converts a northern hemisphere 2.5
degree lat.,lon. 145 by 37 grid to a polar
stereographic 65 by 65 grid. The polar
stereographic map projection is true at 60 deg. North,
the mesh length is 381 km. and the oriention is
80 deg. west.
Input argument list:
ALOLA 145*37 GRID 2.5 LAT,LON GRID N. HEMI.
5365 POINT GRID IS TYPE 29 OR 1D HEX O.N. 84
LINEAR 1 LINEAR INTERPOLATION , NE.1 BIQUADRATIC
Output argument list:
APOLA 65*65 GRID OF NORTHERN HEMI.
4225 POINT GRID IS TYPE 27 OR 1B HEX O.N. 84
W1 65*65 SCRATCH FIELD
W2 65*65 SCRATCH FIELD
REMARKS:
1. W1 AND W2 ARE USED TO STORE SETS OF CONSTANTS WHICH ARE
REUSABLE FOR REPEATED CALLS TO THE SUBROUTINE. IF THEY ARE
OVER WRITTEN BY THE USER, A WARNING MESSAGE WILL BE PRINTED
AND W1 AND W2 WILL BE RECOMPUTED.
2. WIND COMPONENTS ARE NOT ROTATED TO THE 65*65 GRID ORIENTATION
AFTER INTERPOLATION. YOU MAY USE W3FC08 TO DO THIS.
3. THE GRID POINTS VALUES ON THE EQUATOR HAVE BEEN EXTRAPOLATED
OUTWARD TO ALL THE GRID POINTS OUTSIDE THE EQUATOR ON THE 65*65
GRID (ABOUT 1100 POINTS).
4. YOU SHOULD USE THE CRAY VECTORIZED VERSION W3FT05V ON THE CRAY
IT HAS 3 PARAMETERS IN THE CALL, RUNS ABOUT 10 TIMES FASTER. USES
MORE MEMORY.
W3lib.tar |
Library contains Fortran 90 decoder/encoder
routines for GRIB edition 1. (Fortran90)
Date posted: 2/22/2007 |