Changeset 494 for trunk/src/proj_field.m

Timestamp:

Jul 19, 2012, 11:50:39 AM (12 years ago)

Author:

sommeria

Message:

various bugs corrected after testing in Windows OS. Introduction
of filter tps

File:

• trunk/src/proj_field.m (modified) (3 diffs)

trunk/src/proj_field.m

@@ -1238 +1238 @@
                         ProjData.VarAttribute{ivar_new+1+nbcoord}.Role='errorflag';
                     end
+%                     case 'filter'%interpolate data on a regular grid
+%                         errormsg='tps required for filter option'
+                        
             end
             
@@ -2076 +2079 @@
                 eval(['ProjData.' AXName '=[Xbound(1) Xbound(2)];']) %record the new (projected ) x coordinates
             end
-        else       % case with rotation and/or interpolation
+        elseif isfield(FieldData,'A') %TO GENERALISE       % case with rotation and/or interpolation
             if NbDim==2 %2D case
                 [X,Y]=meshgrid(coord_x_proj,coord_y_proj);%grid in the new coordinates
@@ -2104 +2107 @@
                     if test_filter  
                          Aclass=class(FieldData.A);
-                         eval(['ProjData.' VarName '=filter2(Mfilter,FieldData.' VarName ',''valid'');'])
+                         ProjData.(VarName)=filter2(Mfilter,FieldData.(VarName),'valid');
                          if ~isequal(Aclass,'double')
-                             eval(['ProjData.' VarName '=' Aclass '(FieldData.' VarName ');'])%revert to integer values
+                             ProjData.(VarName)=Aclass(FieldData.(VarName));%revert to integer values
                          end
                     end