NCCOOS Trac Projects: Top | Web | Platforms | Processing | Viz | Sprints | Sandbox | (Wind)

Changeset 518

Show
Ignore:
Timestamp:
01/23/14 14:45:08
Author:
cbc
Message:

Add lat/lon and include non-nan final_vx.

Files:

Legend:

Unmodified
Added
Removed
Modified
Copied
Moved
  • gliderproc/trunk/gliderCurrents_Generate_L1_Data.m

    r517 r518  
    6969        vy = []; 
    7070        prior_seg_vy = []; 
     71        lat = []; 
     72        lon = []; 
     73        gps_lat = []; 
     74        gps_lon = []; 
    7175        ptime_dbd = []; 
    7276        ptime_ebd = []; 
     
    123127                        data.data(:,strmatch('x_prior_seg_water_vy',... 
    124128                                             data.vars, 'exact'))]; 
     129                    lat = [lat; ... 
     130                        data.data(:,strmatch('m_lat',... 
     131                                             data.vars, 'exact'))]; 
     132                    lon = [lon; ... 
     133                        data.data(:,strmatch('m_lon',... 
     134                                             data.vars, 'exact'))]; 
     135                    gps_lat = [gps_lat; ... 
     136                        data.data(:,strmatch('m_gps_lat',... 
     137                                             data.vars, 'exact'))]; 
     138                    gps_lon = [gps_lon; ... 
     139                        data.data(:,strmatch('m_gps_lon',... 
     140                                             data.vars, 'exact'))]; 
    125141                    ptime_dbd = [ptime_dbd; ... 
    126142                        data.data(:,strmatch('m_present_time',... 
     
    136152        %*** END READ IN DBD DATA **** 
    137153         
    138         % remove nans from DBD data... 
    139         i = find(~isnan(vx)); 
     154        % include only not NaNs from DBD data for either vx or final_vx... 
     155        i = find(~(isnan(vx) & isnan(final_vx))); 
    140156        vx = vx(i); 
    141157        final_vx = final_vx(i); 
     
    148164        delta_vy = delta_vy(i); 
    149165        prior_seg_vy = prior_seg_vy(i); 
     166        lat = lat(i); 
     167        lon = lon(i); 
     168        gps_lat = gps_lat(i); 
     169        gps_lon = gps_lon(i); 
    150170        ptime_dbd = ptime_dbd(i); 
    151171        ptime_ebd = ptime_ebd(i); 
     
    165185        delta_vy = delta_vy(I); 
    166186        prior_seg_vy = prior_seg_vy(I); 
     187        lat = lat(I); 
     188        lon = lon(I); 
     189        gps_lat = gps_lat(I); 
     190        gps_lon = gps_lon(I); 
    167191        ptime_ebd = ptime_ebd(I); 
    168192 
     
    183207            'delta_vy', 'm/s', ... 
    184208            'prior_seg_vy', 'm/s', ... 
     209            'lat', 'decimal degrees', ... 
     210            'lon', 'decimal degrees', ... 
     211            'gps_lat', 'decimal degrees', ... 
     212            'gps_lon', 'decimal degrees', ... 
    185213            'ptime_dbd', 'seconds since 0000-01-01T00:00', ... 
    186214            'ptime_dbd_datenum', 'days since 1970-01-01T00:00', ... 
     
    199227            'delta_vy', 'Change in northward water velocity', ... 
    200228            'prior_seg_vy', 'Northward water velocity of prior segment', ... 
     229            'lat', 'north/south vehicle position', ... 
     230            'lon', 'east/west vehicle position', ... 
     231            'gps_lat', 'GPS fix latitude', ... 
     232            'gps_lon', 'GPS fix longitude', ... 
    201233            'ptime_dbd', 'flight computer time', ... 
    202234            'ptime_dbd_datenum', 'flight computer date', ... 
     
    227259             'delta_vy', ... 
    228260             'prior_seg_vy', ... 
     261             'lat', ... 
     262             'lon', ... 
     263             'gps_lat', ... 
     264             'gps_lon', ... 
    229265             'ptime_dbd', ... 
    230266             'ptime_dbd_datenum', ...