My Ocean Ferrybox with Matlab Pierre Jaccard 10

My. Ocean Ferrybox with Matlab Pierre Jaccard 10. 06. 2013 - My. Ocean Tutorial, NERSC, Bergen 1

Requirements • Matlab (version? ) • Must have ftp and ncread Pierre Jaccard 10. 06. 2013 - My. Ocean Tutorial, NERSC, Bergen 2

Read net. CDF ncf = '. . DataGO_LATEST_TS_FB_Color. Fantasy_20130607. nc'; nci = ncinfo(ncf); arrayfun(@(x)(x. Name), nci(1). Variables, 'uniformoutput', false) t = ncread(ncf, 'TIME'); T = ncread(ncf, 'TEMP'); Pierre Jaccard 10. 06. 2013 - My. Ocean Tutorial, NERSC, Bergen 3

Calculate Time t = ncread(ncf, 'TIME'); t 0 = datenum(1950, 1, 1); t 1 = t 0 + t; datestr(t 1(1)) datestr(t 1(end)) Pierre Jaccard 10. 06. 2013 - My. Ocean Tutorial, NERSC, Bergen 4

Plot Temperature I T = ncread(ncf, 'TEMP'); plot(t 1, T); grid on datetick Pierre Jaccard 10. 06. 2013 - My. Ocean Tutorial, NERSC, Bergen 5

Quality Flags qc. T = ncread(ncf, qct = qct'; qcp = ncread(ncf, qcp = qcp'; qcm = (qc. T > 0) & any(qcm) all(qcm) 'TEMP_QC'); 'TIME_QC'); 'POSITION_QC'); (qc. T < 3) & (qct > 0) & (qct < 3) & (qcp > 0) & (qcp < 3); Pierre Jaccard 10. 06. 2013 - My. Ocean Tutorial, NERSC, Bergen 6

Plot Temperature II figure; hold on plot(t 1(qcm), T(qcm), 'g. ', 'markersize', 2); plot(t 1(~qcm), T(~qcm), 'r. ', 'markersize', 2); grid on datetick Pierre Jaccard 10. 06. 2013 - My. Ocean Tutorial, NERSC, Bergen 7

FTP host = 'vftp 1. ifremer. fr'; path = '/INSITU_GLO_NRT_OBSERVATIONS_013_030/'; f = ftp(host, user, pass); cd(f, path); list = dir(f); mget(ftp, remote, local) close(f); Pierre Jaccard 10. 06. 2013 - My. Ocean Tutorial, NERSC, Bergen 8