Skip to content

BUG: cdo processed data read by loadGridData are all NaN, how to solve this? #75

@louwangzhiyuwhy

Description

@louwangzhiyuwhy

ERA-interim surface variables and pressure level variables are different grid so if I want merge the grid into one file:

cdo -P 8 remapbil,/home/inspur/working/climate4r/ERA-I/pressure_1979_2018_00.nc surface_1979_2018_00_new.nc surface_1979_2018_00_new_bilinear.nc

cdo merge pressure_1979_2018_00.nc surface_1979_2018_00_new_bilinear.nc merget_pressure_surface.nc

and then I read the merged file with loadR:
print(varnames[[z]])
xp <- lapply(varnames[[z]], function(x) {
loadGridData(dataset = "/home/inspur/working/climate4r/ERA-I/merget_pressure_surface.nc",
var = x,
years = 1979:2018)
}) %>% makeMultiGrid()

the str(xp) is all NaN:
[1] "uas" "vas"
[2021-10-29 16:11:07] Defining geo-location parameters
[2021-10-29 16:11:07] Defining time selection parameters
[2021-10-29 16:11:07] Retrieving data subset ...
[2021-10-29 16:11:14] Done
[2021-10-29 16:11:14] Defining geo-location parameters
[2021-10-29 16:11:14] Defining time selection parameters
[2021-10-29 16:11:15] Retrieving data subset ...
[2021-10-29 16:11:21] Done
List of 4
$ Variable:List of 2
..$ varName: chr [1:2] "uas" "vas"
..$ level : logi [1:2] NA NA
..- attr(, "use_dictionary")= chr [1:2] "FALSE" "FALSE"
..- attr(
, "units")= chr [1:2] "" ""
..- attr(, "longname")= chr [1:2] "uas" "vas"
..- attr(
, "daily_agg_cellfun")= chr [1:2] "none" "none"
..- attr(, "monthly_agg_cellfun")= chr [1:2] "none" "none"
..- attr(
, "verification_time")= chr [1:2] "none" "none"
$ Data : num [1:2, 1, 1:14610, 1:54, 1:41] NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN ...
..- attr(, "dimensions")= chr [1:5] "var" "member" "time" "lat" ...
$ xyCoords:List of 2
..$ x: num [1:41] 100 101 102 102 103 ...
..$ y: num [1:54] 10 10.8 11.5 12.2 13 ...
..- attr(
, "projection")= chr "LatLonProjection"
..- attr(, "resX")= num 0.75
..- attr(
, "resY")= num 0.75
$ Dates :List of 2
..$ :List of 2
.. ..$ start: chr [1:14610] "1979-01-01 00:00:00 GMT" "1979-01-02 00:00:00 GMT" "1979-01-03 00:00:00 GMT" "1979-01-04 00:00:00 GMT" ...
.. ..$ end : chr [1:14610] "1979-01-01 00:00:00 GMT" "1979-01-02 00:00:00 GMT" "1979-01-03 00:00:00 GMT" "1979-01-04 00:00:00 GMT" ...
..$ :List of 2
.. ..$ start: chr [1:14610] "1979-01-01 00:00:00 GMT" "1979-01-02 00:00:00 GMT" "1979-01-03 00:00:00 GMT" "1979-01-04 00:00:00 GMT" ...
.. ..$ end : chr [1:14610] "1979-01-01 00:00:00 GMT" "1979-01-02 00:00:00 GMT" "1979-01-03 00:00:00 GMT" "1979-01-04 00:00:00 GMT" ...

How to walkaroud this bug?

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions