Procedures

add ival

add ival

continuous running of interpo_search_1d

continuous running of auto_interpo_search_1d

Automatic interpolation for 2d data (continuous running of interpolation_2d)

convert char to float

Calculate the nonlinear angle psid at a centain radius
From the geometry,
$\sin{\psi_d} = (R_{TC}/r) \sin{\theta _d}.$ --(1)
$R^2_{TC} = r^2 + D^2 - 2rD\cos{\psi_d} , D=r\dfrac{\sin{\theta}}{\sin{\theta _d}}.$ --(2)
From Eq. (2),
$\cos{\psi_d} = (r^2 + D^2 - R_{TC}^2) / (2rD).$ --(3)
From Eqs. (1) and (3),
$\tan{\psi_d} = 2\rho \dfrac{\sin{\theta}}{\rho ^2+a\rho ^2-1},$
$\rho =r/R_{TC}, a=\dfrac{\sin{\theta}}{\sin{\theta _d}} .$

Convert polar grid with the origin of the storm center to the Cartesian grid (from STPK) The procedure: (1) Define the given Cartesian grid (x-y) on the polar grid (2) Search the 4 nearest points on the polar grid for each Cartesian grid point (3) Performing the bilinear interpolation of the 4 values defined on the polar grid to the Cartesian grid.

Check the innermost radius of data with no undef If the innermost radius has not been found, the value of zero is returned.

Check maximum value in each element of the matrix "a"

check missing grids

Convert double to real

Convert double to real

convert double to real

convert real to double

Convert Vt and Vr to Vx and Vy on R-T coordinates

Convert Vx and Vy to Vr and Vt on R-T coordinates

Convert Vx and Vy to Vr and Vt on X-Y coordinates

decomposition of VD to the components parallel and normal to the direction of thetad = 0

Display the maximum of the array val

Display the maximum of the difference between val1 and val2

Calculation of rotation and divergence from radial and tangential winds
divr = $\dfrac{\partial ru_r}{r\partial r} + \dfrac{\partial v_t}{r\partial \theta}$
curl = $\dfrac{\partial rv_t}{r\partial r} - \dfrac{\partial u_r}{r\partial \theta}$

Calculate inner product of two vectors

Gauss-Jordan method with fully pivotting (from STPK)

Calculate the inverse "xx" for the matrix "ax" (from STPK)

Calculate gradient of "u" in one dimension based on the 2nd order central differential approximation:
$\frac{\partial u}{\partial x}\approx \frac{u_{i+1}-u_{i-1}}{2dx}$

Calculate gradient of "u" in two dimensions based on the 2nd order central differential approximation:
$\frac{\partial u}{\partial x}\approx \frac{u_{i+1}-u_{i-1}}{2dx}$ $\nabla u\approx \left(\frac{p_{i+1,j}-p_{i-1,j}}{2dx} ,\; \frac{p_{i,j+1}-p_{i,j-1}}{2dy} \right)$

Calculation of the Green function: $,\;$
$G_n(r_c;r)=-(2^{-n})(r/r_c)^n, r<r_c$
$G_n(r_c;r)=-(2^{-n})(r_c/r)^n, r\geq r_c$
$n=$nval

Check positive value of the innermost radius. If the innermost radius is negative, the element number located at the innermost positive radius is returned.

Floor function for the real grid points (from STPK)

Floor function for the real grid points (from STPK)

Perform linear interpolation to the "point" on the Cartesian (x) grid (from STPK)

Perform bilinear interpolation to the "point" on the Cartesian (x-y) grid (from STPK)

Calculate a line integral (actually, sum for arguments)

count the line in the fname

calculate arc distance between two points on the sphere

calculate radial and azimuthal location at (lon1,lat1) on the polar coordinate with the center of (lon0,lat0)

Calculate product for a component in a matrix

Calculate product for a component in a matrix

Calculate product for a component in a matrix

Get the max value (from STPK)

Average x

Producing vortex structure with rmax, vmax, and umax

Producing vortex structure with rmax, vmax, and umax in Lee et al. (2006)

projection of a uniform wind [(ux,vy) on Cartesian grids] to the radial and tangential components on the polar coordinate

projection of a uniform wind [(ux,vy) on Cartesian grids] to the line-of-sight component along the radar beam

Convert Vt and Vr to Vx and Vy on R-T coordinates

Calculate Vx and Vy to Vd along with radar beams on X-Y coodinates

read float data from 4-byte unformatted binary

read float data from 4-byte unformatted binary

read ASCII file

Rearrange 2d variable to 1d variable (i,j -> i*j)

Rearrange 3d variable to 2d variable (k,i,j -> i*j,k)

rearrange the radial grids: 1. skip radii with few azimuthal sampling (less than "skip_thres") 2. rearrange the radial coordinate. [Note1]: r, rh, thetad, and Vra must be given by the same as the array number in the original. (that is, size(r_org) == size(r).) [Note2]: nr_out_org >= nr_out. Undef value is given in r, rh, thetad, and Vra outside nr_out.

recover the retrieved values from the rearranged to original radii.

replace ioval2d with undefv at the point undeflag == .true.

replace orgval with repval in ival2d array

Solve unknown variables and return wind velocity on R-T coordinates based on the GBVTD technique.
------------------------------------------------------
[relationship between r and rh] --
------------------------------------------------------
i-1 i i+1
...|-- --|-- --|... : r(1:size(r)) = velocity radii

Solve unknown variables and return wind velocity on R-T coordinates based on the GVTD technique.
------------------------------------------------------
[relationship between r and rh] --
------------------------------------------------------
i-1 i i+1
...|-- --|-- --|... : r(1:size(r)) = velocity radii

Solve unknown variables and return wind velocity on R-T coordinates.
-------------------------------------------------------
-- [relationship between r and rh] --
-------------------------------------------------------
-- i-1 i i+1
-- ...|-- --|-- --|... : r(1:size(r)) = velocity radii
-- |-- --|-- --|-- --| : rh(1,size(r)+1) = potential radii
--i-1 i i+1 i+2

Rotate horizontal winds in the east-west and north-south directions to the storm-relative direction.

calculate (lon,lat) from radial and azimuthal grid (r,theta) on the polar coordinate with the origin of (lon0,lat0)

Convert the polar grid (r-t) to the Cartesian (x-y) grid (from STPK)

Calculate RMSE of x for the "true_val" (from STPK)
Definition: RMSE $= \sqrt{\sum^{N}_{i=1}{(1/N)(x(i)-x_t)^2}},$
$x_t=\mathrm{true\_val}$

Standard output for message

subtract ival

subtract ival

Calculation of sum for 1D variable

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.

Check missing value in "val"

Check missing value in "val"

write float data from 4-byte unformatted binary

write float data from 4-byte unformatted binary

Convert the Cartesian (x-y) grid to the polar (r-t) grid (from STPK)