function a=add_e38(a)
% if one matrix includes the nan and has the trouble to contour, this is for
% changing the NaN to 1.e+38

b=find(isnan(a));
a(b)=1e+38*ones(size(b));
