Last modified: Tue, 03/14/2017 - 00:19
If you have an ensemble constructed such that the Java CDM library can understand the grid incluing the ensemble axis, it's very easy to use the automated tools to include the data set in your LAS server.
If you have carefully constructed your aggregation as described, the the process of configuring the result into LAS is as simple as with any other netCDF data source. You just run the addXML.sh command with the ensemble aggreation as the target data set.
[rhs@dunkel ~/baker]$ ./bin/addXML.sh -n http://dunkel.pmel.noaa.gov:8930/thredds/dodsC/data/atmos3.ncml
Axis has no units. Setting units to 'none'. Check output.
You will receive a warning because the ensemble axis has no units attribute.
With that you should have an output file (called las_from_addXML.xml by default) which you can use to configure your LAS with the ensemble data set. The configuration for our example data set is shown below.
<datasets>
<id-1b5c4759a1 name="model output prepared for GFDL Seasonal-Interannual Experimental Forecasts CM2.1U_CDAef_v1.0" url="http://dunkel.pmel.noaa.gov:8930/thredds/dodsC/data/atmos3.ncml">
<variables>
<zg-id-1b5c4759a1 name="Geopotential Height" units="m" url="#zg">
<link match="/lasdata/grids/grid-lon-lat-plev-time-ensemble-id-1b5c4759a1" />
</zg-id-1b5c4759a1>
<tas-id-1b5c4759a1 name="Surface Air Temperature" units="K" url="#tas">
<link match="/lasdata/grids/grid-lon-lat-time-ensemble-id-1b5c4759a1" />
</tas-id-1b5c4759a1>
</variables>
</id-1b5c4759a1>
</datasets>
<grids>
<grid-lon-lat-time-ensemble-id-1b5c4759a1>
<link match="/lasdata/axes/lon-x-id-1b5c4759a1" />
<link match="/lasdata/axes/lat-y-id-1b5c4759a1" />
<link match="/lasdata/axes/time-t-id-1b5c4759a1" />
<link match="/lasdata/axes/ensemble-e-id-1b5c4759a1" />
</grid-lon-lat-time-ensemble-id-1b5c4759a1>
<grid-lon-lat-plev-time-ensemble-id-1b5c4759a1>
<link match="/lasdata/axes/lon-x-id-1b5c4759a1" />
<link match="/lasdata/axes/lat-y-id-1b5c4759a1" />
<link match="/lasdata/axes/plev-z-id-1b5c4759a1" />
<link match="/lasdata/axes/time-t-id-1b5c4759a1" />
<link match="/lasdata/axes/ensemble-e-id-1b5c4759a1" />
</grid-lon-lat-plev-time-ensemble-id-1b5c4759a1>
</grids>
<axes>
<lon-x-id-1b5c4759a1 type="x" units="degrees_east">
<arange start="1.25" size="144" step="2.5" />
</lon-x-id-1b5c4759a1>
<lat-y-id-1b5c4759a1 type="y" units="degrees_north">
<v>-89.4944</v>
<v>-87.9775</v>
<v>-85.9551</v>
<v>-83.9326</v>
<v>-81.9101</v>
<v>-79.8876</v>
<v>-77.8652</v>
<v>-75.8427</v>
<v>-73.8202</v>
<v>-71.7978</v>
<v>-69.7753</v>
<v>-67.7528</v>
<v>-65.7303</v>
<v>-63.7079</v>
<v>-61.6854</v>
<v>-59.6629</v>
<v>-57.6404</v>
<v>-55.618</v>
<v>-53.5955</v>
<v>-51.573</v>
<v>-49.5506</v>
<v>-47.5281</v>
<v>-45.5056</v>
<v>-43.4831</v>
<v>-41.4607</v>
<v>-39.4382</v>
<v>-37.4157</v>
<v>-35.3933</v>
<v>-33.3708</v>
<v>-31.3483</v>
<v>-29.3258</v>
<v>-27.3034</v>
<v>-25.2809</v>
<v>-23.2584</v>
<v>-21.236</v>
<v>-19.2135</v>
<v>-17.191</v>
<v>-15.1685</v>
<v>-13.1461</v>
<v>-11.1236</v>
<v>-9.1011</v>
<v>-7.0787</v>
<v>-5.0562</v>
<v>-3.0337</v>
<v>-1.0112</v>
<v>1.0112</v>
<v>3.0337</v>
<v>5.0562</v>
<v>7.0787</v>
<v>9.1011</v>
<v>11.1236</v>
<v>13.1461</v>
<v>15.1685</v>
<v>17.191</v>
<v>19.2135</v>
<v>21.236</v>
<v>23.2584</v>
<v>25.2809</v>
<v>27.3034</v>
<v>29.3258</v>
<v>31.3483</v>
<v>33.3708</v>
<v>35.3933</v>
<v>37.4157</v>
<v>39.4382</v>
<v>41.4607</v>
<v>43.4831</v>
<v>45.5056</v>
<v>47.5281</v>
<v>49.5506</v>
<v>51.573</v>
<v>53.5955</v>
<v>55.618</v>
<v>57.6404</v>
<v>59.6629</v>
<v>61.6854</v>
<v>63.7079</v>
<v>65.7303</v>
<v>67.7528</v>
<v>69.7753</v>
<v>71.7978</v>
<v>73.8202</v>
<v>75.8427</v>
<v>77.8652</v>
<v>79.8876</v>
<v>81.9101</v>
<v>83.9326</v>
<v>85.9551</v>
<v>87.9775</v>
<v>89.4944</v>
</lat-y-id-1b5c4759a1>
<time-t-id-1b5c4759a1 type="t" units="month">
<arange start="1982-04-16" size="36" step="1" />
</time-t-id-1b5c4759a1>
<ensemble-e-id-1b5c4759a1 label="Realizations" type="e" units="none">
<v label="Realization01">1</v>
<v label="Realization02">2</v>
<v label="Realization03">3</v>
</ensemble-e-id-1b5c4759a1>
<plev-z-id-1b5c4759a1 type="z" units="Pa">
<v>100000</v>
<v>92500</v>
<v>85000</v>
<v>70000</v>
<v>60000</v>
<v>50000</v>
<v>40000</v>
<v>30000</v>
<v>25000</v>
<v>20000</v>
<v>15000</v>
<v>10000</v>
<v>7000</v>
<v>5000</v>
<v>3000</v>
<v>2000</v>
<v>1000</v>
</plev-z-id-1b5c4759a1>
</axes>
Now we get to the fun part of using the ensemble in LAS.