Changes between Version 5 and Version 6 of udg/ecoms/RPackage/examples/drift


Ignore:
Timestamp:
Sep 4, 2014 12:24:00 PM (7 years ago)
Author:
juaco
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • udg/ecoms/RPackage/examples/drift

    v5 v6  
    66
    77{{{#!text/R
    8 > # The selection of leadMonth = 0 will give a NOTE on screen:
    9 > ref <- loadECOMS(dataset = "CFSv2_seasonal_16", var = "tas", members = 1, lonLim = c(-10,-1), latLim = c(36,40), season = 7, years = 2006, leadMonth = 0, time = "DD")
     8ref <- loadECOMS(dataset = "CFSv2_seasonal_16", var = "tas", members = 1, lonLim = c(-10,-1), latLim = c(36,40), season = 7, years = 2006, leadMonth = 0, time = "DD")
     9}}}
     10
     11Messages on-screen inform about the loading process. Note that the selection of `leadMonth = 0` will give a NOTE on screen:
     12
     13{{{
    1014[2014-09-03 17:50:20] Defining homogeneization parameters for variable "tas"
    1115NOTE: daily mean will be calculated from the 6-h model output
     
    1519[2014-09-03 17:50:26] Retrieving data subset ...
    1620[2014-09-03 17:50:31] Done
    17 > plotMeanField(ref)
    18 > title(main = "Lead month 0 forecast of July 2001")
    19 > # This is the spatial mean of the reference field
    20 > ref.field <- apply(ref$Data, MARGIN = c(3,2), FUN = mean, na.rm = TRUE)
     21}}}
     22
     23{{{#!text/R
     24plotMeanField(ref)
     25title(main = "Lead month 0 forecast of July 2001")
     26# This is the spatial mean of the reference field
     27ref.field <- apply(ref$Data, MARGIN = c(3,2), FUN = mean, na.rm = TRUE)
    2128}}}
    2229
     
    2835
    2936{{{#!text/R
    30 > cfs.list <- lapply(1:6, function(lead.month) {
     37cfs.list <- lapply(1:6, function(lead.month) {
    3138      loadECOMS(dataset = "CFSv2_seasonal_16", var = "tas", members = 1, lonLim = c(-10,-1), latLim = c(36,40), season = 7, years = 2006, leadMonth = lead.month, time = "DD")
    3239})