Changeset 46 for trunk/src/proj_field.m


Ignore:
Timestamp:
Mar 15, 2010, 5:42:11 PM (11 years ago)
Author:
gostiaux
Message:

various bug fixes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/proj_field.m

    r42 r46  
    9090    return
    9191end
     92FieldData
    9293%introduce default field properties (reading old standards)
    9394if ~isfield(ObjectData,'Style')||~isfield(ObjectData,'Coord')||~isfield(ObjectData,'ProjMode')
     
    885886            ProjData.VarAttribute{nbvar+1}.long_name='abscissa along line';
    886887            if nbcolor==3
    887                 ProjData.VarDimName{end}={XName,'rgb'};
     888                ProjData.VarDimName{end}={AXName,'rgb'};
    888889            end
    889890        end     
     
    11771178            for ivar=VarIndex
    11781179                VarName=FieldData.ListVarName{ivar};
    1179                 if ~( ivar==ivar_X | ivar==ivar_Y | ivar==ivar_Z | ivar==ivar_F | ivar==ivar_FF | test_anc(ivar)==1)                 
     1180                if ~( ivar==ivar_X || ivar==ivar_Y || ivar==ivar_Z || ivar==ivar_F || ivar==ivar_FF || test_anc(ivar)==1)                 
    11801181                    ivar_new=ivar_new+1;
    11811182                    ProjData.ListVarName=[ProjData.ListVarName {VarName}];
    11821183                    ProjData.VarDimName=[ProjData.VarDimName {DimCell}];
    11831184                    %ProjData.VarDimIndex=[ProjData.VarDimIndex {[1 2]}];
    1184                     if isfield(FieldData,'VarAttribute') & length(FieldData.VarAttribute) >=ivar
     1185                    if isfield(FieldData,'VarAttribute') && length(FieldData.VarAttribute) >=ivar
    11851186                        ProjData.VarAttribute{ivar_new+nbcoord}=FieldData.VarAttribute{ivar};
    11861187                    end
     
    11901191                        eval(['FieldData.' VarName '=FieldData.' VarName '(indsel);'])
    11911192                    end
    1192                     eval(['ProjData.' VarName '=griddata_uvmat(coord_X,coord_Y,FieldData.' VarName ',coord_x_proj,coord_y_proj'');'])
     1193                    eval(['ProjData.' VarName '=griddata_uvmat(double(coord_X),double(coord_Y),double(FieldData.' VarName '),coord_x_proj,coord_y_proj'');'])
    11931194                    eval(['varline=reshape(ProjData.' VarName ',1,length(coord_y_proj)*length(coord_x_proj));'])
    11941195                    FFlag= isnan(varline); %detect undefined values NaN
     
    15171518    end
    15181519    %projection of  velocity components in the rotated coordinates
    1519     if ~isequal(Phi,0) & length(ivar_U)==1
     1520    if ~isequal(Phi,0) && length(ivar_U)==1
    15201521        if isempty(ivar_V)
    15211522            msgbox_uvmat('ERROR','v velocity component missing in proj_field.m')
Note: See TracChangeset for help on using the changeset viewer.