Skip to content

Dataset problem: Issues with hybrid pressure coordinate of GFDL-CM4 (CMIP6) #541

@schlunma

Description

@schlunma

Describe the dataset issue
The necessary attribute forumula_terms is not given (only formula_term). Bounds of ap and b not set (see #536).

  • Data file has been changed by you in any way
  • Project (CMIP5/CMIP6/OBS/obs4mips/etc): CMIP6
  • Full dataset description (fill out as many as you know, please):
    • dataset: GFDL-CM4
    • experiment: historical
    • mip: Amon
    • grid: gr1
    • version: v20180701
    • variable used: cl
  • Pointer to existing copy of data on ESGF node (it would be very useful if you could provide a physical fullpath to the file(s) that are causing the problem, e.g. on CEDA Jasmin or DKRZ):
    /mnt/lustre02/work/bd0854/DATA/ESMValTool2/CMIP6_DKRZ/CMIP/NOAA-GFDL/GFDL-CM4/historical/r1i1p1f1/Amon/cl/gr1/v20180701/cl_Amon_GFDL-CM4_historical_r1i1p1f1_gr1_195001-201412.nc

ncdump (without global attributes):

netcdf cl_Amon_GFDL-CM4_historical_r1i1p1f1_gr1_195001-201412 {                                                                                                                                                                                                       [108/967]
dimensions:
        lev = 33 ;
        bnds = 2 ;
        time = UNLIMITED ; // (780 currently)
        lat = 180 ;
        lon = 288 ;
variables:
        double ap(lev) ;
                ap:long_name = "vertical coordinate formula term: ap(k)" ;
                ap:units = "Pa" ;
                ap:missing_value = 1.e+20 ;
                ap:_FillValue = 1.e+20 ;
                ap:cell_methods = "time: point" ;
        double ap_bnds(lev, bnds) ;
                ap_bnds:long_name = "vertical coordinate formula term: ap(k+1/2)" ;
                ap_bnds:units = "Pa" ;
                ap_bnds:missing_value = 1.e+20 ;
                ap_bnds:_FillValue = 1.e+20 ;
                ap_bnds:cell_methods = "time: point" ;
        double b(lev) ;
                b:long_name = "vertical coordinate formula term: b(k)" ;
                b:units = "1.0" ;
                b:missing_value = 1.e+20 ;
                b:_FillValue = 1.e+20 ;
                b:cell_methods = "time: point" ;
        double b_bnds(lev, bnds) ;
                b_bnds:long_name = "vertical coordinate formula term: b(k+1/2)" ;
                b_bnds:units = "1.0" ;
                b_bnds:missing_value = 1.e+20 ;
                b_bnds:_FillValue = 1.e+20 ;
                b_bnds:cell_methods = "time: point" ;
        double bnds(bnds) ;
                bnds:long_name = "vertex number" ;
        float cl(time, lev, lat, lon) ;
                cl:long_name = "Cloud Area Fraction" ;
                cl:units = "%" ;
                cl:missing_value = 1.e+20f ;
                cl:_FillValue = 1.e+20f ;
                cl:cell_methods = "area: time: mean" ;
                cl:cell_measures = "area: areacella" ;
                cl:standard_name = "cloud_area_fraction_in_atmosphere_layer" ;
                cl:interp_method = "conserve_order1" ;
                cl:original_name = "cl" ;
        double lat(lat) ;
                lat:long_name = "latitude" ;
                lat:units = "degrees_north" ;
                lat:axis = "Y" ;
                lat:bounds = "lat_bnds" ;
                lat:standard_name = "latitude" ;
                lat:cell_methods = "time: point" ;
        double lat_bnds(lat, bnds) ;
                lat_bnds:long_name = "latitude bounds" ;
                lat_bnds:units = "degrees_north" ;
                lat_bnds:axis = "Y" ;
        double lev(lev) ;
                lev:long_name = "atmospheric model level" ;
                lev:units = "1.0" ;
                lev:axis = "Z" ;
                lev:positive = "down" ;
                lev:formula = "p(n,k,j,i) = ap(k) + b(k)*ps(n,j,i)" ;
                lev:bounds = "lev_bnds" ;
                lev:standard_name = "atmosphere_hybrid_sigma_pressure_coordinate" ;
                lev:description = "Generic atmospheric model vertical coordinate (nondimensional or dimensional) [use CF standard name apporpriate for model vertical coordinate, e.g. model_level_number or atmosphere_sigma_coordinate]" ;
                lev:formula_term = "ap: ap b: b ps: ps" ;
        double lev_bnds(lev, bnds) ;
                lev_bnds:long_name = "atmospheric model level bounds" ;
                lev_bnds:units = "1.0" ;
                lev_bnds:missing_value = 1.e+20 ;
                lev_bnds:_FillValue = 1.e+20 ;
                lev_bnds:cell_methods = "time: point" ;
                lev_bnds:formula = "p(n,k+1/2,j,i) = ap(k+1/2) + b(k+1/2)*ps(n,j,i)" ;
                lev_bnds:formula_terms = "ap: ap_bnds b: b_bnds ps: ps" ;
                lev_bnds:standard_name = "atmosphere_hybrid_sigma_pressure_coordinate" ;
                lev_bnds:axis = "Z" ;
                lev_bnds:positive = "down" ;
        double lon(lon) ;
                lon:long_name = "longitude" ;
                lon:units = "degrees_east" ;
                lon:axis = "X" ;
                lon:bounds = "lon_bnds" ;
                lon:standard_name = "longitude" ;
                lon:cell_methods = "time: point" ;
        double lon_bnds(lon, bnds) ;
                lon_bnds:long_name = "longitude bounds" ;
                lon_bnds:units = "degrees_east" ;
                lon_bnds:axis = "X" ;
        float ps(time, lat, lon) ;
                ps:long_name = "Surface Air Pressure" ;
                ps:units = "Pa" ;
                ps:missing_value = 1.e+20f ;
                ps:_FillValue = 1.e+20f ;
                ps:cell_methods = "time: mean" ;
                ps:cell_measures = "area: areacella" ;
                ps:time_avg_info = "average_T1,average_T2,average_DT" ;
                ps:standard_name = "surface_air_pressure" ;
                ps:interp_method = "conserve_order2" ;
        double time(time) ;
                time:long_name = "time" ;
                time:units = "days since 1850-01-01 00:00:00" ;
                time:axis = "T" ;
                time:calendar_type = "noleap" ;
                time:calendar = "noleap" ;
                time:bounds = "time_bnds" ;
                time:standard_name = "time" ;
                time:description = "Temporal mean" ;
        double time_bnds(time, bnds) ;
                time_bnds:long_name = "time axis boundaries" ;
                time_bnds:units = "days since 1850-01-01 00:00:00" ;

Fix provided in #529.

Metadata

Metadata

Assignees

Type

No type

Projects

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions