Describe the dataset issue
I've recently found an issue in the GFDL-CM4 surface wind (sfcwind) height coordinate. When I was preprocessing the data for analysis, I found that I kept getting an error saying:
ValueError: Can not concatenate cubes: failed to concatenate into a single cube.
Scalar coordinates values or metadata differ: height != height
I was checking the cmor fixes in the ESMValCore directories and it seems there was a height fix where 2m of height was added. From my understanding sfcWind uses 10m height which could maybe causing this issue?
Project (CMIP5/CMIP6/OBS/obs4MIPs/etc): CMIP6
Full dataset description (fill out as many as you know, please):
dataset: GFDL-CM4
experiment: historical and ssp245
mip: day
grid: gr1
version: 20180701
variable used: sfcWind
The full issue is a little hard to read because I was using dask but this is the last call that I can see before the big error message at the bottom of main_log_debug.txt
2024-03-21 21:54:45,701 UTC [1957677] ERROR esmvalcore.preprocessor._io:326 Can not concatenate cubes into a single one: failed to concatenate into a single cube.
Scalar coordinates values or metadata differ: height != height
2024-03-21 21:54:45,701 UTC [1957677] ERROR esmvalcore.preprocessor._io:327 Resulting cubes:
2024-03-21 21:54:45,701 UTC [1957677] ERROR esmvalcore.preprocessor._io:329 wind_speed / (m s-1) (time: 21900; latitude: 180; longitude: 288)
Dimension coordinates:
time x - -
latitude - x -
longitude - - x
Scalar coordinates:
height 10.0 m
Cell methods:
0 area: time: mean
Attributes:
Conventions 'CF-1.7 CMIP-6.0 UGRID-1.0'
activity_id 'CMIP ScenarioMIP'
branch_method 'standard'
branch_time_in_child 0.0
branch_time_in_parent '-638385 60225.0'
contact 'gfdl.climate.model.info@noaa.gov'
data_specs_version '01.00.27'
experiment 'historical update of RCP4.5 based on SSP2'
experiment_id 'historical ssp245'
external_variables 'areacella'
forcing_index 1
frequency 'day'
further_info_url 'https://furtherinfo.es-doc.org/CMIP6.NOAA-GFDL.GFDL-CM4.historical.none.r1i1p1f1 ...'
grid 'atmos data regridded from Cubed-sphere (c96) to 180,288; interpolation ...'
grid_label 'gr1'
initialization_index 1
institution 'National Oceanic and Atmospheric Administration, Geophysical Fluid Dynamics ...'
institution_id 'NOAA-GFDL'
interp_method 'conserve_order2'
license 'CMIP6 model data produced by NOAA-GFDL is licensed under a Creative Commons ...'
mip_era 'CMIP6'
nominal_resolution '100 km'
original_name 'sfcWind'
parent_activity_id 'CMIP'
parent_experiment_id 'historical piControl'
parent_mip_era 'CMIP6'
parent_source_id 'GFDL-CM4'
parent_time_units 'days since 1850-1-1 00:00:00'
parent_variant_label 'r1i1p1f1'
physics_index 1
product 'model-output'
realization_index 1
realm 'atmos'
references 'see further_info_url attribute'
source 'GFDL-CM4 (2018): \naerosol: interactive\natmos: GFDL-AM4.0.1 (Cubed-sphere ...'
source_id 'GFDL-CM4'
source_type 'AOGCM'
sub_experiment 'none'
sub_experiment_id 'none'
table_id 'day'
title 'NOAA GFDL GFDL-CM4 model output prepared for CMIP6 historical NOAA GFDL ...'
variable_id 'sfcWind'
variant_info 'N/A'
variant_label 'r1i1p1f1'
2024-03-21 21:54:45,705 UTC [1957677] ERROR esmvalcore.preprocessor._io:331 From Cell(point=cftime.DatetimeNoLeap(1950, 1, 1, 12, 0, 0, 0, has_year_zero=True), bound=(cftime.DatetimeNoLeap(1950, 1, 1, 0, 0, 0, 0, has_year_zero=True), cftime.DatetimeNoLeap(1950, 1, 2, 0, 0, 0, 0, has_year_zero=True))) to Cell(point=cftime.DatetimeNoLeap(2009, 12, 31, 12, 0, 0, 0, has_year_zero=True), bound=(cftime.DatetimeNoLeap(2009, 12, 31, 0, 0, 0, 0, has_year_zero=True), cftime.DatetimeNoLeap(2010, 1, 1, 0, 0, 0, 0, has_year_zero=True)))
Describe the dataset issue I've recently found an issue in the GFDL-CM4 surface wind (sfcwind) height coordinate. When I was preprocessing the data for analysis, I found that I kept getting an error saying:
I was checking the cmor fixes in the ESMValCore directories and it seems there was a height fix where 2m of height was added. From my understanding sfcWind uses 10m height which could maybe causing this issue?
The full issue is a little hard to read because I was using dask but this is the last call that I can see before the big error message at the bottom of main_log_debug.txt
Cheers :beer: