Merge pull request #31 from NCAR/develop-21-no2-cross-sections

Add XNO2 photolysis configurations for TS mechanisms
NCAR · Jan 19, 2024 · f371af0 · f371af0
2 parents 2d9db9c + 15fd74e
commit f371af0
Show file tree

Hide file tree

Showing 27 changed files with 2,375 additions and 591 deletions.
diff --git a/data/cross_sections/BRONO2_JPL06.nc b/data/cross_sections/BRONO2_JPL06.nc
diff --git a/data/cross_sections/HO2NO2_JPL06.nc b/data/cross_sections/HO2NO2_JPL06.nc
diff --git a/data/cross_sections/HO2NO2_temp_JPL06.nc b/data/cross_sections/HO2NO2_temp_JPL06.nc
diff --git a/examples/tuv_5_4.json b/examples/tuv_5_4.json
@@ -195,15 +195,82 @@
          {
             "name": "HNO4+hv->HO2+NO2",
             "cross section": {
-               "netcdf files": [
-                 { "file path": "data/cross_sections/HNO4_1.nc" }
-               ],
-               "type": "base"
-            },
-            "quantum yield": {
+               "type": "temperature based",
+               "netcdf file": "data/cross_sections/HO2NO2_JPL06.nc",
+               "parameterization": {
+                 "type": "BURKHOLDER",
+                 "netcdf file": {
+                   "file path": "data/cross_sections/HO2NO2_temp_JPL06.nc"
+                 },
+                 "A": -988.0,
+                 "B": 0.69,
+                 "temperature ranges": [
+                   {
+                     "maximum": 279.999999999999,
+                     "fixed value": 280.0
+                   },
+                   {
+                     "minimum": 280.0,
+                     "maximum": 350.0
+                   },
+                   {
+                     "minimum": 350.00000000001,
+                     "fixed value": 350.0
+                   }
+                 ]
+               }
+             },
+             "quantum yield": {
                "type": "base",
-               "constant value": 1.0
-            }
+               "constant value": 0.70,
+               "override bands": [
+                 {
+                   "band": "range",
+                   "minimum wavelength": 200.0,
+                   "value": 0.80
+                 }
+               ]
+             }
+         },
+         {
+            "name": "HNO4+hv->OH+NO3",
+            "cross section": {
+               "type": "temperature based",
+               "netcdf file": "data/cross_sections/HO2NO2_JPL06.nc",
+               "parameterization": {
+                 "type": "BURKHOLDER",
+                 "netcdf file": {
+                   "file path": "data/cross_sections/HO2NO2_temp_JPL06.nc"
+                 },
+                 "A": -988.0,
+                 "B": 0.69,
+                 "temperature ranges": [
+                   {
+                     "maximum": 279.999999999999,
+                     "fixed value": 280.0
+                   },
+                   {
+                     "minimum": 280.0,
+                     "maximum": 350.0
+                   },
+                   {
+                     "minimum": 350.00000000001,
+                     "fixed value": 350.0
+                   }
+                 ]
+               }
+             },
+             "quantum yield": {
+               "type": "base",
+               "constant value": 0.30,
+               "override bands": [
+                 {
+                   "band": "range",
+                   "minimum wavelength": 200.0,
+                   "value": 0.20
+                 }
+               ]
+             }
          },
          {
             "name": "O3+hv->O2+O(1D)",
@@ -957,11 +1024,29 @@
          {
             "name": "BrONO2+hv->BrO+NO2",
             "cross section": {
-               "netcdf files": [
-                 { "file path": "data/cross_sections/BrONO2_1.nc" }
-               ],
-               "type": "base"
-            },
+               "type": "temperature based",
+               "parameterization": {
+                 "type": "TAYLOR_SERIES",
+                 "netcdf file": {
+                   "file path": "data/cross_sections/BRONO2_JPL06.nc"
+                 },
+                 "base temperature": 296.0,
+                 "temperature ranges": [
+                   {
+                     "maximum": 199.999999999999,
+                     "fixed value": 200.0
+                   },
+                   {
+                     "minimum": 200.0,
+                     "maximum": 296.0
+                   },
+                   {
+                     "minimum": 296.00000000001,
+                     "fixed value": 296.0
+                   }
+                 ]
+               }
+             },
             "quantum yield": {
                "type": "base",
                "constant value": 0.15
@@ -970,11 +1055,29 @@
          {
             "name": "BrONO2+hv->Br+NO3",
             "cross section": {
-               "netcdf files": [
-                 { "file path": "data/cross_sections/BrONO2_1.nc" }
-               ],
-               "type": "base"
-            },
+               "type": "temperature based",
+               "parameterization": {
+                 "type": "TAYLOR_SERIES",
+                 "netcdf file": {
+                   "file path": "data/cross_sections/BRONO2_JPL06.nc"
+                 },
+                 "base temperature": 296.0,
+                 "temperature ranges": [
+                   {
+                     "maximum": 199.999999999999,
+                     "fixed value": 200.0
+                   },
+                   {
+                     "minimum": 200.0,
+                     "maximum": 296.0
+                   },
+                   {
+                     "minimum": 296.00000000001,
+                     "fixed value": 296.0
+                   }
+                 ]
+               }
+             },
             "quantum yield": {
                "type": "base",
                "constant value": 0.85

diff --git a/src/cross_sections/CMakeLists.txt b/src/cross_sections/CMakeLists.txt
@@ -31,4 +31,6 @@ target_sources(tuvx_object
     rayliegh.F90
 )
 
+add_subdirectory(util)
+
 ################################################################################