Replace the command line options by .ini files

README.md

@@ -10,7 +10,7 @@ The original code is found at https://github.com/earth-system-radiation/rte-rrtm
 Contacts: Robert Pincus and Eli Mlawer
 email: rrtmgp@aer.com
 
-This C++ implementation can be downloaded from https://github.com/earth-system-radiation/rte-rrtmgp-cpp
+This C++ implementation can be downloaded from https://github.com/microhh/rte-rrtmgp-cpp
 
 Contacts: Chiel van Heerwaarden and Menno Veerman
 
@@ -25,9 +25,11 @@ The source code of the testing executable in the `src_test` and
 `include_test` directory is released under the GPLv3 license,
 see https://www.gnu.org/licenses/gpl-3.0.en.html
 
+The .ini reader is based on TOML and uses toml11: https://github.com/ToruNiina/toml11.
+
 In order to check out the code including the `rte-rrtmgp` submodule, use:
 
-    git clone --recurse-submodules https://github.com/earth-system-radiation/rte-rrtmgp-cpp
+    git clone --recurse-submodules https://github.com/microhh/rte-rrtmgp-cpp
 
 In case you had already checked out the repository, use:
 

allsky/allsky.ini

@@ -0,0 +1,35 @@
+[switches]
+shortwave = true
+longwave = true
+fluxes = true
+cloud-optics = true
+aerosol-optics = false
+output-optical = false
+output-bnd-fluxes = false
+delta-cloud = true
+delta-aerosol = false
+
+# gpu-only (test_rte_rrtmgp_gpu)
+timings = false
+
+# raytracer (test_rte_rrtmgp_rt_gpu)
+sw-two-stream = false
+sw-raytracing = true
+lw-raytracing = true
+independent-column = false
+liq-cloud-optics = false
+ice-cloud-optics = false
+lw-scattering = false
+cloud-mie = false
+single-gpt = false
+profiling = false
+min-mfp-grid-ratio = true
+tica = false
+
+[ints]
+sw-raytracing = 256
+lw-raytracing = 22
+single-gpt = 1
+
+[floats]
+min-mfp-grid-ratio = 1.0

allsky/allsky_init.py

@@ -8,7 +8,7 @@
 n_bnd_lw = 16
 n_bnd_sw = 14
 
-nc_file = nc.Dataset("rte_rrtmgp_input.nc", mode="w", datamodel="NETCDF4", clobber=True)
+nc_file = nc.Dataset("allsky_input.nc", mode="w", datamodel="NETCDF4", clobber=True)
 
 # Radiation profiles.
 z_top = 70.e3

allsky/allsky_run.py

@@ -8,13 +8,13 @@ def remove(filename):
     except OSError:
         pass
 
-remove('rte_rrtmgp_output.nc')
-subprocess.run(['./test_rte_rrtmgp', '--cloud-optics'])
+remove('allsky_output.nc')
+subprocess.run(['./test_rte_rrtmgp', 'allsky'])
 
 remove('rrtmgp-allsky.nc')
 
 
-in_file = nc.Dataset("rte_rrtmgp_output.nc") 
+in_file = nc.Dataset("allsky_output.nc")
 out_file = nc.Dataset("rrtmgp-allsky.nc", "w")
 
 for dim_name in ["band_sw", "band_lw", "gpt_sw", "gpt_lw", "pair", "lay", "lev"]: