Commit f1d33131da88b1a2963556c2fc7bf481cc8f24ce

Authored by Marion Bandet
1 parent 06df0a3d
Exists in master

update

source/Comparison_WERA-ADCP_Sept2013.rst 0 → 100644
... ... @@ -0,0 +1,77 @@
  1 +Comparison WERA and ADCP radial velocities September 2013
  2 +===============================================================
  3 +
  4 +Matlab script used to do the analysis is **RadialVeloCompa_Sept2013_HFR_ADCP_WERA_ADCPcorrected_v2.m**
  5 +
  6 +Steps of the algorithm:
  7 +
  8 +1- we choose a site to analyze, either PAO (site 1) or PAB (site 2)
  9 +
  10 +2- we load the bathymetry and the HFR sites'coordinates
  11 +
  12 +3- we specify if we use the "old" or "new" calibrations. Usually "New".
  13 +
  14 +4- we load IML-4 ADCP data (2013) and gps data
  15 +
  16 +5- we make a list of all the HFR files available for the site of interest.
  17 +
  18 +6- we analyze the ADCP data (bin 2 only to get near-surface data):
  19 +
  20 + We look for the indices of the ADCP time vector corresponding to the time period of interest (here 01 to 30 September 2013).
  21 +
  22 + For each index, we look for the ADCP position (xadcp,yadcp)(via the gps coordinates).
  23 +
  24 + We then look for the positions 15min before and 15min after.
  25 +
  26 + We estimate the mean drift of the ADCP within these 30 minutes (Umean_adcp_u, Umean_adcp_v and U_mean_adcp (total)).
  27 +
  28 + We find the radial orientation at the ADCP position given by xadcp,yadcp: azimuth in degrees from HFR site to adcp position
  29 +
  30 + We project U_mean_adcp onto the radial direction to get the radial current components (ur_adcp_proj, vr_adcp_proj).
  31 +
  32 + We store all the meaningful parameters into a structure named adcp.
  33 +
  34 +7- We analyze the HFR data:
  35 +
  36 + We look for the indices of the HFR time vector corresponding to the time period of interest (here 01 to 30 September 2013)
  37 +
  38 + For each index, we load the corresponding HFR file (i.e. 20132440000_pao.mat). Each file contains a grid (X,Y) with radial currents at each grid point (Ur), a corresponding time vector (t) and the angle and distance from the site (theta,r).
  39 +
  40 +[[Fichier:Exple_allgridpt.png]] [[Fichier:Gridwithnan.png]]
  41 +
  42 +
  43 + We then look for the HFR grid point nearest the ADCP position at the time of interest.
  44 +
  45 +[[Fichier:gridwithnan_adcp.png]]
  46 +
  47 + We find the HFR Radial velocity at this position and project it onto zonal and meridional components.
  48 +
  49 + We store all the meaningful parameters into a structure named radar.
  50 +
  51 +
  52 +This analysis allows us to compare radial velocity components from each HFR site with the ADCP components.
  53 +
  54 +
  55 +
  56 +To verify that the HFR look in the right direction:
  57 +
  58 +For each time step of the time period of interest, we compute the distance from the HFR site to each grid point. We select the grid points that are at the same range from the HFR site than the grid point closest to the ADCP.
  59 +
  60 +[[Fichier:rangeADCP.png]]
  61 +
  62 +For each of these grid point (in cyan in the figure), we compute the azimuth the the site coordinates and we extract the radial velocity (total) and its components (zonal and meridional).
  63 +
  64 +We store all the meaningful parameters into a structure named range.
  65 +
  66 +Results from the whole analysis are stored, for example, under:
  67 +/run/media/bandma01/SeagateWireless/MEOPAR/Matlab/Analysis/Radial_Comparison/201309/ADCP-PAO/ADCP_PAO_version5adcpCorrected_newcalibrations_v2.mat
  68 +
  69 +
  70 +[[Fichier:CompsNewCalib_RadialVelo_ADCP_PAO_Sept2013.png|200px|thumb|left|Comparison of radial velocities from ADCP bin 2 and HFR WERA at PAO site September 2013]]
  71 +
  72 +
  73 +[[Fichier:orientationPAOSept2013.png|200px|thumb|left|RMS diff and correlation from radial velocities compared betwen ADCP bin 2 and HFR WERA at PAO site September 2013]]
  74 +
  75 +
  76 +[[Fichier:ADCPbin2vsPAOSept2013.png|200px|thumb|left|radial velocities from ADCP bin 2 vs HFR WERA at PAO site September 2013 - in red: 48h averages]]]]
  77 +
... ...
source/OGSL.rst deleted
... ... @@ -1,3 +0,0 @@
1   -OGSL quasi-realtime data feed
2   -=============================
3   -
source/WERAprocessing.rst
1   -WERA POST-PROCESSING
2   -====================
  1 +WERA radial currents processing
  2 +==================================
3 3  
4 4 Post-processing of WERA radial currents is much simpler than for the CODAR system.
5 5  
... ...
source/index.rst
... ... @@ -20,8 +20,8 @@ Documentation about the steps of the analysis
20 20  
21 21 CODARprocessing
22 22 WERAprocessing
  23 + Comparison_WERA-ADCP_Sept2013
23 24 CurrentsCombination
24   - OGSL
25 25  
26 26  
27 27  
... ...