santandermetgroup / climate4r.climdex Goto Github PK

As a user of RClimDex (https://github.com/rodneychan-ec/RClimDex), it would be easier to shift to this package if there was one function call to calculate all 27 indices at once.

I am trying to calculate climate change indices for a specific crossing season (October to April) from 1981 to 2020 but I get results as months instead of the season year. For the given period 1981 to 2020, I expected the output to have 39 seasons layers but it returns 273 monthly layers. Is there a way for me to specify that a year is defined after 212 days(days between October and April no leap)?

Please find my code below and attached data

#Calculating the indices
ds1<-"/Volumes/Untitled 2/Daily kk/Daily2_pr19812020.nc"

ppt.eobs <- loadGridData(dataset = ds1,
var = "ppt",
years = 1982:2020,
season = c(10,11,12,1,2,3,4))

PRCPTOT.indices=climdexGrid('PRCPTOT',pr =ppt.eobs,input.arg.list = list(base.range=c(1982,2020),northern.hemisphere=FALSE),cal = "365_day")
Daily2_pr19812020.nc.zip

Hello!
According to the definition of the climate index "Rnnmm", nn is a user defined threshold. When I tried to calculate the "R1mm" by using "climdexGrid()", it told me error in match.arg(index.code...). How can I define the nn threshold?
The code is:

R1mm <- climdexGrid(pr = prdata, index.code = "R1mm")
Error in match.arg(index.code, choices = c("FD", "SU", "ID", "TR", "GSL", :
'arg' should be one of “FD”, “SU”, “ID”, “TR”, “GSL”, “TXx”, “TNx”, “TXn”, “TNn”, “TN10p”, “TX10p”, “TN90p”, “TX90p”, “WSDI”, “CSDI”, “DTR”, “Rx1day”, “Rx5day”, “SDII”, “R10mm”, “R20mm”, “Rnnmm”, “CDD”, “CWD”, “R95pTOT”, “R99pTOT”, “PRCPTOT”

Thank you in advance!

The RDA files in data and inst are binary files that could be considered to move as external package or merge into climate4R.datasets

I have used climdexGrid to calculate growth season length using Cordex data:

GSL.grid <- climdexGrid(tn = tasmin,tx=tasmax, index.code = "GSL")

It works perfectly, but I would like to know the date that the growth season is starting to be calculated (index when the mean temperature above 5 has been achieved for 6 days and the 'counting' of the GSL is initiated ) and the date that it ends. Is that something achievable with your package?

I was wondering if it is possible to change the 0 degree threshold for the calculation of frost days? In the calculation of certain frost events, particularly in relation to plant damage, other thresholds are commonly used. It would be great to have the option to set the threshold for the "FD" index. Muchas gracias por tu ayuda!

while installing this package with:

/usr/bin/R -e "devtools::install_github(c('SantanderMetGroup/loadeR.java', 'SantanderMetGroup/transformeR', 'SantanderMetGroup/climate4R.climdex','SantanderMetGroup/climate4R.UDG','SantanderMetGroup/loadeR', 'SantanderMetGroup/visualizeR','SantanderMetGroup/downscaleR', 'SantanderMetGroup/loadeR.2nc'))"

on debian 12 stable (which has R 4.2.2) we receive this warning:

Warning messages:
1: package ‘climdex.pcic’ is not available for this version of R

A version of this package for your version of R might be available elsewhere,
see the ideas at
https://cran.r-project.org/doc/manuals/r-patched/R-admin.html#Installing-packages 
2: In i.p(...) :
  installation of package ‘/tmp/RtmpHvShJw/file921406163fc/climate4R.climdex_0.2.3.tar.gz’ had non-zero exit status

the installation works fine on debian 11 oldstable (which has R 4.0.4)

Would be great if there were functions for performing trend analysis for ETCCDI indicators using the climate4R data structures so that multiple ensembles can be analyzed in one go.

Does climdexGrid not take climate4R station data type? I tried the following:

download.file("http://meteo.unican.es/work/loadeR/data/VALUE_ECA_86_v2.tar.gz", destfile = "VALUE_ECA_86_v2.tar.gz")
untar("VALUE_ECA_86_v2.tar.gz")
value <- "VALUE_ECA_86_v2"
example1 <- loadStationData(dataset = value, var="tmax", stationID = c("000234", "003946"), years = 1981:2000)
SU1<-climdexGrid(tx=example1, index.code="SU")

which gives the following error message:
Error in aperm.default(array3D, match(c("time", "lat", "lon"), dimNames)) : perm is of wrong length 3 (!= 2)