Changes between Version 13 and Version 14 of udg/ecoms/RPackage/examples/continentalSelection

Timestamp:

Sep 2, 2014 5:27:37 PM (7 years ago)

Author:

juaco

Comment:

--

udg/ecoms/RPackage/examples/continentalSelection

@@ -65 +65 @@
 #!text/R
 > obs.regridded <- interpGridData(gridData = ex2.obs, new.grid = getGrid(ex2), method = "bilinear")
+[2014-09-02 17:22:11] Performing bilinear interpolation... may take a while
+[2014-09-02 17:22:30] Done
+Warning messages:
+1: In interpGridData(gridData = ex2.obs, new.grid = getGrid(ex2), method = "bilinear") :
+  The new longitudes are outside the data extent
+2: In interpGridData(gridData = ex2.obs, new.grid = getGrid(ex2), method = "bilinear") :
+  The new latitudes are outside the data extent
+}}}
+
+Note the warnings reminding us that the extent of the input grid is wider that that of CFS. However, in this case we can safely ignore this warnings, since all the land areas we are interest in are within the CFS domain.
+
+{{{
+#!text/R
 > plotMeanField(obs.regridded)
 }}}
 
+[[Image(image-20140902-172548.png)]]
 
 
+Now that both model data and observations are in the same grid, we can compute the bias. First, we calculate the spatial mean of WFDEI, which is the reference against which to compute the biases:
+
+{{{
+#!text/R
+ref <- apply(obs.regridded$Data, MARGIN = c(3,2), mean, na.rm = TRUE)
+}}}
 
 
-
-