| 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_hfac=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_hfac); mskC=reshape(mskC,nPg,nr_hfac); |
%hFacC=reshape(hFacC,nPg,nr_hfac); mskC=reshape(mskC,nPg,nr_hfac); |
| 119 |
|
if nr == 1, hFacC = hFacC(:,:,1); mskC = mskC(:,:,1); end % Cludge by Daniel. |
| 120 |
|
hFacC=reshape(hFacC,nPg,nr); mskC=reshape(mskC,nPg,nr); |
| 121 |
a6c=reshape(arc,nPg,1); |
a6c=reshape(arc,nPg,1); |
| 122 |
if nBas > 0, |
if nBas > 0, |
| 123 |
%- standard 3. Sector definition: |
%- standard 3. Sector definition: |