Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

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

Open
louwangzhiyuwhy opened this issue Oct 29, 2021 · 0 comments

Comments

@louwangzhiyuwhy
Copy link

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?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant