processed PIV data for wave-current flow over vortex ripples

Processed PIV data: Coordinates:Vertical coordinate, the 0 datum is set at ripple trough. Data structure For wave current cases: 1DV profile:Ripple_ave.u: horizontal velocity profile @ 32 phases [m/s]Ripple_ave.w: vertical velocity profile @ 32 phases [m/s]Ripple_ave.z: vertical coordinates [mm]Ripple_ave.TKE: Turbulence kinetic energy profile @ 32 phasesRipple_ave.tau_conv: convective stress profile @ 32 phasesRipple_ave.tau_Re: Reynolds-stress profile @ 32 phasesRipple_ave.theta: phase (32 phases)Ripple_ave.u0: time-averaged velocity profileRipple_ave.TKE0: time-averaged Turbulence kinetic energy profileRipple_ave.tau_conv0: time-averaged convective stress profileRipple_ave.tau_Re0: time-averaged Reynolds-stress profile 2D fields:The size of gridded data is I- vertical direction, J-horizontal direction, M-phasesThe 2D fields cover one ripple length. The ripple trough is at the middle, while the ripple crests are at the two lateral ends.The 2D fields are stored in a structure called "Fields_g". The elements are:x: horizontal coordinate, size=J, unit= [mm]z: vertical coordinate, size=I, unit =[mm]theta: phase coordinate, size=M, unit= [degree]u_pa: phase-averaged horizontal velocity. size=[I,J,M], unit=[m/s]w_pa: phase-averaged horizontal velocity. size=[I,J,M], unit=[m/s]TKE: turbulence kinetic energy, size=[I,J,M], unit=[m^2/s^2]tau_Re: turbulence Reynolds stress, size=[I,J,M], unit=[m^2/s^2]w_rms: root-mean-square value of vertical turbulence component. size=[I,J,M], unit=[m/s]u_rms: root-mean-square value of horizontal turbulence component. size=[I,J,M], unit=[m/s]vorticity: vorticity based on (u_pa,w_pa), size=[I,J,M], unit=[1/s] Form drag:The form drag, which is given by integrating pressure over the bottom surface, is stored in structure Form_drag, the elements are:x: horizontal coordinate, size=J, unit= [mm]zb: vertical coordinate of the bottom surface @x, , size=J, unit= [mm]pb: pressure on the bottom surface, size=[I,J,M], unit=[m^2/s^2]f: normalized form_drag , size=Mtheta: : phase coordinate, size=M, unit= [degree] Bottom_g: a structure stores the bottom profilexb: x-coordinate;zb: vertical coordinate of bottom surface @xb. example: plot(xb,zb) gives the bottom profile