Skip to end of metadata
Go to start of metadata

You are viewing an old version of this page. View the current version.

Compare with Current View Page History

« Previous Version 49 Next »

This is a working document showing the conversion status for publishing E3SM data to CMIP6. The variables listed are unordered, and pulled from the priority 1 request list as well as scientist input. If you think that an important variable is missing please add it.

The conversion formulas are for the most part compiled by a NON DOMAIN EXPERT (me) using the following resources. There is the very real possibility they may be incorrect, please verify the conversions in your area of expertise and mark the "verified" column with your name.

Much of this information is pulled from the previous discussion page E3SM data conversion for CMIP6, but I wanted to distill the information down to the bare necessities. 

Resources:

ATM

Amon (monthly atmospheric data)

Notes:



CMIP6 nameCMIP6 descriptionE3SM variable(s)conversion formulaCMOR handler completeconversion formula verified by
1tasNear-Surface Air TemperatureTREFHTtas = TREFHT, no changeyesJ.Zhang
2tsSurface TemperatureTSts = TS, no changeyesJ.Zhang
3pslSea Level PressurePSLpsl = PSL, no changeyesJ.Zhang
4psSurface Air PressurePSps = PS, no changeyesJ.Zhang
5sfcWindNear-Surface Wind SpeedU10sfcWind = U10, no changeyesJ.Zhang
6hussNear-Surface Specific HumidityQREFHThuss = QREFHT, no changeyesJ.Zhang
7prPrecipitationPRECC PRECLpr = (PRECC  + PRECL) * 1000yesJ.Zhang
8prcConvective PrecipitationPRECCprc = PRECC * 1000
J.Zhang
9prsnSnowfall FluxPRECSC, PRECSLprsn = (PRECSC + PRECSL) * 1000yesJ.Zhang
10evspsblEvaporation Including Sublimation and TranspirationQFLXevspsbl = QFLX, no changeyesJ.Zhang
11tauuSurface Downward Eastward Wind StressTAUXtauu = TAUX, no changeyesJ.Zhang
12tauvSurface Downward Northward Wind StressTAUYtauv = TAUY, no changeyesJ.Zhang
13hflsSurface Upward Latent Heat FluxLHFLXhfls = LHFLX, no changeyesJ.Zhang
14cltTotal Cloud Cover PercentageCLDTOTclt = CLDTOT * 100., unit conversion
yesChris Golaz
15rldsSurface Downwelling Longwave RadiationFLDSrlds = FLDS, no changeyes

J.Zhang

Chris Golaz

16rlusSurface Upwelling Longwave RadiationFLDS, FLNSrlus = FLDS + FLNSyes

J.Zhang

Chris Golaz

17rsdsSurface Downwelling Shortwave RadiationFSDSrsds = FSDSyes

J.Zhang

Chris Golaz

18rsusSurface Upwelling Shortwave RadiationFSDS, FSNSrsus = FSDS - FSNSyes

J.Zhang

Chris Golaz

19hfssSurface Upward Sensible Heat FluxSHFLXhfss = SHFLXyesJ.Zhang
20clPercentage Cloud CoverCLOUDcl = CLOUD *100., on model levels, unit conversion

Chris Golaz
21cliMass Fraction of Cloud IceCLDICEcli = CLDICE on model levels

Kai Zhang

Note: for consistency with clivi, don't include snow in cli.

22cliviIce Water PathTGCLDIWPclivi = TGCLDIWP

Kai Zhang

Note: TGCLDIWP doesn't include snow water path.

23clwMass Fraction of Cloud Liquid WaterCLDLIQclw = CLDLIQ on model levels

Kai Zhang

Note: for consistency with clwvi, don't include rain in clw. 

24clwviCondensed Water PathTGCLDCWPclwvi = TGCLDCWP

Kai Zhang

Note: TGCLDCWP doesn't include rain and snow water path.

25hurRelative HumidityRELHUMhur = RELUM interpolated to 19 pressure levels
Chris Golaz
26hursNear-Surface Relative HumidityRHREFHThurs = RHREFHT

27husSpecific HumidityQhus = Q interpolated to 19 pressure levels
Chris Golaz
28mcConvective Mass Flux



29o3Mole Fraction of O3O3To verify, interpolated to 19 pressure levelsNeed input from Philip Cameron-Smith (Unlicensed)
30pfullPressure at Model Full-LevelsP0, PS, hyam, hybmpfull = P0*hyam + PS*hybm
Chris Golaz
31phalfPressure on Model Half-LevelsP0, PS, hyai, hybiphalf = P0*hyai + PS*hybi
Chris Golaz
32prwWater Vapor PathTMQprw = TMQ
Chris Golaz
33rldscsSurface Downwelling Clear-Sky Longwave RadiationFLDS, FLNS, FLNSC
rldscs = FLDS + FLNS - FLNSC

Chris Golaz
34rlutTOA Outgoing Longwave RadiationFSNTOA, FSNT, FLNT

rlut = FSNTOA - FSNT + FLNT


Note: equation originally from NCAR. Using this ensures that individual TOA terms (SW, LW, up/down) are consistent with TOM net flux (rtmt):

rtmt = rsdt - rsut - rlut

Chris Golaz

35rlutcsTOA Outgoing Clear-Sky Longwave RadiationFLUTCrlutcs = FLUTC

Note: slight approximation here since we are using TOM, but NCAR does the same.

Chris Golaz

36rsdscsSurface Downwelling Clear-Sky Shortwave RadiationFSDSCrsdscs = FSDSC
Chris Golaz
37rsdtTOA Incident Shortwave RadiationSOLINrsdt = SOLIN
Chris Golaz
38rsuscsSurface Upwelling Clear-Sky Shortwave RadiationFSDSC, FSNSCrsuscs = FSDSC - FSNSC
Chris Golaz
39rsutTOA Outgoing Shortwave RadiationFSUTOArsut = FSUTOA
Chris Golaz
40rsutcsTOA Outgoing Clear-Sky Shortwave RadiationFSUTOACrsutcs = FSUTOAC
Chris Golaz
41rtmtNet Downward Radiative Flux at Top of ModelFSNT, FLNTrtmt = FSNT - FLNT
Chris Golaz
42taAir TemperatureTta = T, interpolated to 19 pressure levels
Chris Golaz
43tasmaxDaily Maximum Near-Surface Air TemperatureTREFHTMXWe cannot use monthly h0 output for this. But the necessary data is in daily h1 files. We would have to create time series of monthly averages from these daily files.

44tasminDaily Minimum Near-Surface Air TemperatureTREFHTMNSame as above

45uaEastward WindUua = U, interpolated to 19 pressure levels
Chris Golaz
46uasEastward Near-Surface Wind



47vaNorthward WindVva = V, interpolated to 19 pressure levels
Chris Golaz
48vasNorthward Near-Surface Wind



49wapOmega (=dp/dt)OMEGAwap = OMEGA, interpolated to 19 pressure levels
Chris Golaz
50zgGeopotential HeightZ3zg = Z3, interpolated to 19 pressure levels
Chris Golaz

AERmon
51abs550aerAmbient Aerosol Absorption Optical Thickness at 550nmAODABSabs550aer = AODABS

52od550aerAmbient Aerosol Optical Thickness at 550nmAODVISod550aer = AODVIS

53reffclwtopCloud-Top Effective Droplet RadiusARELreffclwtop = AREL

Additional tables to consider


CFmon


CMIP6 nameCMIP6 descriptionE3SM variableconversion formulaunit and dimensionconversion formula verified by
1clisccp

ISCCP Cloud Area Percentage

FISCCP1_COSP

clisccp = FISCCP1_COSP

%

(time, plev7c, tau, lat, lon)

2cltisccp

ISCCP Total Cloud Cover Percentage

CLDTOT_ISCCP

cltisccp = CLDTOT_ISCCP%

(time, lat, lon)

3albisccp

ISCCP Mean Cloud Albedo; (Time-means are weighted by the ISCCP Total Cloud Fraction {:cltisccp})

MEANCLDALB_ISCCP

albisccp = MEANCLDALB_ISCCP1

(time, lat, lon)

4pctisccp

ISCCP Mean Cloud Top Pressure; (Time-means are weighted by the ISCCP Total Cloud Fraction {:cltisccp})

MEANPTOP_ISCCP

pctisccp = MEANPTOP_ISCCPPa

(time, lat, lon)

5clcalipso

CALIPSO Percentage Cloud Cover

CLD_CAL

clcalipso = CLD_CAL

%

(time, alt40, lat, lon)

6cltcalipso

CALIPSO Total Cloud Cover Percentage

CLDTOT_CAL

cltcalipso = CLDTOT_CAL%

(time, lat, lon)

7cllcalipso

CALIPSO Low Level Cloud Cover Percentage

CLDLOW_CAL

cllcalipso = CLDLOW_CAL%

(time, lat, lon)

8clmcalipso

CALIPSO Mid Level Cloud Cover Percentage

CLDMED_CAL

clmcalipso = CLDMED_CAL%

(time, lat, lon)

9clhcalipso

CALIPSO High Level Cloud Area Percentage

CLDHGH_CAL

clhcalipso = CLDHGH_CAL%

(time, lat, lon)



First set of atmospheric data has been cmorized provided by Sterling:

Data access:

  • acme1 -> /p/user_pub/work/E3SM/1_0/cmip6_variables/piControl/atm/CMIP6
  • Edison -> /global/cscratch1/sd/sbaldwin/CMIP6

{

    dimensions:

       time = UNLIMITED ; // (6000 currently)

       lat = 180 ;

       lon = 360 ;

       bnds = 2 ;

    variables:

       double time(time) ;

             time:bounds = "time_bnds" ;

             time:units = "days since 0001-01-01 00:00:00" ;

             time:calendar = "noleap" ;

             time:axis = "T" ;

             time:long_name = "time" ;

             time:standard_name = "time" ;

       double time_bnds(time, bnds) ;

       double lat(lat) ;

             lat:bounds = "lat_bnds" ;

             lat:units = "degrees_north" ;

             lat:axis = "Y" ;

             lat:long_name = "latitude" ;

             lat:standard_name = "latitude" ;

       double lat_bnds(lat, bnds) ;

       double lon(lon) ;

             lon:bounds = "lon_bnds" ;

             lon:units = "degrees_east" ;

             lon:axis = "X" ;

             lon:long_name = "longitude" ;

             lon:standard_name = "longitude" ;

       double lon_bnds(lon, bnds) ;

       double height ;

             height:units = "m" ;

             height:axis = "Z" ;

             height:positive = "up" ;

             height:long_name = "height" ;

             height:standard_name = "height" ;

       float tas(time, lat, lon) ;

             tas:standard_name = "air_temperature" ;

             tas:long_name = "Near-Surface Air Temperature" ;

             tas:comment = "near-surface (usually, 2 meter) air temperature" ;

             tas:units = "K" ;

             tas:cell_methods = "area: time: mean" ;

             tas:cell_measures = "area: areacella" ;

             tas:history = "2019-01-26T00:52:57Z altered by CMOR: Treated scalar dimension: \'height\'." ;

             tas:coordinates = "height" ;

             tas:missing_value = 1.e+20f ;

             tas:_FillValue = 1.e+20f ;

// global attributes:
:Conventions = "CF-1.7 CMIP-6.2" ;
:activity_id = "CMIP" ;
:branch_method = "Spin-up documentation" ;
:branch_time_in_child = 0. ;
:branch_time_in_parent = 0. ;
:contact = "Dave Bader (bader2@llnl.gov)" ;
:creation_date = "2019-01-26T00:52:57Z" ;
:data_specs_version = "01.00.23" ;
:experiment = "pre-industrial control" ;
:experiment_id = "piControl" ;
:external_variables = "areacella" ;
:forcing_index = 1 ;
:frequency = "mon" ;
:further_info_url = "https://furtherinfo.es-doc.org/CMIP6.E3SM-Project.E3SM-1-0.piControl.none.r1i1p1f1" ;
:grid = "gs1x1" ;
:grid_label = "gr" ;
:history = "2019-01-26T00:52:57Z ;rewrote data to be consistent with CMIP for variable tas found in table Amon.;\n",
"Output from 20180129.DECKv1b_piControl.ne30_oEC.edison" ;
:initialization_index = 1 ;
:institution = "LLNL (Lawrence Livermore National Laboratory, Livermore, CA 94550, USA); ANL (Argonne National Laboratory, Argonne, IL 60439, USA); BNL (Brookhaven National Laboratory, Upton, NY 11973, USA); LANL (Los Alamos National Laboratory, Los Alamos, NM 87545, USA); LBNL (Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA); ORNL (Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA); PNNL (Pacific Northwest National Laboratory, Richland, WA 99352, USA); SNL (Sandia National Laboratories, Albuquerque, NM 87185, USA). Mailing address: LLNL Climate Program, c/o David C. Bader, Principal Investigator, L-103, 7000 East Avenue, Livermore, CA 94550, USA" ;
:institution_id = "E3SM-Project" ;
:mip_era = "CMIP6" ;
:nominal_resolution = "100 km" ;
:parent_activity_id = "CMIP" ;
:parent_experiment_id = "piControl-spinup" ;
:parent_mip_era = "CMIP6" ;
:parent_source_id = "E3SM-1-0" ;
:parent_time_units = "days since 0001-01-01" ;
:parent_variant_label = "r1i1p1f1" ;
:physics_index = 1 ;
:product = "model-output" ;
:realization_index = 1 ;
:realm = "atmos" ;
:references = "Golaz, J.-C., P. M. Caldwell, L. P. Van Roekel and co-authors, 2019: The DOE E3SM coupled model version 1: Overviewand evaluation at standard resolution. JAMES, submitted; http://e3sm.org" ;
:source = "E3SM 1.0 (2018): \n",
       "aerosol: MAM4 with resuspension, marine organics, and secondary organics (same grid as atmos)\n",
       "atmos: EAM (v1.0, cubed sphere spectral-element grid; 5400 elements with p=3; 1 deg average grid spacing; 90 x 90 x 6 longitude/latitude/cubeface; 72 levels; top level 0.1 hPa)\n",
       "atmosChem: Troposphere specified oxidants for aerosols. Stratosphere linearized interactive ozone (LINOZ v2) (same grid as atmos)\n",
       "land: ELM (v1.0, cubed sphere spectral-element grid; 5400 elements with p=3; 1 deg average grid spacing; 90 x 90 x 6 longitude/latitude/cubeface; satellite phenology mode), MOSART (v1.0, 0.5 degree latitude/longitude grid)\n",
       "landIce: none\n",
       "ocean: MPAS-Ocean (v6.0, oEC60to30 unstructured SVTs mesh with 235160 cells and 714274 edges, variable resolution 60 km to 30 km; 60 levels; top grid cell 0-10 m)\n",
       "ocnBgchem: none\n",
       "seaIce: MPAS-Seaice (v6.0, same grid as ocean)" ;
:source_id = "E3SM-1-0" ;
:source_type = "AOGCM AER" ;
:sub_experiment = "none" ;
:sub_experiment_id = "none" ;
:table_id = "Amon" ;
:table_info = "Creation Date:(21 March 2018) MD5:528cdc7b6fa35bf79fa0b69273987051" ;
:title = "E3SM-1-0 output prepared for CMIP6" ;
:variable_id = "tas" ;
:variant_label = "r1i1p1f1" ;
:license = "CMIP6 model data produced by E3SM 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." ;
:cmor_version = "3.4.0" ;
:tracking_id = "hdl:21.14100/c1e2e21a-5d05-4a26-8b86-21f947ab44a1" ;
:e3sm_source_code_doi = "10.11578/E3SM/dc.20180418.36" ;
:e3sm_source_code_reference = "https://github.com/E3SM-Project/E3SM/releases/tag/v1.0.0" ;
:doe_acknowledgement = "This research was supported as part of the Energy Exascale Earth System Model (E3SM) project, funded by the U.S. Department of Energy, Office of Science, Office of Biological and Environmental Research." ;
:computational_acknowledgement = "The data were produced using resources of the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231." ;
:ncclimo_generation_command = "drc_in=\'/p/user_pub/work/E3SM/1_0/piControl/1deg_atm_60-30km_ocean/atmos/native/model-output/mon/ens1/v1\' # Input directory\n",
       "drc_out=\"/p/user_pub/work/E3SM/1_0/cmip6_variables/piControl/atm/ne30/climo\" # Native grid output directory\n",
       "drc_rgr=\"/p/user_pub/work/E3SM/1_0/cmip6_variables/piControl/atm/ne30/rgr\" # Regridded output directory\n",
       "drc_tmp=\'/p/cscratch/acme/baldwin32/tmp\' # Temporary/intermediate-file directory\n",
       "map=\"/export/zender1/data/maps/map_ne30np4_to_cmip6_180x360_aave.20181001.nc\" # Regridding map-file\n",
       "cmip6_opt=\'-7 --dfl_lvl=1 --no_cll_msr --no_frm_trm --no_stg_grd\' # CMIP6-specific options\n",
       "spl_opt=\'--yr_srt=1 --yr_end=500 --ypf=500\' # Splitter options\n",
       "vars=(\'TREFHT\' \'TS\' \'PSL\' \'PS\' \'QREFHT\' \'PRECC\' \'PRECL\' \'PRECSC\' \'PRECSL\' \'QFLX\' \'TAUX\' \'TAUY\' \'LHFLX\' \'CLDTOT\' \'FLDS\' \'FLNS\' \'FSDS\' \'SHFLX\' \'CLOUD\' \'CLDICE\' \'TGCLDIWP\' \'CLDLIQ\' \'RELHUM\' \'Q\' \'O3\' \'P0\' \'hyam\' \'hybi\' \'TMQ\' \'FLDSC\' \'FLUT\' \'FSDSC\' \'SOLIN\' \'FSNSC\' \'FSUTOA\' \'FSUTOAC\' \'FSNT\' \'FLNT\' \'T\' \'U\' \'V\' \'OMEGA\' \'Z3\')\n",
       "export TMPDIR=${drc_tmp}\n",
       "cd ${drc_in}\n",
       "\n",
       "for i in \"${vars[@]}\"\n",
       "do\n",
             " ls 20180129.DECKv1b_piControl.ne30_oEC.edison.cam.h0.0???-*.nc | ncclimo --var=${i} ${cmip6_opt} ${spl_opt} --map=${map} --drc_out=${drc_out} --drc_rgr=${drc_rgr} &> ${drc_out}/ncclimo_${i}.out &\n",
       "done" ;
:ncclimo_version = "4.7.9-alpha04" ;
}

LND


CMIP6 nameE3SM variable(s)conversion formulaCMOR handler completeverified
1mrsosSOILWATER_10CMmrsos = SOILWATER_10CM, no changeyes
2mrsoSOILICE, SOILLIQmrso = sum_over_z(SOILICE + SOILLIQ, capped_at=5000)yes
3mrfsoSOILICEmrfso = sum_over_z(SOILICE, capped_at=5000)yes
4mrrosQOVERmrros = QOVER, no changeyes
5mrroQRUNOFF


6prvegQINTRprveg = QINTR, no changeyes
7evspsblvegQVEGEevspsblveg = QVEGE, no changeyes
8evspsblsoiQSOILevspsblsoi = QSOIL, no changeyes
9tranQSOIL, QVEGTtran = QSOIL + QVEGTyes
10tslTSOItsl = TSOI, no changeyes
11laiLAISHA, LAISUNlai = LAISHA + LAISUNyes
12cLitterTOTLITC, CWDCcLitter = TOTLITC + CWDC

13cProductTOTPRODCcProduct = TOTPRODC, no change

14cSoilTOTSOMCcSoil = TOTSOMC, no change

15fFireCOL_FIRE_CLOSSfFire = COL_FIRE_CLOSS, no change

16fHarvest

WOOD_HARVESTC

fHarvest = WOOD_HARVESTC, no change

17fVegLitterLITFALL, gap_mortality_c_to_cwdc_col, fire_mortality_c_to_cwdc_col, harvest_c_to_cwdc_col, dwt_livecrootc_to_cwdc_col, dwt_deadcrootc_to_cwdc_col fVegLitter = LITFALL + gap_mortality_c_to_cwdc_col + fire_mortality_c_to_cwdc_col + harvest_c_to_cwdc_col + dwt_livecrootc_to_cwdc_col  + dwt_deadcrootc_to_cwdc_col 

18fLitterSoilLITTERC_LOSS, LITTERC_HR, M_LITR1C_TO_FIRE, M_LITR2C_TO_FIRE, M_LITR3C_TO_FIREfLitterSoil = LITTERC_LOSS - LITTERC_HR - M_LITR1C_TO_FIRE - M_LITR2C_TO_FIRE - M_LITR3C_TO_FIRE

19cVegTOTVEGCcVeg = TOTVEGC, no change

20nbpNBPnbp = NBP, no change

21gppGPPgpp = GPP, no change

22raARra = AR, no change

23rhHRrh = HR, no change

OCN


CMIP6 nameE3SM variable(s)conversion formulaCMOR handler completeverified
1masso



2sf6



3cfc11



4pbo



5cfc12



6pso



7dissicnatos



8fgsf6



9volo



10dissicabioos



11fgcfc11



12zos



13dissi14cabioos



14fgcfc12



15zossq



16dissi13cos



17zostoga



18masscello



ICE. – Notz et al. 2016 – https://www.geosci-model-dev.net/9/3427/2016/gmd-9-3427-2016.pdf


CMIP6 name

E3SM variable(s)


conversion formulaCMOR handler completeverified
1

Fraction of time steps with sea ice (sitimefrac)





2

Sea-ice area fraction (siconc)





3

Sea-ice mass per area (simass)





4

Sea-ice thickness (sithick)





5

Snow area fraction (sisnconc)





6

Snow mass per area (sisnmass)





7

Snow thickness (sisnthick)





8

Surface temperature (sitemptop)





9

Sea-ice volume per area (sivol)





  • No labels