Ignore:
Timestamp:
Jan 27, 2015, 8:53:50 PM (6 years ago)
Author:
sommeria
Message:

Testpatch_fixed

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/series/civ_series.m

    r859 r860  
    737737        Data.Civ2_C=reshape(ctable,[],1);
    738738        Data.Civ2_F=reshape(F,[],1);
     739    elseif ~isfield(Data,'ListVarName') % we start there, using existing Civ2 data
     740        if exist('ncfile','var')
     741            CivFile=ncfile;
     742            [Data,tild,tild,errormsg]=nc2struct(CivFile);%read civ1 and fix1 data in the existing netcdf file
     743            if ~isempty(errormsg)
     744                disp_uvmat('ERROR',errormsg,checkrun)
     745                return
     746            end         
     747        elseif isfield(Param,'Civ2_X')% use Civ2 data as input in Param (test mode)
     748            Data.ListGlobalAttribute={};
     749            Data.ListVarName={};
     750            Data.VarDimName={};
     751            Data.Civ2_X=Param.Civ2_X;
     752            Data.Civ2_Y=Param.Civ2_Y;
     753            Data.Civ2_U=Param.Civ2_U;
     754            Data.Civ2_V=Param.Civ2_V;
     755            Data.Civ2_FF=Param.Civ2_FF;
     756        end
    739757    end
    740758   
     
    786804        Data.ListGlobalAttribute=[Data.ListGlobalAttribute Patch2_param];
    787805       
    788        
    789         %         Data.ListGlobalAttribute=[Data.ListGlobalAttribute {'Patch2_FieldSmooth','Patch2_MaxDiff','Patch2_SubDomainSize'}];
    790         %         Data.Patch2_FieldSmooth=Param.ActionInput.Patch2.FieldSmooth;
    791         %         Data.Patch2_MaxDiff=Param.ActionInput.Patch2.MaxDiff;
    792         %         Data.Patch2_SubDomainSize=Param.ActionInput.Patch2.SubDomainSize;
    793806        nbvar=length(Data.ListVarName);
    794807        Data.ListVarName=[Data.ListVarName {'Civ2_U_smooth','Civ2_V_smooth','Civ2_SubRange','Civ2_NbCentres','Civ2_Coord_tps','Civ2_U_tps','Civ2_V_tps'}];
     
    801814        Data.VarAttribute{nbvar+6}.Role='vector_x';
    802815        Data.VarAttribute{nbvar+7}.Role='vector_y';
    803         Data.Civ2_U_smooth=zeros(size(Data.Civ2_X));
    804         Data.Civ2_V_smooth=zeros(size(Data.Civ2_X));
     816        Data.Civ2_U_smooth=Data.Civ2_U; % zeros(size(Data.Civ2_X));
     817        Data.Civ2_V_smooth=Data.Civ2_V; %zeros(size(Data.Civ2_X));
    805818        if isfield(Data,'Civ2_FF')
    806819            ind_good=find(Data.Civ2_FF==0);
Note: See TracChangeset for help on using the changeset viewer.