From 324b6d67517478942436290b73b85dbcc7f46780 Mon Sep 17 00:00:00 2001
From: Cheng Li <69489965+chengcli@users.noreply.github.com>
Date: Sun, 24 Mar 2024 17:23:07 -0400
Subject: [PATCH] Add RRTMG-SW (#135)

- Bind RRTMG-SW to C
---
 CMakeLists.txt                          |   5 +-
 cmake/{ => examples}/bryan.cmake        |   0
 cmake/{ => examples}/cart3.cmake        |   0
 cmake/{ => examples}/dam.cmake          |   0
 cmake/{ => examples}/default.cmake      |   0
 cmake/{ => examples}/exo2.cmake         |   0
 cmake/{ => examples}/exo3.cmake         |   0
 cmake/{ => examples}/hjupiter.cmake     |   0
 cmake/{ => examples}/hywater.cmake      |   0
 cmake/{ => examples}/juno.cmake         |   0
 cmake/{ => examples}/jupiter_crm.cmake  |   0
 cmake/{ => examples}/ktable_earth.cmake |   0
 cmake/{ => examples}/ktable_jup.cmake   |   0
 cmake/{ => examples}/mars_far_ir.cmake  |   0
 cmake/{ => examples}/polar3.cmake       |   0
 cmake/examples/rcemip.cmake             |   9 ++
 cmake/{ => examples}/robert.cmake       |   0
 cmake/{ => examples}/saturn_radio.cmake |   0
 cmake/{ => examples}/straka.cmake       |   0
 cmake/{ => examples}/uranus_ir.cmake    |   0
 cmake/{ => examples}/uranus_mwr.cmake   |   0
 cmake/{ => examples}/vortex.cmake       |   0
 cmake/parameters.cmake                  |  12 +++
 cmake/rrtmg.cmake                       |  29 ++++++
 examples/2024-CLi-RCEMIP/CMakeLists.txt |  11 +++
 examples/2024-CLi-RCEMIP/rce_large.inp  |   0
 examples/2024-CLi-RCEMIP/rce_small.inp  | 107 +++++++++++++++++++++
 examples/2024-CLi-RCEMIP/rcemip.cpp     | 120 ++++++++++++++++++++++++
 examples/2024-CLi-RCEMIP/rcemip.yaml    |   9 ++
 examples/CMakeLists.txt                 |   4 +
 src/outputs/output_utils.cpp            |  89 ++++++++++++++++++
 31 files changed, 393 insertions(+), 2 deletions(-)
 rename cmake/{ => examples}/bryan.cmake (100%)
 rename cmake/{ => examples}/cart3.cmake (100%)
 rename cmake/{ => examples}/dam.cmake (100%)
 rename cmake/{ => examples}/default.cmake (100%)
 rename cmake/{ => examples}/exo2.cmake (100%)
 rename cmake/{ => examples}/exo3.cmake (100%)
 rename cmake/{ => examples}/hjupiter.cmake (100%)
 rename cmake/{ => examples}/hywater.cmake (100%)
 rename cmake/{ => examples}/juno.cmake (100%)
 rename cmake/{ => examples}/jupiter_crm.cmake (100%)
 rename cmake/{ => examples}/ktable_earth.cmake (100%)
 rename cmake/{ => examples}/ktable_jup.cmake (100%)
 rename cmake/{ => examples}/mars_far_ir.cmake (100%)
 rename cmake/{ => examples}/polar3.cmake (100%)
 create mode 100644 cmake/examples/rcemip.cmake
 rename cmake/{ => examples}/robert.cmake (100%)
 rename cmake/{ => examples}/saturn_radio.cmake (100%)
 rename cmake/{ => examples}/straka.cmake (100%)
 rename cmake/{ => examples}/uranus_ir.cmake (100%)
 rename cmake/{ => examples}/uranus_mwr.cmake (100%)
 rename cmake/{ => examples}/vortex.cmake (100%)
 create mode 100644 cmake/rrtmg.cmake
 create mode 100644 examples/2024-CLi-RCEMIP/CMakeLists.txt
 create mode 100644 examples/2024-CLi-RCEMIP/rce_large.inp
 create mode 100644 examples/2024-CLi-RCEMIP/rce_small.inp
 create mode 100644 examples/2024-CLi-RCEMIP/rcemip.cpp
 create mode 100644 examples/2024-CLi-RCEMIP/rcemip.yaml

diff --git a/CMakeLists.txt b/CMakeLists.txt
index ae5160f7..540f1bb3 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -33,10 +33,10 @@ message(STATUS "2. Set up project parameters ...")
 # load custom task
 message(STATUS "Load custom task = ${TASK}")
 if (NOT ${TASK} STREQUAL "")
-  if (EXISTS ${CMAKE_CURRENT_SOURCE_DIR}/cmake/${TASK}.cmake)
+  if (EXISTS ${CMAKE_CURRENT_SOURCE_DIR}/cmake/examples/${TASK}.cmake)
     message(STATUS "Load custom setting - ${TASK} -")
     message(STATUS "Include ${CMAKE_CURRENT_SOURCE_DIR}/cmake/${TASK}.cmake")
-    include(${CMAKE_CURRENT_SOURCE_DIR}/cmake/${TASK}.cmake)
+    include(${CMAKE_CURRENT_SOURCE_DIR}/cmake/examples/${TASK}.cmake)
   else()
     message(FATAL_ERROR "Cannot find custom file cmake/${TASK}.cmake")
   endif()
@@ -59,6 +59,7 @@ include(${CMAKE_SOURCE_DIR}/cmake/athena.cmake)
 include(${CMAKE_SOURCE_DIR}/cmake/application.cmake)
 include(${CMAKE_SOURCE_DIR}/cmake/rfm.cmake)
 include(${CMAKE_SOURCE_DIR}/cmake/disort.cmake)
+include(${CMAKE_SOURCE_DIR}/cmake/rrtmg.cmake)
 include(${CMAKE_SOURCE_DIR}/cmake/pvfmm.cmake)
 include(${CMAKE_SOURCE_DIR}/cmake/pybind11.cmake)
 
diff --git a/cmake/bryan.cmake b/cmake/examples/bryan.cmake
similarity index 100%
rename from cmake/bryan.cmake
rename to cmake/examples/bryan.cmake
diff --git a/cmake/cart3.cmake b/cmake/examples/cart3.cmake
similarity index 100%
rename from cmake/cart3.cmake
rename to cmake/examples/cart3.cmake
diff --git a/cmake/dam.cmake b/cmake/examples/dam.cmake
similarity index 100%
rename from cmake/dam.cmake
rename to cmake/examples/dam.cmake
diff --git a/cmake/default.cmake b/cmake/examples/default.cmake
similarity index 100%
rename from cmake/default.cmake
rename to cmake/examples/default.cmake
diff --git a/cmake/exo2.cmake b/cmake/examples/exo2.cmake
similarity index 100%
rename from cmake/exo2.cmake
rename to cmake/examples/exo2.cmake
diff --git a/cmake/exo3.cmake b/cmake/examples/exo3.cmake
similarity index 100%
rename from cmake/exo3.cmake
rename to cmake/examples/exo3.cmake
diff --git a/cmake/hjupiter.cmake b/cmake/examples/hjupiter.cmake
similarity index 100%
rename from cmake/hjupiter.cmake
rename to cmake/examples/hjupiter.cmake
diff --git a/cmake/hywater.cmake b/cmake/examples/hywater.cmake
similarity index 100%
rename from cmake/hywater.cmake
rename to cmake/examples/hywater.cmake
diff --git a/cmake/juno.cmake b/cmake/examples/juno.cmake
similarity index 100%
rename from cmake/juno.cmake
rename to cmake/examples/juno.cmake
diff --git a/cmake/jupiter_crm.cmake b/cmake/examples/jupiter_crm.cmake
similarity index 100%
rename from cmake/jupiter_crm.cmake
rename to cmake/examples/jupiter_crm.cmake
diff --git a/cmake/ktable_earth.cmake b/cmake/examples/ktable_earth.cmake
similarity index 100%
rename from cmake/ktable_earth.cmake
rename to cmake/examples/ktable_earth.cmake
diff --git a/cmake/ktable_jup.cmake b/cmake/examples/ktable_jup.cmake
similarity index 100%
rename from cmake/ktable_jup.cmake
rename to cmake/examples/ktable_jup.cmake
diff --git a/cmake/mars_far_ir.cmake b/cmake/examples/mars_far_ir.cmake
similarity index 100%
rename from cmake/mars_far_ir.cmake
rename to cmake/examples/mars_far_ir.cmake
diff --git a/cmake/polar3.cmake b/cmake/examples/polar3.cmake
similarity index 100%
rename from cmake/polar3.cmake
rename to cmake/examples/polar3.cmake
diff --git a/cmake/examples/rcemip.cmake b/cmake/examples/rcemip.cmake
new file mode 100644
index 00000000..458d8f55
--- /dev/null
+++ b/cmake/examples/rcemip.cmake
@@ -0,0 +1,9 @@
+# configure RCEMIP case
+set(NUMBER_GHOST_CELLS 3)
+set(NVAPOR 1)
+set(NCLOUD 2)
+set(NETCDF ON)
+set(PNETCDF ON)
+set(MPI ON)
+set(RRTMG_SW ON)
+set(RSOLVER lmars)
diff --git a/cmake/robert.cmake b/cmake/examples/robert.cmake
similarity index 100%
rename from cmake/robert.cmake
rename to cmake/examples/robert.cmake
diff --git a/cmake/saturn_radio.cmake b/cmake/examples/saturn_radio.cmake
similarity index 100%
rename from cmake/saturn_radio.cmake
rename to cmake/examples/saturn_radio.cmake
diff --git a/cmake/straka.cmake b/cmake/examples/straka.cmake
similarity index 100%
rename from cmake/straka.cmake
rename to cmake/examples/straka.cmake
diff --git a/cmake/uranus_ir.cmake b/cmake/examples/uranus_ir.cmake
similarity index 100%
rename from cmake/uranus_ir.cmake
rename to cmake/examples/uranus_ir.cmake
diff --git a/cmake/uranus_mwr.cmake b/cmake/examples/uranus_mwr.cmake
similarity index 100%
rename from cmake/uranus_mwr.cmake
rename to cmake/examples/uranus_mwr.cmake
diff --git a/cmake/vortex.cmake b/cmake/examples/vortex.cmake
similarity index 100%
rename from cmake/vortex.cmake
rename to cmake/examples/vortex.cmake
diff --git a/cmake/parameters.cmake b/cmake/parameters.cmake
index 22ed030b..d41dad8a 100644
--- a/cmake/parameters.cmake
+++ b/cmake/parameters.cmake
@@ -95,6 +95,18 @@ else()
   set(DISORT_OPTION "RT_DISORT")
 endif()
 
+if(NOT RRTMG_SW OR NOT DEFINED RRTMG_SW)
+  set(RRTMG_SW_OPTION "NOT_RT_RRTMG_SW")
+else()
+  set(RRTMG_SW_OPTION "RT_RRTMG_SW")
+endif()
+
+if(NOT RRTMG_LW OR NOT DEFINED RRTMG_LW)
+  set(RRTMG_LW_OPTION "NOT_RT_RRTMG_LW")
+else()
+  set(RRTMG_LW_OPTION "RT_RRTMG_LW")
+endif()
+
 if(NOT PVFMM OR NOT DEFINED PVFMM)
   set(PVFMM_OPTION "DISABLE_PVFMM")
 else()
diff --git a/cmake/rrtmg.cmake b/cmake/rrtmg.cmake
new file mode 100644
index 00000000..737a23e0
--- /dev/null
+++ b/cmake/rrtmg.cmake
@@ -0,0 +1,29 @@
+include(FetchContent)
+
+set(FETCHCONTENT_QUIET FALSE)
+
+macro(set_if_empty _variable)
+  if("${${_variable}}" STREQUAL "")
+    set(${_variable} ${ARGN})
+  endif()
+endmacro()
+
+option(RRTMG_SW "Build RRTMG_SW" OFF)
+
+if(RRTMG_SW)
+  FetchContent_Declare(
+    RRTMG_SW
+    GIT_REPOSITORY https://github.com/chengcli/RRTMG_SW
+    GIT_TAG master)
+  FetchContent_MakeAvailable(RRTMG_SW)
+endif()
+
+option(RRTMG_LW "Build RRTMG_LW" OFF)
+
+if(RRTMG_LW)
+  FetchContent_Declare(
+    RRTMG_LW
+    GIT_REPOSITORY https://github.com/chengcli/RRTMG_LW
+    GIT_TAG master)
+  FetchContent_MakeAvailable(RRTMG_LW)
+endif()
diff --git a/examples/2024-CLi-RCEMIP/CMakeLists.txt b/examples/2024-CLi-RCEMIP/CMakeLists.txt
new file mode 100644
index 00000000..62e1fb22
--- /dev/null
+++ b/examples/2024-CLi-RCEMIP/CMakeLists.txt
@@ -0,0 +1,11 @@
+
+# 1. Compile rcemip
+
+# 2. Compile hjupiter
+setup_problem(rcemip)
+
+# 2. Copy input files to run directory
+file(GLOB inputs *.py *.inp *.yaml)
+foreach(input ${inputs})
+  execute_process(COMMAND ln -sf ${input} ${CMAKE_BINARY_DIR}/bin/${inp})
+endforeach()
diff --git a/examples/2024-CLi-RCEMIP/rce_large.inp b/examples/2024-CLi-RCEMIP/rce_large.inp
new file mode 100644
index 00000000..e69de29b
diff --git a/examples/2024-CLi-RCEMIP/rce_small.inp b/examples/2024-CLi-RCEMIP/rce_small.inp
new file mode 100644
index 00000000..f66e5de3
--- /dev/null
+++ b/examples/2024-CLi-RCEMIP/rce_small.inp
@@ -0,0 +1,107 @@
+<comment>
+problem     = RCEMPI-1 Small Domain
+
+<job>
+problem_id  = rce1_small    # problem ID: basename of output filenames
+
+<output0>
+file_type   = rst
+dt          = 10.E5
+
+<output1>
+file_type   = hst           # History data dump
+dt          = 2.E4          # time increment between outputs
+
+<output2>
+file_type   = pnetcdf        # Netcdf format output
+variable    = prim          # variables to be output
+dt          = 1.E4          # time increment between outputs
+
+<output3>
+file_type   = pnetcdf        # Netcdf format output
+variable    = uov           # diagnostic variables
+dt          = 1.E4          # time increment
+
+<output3>
+file_type   = pnetcdf        # Netcdf format output
+variable    = diag           # diagnostic variables
+dt          = 1.E4          # time increment
+
+<time>
+cfl_number  = 0.9
+nlim        = -1            # cycle limit
+tlim        = 100.E5
+xorder      = 5             # horizontal reconstruction order
+integrator  = rk3           # integration method
+
+<mesh>
+nx1         = 80          # Number of zones in X1-direction
+x1min       = 0.E3     # minimum value of X1
+x1max       = 34.E3       # maximum value of X1
+ix1_bc      = reflecting  # Inner-X1 boundary condition flag
+ox1_bc      = reflecting  # Outer-X1 boundary condition flag
+
+nx2         = 100         # Number of zones in X2-direction
+x2min       = 0.          # minimum value of X2
+x2max       = 100.E3      # maximum value of X2
+ix2_bc      = periodic    # Inner-X2 boundary condition flag
+ox2_bc      = periodic    # Outer-X2 boundary condition flag
+
+nx3         = 100         # Number of zones in X3-direction
+x3min       = 0.          # minimum value of X3
+x3max       = 100.E3      # maximum value of X3
+ix3_bc      = periodic    # Inner-X3 boundary condition flag
+ox3_bc      = periodic    # Outer-X3 boundary condition flag
+
+<meshblock>
+nx1         = 80
+nx2         = 50
+nx3         = 50
+
+<hydro>
+gamma       = 1.42
+grav_acc1   = -9.79764
+implicit_flag = 1
+
+<species>
+vapor = H2O
+cloud = H2O(c), H2O(p)
+
+<chemistry>
+microphysics_config = rcemip.yaml
+Nc = 1.e8
+Ni = 1.e5
+
+<thermodynamics>
+Rd          = 287.04   # mu = 2.3175 g/mol
+
+eps1        = 8.18    8.18    8.18
+beta1       = 0.      24.845  24.845
+rcp1        = 0.15    0.33    0.33
+Ttriple1    = 273.16
+Ptriple1    = 611.7
+
+sa.relax    = 1.0
+sa.max_iter = 4
+sa.ftol     = 1.e-2
+
+#<forcing>
+#packages = fix_bot_temp
+
+<problem>
+radius      = 6471.e3
+xCO2        = 348.e-6
+xCH4        = 1650.e-12
+xN2O        = 306.e-12
+
+zt          = 15.e3
+zq1         = 4.e3
+zq2         = 7.5e3
+
+qt          = 1.e-14
+q0          = 12.e-3
+
+P0          = 1.0148e5
+T0          = 295.
+gamma       = 0.0067
+diagnostics = div, curl, b, mean, div_h, tempa, presa, hydroflux
diff --git a/examples/2024-CLi-RCEMIP/rcemip.cpp b/examples/2024-CLi-RCEMIP/rcemip.cpp
new file mode 100644
index 00000000..8e16737c
--- /dev/null
+++ b/examples/2024-CLi-RCEMIP/rcemip.cpp
@@ -0,0 +1,120 @@
+// athena
+#include <athena/athena.hpp>
+#include <athena/coordinates/coordinates.hpp>
+#include <athena/eos/eos.hpp>
+#include <athena/field/field.hpp>
+#include <athena/hydro/hydro.hpp>
+#include <athena/mesh/mesh.hpp>
+#include <athena/parameter_input.hpp>
+
+// climath
+#include <climath/core.h>
+
+// snap
+#include <snap/thermodynamics/thermodynamics.hpp>
+
+//! variables used in initial condition
+Real grav, P0, T0, gamma;
+Real xCO2, xCH4, xN2O;
+Real zt, zq1, zq2;
+Real qt, q0;
+int iH2O = 1, iCO2 = 0, iCH4 = 1, iN2O = 2, iO3 = 3;
+
+void MeshBlock::InitUserMeshBlockData(ParameterInput *pin) {
+  AllocateUserOutputVariables(5);
+  SetUserOutputVariableName(0, "temp");
+  SetUserOutputVariableName(1, "theta");
+  SetUserOutputVariableName(2, "thetav");
+  SetUserOutputVariableName(3, "mse");
+  SetUserOutputVariableName(4, "rh_H2O");
+
+  AllocateRealUserMeshBlockDataField(1);
+
+  // CO2, CH4, N2O, O3 in mole fraction
+  ruser_meshblock_data[0].NewAthenaArray(4, ncells1);
+  for (int i = is; i <= ie; ++i) {
+    ruser_meshblock_data[0](iCO2, i) = xCO2;
+    ruser_meshblock_data[0](iCH4, i) = xCH4;
+    ruser_meshblock_data[0](iN2O, i) = xN2O;
+  }
+}
+
+void MeshBlock::UserWorkBeforeOutput(ParameterInput *pin) {
+  auto pthermo = Thermodynamics::GetInstance();
+
+  for (int k = ks; k <= ke; ++k)
+    for (int j = js; j <= je; ++j)
+      for (int i = is; i <= ie; ++i) {
+        user_out_var(0, k, j, i) = pthermo->GetTemp(this, k, j, i);
+        user_out_var(1, k, j, i) = pthermo->PotentialTemp(this, P0, k, j, i);
+        // theta_v
+        user_out_var(2, k, j, i) =
+            user_out_var(1, k, j, i) * pthermo->RovRd(this, k, j, i);
+        // mse
+        user_out_var(3, k, j, i) =
+            pthermo->MoistStaticEnergy(this, grav * pcoord->x1v(i), k, j, i);
+        user_out_var(4, k, j, i) = pthermo->RelativeHumidity(this, 1, k, j, i);
+      }
+}
+
+void Mesh::InitUserMeshData(ParameterInput *pin) {
+  grav = -pin->GetReal("hydro", "grav_acc1");
+
+  xCO2 = pin->GetReal("problem", "xCO2");
+  xCH4 = pin->GetReal("problem", "xCH4");
+  xN2O = pin->GetReal("problem", "xN2O");
+
+  zt = pin->GetReal("problem", "zt");
+  zq1 = pin->GetReal("problem", "zq1");
+  zq2 = pin->GetReal("problem", "zq2");
+
+  qt = pin->GetReal("problem", "qt");
+  q0 = pin->GetReal("problem", "q0");
+
+  P0 = pin->GetReal("problem", "P0");
+  T0 = pin->GetReal("problem", "T0");
+  gamma = pin->GetReal("problem", "gamma");
+}
+
+void MeshBlock::ProblemGenerator(ParameterInput *pin) {
+  auto pthermo = Thermodynamics::GetInstance();
+
+  Real Rd = pthermo->GetRd();
+  Real Tv0 = T0 * (1. + 0.608 * q0);
+  Real Tvt = Tv0 - gamma * zt;
+  Real Pt = P0 * pow(Tvt / Tv0, grav / (Rd * gamma));
+
+  // O3 parameters
+  Real g1 = 3.6478;  // ppmv hPa−g2
+  Real g2 = 0.83209;
+  Real g3 = 11.3515;  // hPa
+
+  for (int k = ks; k <= ke; ++k)
+    for (int j = js; j <= je; ++j)
+      for (int i = is; i <= ie; ++i) {
+        Real z = pcoord->x1v(i);
+        Real Tv;
+
+        if (z < zt) {
+          phydro->w(iH2O, k, j, i) = q0 * exp(-z / zq1) * exp(-sqr(z / zq2));
+          phydro->w(IPR, k, j, i) =
+              P0 * pow((Tv0 - gamma * z) / Tv0, grav / (Rd * gamma));
+          Tv = Tv0 - gamma * z;
+        } else {
+          phydro->w(iH2O, k, j, i) = qt;
+          phydro->w(IPR, k, j, i) = Pt * exp(-grav * (z - zt) / (Rd * Tvt));
+          Tv = Tvt;
+        }
+
+        phydro->w(IDN, k, j, i) = phydro->w(IPR, k, j, i) / (Rd * Tv);
+
+        // pa -> hpa
+        Real pre = phydro->w(IPR, k, j, i) / 100.;
+        // ppmv -> mole fraction
+        ruser_meshblock_data[0](iO3, i) =
+            1.e-6 * g1 * pow(pre, g2) * exp(-pre / g3);
+      }
+
+  peos->PrimitiveToConserved(phydro->w, pfield->bcc, phydro->u, pcoord, is, ie,
+                             js, je, ks, ke);
+}
diff --git a/examples/2024-CLi-RCEMIP/rcemip.yaml b/examples/2024-CLi-RCEMIP/rcemip.yaml
new file mode 100644
index 00000000..efa67ff9
--- /dev/null
+++ b/examples/2024-CLi-RCEMIP/rcemip.yaml
@@ -0,0 +1,9 @@
+name: giant planet chemistry model
+
+microphysics:
+  - water-system
+
+water-system:
+  scheme: Kessler94
+  dependent-species: [vapor.H2O, cloud.H2O(c), cloud.H2O(p)]
+  parameters: {autoconversion: 1.e-3, accretion: 0.0, evaporation: 1.e-3, sedimentation: -20.}
diff --git a/examples/CMakeLists.txt b/examples/CMakeLists.txt
index 9394182c..77a6ad24 100644
--- a/examples/CMakeLists.txt
+++ b/examples/CMakeLists.txt
@@ -44,3 +44,7 @@ endif()
 if (${TASK} STREQUAL "hjupiter")
   add_subdirectory(2024-XZhang-cloud-rt)
 endif()
+
+if (${TASK} STREQUAL "rcemip")
+  add_subdirectory(2024-CLi-RCEMIP)
+endif()
diff --git a/src/outputs/output_utils.cpp b/src/outputs/output_utils.cpp
index 5bb8be09..1dcfbcd8 100644
--- a/src/outputs/output_utils.cpp
+++ b/src/outputs/output_utils.cpp
@@ -72,6 +72,95 @@ __attribute__((weak)) MetadataTable::MetadataTable() {
       {"v1v2", "total upward velocity 2 flux", "(m/s)/(m^2.s)", "--C"},
       {"v1v3", "total upward velocity 3 flux", "(m/s)/(m^2.s)", "--C"},
       {"v1T", "total upward temperature flux", "K/(m^2.s)", "--C"},
+      // CMPI6 1D variables description
+      {"ta_avg", "domain avg. air temperature profile", "K", "--C"},
+      {"ua_avg", "domain avg. zonal wind profile", "m/s", "--C"},
+      {"va_avg", "domain avg. meridional wind profile", "m/s", "--C"},
+      {"hus_avg", "domain avg. specific humidity profile profile", "kg/kg",
+       "--C"},
+      {"hur_avg", "domain avg. relative humidity profile", "%", "--C"},
+      {"clw_avg", "domain avg. mass fraction of cloud liquid water profile",
+       "kg/kg", "--C"},
+      {"cli_avg", "domain avg. mass fraction of cloud ice profile", "kg/kg",
+       "--C"},
+      {"plw_avg",
+       "domain avg. mass fraction of precipitating liquid water profile",
+       "kg/kg", "--C"},
+      {"pli_avg", "domain avg. mass fraction of precipitating ice profile",
+       "kg/kg", "--C"},
+      {"theta_avg", "domain avg. potential temperature profile", "K", "--C"},
+      {"thetae_avg", "domain avg. equivalent potential temperature profile",
+       "K", "--C"},
+      {"tntrs_avg", "domain avg. shortwave radiative heating rate profile",
+       "K/s", "--C"},
+      {"tntrl_avg", "domain avg. longwave radiative heating rate profile",
+       "K/s", "--C"},
+      {"tntrscs_avg",
+       "domain avg. shortwave radiative heating rate profile – clear sky",
+       "K/s", "--C"},
+      {"tntrlcs_avg",
+       "domain avg. longwave radiative heating rate profile – clear sky", "K/s",
+       "--C"},
+      {"cldfrac_avg", "global cloud fraction profile", "%", "--C"},
+      // CMPI6 0D variables description
+      {"pr_avg", "domain avg. surface precipitation rate", "kg/(m^2 s)", "---"},
+      {"hfls_avg", "domain avg. surface upward latent heat flux", "W/m^2",
+       "---"},
+      {"hfss_avg", "domain avg. surface upward sensible heat flux", "W/m^2",
+       "---"},
+      {"prw_avg", "domain avg. water vapor path", "kg/m^2", "---"},
+      {"clwvi_avg", "domain avg. condensed water path", "kg/m^2", "---"},
+      {"clivi_avg", "domain avg. condensed ice water path", "kg/m^2", "---"},
+      {"spwr_avg", "domain avg. saturated water vapor path", "kg/m^2", "---"},
+      {"rlds_avg", "domain avg. surface downwelling longwave flux", "W/m^2",
+       "---"},
+      {"rlus_avg", "domain avg. surface upwelling longwave flux", "W/m^2",
+       "---"},
+      {"rsds_avg", "domain avg. surface downwelling shortwave flux", "W/m^2",
+       "---"},
+      {"rsus_avg", "domain avg. surface upward shortwave flux", "W/m^2", "---"},
+      {"rsdscs_avg",
+       "domain avg. surface downwelling shortwave flux – clear sky", "W/m^2",
+       "---"},
+      {"rsuscs_avg", "domain avg. surface upward shortwave flux – clear sky",
+       "W/m^2", "---"},
+      {"rldscs_avg",
+       "domain avg. surface downwelling longwave flux – clear sky", "W/m^2",
+       "---"},
+      {"rluscs_avg", "domain avg. surface upwelling longwave flux – clear sky",
+       "W/m^2", "---"},
+      {"rsdt_avg", "domain avg. TOA incoming shortwave flux", "W/m^2", "---"},
+      {"rsut_avg", "domain avg. TOA outgoing shortwave flux", "W/m^2", "---"},
+      {"rlut_avg", "domain avg. TOA outgoing longwave flux", "W/m^2", "---"},
+      {"rsutcs_avg", "domain avg. TOA outgoing shortwave flux – clear sky",
+       "W/m^2", "---"},
+      {"rlutcs_avg", "domain avg. TOA outgoing longwave flux – clear sky",
+       "W/m^2", "---"},
+      // CMPI6 2D variable description
+      {"pr", "surface precipitation rate", "kg/(m^2 s)", "CCC"},
+      {"pr_conv", "surface convective precipitation rate", "kg/(m^2 s)", "CCC"},
+      {"evspsbl", "evaporation flux", "kg/(m^2 s)", "CCC"},
+      {"hfls", "surface upward latent heat flux", "W/m^2", "CCC"},
+      {"hfss", "surface upward sensible heat flux", "W/m^2", "CCC"},
+      {"prw", "water vapor path", "kg/m^2", "CCC"},
+      {"clwvi", "condensed water path", "kg/m^2", "CCC"},
+      {"clivi", "condensed ice water path", "kg/m^2", "CCC"},
+      {"spwr", "saturated water vapor path", "kg/m^2", "CCC"},
+      {"rlds", "surface downwelling longwave flux", "W/m^2", "CCC"},
+      {"rlus", "surface upwelling longwave flux", "W/m^2", "CCC"},
+      {"rsds", "surface downwelling shortwave flux", "W/m^2", "CCC"},
+      {"rsus", "surface upward shortwave flux", "W/m^2", "CCC"},
+      {"rsdscs", "surface downwelling shortwave flux – clear sky", "W/m^2",
+       "CCC"},
+      {"rsuscs", "surface upward shortwave flux – clear sky", "W/m^2", "CCC"},
+      {"rldscs", "surface downwelling longwave flux – clear sky", "W/m^2",
+       "CCC"},
+      {"rluscs", "surface upwelling longwave flux – clear sky", "W/m^2", "CCC"},
+      {"rsdt", "TOA incoming shortwave flux", "W/m^2", "CCC"},
+      {"rsut", "TOA outgoing shortwave flux", "W/m^2", "CCC"},
+      {"rlut", "TOA outgoing longwave flux", "W/m^2", "CCC"},
+      {"rsutcs", "TOA outgoing shortwave flux – clear sky", "W/m^2", "CCC"},
+      {"rlutcs", "TOA outgoing longwave flux – clear sky", "W/m^2", "CCC"},
   };
 }