function extract_ECCO_frcst(FRCST_dir,FRCST_prefix,prefix,suffix,tindex,missval,... lon,lon_u,lat,lat_v,depth,... krange,jrange,jrange_v,... i1min,i1max,i2min,i2max,i3min,i3max,... i1min_u,i1max_u,i2min_u,i2max_u,i3min_u,i3max_u,... time,Yorig) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % Extract a subset from ECCO using DODS % Write it in a local file (keeping the classic % SODA netcdf format) % % % Further Information: % http://www.croco-ocean.org % % This file is part of CROCOTOOLS % % CROCOTOOLS is free software; you can redistribute it and/or modify % it under the terms of the GNU General Public License as published % by the Free Software Foundation; either version 2 of the License, % or (at your option) any later version. % % CROCOTOOLS is distributed in the hope that it will be useful, but % WITHOUT ANY WARRANTY; without even the implied warranty of % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the % GNU General Public License for more details. % % You should have received a copy of the GNU General Public License % along with this program; if not, write to the Free Software % Foundation, Inc., 59 Temple Place, Suite 330, Boston, % MA 02111-1307 USA % % Copyright (c) 2006 by Pierrick Penven % e-mail:Pierrick.Penven@ird.fr % % Updated 9-Sep-2006 by Pierrick Penven % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % disp([' Download ECCO']) trange=['[',num2str(tindex-1),':',num2str(tindex-1),']']; % % Get the time % %vname='Have'; %fname=[prefix,vname,suffix]; %time=readdap(fname,'time',trange) %time=floor(time/24+datenum(1970,1,1)); %disp([' Date: ',datestr(time)]) %time=time-datenum(Yorig,1,1) % % Get SSH % disp(' ...SSH') vname='Have'; fname=[vname,suffix]; ssh=getdap(prefix,fname,vname,... trange,'',jrange,... i1min,i1max,i2min,i2max,i3min,i3max); ssh=shiftdim(ssh,2); ssh(ssh<=missval)=NaN; % % Get U % disp(' ...U') vname='Uave'; fname=[vname,suffix]; u=getdap(prefix,fname,vname,... trange,krange,jrange,... i1min_u,i1max_u,i2min_u,i2max_u,i3min_u,i3max_u); u=permute(u,[4 3 1 2]); u(u<=missval)=NaN; % % Get V % disp(' ...V') vname='Vave'; fname=[vname,suffix]; v=getdap(prefix,fname,vname,... trange,krange,jrange_v,... i1min,i1max,i2min,i2max,i3min,i3max); v=permute(v,[4 3 1 2]); v(v<=missval)=NaN; % % Get TEMP % disp(' ...TEMP') vname='Tave'; fname=[vname,suffix]; temp=getdap(prefix,fname,vname,... trange,krange,jrange,... i1min,i1max,i2min,i2max,i3min,i3max); temp=permute(temp,[4 3 1 2]); temp(temp<=missval)=NaN; % % Get SALT % disp(' ...SALT') vname='Save'; fname=[vname,suffix]; salt=getdap(prefix,fname,vname,... trange,krange,jrange,... i1min,i1max,i2min,i2max,i3min,i3max); salt=permute(salt,[4 3 1 2]); salt(salt<=missval)=NaN; % % Create the ECCO file % create_OGCM([FRCST_dir,FRCST_prefix,num2str(time),'.cdf'],... lon,lat,lon_u,lat,lon,lat_v,depth,time,... squeeze(temp),squeeze(salt),squeeze(u),... squeeze(v),squeeze(ssh),Yorig) % return