Home > atmlab > retrieval > spareice > CollocatedNNIWP.m

CollocatedNNIWP

PURPOSE ^

SYNOPSIS ^

This is a script file.

DESCRIPTION ^

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SUBFUNCTIONS ^

DOWNLOAD ^

CollocatedNNIWP.m

SOURCE CODE ^

0001 classdef CollocatedNNIWP < NNTrainedProduct & InstrumentVisualiser & CollocationRetrievalDatabaseProduct
0002     
0003     % To be documented.
0004     %
0005     % See also: NNTrainedProduct, InstrumentVisualiser, CollocationRetrievalDatabaseProduct
0006 
0007     properties
0008        maxsize = 200000;
0009        
0010        inputs = struct();
0011 %             'MEAN_AVHRR_Y', struct(...
0012 %                 'transform', @(x)x, ...
0013 %                 'invtransform', @(x)x, ...
0014 %                 'lims', [100, 400], ...
0015 %                 'chans', [4, 5]));
0016 %                'chans', [1, 2, 4, 5]), ...
0017 %            'MEAN_ECMWF_Skin_temperature', struct(...
0018 %                'transform', @(x)x, ...
0019 %                'invtransform', @(x)x, ...
0020 %                'lims', [100 400]), ...
0021 %            'B_SZA', struct(...
0022 %                'transform', @sind, ...
0023 %                'invtransform', @asind, ...
0024 %                'lims', [0 80]));
0025 %        inputs = struct(...
0026 %            'B_BT', struct(...
0027 %                 'transform', @(x)x, ...
0028 %                 'invtransform', @(x)x, ...
0029 %                 'lims', [100 400], ...
0030 %                 'chans', [3, 4, 5]));
0031             
0032        targets = struct(...
0033             'MEAN_IWP_2C', struct(...
0034                 'transform', @log10, ...
0035                 'invtransform', @(x)10.^x, ...
0036                 'classify', 10));
0037                  %'lims', [realmin realmax], ...
0038            
0039         dbname = 'Collocated_cacmah_avhrr';
0040         date1 = [2007, 1, 1];
0041         date2 = [2008, 1, 0];
0042         sat1 = 'noaa18';
0043         sat2 = 'noaa18';
0044         fields = {'POS2', 'LAT1', 'LON1', 'TIME1', 'LAT2', 'LON2', ...
0045                   'MHS_LAT', 'MHS_LON', 'A_LAT', 'A_LONG', ...
0046                   'H_LAT', 'H_LONG', 'MEAN_ROIWP', 'MEAN_IOROIWP', ...
0047                   'MEAN_IWP_2C', 'MEAN_dardar_IWP', ...
0048                   'MEAN_ECMWF_Skin_temperature', 'FRAC40_dardar_IWP', ...
0049                   'NO_ROIWP', 'NO_IOROIWP', 'NO_IWP_2C', ...
0050                   'NO_dardar_IWP', 'MEAN_AVHRR_Y', ...
0051                   'MEAN_AVHRR_FLAG_3AB', 'H_BT', 'B_BT', ...
0052                   'NO_AVHRR_Y', 'B_SZA', 'B_LZA', 'B_LAA', 'B_SAA'};
0053 
0054         collocation_limits = struct();
0055 %         struct('B_BT', [100 400], ...
0056 %             'MEAN_AVHRR_Y', [0 400], ...
0057 %             'LAT1', [-40, 40], ...
0058 %             'B_SZA', [0 80]);
0059 
0060         collocation_filters = {};
0061 %         ...
0062 %             {@(avhrr, surf) avhrr(:, 5) - surf < -20, ...
0063 %              {'MEAN_AVHRR_Y', 'MEAN_ECMWF_Skin_temperature'}}, ...
0064 %             {@(avhrr) avhrr(:, 5) - avhrr(:, 4) > -2, ...
0065 %              {'MEAN_AVHRR_Y'}}, ...
0066 %             };
0067     end 
0068     
0069     properties (Constant)
0070        chans = [1, 2, 4, 5];
0071        freqs = {'0.58-0.68', '0.725-1.0', '3.55-3.93', '10.3-11.3', '11.5-12.5'};
0072         
0073     end
0074 
0075     methods
0076         function self = CollocatedNNIWP(varargin)
0077             self = self@NNTrainedProduct(varargin{:});
0078         end
0079         
0080         function cc = fix_cols_in_postprocessing(self, cc)
0081             % when reprocessing, don't want to recollocate but use stored
0082             % fields.  However, they may have different names in the
0083             % cols-structure.  Re-arrange so that this can be used.
0084             if isfield(cc, 'MHS')
0085                 %cc.B_BT = [0, 0, cc.MHS];
0086                 cc.B_BT = cc.MHS;
0087             end
0088             cc.LAT1 = cc.LAT;
0089             cc.LON1 = cc.LON;
0090             cc.TIME1 = cc.TIME;
0091         end
0092     end
0093 
0094 end
0095 %[M, c] = D.Collocated_cacmah_avhrr.read([2007, 1, 1], [2007, 1, 3], {'noaa18', 'noaa18'}, {'POS2', 'LAT1', 'LON1', 'LAT2', 'LON2', 'MHS_LAT', 'MHS_LON', 'A_LAT', 'A_LONG', 'H_LAT', 'H_LONG', 'MEAN_ROIWP', 'MEAN_IOROIWP', 'MEAN_IWP_2C', 'MEAN_dardar_IWP', 'MEAN_ECMWF_Skin_temperature', 'FRAC100_dardar_IWP', 'NO_ROIWP', 'NO_IOROIWP', 'NO_IWP_2C', 'NO_dardar_IWP', 'MEAN_AVHRR_Y', 'MEAN_AVHRR_FLAG_3AB', 'H_BT', 'B_BT', 'NO_AVHRR_Y', 'B_SZA', 'B_LZA', 'B_LAA'}, struct('B_BT', [100 400], 'MEAN_AVHRR_Y', [0 400], 'LAT1', [-40, 40]));
0096

Generated on Mon 15-Sep-2014 13:31:28 by m2html © 2005