Error in ncdf4::nc_create(paste(out.dir, cdx.dat$filename[x], sep = "/")

[Lotus-K]

Hello, when I use CMIP6 Data(pr_day_CESM2_historical_r4i1p1f1_gn_18500101-18591231.nc), there shows a problem

Error in ncdf4::nc_create(paste(out.dir, cdx.dat$filename[x], sep = "/"), : Error, found an element of the vars list that is NOT an object created by a call to ncvar_def...element #3 Calls: create.indices.from.files ... create.ncdf.output.files -> lapply -> FUN -> Execution halted

Can you please help me?

about the data(pr_day_CESM2_historical_r4i1p1f1_gn_18500101-18591231.nc): dimensions: lat = 192 ; nbnd = 2 ; lon = 288 ; time = UNLIMITED ; // (3650 currently) variables: float pr(time, lat, lon) ; pr:_FillValue = 1.e+20f ; pr:cell_measures = "area: areacella" ; pr:cell_methods = "area: time: mean" ; pr:comment = "at surface; includes both liquid and solid phases from all types of clouds (both large-scale and convective)" ; pr:coordinates = "time lat lon" ; pr:description = "at surface; includes both liquid and solid phases from all types of clouds (both large-scale and convective)" ; pr:frequency = "day" ; pr:id = "pr" ; pr:long_name = "Precipitation" ; pr:mipTable = "day" ; pr:missing_value = 1.e+20 ; pr:out_name = "pr" ; pr:prov = "day ((isd.003))" ; pr:realm = "atmos" ; pr:standard_name = "precipitation_flux" ; pr:time = "time" ; pr:time_label = "time-mean" ; pr:time_title = "Temporal mean" ; pr:title = "Precipitation" ; pr:type = "real" ; pr:units = "kg m-2 s-1" ; pr:variable_id = "pr" ; double lat(lat) ; lat:axis = "Y" ; lat:bounds = "lat_bnds" ; lat:standard_name = "latitude" ; lat:title = "Latitude" ; lat:type = "double" ; lat:units = "degrees_north" ; lat:valid_max = 90. ; lat:valid_min = -90. ; double lon(lon) ; lon:axis = "X" ; lon:bounds = "lon_bnds" ; lon:standard_name = "longitude" ; lon:title = "Longitude" ; lon:type = "double" ; lon:units = "degrees_east" ; lon:valid_max = 360. ; lon:valid_min = 0. ; double time(time) ; time:axis = "T" ; time:bounds = "time_bnds" ; time:standard_name = "time" ; time:title = "time" ; time:type = "double" ; time:units = "days since 0001-01-01 00:00:00" ; time:calendar = "noleap" ; double time_bnds(time, nbnd) ; time_bnds:calendar = "noleap" ; time_bnds:units = "days since 0001-01-01 00:00:00" ; double lat_bnds(lat, nbnd) ; lat_bnds:units = "degrees_north" ; double lon_bnds(lon, nbnd) ; lon_bnds:units = "degrees_east" ;

// global attributes: :Conventions = "CF-1.7 CMIP-6.2" ; :activity_id = "CMIP" ; :branch_method = "standard" ; :branch_time_in_child = 674885. ; :branch_time_in_parent = 182500. ; :case_id = "18" ; :cesm_casename = "b.e21.BHIST.f09_g17.CMIP6-historical.004" ; :contact = "cesm_cmip6@ucar.edu" ; :creation_date = "2019-03-13T15:05:58Z" ; :data_specs_version = "01.00.29" ; :experiment = "all-forcing simulation of the recent past" ; :experiment_id = "historical" ; :external_variables = "areacella" ; :forcing_index = 1LL ; :frequency = "day" ; :further_info_url = "https://furtherinfo.es-doc.org/CMIP6.NCAR.CESM2.historical.none.r4i1p1f1" ; :grid = "native 0.9x1.25 finite volume grid (192x288 latxlon)" ; :grid_label = "gn" ; :initialization_index = 1LL ; :institution = "National Center for Atmospheric Research, Climate and Global Dynamics Laboratory, 1850 Table Mesa Drive, Boulder, CO 80305, USA" ; :institution_id = "NCAR" ; :license = "CMIP6 model data produced by is licensed under a Creative Commons Attribution-[]ShareAlike 4.0 International License (https://creativecommons.org/licenses/). Consult https://pcmdi.llnl.gov/CMIP6/TermsOfUse for terms of use governing CMIP6 output, including citation requirements and proper acknowledgment. Further information about this data, including some limitations, can be found via the further_info_url (recorded as a global attribute in this file)[]. The data producers and data providers make no warranty, either express or implied, including, but not limited to, warranties of merchantability and fitness for a particular purpose. All liabilities arising from the supply of the information (including any liability arising in negligence) are excluded to the fullest extent permitted by law." ; :mip_era = "CMIP6" ; :model_doi_url = "https://doi.org/10.5065/D67H1H0V" ; :nominal_resolution = "100 km" ; :parent_activity_id = "CMIP" ; :parent_experiment_id = "piControl" ; :parent_mip_era = "CMIP6" ; :parent_source_id = "CESM2" ; :parent_time_units = "days since 0001-01-01 00:00:00" ; :parent_variant_label = "r1i1p1f1" ; :physics_index = 1LL ; :product = "model-output" ; :realization_index = 4LL ; :realm = "atmos" ; :source = "CESM2 (2017): atmosphere: CAM6 (0.9x1.25 finite volume grid; 288 x 192 longitude/latitude; 32 levels; top level 2.25 mb); ocean: POP2 (320x384 longitude/latitude; 60 levels; top grid cell 0-10 m); sea_ice: CICE5.1 (same grid as ocean); land: CLM5 0.9x1.25 finite volume grid; 288 x 192 longitude/latitude; 32 levels; top level 2.25 mb); aerosol: MAM4 (0.9x1.25 finite volume grid; 288 x 192 longitude/latitude; 32 levels; top level 2.25 mb); atmoschem: MAM4 (0.9x1.25 finite volume grid; 288 x 192 longitude/latitude; 32 levels; top level 2.25 mb); landIce: CISM2.1; ocnBgchem: MARBL (320x384 longitude/latitude; 60 levels; top grid cell 0-10 m)" ; :source_id = "CESM2" ; :source_type = "AOGCM BGC" ; :sub_experiment = "none" ; :sub_experiment_id = "none" ; :table_id = "day" ; :tracking_id = "hdl:21.14100/46c92cff-a471-4f1e-aaf9-d2bda7c3a95e" ; :variable_id = "pr" ; :variant_info = "CMIP6 20th century experiments (1850-2014) with CAM6, interactive land (CLM5), coupled ocean (POP2) with biogeochemistry (MARBL), interactive sea ice (CICE5.1), and non-evolving land ice (CISM2.1)\r\n" ; :variant_label = "r4i1p1f1" ; }

[GeoFelpave]

Hi Lotus-K,

Have you had a response to the issue?. I am having the same issue when trying to run the tool with data taken from GPM.

Thanks

[Lotus-K]

Hi Lotus-K,

Have you had a response to the issue?. I am having the same issue when trying to run the tool with data taken from GPM.

Thanks

I did not receive a response. Climpact worked fine when calculated with the previous CMIP5 data(CCSM4). I am now trying to calculate the extreme indices using Python-xclim but the results are a little different from those at Climpact with CCSM4 data.

[GeoFelpave]

Hi Lotus-K, Have you had a response to the issue?. I am having the same issue when trying to run the tool with data taken from GPM. Thanks

I did not receive a response. Climpact worked fine when calculated with the previous CMIP5 data(CCSM4). I am now trying to calculate the extreme indices using Python-xclim but the results are a little different from those at Climpact with CCSM4 data.

Ok. Thanks for the reply. I will check Python-xclim. I was trying to run Climpact with data from GPM-IMERG but still not luck.

[heroldn]

I think someone emailed me and the issue was non-daily data(?). If this is still an issue please upload a couple of years of data and I'll take a look when I get time.