Changeset 93 for trunk/src/get_field


Ignore:
Timestamp:
May 27, 2010, 7:25:58 PM (15 years ago)
Author:
sommeria
Message:

FFT: improved to deal with NaN data
merge_proj: corrected to reproduce dt if unique
uvmat.fig: tooltip corrected
plot_field: bug for isocontour corrected,
im_filter: cleaning
phys_polar: spatial derivative included (still to check)
set_obeject.fig: minor correction
struct2nc: comments improved
uvmat: button NB implemented
read_civxdata: error message improved

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/get_field/FFT.m

    r89 r93  
    66    return %exit the function
    77end
    8 
     8GUI_input=[];
    99%initiation
    1010hhget_field=guidata(hget_field);
     
    2929% get variable
    3030eval(['Var= Field.' ordinate_name ';']);
    31 np=size(Var);
     31np=size(Var)
    3232np_freq=floor(np(1)/2);
    3333dx=1;%default
    3434dfreq=1/np(1);%default frequency interval (abscissa= array index)
     35sum_data=sum(Var,2);
    3536if ~isequal(abscissa_name,'')
    3637    eval(['Coord_x= Field.' abscissa_name ';']);
    37     ind_select=find(~isinf(Coord_x));%detect infinite values
     38    ind_select=find(~isinf(Coord_x)&~isnan(sum_data));%detect infinite values
    3839    Coord_x=Coord_x(ind_select);
    3940    Var=Var(ind_select,:);
Note: See TracChangeset for help on using the changeset viewer.