112 |
if kpr > 0 & krd > 0, |
if kpr > 0 & krd > 0, |
113 |
rac=dir; |
rac=dir; |
114 |
%if nBas < 0, hFacC=ones(ncx,nc,dim3); else hFacC=rdmds([rac,'hFacC']); end |
%if nBas < 0, hFacC=ones(ncx,nc,dim3); else hFacC=rdmds([rac,'hFacC']); end |
115 |
%nr=size(hFacC,3); |
nr_hfac=size(hFacC,3); |
116 |
if nBas < 0, hFacC=ones(ncx,nc,nr); end%else hFacC=rdmds([rac,'hFacC']); end |
if nBas < 0, hFacC=ones(ncx,nc,nr); end%else hFacC=rdmds([rac,'hFacC']); end |
117 |
mskC=min(1,hFacC); mskC=ceil(mskC); |
mskC=min(1,hFacC); mskC=ceil(mskC); |
118 |
hFacC=reshape(hFacC,nPg,nr); mskC=reshape(mskC,nPg,nr); |
hFacC=reshape(hFacC,nPg,nr_hfac); mskC=reshape(mskC,nPg,nr_hfac); |
119 |
a6c=reshape(arc,nPg,1); |
a6c=reshape(arc,nPg,1); |
120 |
if nBas > 0, |
if nBas > 0, |
121 |
%- standard 3. Sector definition: |
%- standard 3. Sector definition: |