ideas-lab-nus / epwshiftr

Create future EnergyPlus Weather files using CMIP6 data
https://ideas-lab-nus.github.io/epwshiftr/
Other
30 stars 7 forks source link

Can't find MIROC6 model results #17

Closed dvictori closed 3 years ago

dvictori commented 3 years ago

I'm trying to download daily data from the MIROC6 model. I've checked that the result in question is available in the ESGF system. Looking at the idx response attributes, the shards looks funny.

Is there a limit as to which model we can search?


library(epwshiftr)

idx <- init_cmip6_index(
    # only consider ScenarioMIP activity
    activity = "ScenarioMIP",
    # specify variables
    variable = c("pr"),
    # specify report frequent
    frequency = "day",
    # specify experiment name
    experiment = c("ssp245"),
    # specify GCM name
    source = c("MIROC6"),
    # specify variant,
    variant = "r1i1p1f1"
    )
#> No matched data. Please examine the actual response using 'attr(x, "response")'.

# This query should return one result
# looking at the response attribute
attr(idx, 'response')
#> $responseHeader
#> $responseHeader$status
#> [1] 0
#> 
#> $responseHeader$QTime
#> [1] 28
#> 
#> $responseHeader$params
#> $responseHeader$params$df
#> [1] "text"
#> 
#> $responseHeader$params$q.alt
#> [1] "*:*"
#> 
#> $responseHeader$params$indent
#> [1] "true"
#> 
#> $responseHeader$params$echoParams
#> [1] "all"
#> 
#> $responseHeader$params$fl
#> [1] "*,score"
#> 
#> $responseHeader$params$start
#> [1] "0"
#> 
#> $responseHeader$params$fq
#> $responseHeader$params$fq[[1]]
#> [1] "type:Dataset"
#> 
#> $responseHeader$params$fq[[2]]
#> [1] "project:\"CMIP6\""
#> 
#> $responseHeader$params$fq[[3]]
#> [1] "activity_id:\"ScenarioMIP\""
#> 
#> $responseHeader$params$fq[[4]]
#> [1] "experiment_id:\"ssp245\""
#> 
#> $responseHeader$params$fq[[5]]
#> [1] "source_id:\"MIROC6\""
#> 
#> $responseHeader$params$fq[[6]]
#> [1] "variable_id:\"pr\""
#> 
#> $responseHeader$params$fq[[7]]
#> [1] "nominal_resolution:\"100km\" || nominal_resolution:\"50km\" || nominal_resolution:\"100 km\" || nominal_resolution:\"50 km\""
#> 
#> $responseHeader$params$fq[[8]]
#> [1] "variant_label:\"r1i1p1f1\""
#> 
#> $responseHeader$params$fq[[9]]
#> [1] "frequency:\"day\""
#> 
#> $responseHeader$params$fq[[10]]
#> [1] "replica:false"
#> 
#> $responseHeader$params$fq[[11]]
#> [1] "latest:true"
#> 
#> 
#> $responseHeader$params$rows
#> [1] "10000"
#> 
#> $responseHeader$params$q
#> [1] "*:*"
#> 
#> $responseHeader$params$shards
#> [1] "localhost:8983/solr/datasets,localhost:8985/solr/datasets,localhost:8987/solr/datasets,localhost:8988/solr/datasets,localhost:8990/solr/datasets,localhost:8993/solr/datasets,localhost:8994/solr/datasets,localhost:8995/solr/datasets,localhost:8996/solr/datasets,localhost:8997/solr/datasets"
#> 
#> $responseHeader$params$tie
#> [1] "0.01"
#> 
#> $responseHeader$params$facet.limit
#> [1] "2048"
#> 
#> $responseHeader$params$qf
#> [1] "text"
#> 
#> $responseHeader$params$facet.method
#> [1] "fc"
#> 
#> $responseHeader$params$facet.mincount
#> [1] "1"
#> 
#> $responseHeader$params$wt
#> [1] "json"
#> 
#> $responseHeader$params$facet.sort
#> [1] "lex"
#> 
#> 
#> 
#> $response
#> $response$numFound
#> [1] 0
#> 
#> $response$start
#> [1] 0
#> 
#> $response$maxScore
#> [1] 0
#> 
#> $response$docs
#> list()

Created on 2021-08-12 by the reprex package (v2.0.0)

This is the expected result

image

dvictori commented 3 years ago

It turns out that init_cmip6_index default setting looks only for models with 100km spatial resolution. However, MIROC6 model has 250km resolution. By explicitly setting resolution = NULL in the function parameters, it worked.

hongyuanjia commented 3 years ago

That's what I exactly tried to point out. Glad that you find it by yourself.