Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

Points precipitation to grid example

noformat
Code Block
xml
xml
<?xml version="1.0" encoding="UTF-8"?>
<pcrTransformationSets xmlns="http://www.wldelft.nl/fews" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://www.wldelft.nl/fews
http://fews.wldelft.nl/schemas/version1.0/pcrTransformationSets.xsd" version="1.1">
	<logLevel>WARN</logLevel>
	<pcrTransformationSet id="Thiessen">
		<areaMap>
			<locationId>FineGrid</locationId>
		</areaMap>
		<definitions>
			<dataExchange>memory</dataExchange>
			<inputVariable variableId="P" dataType="scalar" convertDatum="false">
				<timeSeriesSet>
					<moduleInstanceId>ImportPubRts</moduleInstanceId>
					<valueType>scalar</valueType>
					<parameterId>P.obs</parameterId>
					<locationSetId>MetGauges_P.obs</locationSetId>
					<timeSeriesType>external historical</timeSeriesType>
					<timeStep unit="minute" multiplier="10"/>
					<relativeViewPeriod unit="hour" start="-96" startOverrulable="true" end="0"/>
					<readWriteMode>add originals</readWriteMode>
				</timeSeriesSet>
			</inputVariable>
			<outputVariable variableId="MeasMap" dataType="scalar" convertDatum="false">
				<timeSeriesSet>
					<moduleInstanceId>PrecipitationGaugeToGrid_Historical</moduleInstanceId>
					<valueType>grid</valueType>
					<parameterId>P.obs</parameterId>
					<locationId>FineGrid</locationId>
					<timeSeriesType>simulated forecasting</timeSeriesType>
					<timeStep unit="minute" multiplier="10"/>
					<relativeViewPeriod unit="hour" start="-96" end="0" startOverrulable="true"/>
					<readWriteMode>add originals</readWriteMode>
				</timeSeriesSet>
			</outputVariable>
		</definitions>
		<pcrModel id="String">
			<text><![CDATA[#! --unittrue --degrees
dynamic
# Simple Thiesen polygons to get spatial average precipitation on a grid

# Creat unique Id's for input stations
Unq = uniqueid(boolean(P));
# Now generate polygons and fill those
GaugeArea = spreadzone(ordinal(cover(Unq,0)),0,1);
MeasMap = areaaverage(P,GaugeArea);	
]]></text>
		</pcrModel>
	</pcrTransformationSet>
</pcrTransformationSets>