Ignore:
Timestamp:
Apr 1, 2016, 7:22:25 PM (5 years ago)
Author:
sommeria
Message:

group writing access set for result files

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/calc_field_tps.m

    r924 r938  
    122122    end
    123123    for ilist=1:length(FieldName)
    124         %Operator{ilist}='';%default empty operator (vec, norm,...)
    125         %r=regexp(FieldName{ilist},'(?<Operator>(^vec|^norm|^curl|^div|^strain))\((?<UName>.+),(?<VName>.+)\)$','names');% TODO, replace U, V
    126124        switch FieldName{ilist}
    127125            case 'vec(U,V)'
    128 %                 ListVar=[ListVar {'U', 'V'}];
    129 %                 VarAttribute{var_count+1}.Role='vector_x';
    130 %                 VarAttribute{var_count+2}.Role='vector_y';
    131126                DataOut.U(ind_sel)=DataOut.U(ind_sel)+EM *FieldVar(1:nbvec_sub+3,isub,1);
    132127                DataOut.V(ind_sel)=DataOut.V(ind_sel)+EM *FieldVar(1:nbvec_sub+3,isub,2);
    133128            case 'U'
    134 %                 ListVar=[ListVar {'U'}];
    135 %                 VarAttribute{var_count+1}.Role='scalar';
    136129                DataOut.U(ind_sel)=DataOut.U(ind_sel)+EM *FieldVar(1:nbvec_sub+3,isub,1);
    137130            case 'V'
    138 %                 ListVar=[ListVar {'V'}];
    139 %                 VarAttribute{var_count+1}.Role='scalar';
    140131                DataOut.V(ind_sel)=DataOut.V(ind_sel)+EM *FieldVar(1:nbvec_sub+3,isub,2);
    141132            case 'norm(U,V)'
    142 %                 ListVar=[ListVar {'norm'}];
    143 %                 VarAttribute{var_count+1}.Role='scalar';
    144133                U=DataOut.U(ind_sel)+EM *FieldVar(1:nbvec_sub+3,isub,1);
    145134                V=DataOut.V(ind_sel)+EM *FieldVar(1:nbvec_sub+3,isub,2);
    146135                DataOut.norm(ind_sel)=sqrt(U.*U+V.*V);
    147136            case 'curl(U,V)'
    148 %                 ListVar=[ListVar {'curl'}];
    149 %                 VarAttribute{var_count+1}.Role='scalar';
    150137                DataOut.curl(ind_sel)=DataOut.curl(ind_sel)-EMDY *FieldVar(1:nbvec_sub+3,isub,1)+EMDX *FieldVar(1:nbvec_sub+3,isub,2);
    151138            case 'div(U,V)'
    152 %                 ListVar=[ListVar {'div'}];
    153 %                 VarAttribute{var_count+1}.Role='scalar';
    154139                DataOut.div(ind_sel)=DataOut.div(ind_sel)+EMDX*FieldVar(1:nbvec_sub+3,isub,1)+EMDY *FieldVar(1:nbvec_sub+3,isub,2);
     140%             case 'curl_polar(U,V)'
     141%                 DataOut.curl(ind_sel)=DataOut.curl(ind_sel)+(-EMDY *FieldVar(1:nbvec_sub+3,isub,1)+EMDX *(Coord_interp(ind_sel,1).*FieldVar(1:nbvec_sub+3,isub,2)))./Coord_interp(ind_sel,1);
    155142            case 'strain(U,V)'
    156 %                 ListVar=[ListVar {'strain'}];
    157 %                 VarAttribute{var_count+1}.Role='scalar';
    158143                DataOut.strain(ind_sel)=DataOut.strain(ind_sel)+EMDY*FieldVar(1:nbvec_sub+3,isub,1)+EMDX *FieldVar(1:nbvec_sub+3,isub,2);
    159144        end
    160145    end
    161146end
    162 %DataOut.FF=nbval==0; %put errorflag to 1 for points outside the interpolation rang
    163147nbval(nbval==0)=NaN;% to avoid division by zero for averaging
    164148ListFieldOut=fieldnames(DataOut);
Note: See TracChangeset for help on using the changeset viewer.