Changeset 247 for trunk/src/calc_field.m
 Timestamp:
 Apr 29, 2011, 6:44:35 PM (13 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

trunk/src/calc_field.m
r246 r247 105 105 EM = tps_eval(epoints,ctrs); 106 106 case{'vort','div'} 107 EMDXY= tps_eval_dxy(epoints,ctrs);%2D matrix of distances between extrapolation points epoints and spline centres (=site points) ctrs107 [EMDX,EMDY] = tps_eval_dxy(epoints,ctrs);%2D matrix of distances between extrapolation points epoints and spline centres (=site points) ctrs 108 108 109 109 end … … 132 132 % DataIn.V_tps(nbvec_sub+2,isub)=0;% X coefficient suppressed for x wise derivatives 133 133 % DataIn.U_tps(nbvec_sub+3,isub)=0;% Y coefficient suppressed for x wise derivatives 134 DataOut.vort(ind_sel)=DataOut.vort(ind_sel)+EMD XY(:,:,2) *DataIn.U_tps(1:nbvec_sub+3,isub)EMDXY(:,:,1)*DataIn.V_tps(1:nbvec_sub+3,isub);134 DataOut.vort(ind_sel)=DataOut.vort(ind_sel)+EMDY *DataIn.U_tps(1:nbvec_sub+3,isub)EMDX *DataIn.V_tps(1:nbvec_sub+3,isub); 135 135 case 'div' 136 136 ListFields={'div'}; … … 142 142 % DataIn.V_tps(nbvec_sub+2,isub)=0;% X coefficient suppressed for x wise derivatives 143 143 % DataIn.U_tps(nbvec_sub+3,isub)=0;% Y coefficient suppressed for x wise derivatives 144 DataOut.div(ind_sel)=DataOut.div(ind_sel)+EMDX Y(:,:,1)*DataIn.U_tps(1:nbvec_sub+3,isub)+EMDXY(:,:,2)*DataIn.V_tps(1:nbvec_sub+3,isub);144 DataOut.div(ind_sel)=DataOut.div(ind_sel)+EMDX*DataIn.U_tps(1:nbvec_sub+3,isub)+EMDY *DataIn.V_tps(1:nbvec_sub+3,isub); 145 145 end 146 146 end
Note: See TracChangeset
for help on using the changeset viewer.