Changeset 748 for trunk/src/proj_field.m
- Timestamp:
- Apr 23, 2014, 4:57:58 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/proj_field.m
r747 r748 1360 1360 VarDimName=[VarDimName {DimCell}]; 1361 1361 VarAttribute{length(ListVarName)}=FieldData.VarAttribute{ivar}; %reproduce the variable attributes 1362 eval(['ProjData.' VarName '=permute(FieldData.' VarName ',ind_new);'])% permute x and z indices for 90 degree rotation1363 eval(['ProjData.' VarName '=squeeze(ProjData.' VarName '(iz,:,:));'])% select the z index iz1362 ProjData.(VarName)=permute(FieldData.(VarName),ind_new);% permute x and z indices for 90 degree rotation 1363 ProjData.(VarName)=squeeze(ProjData.(VarName)(iz,:,:));% select the z index iz 1364 1364 end 1365 1365 ProjData.(AYName)=[Ybound(1) Ybound(2)]; %record the new (projected ) y coordinates … … 1390 1390 end 1391 1391 end 1392 ProjData.(AYName)=Coord{NbDim-1}(1)+DY*(YIndexRange-1); %record the new (projected ) y coordinates 1393 ProjData.(AXName)=Coord{NbDim}(1)+DX*(XIndexRange-1); %record the new (projected ) x coordinates 1392 if testXMax 1393 ProjData.(AXName)=Coord{NbDim}(1)+DX*(XIndexRange-1); %record the new (projected ) x coordinates 1394 else 1395 ProjData.(AXName)=FieldData.(AXName); 1396 end 1397 if testYMax 1398 ProjData.(AYName)=Coord{NbDim-1}(1)+DY*(YIndexRange-1); %record the new (projected ) x coordinates 1399 else 1400 ProjData.(AYName)=FieldData.(AYName); 1401 end 1394 1402 end 1395 1403 end
Note: See TracChangeset
for help on using the changeset viewer.