Convert the Cartesian grid to polar grid with the origin of the storm center (from STPK) The procedure: (1) Define the given polar grid (r-theta) on the Cartesian grid (2) Search the 4 nearest points on the Cartesian grid (x-y) for each polar grid point (3) Performing the bilinear interpolation of the 4 values defined on the Cartesian grid to the polar grid.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| double precision, | intent(in) | :: | x(:) |
X-coordinate on the Cartesian grid |
||
| double precision, | intent(in) | :: | y(:) |
Y-coordinate on the Cartesian grid |
||
| double precision, | intent(in) | :: | xc |
X-component of the storm center |
||
| double precision, | intent(in) | :: | yc |
Y-component of the storm center |
||
| double precision, | intent(in) | :: | u(size(x),size(y)) |
Values defined on the Cartesian grid |
||
| double precision, | intent(in) | :: | r(:) |
R-coordinate on the polar grid with the origin (xc, yc). |
||
| double precision, | intent(in) | :: | theta(:) |
theta-coordinate of the polar grid with the origin (xc, yc) [rad] |
||
| double precision, | intent(out) | :: | v(size(r),size(theta)) |
Values converted to the polar grid |
||
| double precision, | intent(in), | optional | :: | undef |
Missing value for the outside of the polar grid area (default: -999.0) |
|
| double precision, | intent(in), | optional | :: | undefg |
Missing value for the inside of the polar grid area (default: -999.0) |
|
| logical, | intent(in), | optional | :: | stdopt |
Display debug messages. (default: .false. == No display) |