[
https://issues.apache.org/jira/browse/CLIMATE-236?page=com.atlassian.jira.plugin.system.issuetabpanels:comment-tabpanel&focusedCommentId=13725689#comment-13725689
]
Paul Ramirez commented on CLIMATE-236:
--------------------------------------
Instead of "private" what about just have it take in the datasets and return a
dataset. That way this utility function is available to others.
> Add Dataset Ensemble Support to the ocw.dataset_processor module
> ----------------------------------------------------------------
>
> Key: CLIMATE-236
> URL: https://issues.apache.org/jira/browse/CLIMATE-236
> Project: Apache Open Climate Workbench
> Issue Type: Sub-task
> Components: rcmet
> Affects Versions: 0.1-incubating
> Environment: *nix
> Reporter: Cameron Goodale
> Assignee: Cameron Goodale
> Fix For: 0.2-incubating
>
> Original Estimate: 4h
> Remaining Estimate: 4h
>
> Currently within the rcmes code base there is an ability to take in multiple
> datasets and return an 'ensemble' dataset. This ensemble is just the mean of
> all the input datasets.
> The plan is to grab this code block (around line 250) from do_data_prep:
> {code}
> # TODO: Refactor this into a function within the toolkit module
> # compute the simple multi-obs ensemble if multiple obs are used
> if numOBSs > 1:
> print 'numOBSs = ', numOBSs
> oData = obsData
> print 'oData shape = ', oData.shape
> obsData = ma.zeros((numOBSs + 1, nT, ngrdY, ngrdX))
> print 'obsData shape = ', obsData.shape
> avg = ma.zeros((nT, ngrdY, ngrdX))
>
> for i in np.arange(numOBSs):
> obsData[i, :, :, :] = oData[i, :, :, :]
> avg[:, :, :] = avg[:, :, :] + oData[i, :, :, :]
> avg = avg / float(numOBSs)
> obsData[numOBSs, :, :, :] = avg[:, :, :] # store the
> model-ensemble data
> numOBSs = numOBSs + 1 # update the number of obs
> data to include the model ensemble
> obsList.append('ENS-OBS')
> {code}
> Port all of that into a private function in dataset_processor called
> _rcmes_make_dataset_ensemble() and setup the dataset_processor.ensemble()
> function use it initially.
> Once the code move, documentation, unit tests are all complete, then I will
> resolve this issue.
--
This message is automatically generated by JIRA.
If you think it was sent incorrectly, please contact your JIRA administrators
For more information on JIRA, see: http://www.atlassian.com/software/jira
