source: altifloat/exp_twin_aviso/script_prepare_exp_twin.m @ 199

addpath('../matlab_toolbox');
rootdir = './'

%STEP ONE : produce the meshfile

%meshgrid Aviso
meshfile=[rootdir 'meshgrid_aviso.dat'];

%important variables because it is read by script bash
%WARNING : PUT A SPACE AROUND EQUAL SIGN AND BEFORE SEMI-COLON
nlon = 87 ;
nlat = 58 ;
make_meshgrid(meshfile,30.0625,26.1875,37.1875,36.9375,nlon,nlat);

%% STEP TWO : transform and interpolate aviso file to YAO file to make the truth
dir_uvfiles='../../data/AVISO/';
aviso_netcdf={[dir_uvfiles 'dt_med_allsat_madt_uv_20130828_20140704.nc'],...
              [dir_uvfiles 'dt_med_allsat_madt_uv_20130829_20140704.nc'],...
              [dir_uvfiles 'dt_med_allsat_madt_uv_20130830_20140704.nc']};

it=1:24:1+24*(length(aviso_netcdf)-1);

outfile=[rootdir 'uv_truth.dat'];
outmask=[rootdir 'mask.dat'];

dt=1;

interp_bck(aviso_netcdf,it,dt,meshfile,outfile,outmask);

%% STEP THREE : Prepare filter
make_filter(meshfile,outmask,20000,rootdir);

%% STEP THREE BIS : Prepare wind
windfile = [rootdir 'wind_interp.nc'];
infile = '../../data/wind/wind_20130825-20130930.nc';
interp_wind(infile,meshfile,windfile);


%% STEP FOUR: Prepare init
initfile=[rootdir 'init5.dat'];

%% STEP FIVE
% RUN forward.sh to set grid and filter parameters
% Modify if necessary the number of floater and the number of time steps
% make the project (with -DINC_FIL and -DINC_GRID options)
% Make a symbolic link on the executable
% run the project with forward.i

%% STEP SIX : modify obs file (e.g. subsampling)
obsallfile=[rootdir 'rfloat_total.dat'];
obsfile=[rootdir 'obs_jum.dat'];

init=load(initfile);
obs=load(obsallfile);

istep=5;
I=find(obs(:,end)>1e-8);
obs=obs(I,:);

obs=[init;obs];
obs=sortrows(obs,[2 1]);

t=unique(obs(:,1));
tech=t(1:istep:end);

I=find(ismember(obs(:,1),tech));
obs=obs(I,:);
dlmwrite(obsfile,obs,'\t');

%% STEP SEVEN : Create bck
aviso_netcdf={[dir_uvfiles 'dt_med_allsat_madt_uv_20130912_20140704.nc'],...
              [dir_uvfiles 'dt_med_allsat_madt_uv_20130913_20140704.nc'],...
              [dir_uvfiles 'dt_med_allsat_madt_uv_20130914_20140704.nc']};

outfile=[rootdir 'uv_bck.dat'];
interp_bck(aviso_netcdf,it,dt,meshfile,outfile,outmask);


%% STEP EIGHT : ASSIMILATION
% run the project with assim.i