Isothermal Layer Depth Computation

% USAGE: ILD=ra_ild(temp,Z,dT)
%
% DESCRIPTION: This function determines Isothermal Layer Depth (ILD) from profile data
% sets based on subjective method. If you have 3D data sets i.e. level, lat and lon and want
% to compute the ILD, then this function will be very handy. Because this function is specifically
% designed for those cases. However, it can evaluate ILD from profile data too.
%
% INPUTS:
% temp = Temperature profiles over the study region [deg. C], either 3D or vector
% Z = Levels [m], Must be vector
% dT = temperature difference criterion [deg. C], Must be scalar
%
% OUTPUT:
% ild = Isothermal layer depth, spatial output [m]
%
% DISCLAIMER:
% Albeit this function is designed only for academic purpose, it can be implemented in
% research. Nonetheless, author does not guarantee the accuracy.