WACCMX (free-run) Electron Density, winds and Vertical drift data

The netcdf file e.g., "waccmx_cam6_0_032_smin_clim.cam.h1.0004-09-17-00000.nc" contains select outputs from a Free running WACCM-X with constant solar flux of f107=70 sfu and constant kp=1, with Heelis model for high-latitude fields. These a subset of the original model output. As this is a free-run case, the data represents a generative year, e.g., in this case Year='0004'. The times are in UT and counted in days since year=0001.

Below is the netcdf header for this example file:

netcdf waccmx_cam6_0_032_smin_clim.cam.h1.0004-09-17-00000 {dimensions:        time = UNLIMITED ; // (24 currently)        lat = 96 ;        lon = 144 ;        lev = 60 ;        ilev = 127 ;        nbnd = 2 ;variables:        double P0 ;                P0:long_name = "reference pressure" ;                P0:units = "Pa" ;        float PS(time, lat, lon) ;                PS:units = "Pa" ;                PS:long_name = "Surface pressure" ;        float U(time, lev, lat, lon) ;                U:mdims = 1 ;                U:units = "m/s" ;                U:long_name = "Zonal wind" ;        float V(time, lev, lat, lon) ;                V:mdims = 1 ;                V:units = "m/s" ;                V:long_name = "Meridional wind" ;        float WI(time, lev, lat, lon) ;                WI:mdims = 1 ;                WI:units = "m/s" ;                WI:long_name = "WI Vertical ion drift from edynamo" ;        float Z3(time, lev, lat, lon) ;                Z3:mdims = 1 ;                Z3:units = "m" ;                Z3:long_name = "Geopotential Height (above sea level)" ;        float e(time, lev, lat, lon) ;                e:mdims = 1 ;                e:units = "mol/mol" ;                e:long_name = "e concentration" ;        double hyai(ilev) ;                hyai:long_name = "hybrid A coefficient at layer interfaces" ;        double hyam(lev) ;                hyam:long_name = "hybrid A coefficient at layer midpoints" ;        double hybi(ilev) ;                hybi:long_name = "hybrid B coefficient at layer interfaces" ;        double hybm(lev) ;                hybm:long_name = "hybrid B coefficient at layer midpoints" ;        double ilev(ilev) ;                ilev:long_name = "hybrid level at interfaces (1000*(A+B))" ;                ilev:units = "hPa" ;                ilev:positive = "down" ;                ilev:standard_name = "atmosphere_hybrid_sigma_pressure_coordinate" ;                ilev:formula_terms = "a: hyai b: hybi p0: P0 ps: PS" ;        double lat(lat) ;                lat:_FillValue = -900. ;                lat:long_name = "latitude" ;                lat:units = "degrees_north" ;        double lev(lev) ;                lev:long_name = "hybrid level at midpoints (1000*(A+B))" ;                lev:units = "hPa" ;                lev:positive = "down" ;                lev:standard_name = "atmosphere_hybrid_sigma_pressure_coordinate" ;                lev:formula_terms = "a: hyam b: hybm p0: P0 ps: PS" ;        double lon(lon) ;                lon:_FillValue = -900. ;                lon:long_name = "longitude" ;                lon:units = "degrees_east" ;        double time(time) ;                time:long_name = "time" ;                time:units = "days since 0001-01-01 00:00:00" ;                time:calendar = "noleap" ;                time:bounds = "time_bnds" ;        double time_bnds(time, nbnd) ;                time_bnds:long_name = "time interval endpoints" ;