NCAR / esds

Repository for hosting material related to the Earth System Data Science (ESDS) initiative efforts
https://ncar.github.io/esds
Apache License 2.0
24 stars 20 forks source link

Grid/PFT variable to lat lon variable conversion #222

Open mm3975 opened 1 year ago

mm3975 commented 1 year ago

To test a hypothesis, we want to use the University of Texas coupled high and low-resolution outputs for ET. They have saved the ET as column variables (maybe at grid level). Can anyone help us convert the column variables into lat lon variables? Is there a script or function to process column data? You may find the metadata details below:

netcdf evspsbl_Amon_CESM1-CAM5-SE-HR_hist-1950_r1i1p1f1_gn_195001-201412 {
dimensions:
        nbnd = 2 ;
        grid_corners = 5 ;
        ncol = 777602 ;
        time = UNLIMITED ; // (780 currently)
variables:
        float evspsbl(time, ncol) ;
                evspsbl:_FillValue = 1.e+20f ;
                evspsbl:cell_measures = "area: areacella" ;
                evspsbl:cell_methods = "area: time: mean" ;
                evspsbl:comment = "QFLX" ;
                evspsbl:coordinates = "time lat lon" ;
                evspsbl:description = "at surface; flux of water into the atmosphere due to conversion of both liquid and solid phases to vapor (from underlying surface and vegetation)" ;
                evspsbl:frequency = "mon" ;
                evspsbl:id = "evspsbl" ;
                evspsbl:long_name = "Evaporation Including Sublimation and Transpiration" ;
                evspsbl:mipTable = "Amon" ;
                evspsbl:missing_value = 1.e+20 ;
                evspsbl:out_name = "evspsbl" ;
                evspsbl:prov = "Amon ((isd.003))" ;
                evspsbl:realm = "atmos" ;
                evspsbl:standard_name = "water_evapotranspiration_flux" ;
                evspsbl:time = "time" ;
                evspsbl:time_label = "time-mean" ;
                evspsbl:time_title = "Temporal mean" ;
                evspsbl:title = "Evaporation Including Sublimation and Transpiration" ;
                evspsbl:type = "real" ;
                evspsbl:units = "kg m-2 s-1" ;
                evspsbl:variable_id = "evspsbl" ;
        double lat(ncol) ;
                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 lat(ncol) ;
                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(ncol) ;
                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. ;
katiedagon commented 1 year ago

@mm3975 sorry for the very long delay responding to this issue! We have a couple of recent ESDS blog posts that cover working with unstructured model output (i.e., column level) and remapping to latitude x longitude. Hopefully these might help: https://ncar.github.io/esds/posts/2023/cam-se-analysis/ https://ncar.github.io/esds/posts/2022/cam-se-regridding/