Changeset 1142 for trunk/src/series

Timestamp:

May 5, 2024, 12:29:30 PM (8 months ago)

Author:

sommeria

Message:

filter parameter improved in filter_tps

File:

• trunk/src/series/test_patch_tps.m (modified) (4 diffs)

trunk/src/series/test_patch_tps.m

@@ -33 +33 @@
 %              .FieldName: name(s) of the field
 %              .VelType: velocity type
-%              .FieldName_1: name of the second field in case of two input series
-%              .VelType_1: velocity type of the second field in case of two input series
-%              .Coord_y: name of y coordinate variable
-%              .Coord_x: name of x coordinate variable
 
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
@@ -110 +106 @@
             case {4,5}
                 CivStage='civ2';
-                MaxDiff=1.5; SubDomainSize=500; FieldSmooth=2; %default
+                MaxDiff=1.5; SubDomainSize=250; FieldSmooth=5; %default
             case 6
                 CivStage='civ2';
@@ -211 +207 @@
 
 %% Prepare the structure of output netcdf file
-DataOut.ListGlobalAttribute={'CivStage'};
+DataOut.ListGlobalAttribute={'CivStage','SubDomainSize','MaxDiff','CoordUnit'};
 DataOut.CivStage=Param.InputFields.VelType;
-DataOut.ListVarName=[{'FieldSmooth','Diff_rms','NbExclude','FF','X','Y'} Ustr Vstr] ;
-DataOut.VarDimName=[{'FieldSmooth','FieldSmooth','FieldSmooth','NbVec','NbVec','NbVec'} Dimstr Dimstr]; 
+DataOut.SubDomainSize=SubDomainSize;
+DataOut.MaxDiff=MaxDiff;
+DataOut.CoordUnit='pixel';
+DataOut.ListVarName=[{'FieldSmooth','Diff_rms','NbExclude','X','Y'} Ustr Vstr] ;
+DataOut.VarDimName=[{'FieldSmooth','FieldSmooth','FieldSmooth','NbVec','NbVec'} Dimstr Dimstr]; 
 DataOut.VarAttribute{4}.Role='falseflag';
 DataOut.FieldSmooth=FieldSmooth;
@@ -249 +248 @@
         DataOut.Diff_rms(irho)=sqrt(mean(U_Diff.*U_Diff+V_Diff.*V_Diff)/2);
         DataOut.NbExclude(irho)=(NbGood-numel(ind_good))/NbGood;
+        U_smooth(ind_false)=NaN;
+        V_smooth(ind_false)=NaN;
         DataOut.(['U_' str_i{irho}])=U_smooth;
         DataOut.(['V_' str_i{irho}])=V_smooth;
-        DataOut.FF(ind_false)=FieldSmooth(irho);
+        % DataOut.FF(ind_false)=FieldSmooth(irho);
     end
     time=toc