% extract surface U/V/T/S and Eta, theta, salt, oceTauX, oceTauY
% 45 N to 55 N, -25 E to -15 E

% define desired region
nx=4320;
region_name='Osmosis3';
minlat=45;
maxlat=55;
minlon=335;
maxlon=345;
kx=1;

% extract indices for desired region
nx=4320;
prec='real*4';
gdir='/nobackupp8/dmenemen/llc/llc_4320/grid/';
fnam=[gdir 'Depth.data'];
[fld fc ix jx]=quikread_llc(fnam,nx,1,prec,gdir,minlat,maxlat,minlon,maxlon);
quikpcolor(fld')
thincolorbar

% get and save grid information
pin='/nobackupp8/dmenemen/llc/llc_4320/grid/';
pout=['/nobackupp8/dmenemen/llc/llc_4320/regions/' region_name '/grid/'];
eval(['mkdir ' pout])
eval(['cd ' pout])
for fnm={'AngleCS','AngleSN','DXC','DXG','DYC','DYG','Depth', ...
         'RAC','RAS','RAW','RAZ','U2zonDir','V2zonDir', ...
         'XC','XG','YC','YG','hFacC','hFacS','hFacW'}
    fin=[pin fnm{1} '.data'];
    fout=[fnm{1} '_' int2str(length(ix)) 'x' int2str(length(jx))];
    switch fnm{1}
      case{'hFacC','hFacS','hFacW'}
        fld=read_llc_fkij(fin,nx,fc,kx,ix,jx);
        fout=[fout 'x' int2str(length(kx))];
      otherwise
        fld=read_llc_fkij(fin,nx,fc,1,ix,jx);
    end
    writebin(fout,fld);
end

% get and save regional S/T/U/V/W
pin='/nobackupp8/dmenemen/llc/llc_4320/MITgcm/run/';
pout=['/nobackupp8/dmenemen/llc/llc_4320/regions/' region_name '/'];
for fnm={'Salt','Theta','U','V','Eta'}
    %eval(['mkdir ' pout fnm{1}])
    eval(['cd ' pout fnm{1}])
    %for ts=10368:144:279360;
    for ts=10368:144:485568%
        fin=[pin fnm{1} '.' myint2str(ts,10) '.data'];
        dy=ts2dte(ts,25,2011,9,10,30);
        fout=[fnm{1} '_' int2str(length(ix)) 'x' int2str(length(jx)) ...
              'x' int2str(length(kx)) '.' dy];
        for k=1:length(kx);% mydisp(k)
            fld=read_llc_fkij(fin,nx,fc,kx(k),ix,jx);
            writebin(fout,fld,1,'real*4',k-1);
        end
    end
end

% get and save regional surface fluxes
%for fnm={'oceFWflx','oceQnet','oceQsw','oceSflux','oceTAUX','oceTAUY'}
for fnm={'oceTAUX','oceTAUY'}

    %eval(['mkdir ' pout fnm{1}])
    eval(['cd ' pout fnm{1}])
    %for ts=10368:144:279360%, disp(ts)
    for ts=10368:144:485568%
        fin=[pin fnm{1} '.' myint2str(ts,10) '.data'];
        dy=ts2dte(ts,25,2011,9,10,30);
        fout=[fnm{1} '_' int2str(length(ix)) 'x' int2str(length(jx)) '.' dy];
        fld=read_llc_fkij(fin,nx,fc,1,ix,jx);
        writebin(fout,fld);
    end
end