rc-qlc 0.3.24__cp311-cp311-macosx_11_0_arm64.whl → 0.3.26__cp311-cp311-macosx_11_0_arm64.whl

qlc/cli/__init__.py

@@ -2,20 +2,111 @@
 import os
 import subprocess
 import sys
+from datetime import datetime
+from pathlib import Path
 
 def run_shell_driver():
+    """
+    Finds and executes qlc_main.sh, capturing its output for logging.
+    This acts as the entry point for the 'qlc' command.
+    """
     # Correctly locate the 'sh' directory relative to the package installation
     sh_dir = os.path.join(os.path.dirname(__file__), '..', 'sh')
     script = os.path.join(sh_dir, "qlc_main.sh")
-    subprocess.run(["bash", script] + sys.argv[1:])
 
-def run_python_driver():
-    # This function seems to be for a different purpose, ensure qlc_main is correct
-    from qlc.cli.qlc_main import main
-    main()
+    # Determine log directory from QLC_HOME or fall back to a default
+    # This mirrors the logic that the shell scripts would use.
+    qlc_home_str = os.environ.get("QLC_HOME", f"{os.path.expanduser('~')}/qlc")
+    log_dir = os.path.join(qlc_home_str, "log")
+    os.makedirs(log_dir, exist_ok=True)
+
+    # Create a timestamped log file for the shell script's output
+    timestamp = datetime.now().strftime("%Y%m%d_%H%M%S")
+    log_file_path = os.path.join(log_dir, f"qlc_shell_main_{timestamp}.log")
+    print(f"[QLC Wrapper] Logging shell script output to: {log_file_path}")
+
+    try:
+        # Use a list of strings for Popen
+        command = [str(script)] + sys.argv[1:]
+        
+        with open(log_file_path, 'w', encoding='utf-8') as log_file:
+            process = subprocess.Popen(
+                command,
+                stdout=subprocess.PIPE,
+                stderr=subprocess.STDOUT,
+                text=True,
+                bufsize=1, # Line-buffered
+                universal_newlines=True
+            )
+            
+            # Real-time stream processing
+            for line in process.stdout:
+                # Write to file without adding a newline, as 'line' already has one
+                log_file.write(line)
+                # Print to console, stripping the newline to avoid double spacing
+                sys.stdout.write(line)
+            
+            process.wait()
+
+        if process.returncode != 0:
+            print(f"\n[ERROR] Shell script exited with non-zero code: {process.returncode}", file=sys.stderr)
+            sys.exit(process.returncode)
+
+    except FileNotFoundError:
+        print(f"Error: Could not find the qlc_main.sh script at {script}", file=sys.stderr)
+        sys.exit(1)
+    except Exception as e:
+        print(f"\n[ERROR] An unexpected error occurred: {e}", file=sys.stderr)
+        sys.exit(1)
+
 
 def run_batch_driver():
-    # Correctly locate the 'sh' directory relative to the package installation
-    sh_dir = os.path.join(os.path.dirname(__file__), '..', 'sh')
-    script = os.path.join(sh_dir, "qlc_start_batch.sh")
-    subprocess.run(["bash", script] + sys.argv[1:])
+    """
+    Finds and executes qlc_batch.sh, capturing its output for logging.
+    This acts as the entry point for the 'sqlc' command.
+    """
+    try:
+        sh_dir = Path(__file__).resolve().parent.parent / "sh"
+        script_path = sh_dir / "qlc_batch.sh"
+        if not script_path.is_file():
+            print(f"[ERROR] Batch script not found at: {script_path}", file=sys.stderr)
+            sys.exit(1)
+
+        # Ensure the script is executable
+        script_path.chmod(script_path.stat().st_mode | 0o111)
+
+        # Determine log directory from config
+        home = Path.home()
+        log_dir_str = os.environ.get("QLC_LOG_DIR", str(home / "qlc" / "log"))
+        log_dir = Path(log_dir_str)
+        log_dir.mkdir(parents=True, exist_ok=True)
+        
+        timestamp = datetime.now().strftime("%Y%m%d_%H%M%S")
+        log_file_path = log_dir / f"sqlc_shell_main_{timestamp}.log"
+        print(f"[QLC Batch Wrapper] Logging shell script output to: {log_file_path}")
+
+        command = [str(script_path)] + sys.argv[1:]
+        
+        with open(log_file_path, 'w', encoding='utf-8') as log_file:
+            process = subprocess.Popen(
+                command,
+                stdout=subprocess.PIPE,
+                stderr=subprocess.STDOUT,
+                text=True,
+                bufsize=1,
+                universal_newlines=True
+            )
+            
+            for line in process.stdout:
+                log_file.write(line)
+                sys.stdout.write(line)
+            
+            process.wait()
+
+        if process.returncode != 0:
+            print(f"\n[ERROR] Batch script exited with non-zero code: {process.returncode}", file=sys.stderr)
+            sys.exit(process.returncode)
+
+    except Exception as e:
+        print(f"\n[ERROR] An unexpected error occurred: {e}", file=sys.stderr)
+        sys.exit(1)

qlc/cli/installer.py

@@ -3,14 +3,32 @@ from qlc.install import setup
 from qlc.py.version import QLC_VERSION
 
 def main():
-    parser = argparse.ArgumentParser(description="QLC Installer")
+    parser = argparse.ArgumentParser(
+        description="Install the QLC runtime environment.",
+        formatter_class=argparse.RawTextHelpFormatter
+    )
     parser.add_argument(
         "--mode",
-        choices=["cams", "test", "interactive"],
-        default="test",
-        help="Installation mode: test (for tests), cams (for CAMS), or interactive (development)"
+        type=str,
+        choices=['test', 'cams', 'interactive'],
+        required=True,
+        help="The installation mode.\n"
+             "'test': A standalone mode with bundled example data.\n"
+             "'cams': An operational mode for CAMS environments.\n"
+             "'interactive': A mode for developers to use a custom config."
     )
+    parser.add_argument("--version", type=str, help="Override QLC version (for development)")
+    parser.add_argument("--config", type=str, help="[interactive mode only] Path to a custom config file.")
+
+
     args = parser.parse_args()
 
     print(f"[QLC-INSTALL] Installing QLC version {QLC_VERSION} in '{args.mode}' mode")
-    setup(mode=args.mode, version=QLC_VERSION)
+    if args.mode == 'interactive' and not args.config:
+        parser.error("--config is required when using --mode interactive")
+    
+    # Use a dummy version if not provided, as the setup script will find the real one
+    version = args.version if args.version else "0.0.0"
+    
+    setup(mode=args.mode, version=version, config_file=args.config)
+    return 0

qlc/cli/qlc_main.py

@@ -1,3 +1,7 @@
+# --- Set Matplotlib backend before any other imports ---
+import matplotlib
+matplotlib.use('Agg')
+
 import os
 import sys
 import json
@@ -6,16 +10,22 @@ import logging
 import argparse
 import traceback
 import multiprocessing
+import pandas as pd
 import concurrent.futures
 from datetime import datetime
 from multiprocessing import current_process
 from concurrent.futures import ThreadPoolExecutor, as_completed
 from qlc.py.control import run_main_processing, process_multiple_configs_sp, process_single_config
-from qlc.py.utils import expand_paths, get_timestamp, setup_logging, validate_paths, merge_global_attributes
+from qlc.py.utils import expand_paths, get_timestamp, validate_paths, merge_global_attributes
 from qlc.py.version import QLC_VERSION, QLC_RELEASE_DATE, QLC_DISTRIBUTION
 from qlc.py.plugin_loader import try_import_plugin_module
-from qlc.py.logging_utils import log_input_availability
-from qlc.py.logging_utils import log_qlc_banner
+from qlc.py.logging_utils import log_input_availability, log_qlc_banner, setup_logging
+
+# Globally ignore the ChainedAssignmentError FutureWarning by its message.
+# This is the most robust way to handle this persistent, non-critical warning
+# that arises from the complex interaction between pandas and Cython.
+import warnings
+warnings.filterwarnings('ignore', '.*ChainedAssignmentError.*')
 
 # -----------------------------------------------
 # QLC main controller (Master Orchestration)
@@ -94,15 +104,7 @@ def load_config_with_defaults(config):
 
 def run_single_config(config_entry, idx, total_configs):
     try:
-        # Minimal logging setup for subprocesses only
-        if current_process().name != "MainProcess":
-            if not logging.getLogger().hasHandlers():
-                logging.basicConfig(
-                    level=logging.INFO,
-                    format='[%(asctime)s] %(levelname)s: [%(processName)s] %(message)s',
-                    handlers=[logging.StreamHandler(sys.stdout)]
-                )
-
+        # Use the main logger configured in run_with_file, no separate config needed here
         logging.info(f"(Process {idx+1}/{total_configs}) Starting configuration '{config_entry.get('name', f'config_{idx+1}')}'...")
 
         config = load_config_with_defaults(config_entry)
@@ -112,7 +114,11 @@ def run_single_config(config_entry, idx, total_configs):
             logging.warning("⚠️ Nothing to process: no model or observation input available. Skipping execution.")
             return
 
-        run_main_processing(config)
+        # Suppress persistent ChainedAssignmentError FutureWarning from pandas.
+        # This is a known issue with the Cython execution path where standard
+        # fixes (.copy(), .loc) are not sufficient. The underlying code is correct.
+        with pd.option_context('mode.chained_assignment', None):
+            run_main_processing(config)
 
         logging.info(f"(Process {idx+1}/{total_configs}) Finished configuration '{config_entry.get('name', f'config_{idx+1}')}'.")
     except Exception as e:
@@ -121,32 +127,48 @@ def run_single_config(config_entry, idx, total_configs):
 
 def run_with_file(file_path):
     start_time = time.time()
-    print("************************************************************************************************")
-    print(f"Start execution: {time.strftime('%Y-%m-%d %H:%M:%S', time.localtime(start_time))}")
-    if not os.path.isfile(file_path):
-        print(f"ERROR: Input file '{file_path}' not found.")
+    
+    config_data = None
+    try:
+        if file_path == '-':
+            # Read from stdin if the file path is '-'
+            logging.debug("Reading configuration from stdin.")
+            config_data = expand_paths(json.load(sys.stdin))
+        else:
+            if not os.path.isfile(file_path):
+                # Using print because logging may not be set up if config fails to load.
+                print(f"ERROR: Input file '{file_path}' not found.")
+                sys.exit(1)
+            with open(file_path, 'r') as f:
+                config_data = expand_paths(json.load(f))
+    except json.JSONDecodeError as e:
+        print(f"ERROR: Failed to parse JSON from {'stdin' if file_path == '-' else file_path}: {e}")
         sys.exit(1)
-
-    with open(file_path, 'r') as f:
-        try:
-            config_data = expand_paths(json.load(f))
-            validate_paths(config_data if isinstance(config_data, dict) else config_data[0])
-            print(json.dumps(config_data, indent=2))
-        except json.JSONDecodeError as e:
-            print(f"ERROR: Failed to parse JSON: {e}")
-            sys.exit(1)
-
+            
+    # --- Setup logging as early as possible ---
     if isinstance(config_data, list):
         first_config = config_data[0]
     else:
         first_config = config_data
-    
+        
     log_dir = first_config.get("logdir", "~/qlc/log")
-    log_filename = os.path.join(log_dir, f"qlc_{get_timestamp()}.log")
+    log_filename = f"qlc_{get_timestamp()}.log"
     full_log_path = os.path.join(os.path.expanduser(log_dir), log_filename)
-    log_level = logging.DEBUG if first_config.get("debug", False) else logging.INFO
-#   setup_logging(log_dir=os.path.dirname(full_log_path), level=log_level)   
-    setup_logging(log_dir=log_dir, level=log_level)
+    setup_logging(log_file_path=full_log_path, debug=first_config.get("debug", False))
+    
+    # --- Now, use the logger for all output ---
+    logging.info("************************************************************************************************")
+    logging.info(f"Start execution: {time.strftime('%Y-%m-%d %H:%M:%S', time.localtime(start_time))}")
+    
+    try:
+        validate_paths(first_config)
+        # Pretty print JSON to a string, then log it
+        config_str = json.dumps(config_data, indent=2)
+        for line in config_str.splitlines():
+            logging.info(line)
+    except Exception as e:
+        logging.error(f"Configuration error: {e}")
+        sys.exit(1)
 
 #   If the config is a list, loop over multiple configs
     if isinstance(config_data, list) and len(config_data) == 1:

qlc/cli/qlc_py_main.py

@@ -7,67 +7,72 @@ from pathlib import Path
 def main():
     """
     Python entry point for running the main QLC processing script.
+    This script now acts as a silent wrapper. All console output
+    is handled by the logging configuration in the downstream qlc_main module.
     """
-    script_name = sys.argv[0]
-    print("________________________________________________________________________________________")
-    print(f"Start {script_name} at {datetime.now()}")
-    print("----------------------------------------------------------------------------------------")
-
-    # --- Change to QLC home directory ---
+    # --- Ensure QLC home directory exists ---
     qlc_home = Path.home() / "qlc"
+    # The qlc-install script creates this directory.
+    # If it's not here, the user needs to run the installer.
     if not qlc_home.is_dir():
+        # Use print here because logging is not yet configured and this is a pre-flight check.
         print(f"[ERROR] QLC home directory not found at: {qlc_home}")
         print("Please run 'qlc-install' to set up the required structure.")
         sys.exit(1)
-    
-    os.chdir(qlc_home)
-    print(f"Changed working directory to: {os.getcwd()}")
-    # ---
 
-    config_path = sys.argv[1] if len(sys.argv) > 1 else str(Path.home() / "qlc" / "config" / "json" / "qlc_config.json")
-    log_dir = qlc_home / "log"
-    log_dir.mkdir(parents=True, exist_ok=True)
+    os.chdir(qlc_home)
 
-    timestamp = datetime.now().strftime("%Y%m%d_%H%M%S")
-    log_file = log_dir / f"qlc_run_{timestamp}.log"
+    # Determine config path. Default to the standard location if no args are given.
+    # Look for a --config flag.
+    config_arg_present = "--config" in sys.argv
+    stdin_mode = False
+    config_path_str = str(qlc_home / "config" / "json" / "qlc_config.json") # Default
 
-    print("************************************************************************************************")
-    print(f"Running {script_name} with output logged to {log_file}")
+    if config_arg_present:
+        config_idx = sys.argv.index("--config")
+        if len(sys.argv) > config_idx + 1:
+            config_value = sys.argv[config_idx + 1]
+            if config_value == '-':
+                stdin_mode = True
+                config_path_str = '-' # Pass the stdin marker to the subprocess
+            else:
+                config_path_str = config_value # A specific file was passed
+        else:
+            # This case should be caught by argparse in the child, but we can be safe.
+            print("[ERROR] --config flag was provided without a value.", file=sys.stderr)
+            sys.exit(1)
 
-    # The command will use the same Python interpreter that is running this script
     python_executable = sys.executable
     command = [
         python_executable,
         "-m",
         "qlc.cli.qlc_main",
         "--config",
-        config_path,
+        config_path_str,
     ]
     
-    print(f"Executing command: {' '.join(command)}")
-
     try:
-        with open(log_file, "w") as lf:
-            process = subprocess.Popen(
+        # If in stdin mode, we need to pass our stdin to the subprocess.
+        # Otherwise, the subprocess runs without piped data.
+        if stdin_mode:
+            subprocess.run(
                 command,
-                stdout=subprocess.PIPE,
-                stderr=subprocess.STDOUT,
+                check=True,
                 text=True,
-                bufsize=1,
-                universal_newlines=True,
+                stdin=sys.stdin  # Pass stdin through
             )
-            if process.stdout:
-                for line in process.stdout:
-                    print(line, end="")
-                    lf.write(line)
-            process.wait()
-
+        else:
+            subprocess.run(
+                command,
+                check=True,
+                text=True,
+            )
+    except subprocess.CalledProcessError as e:
+        # The error output from the subprocess will have already been printed.
+        sys.exit(e.returncode)
     except Exception as e:
-        print(f"An error occurred: {e}")
-
-    print("________________________________________________________________________________________")
-    print(f"End   {script_name} at {datetime.now()}")
-    print("________________________________________________________________________________________")
+        print(f"A critical error occurred in the QLC wrapper: {e}")
+        sys.exit(1)
 
 if __name__ == "__main__":
     main()

qlc/config/json/qlc_config.json

@@ -1,19 +1,20 @@
 [
   {
-    "name": "CAMS2_35",
+    "name": "CAMS",
     "logdir": "./log",
     "workdir": "./run",
-    "output_base_name": "$HOME/qlc/output/GMD",
+    "output_base_name": "$HOME/qlc/Plots/PY",
     "station_file": "$HOME/qlc/obs/data/ebas_station-locations.csv",
     "obs_path": "$HOME/qlc/obs/data/ver0d",
     "obs_dataset_type": "ebas_daily",
-    "start_date": "2018-01-01",
-    "end_date": "2018-12-31",
-    "variable": "SO2,SO4,HNO3,NO3,NH3,NH4",
-    "station_radius_deg": 10.0,
+    "obs_dataset_version": "latest",
+    "start_date": "2018-12-01",
+    "end_date": "2018-12-21",
+    "variable": "NH3,NH4_as",
+    "station_radius_deg": 0.5,
     "plot_type": "",
     "plot_region": "EU",
-    "time_average": "daily, weekly, monthly",
+    "time_average": "daily",
     "station_plot_group_size": 5,
     "show_stations": false,
     "show_min_max": true,
@@ -26,15 +27,98 @@
     "show_station_timeseries_com": false,
     "save_plot_format": "pdf",
     "save_data_format": "nc",
-    "multiprocessing": true,
+    "multiprocessing": false,
     "n_threads": "20",
     "debug": false,
 
     "global_attributes": {
       "title": "Air pollutants over Europe, SO2,SO4,HNO3,NO3,NH3,NH4",
       "summary": "Custom summary for netCDF output: Ebas daily observations for selected EU stations.",
-      "author": "Plain & Simple, plain.simple@example.com",
-      "history": "Processed for CAMS (test 1.1)",
+      "author": "Swen Metzger, sm@researchconcepts.io",
+      "history": "Processed for CAMS2_35bis (qlc_v0.3.25)",
+      "Conventions": "CF-1.8"
+    }
+  },
+  {
+    "name": "CAMS",
+    "logdir": "./log",
+    "workdir": "./run",
+    "output_base_name": "$HOME/qlc/Plots/PY",
+    "station_file": "$HOME/qlc/obs/data/ebas_station-locations.csv",
+    "mod_path": "$HOME/qlc/Analysis/",
+    "model": "ifs",
+    "experiments": "b2ro,b2rn",
+    "start_date": "2018-12-01",
+    "end_date": "2018-12-21",
+    "variable": "NH3,NH4_as",
+    "station_radius_deg": 0.5,
+    "model_level": 9,
+    "plot_type": "",
+    "plot_region": "EU",
+    "time_average": "daily",
+    "station_plot_group_size": 5,
+    "show_stations": false,
+    "show_min_max": true,
+    "log_y_axis": false,
+    "fix_y_axis": true,
+    "show_station_map": true,
+    "show_station_timeseries_obs": false,
+    "show_station_timeseries_mod": true,
+    "show_station_timeseries_com": false,
+    "save_plot_format": "pdf",
+    "save_data_format": "nc",
+    "multiprocessing": false,
+    "n_threads": "20",
+    "debug": false,
+
+    "global_attributes": {
+      "title": "Air pollutants over EU",
+      "summary": "Output of b2ro,iqi9 / ifs for selected EU stations.",
+      "author": "Swen Metzger, sm@researchconcepts.io",
+      "history": "Processed for CAMS2_35bis (qlc_v0.3.25)",
+      "Conventions": "CF-1.8"
+    }
+  },
+  {
+    "name": "CAMS",
+    "logdir": "./log",
+    "workdir": "./run",
+    "output_base_name": "$HOME/qlc/Plots/PY",
+    "station_file": "$HOME/qlc/obs/data/ebas_station-locations.csv",
+    "obs_path": "$HOME/qlc/obs/data/ver0d",
+    "obs_dataset_type": "ebas_daily",
+    "mod_path": "$HOME/qlc/Analysis/",
+    "model": "ifs",
+    "experiments": "b2ro,b2rn",
+    "exp_labels": "E4C,REF",
+    "start_date": "2018-12-01",
+    "end_date": "2018-12-21",
+    "variable": "NH3,NH4_as",
+    "station_radius_deg": 0.5,
+    "plot_type": "",
+    "plot_region": "EU",
+    "time_average": "daily",
+    "station_plot_group_size": 5,
+    "show_stations": false,
+    "show_min_max": true,
+    "log_y_axis": false,
+    "fix_y_axis": true,
+    "show_station_map": true,
+    "load_station_timeseries_obs": false,
+    "show_station_timeseries_obs": false,
+    "show_station_timeseries_mod": false,
+    "show_station_timeseries_com": true,
+    "save_plot_format": "pdf",
+    "save_data_format": "nc",
+    "multiprocessing": false,
+    "n_threads": "20",
+    "debug": false,
+
+    "global_attributes": {
+      "title": "Air pollutants over Europe, PM2.5,NH3 (CAMAERA-D3.2-V1.3_M7_in_IFS)",
+      "summary": "EBAS daily observations for selected EU stations.",
+      "author": "Swen Metzger, sm@researchconcepts.io",
+      "history": "Processed for CAMS2_35bis (qlc_v0.3.25)",
       "Conventions": "CF-1.8"
     }
   }

qlc/config/nml/mars_A1_sfc.nml

@@ -1,8 +1,8 @@
 # adopt to your needs, except for EXP, SDATE, EDATE, MYPATH, MYFILE, which will be auto-replaced
 # add your param_, ncvar_ and myvar_ to your $HOME/qlc/qlc.conf file
-# param_A=("207.210" "72.210"  "73.210"  "74.210")
-# ncvar_A=("var207"  "var72"   "var73"   "var74")
-# myvar_A=("AOD"     "PM1"     "PM2.5"   "PM10")
+# param_A1_sfc=("73.210")
+# ncvar_A1_sfc=("var73")
+# myvar_A1_sfc=("PM25")
 # surface level
 RETRIEVE,
 CLASS      = XCLASS,
@@ -13,7 +13,6 @@ PARAM      = 73.210,
 LEVTYPE    = sfc,
 GRID       = 1.0/1.0,
 TIME       = 00:00:00,
-STEP       = 0/6/12/18,
-LEVELIST   = 137,
+STEP       = 12,
 DATE       = SDATE/to/EDATE,
 target     = "MYPATH/MYFILE_A1_sfc.grb"

qlc/config/nml/mars_A3_sfc.nml

@@ -14,6 +14,5 @@ LEVTYPE    = sfc,
 GRID       = 1.0/1.0,
 TIME       = 00:00:00,
 STEP       = 0/6/12/18,
-LEVELIST   = 137,
 DATE       = SDATE/to/EDATE,
 target     = "MYPATH/MYFILE_A3_sfc.grb"

qlc/config/nml/mars_B1_pl.nml

@@ -13,7 +13,7 @@ PARAM      = 19.217,
 LEVTYPE    = pl,
 GRID       = 1.0/1.0,
 TIME       = 00:00:00,
-STEP       = 0/6/12/18,
-LEVELIST   = 10/20/30/50/70/100/150/200/250/300/400/500/700/850/925/1000,
+STEP       = 12,
+LEVELIST   = 10/50/100/150/200/400/700/850/925/1000,
 DATE       = SDATE/to/EDATE,
 target     = "MYPATH/MYFILE_B1_pl.grb"

qlc/{examples/cams_case_1/config/nml/mars_A3_sfc.nml → config/nml/mars_B1_sfc.nml}

@@ -1,19 +1,18 @@
 # adopt to your needs, except for EXP, SDATE, EDATE, MYPATH, MYFILE, which will be auto-replaced
-# add your param_, ncvar_ and myvar_ to your $HOME/qlc/config/qlc.conf file
-# param_A3-sfc=("203.210")
-# ncvar_A3-sfc=("go3")
-# myvar_A3-sfc=("O3")
+# add your param_, ncvar_ and myvar_ to your $HOME/qlc/qlc.conf file
+# param_B1_sfc=("19.217")
+# ncvar_B1_sfc=("nh3")
+# myvar_B1_sfc=("NH3")
 # surface level
 RETRIEVE,
 CLASS      = XCLASS,
 TYPE       = fc, 
 STREAM     = oper,
 EXPVER     = EXP,
-PARAM      = 203.210,
+PARAM      = 19.217,
 LEVTYPE    = sfc,
 GRID       = 1.0/1.0,
 TIME       = 00:00:00,
-STEP       = 0/6/12/18,
-LEVELIST   = 137,
+STEP       = 12,
 DATE       = SDATE/to/EDATE,
-target     = "MYPATH/MYFILE_A3_sfc.grb"
+target     = "MYPATH/MYFILE_B1_sfc.grb"

qlc/config/nml/mars_C1_pl.nml

@@ -3,7 +3,7 @@
 # param_C2=("36.212" "35.212" "2.212" "1.212")
 # ncvar_C2=("var36"  "var35"  "var2"  "var1" )
 # myvar_C2=("NO3_as" "NH4_as" "SO4_as" "N_as")
-# model levels
+# pressure levels
 RETRIEVE,
 CLASS      = XCLASS,
 TYPE       = fc, 

qlc/{examples/cams_case_1/config/nml/mars_C1_pl.nml → config/nml/mars_C1_sfc.nml}

@@ -1,19 +1,18 @@
 # adopt to your needs, except for EXP, SDATE, EDATE, MYPATH, MYFILE, which will be auto-replaced
 # add your param_, ncvar_ and myvar_ to your $HOME/qlc/qlc.conf file
-# param_C2=("36.212" "35.212" "2.212" "1.212")
-# ncvar_C2=("var36"  "var35"  "var2"  "var1" )
-# myvar_C2=("NO3_as" "NH4_as" "SO4_as" "N_as")
-# model levels
+# param_C1_sfc=("35.212")
+# ncvar_C1_sfc=("var35")
+# myvar_C1_sfc=("NH4_as")
+# surface level
 RETRIEVE,
 CLASS      = XCLASS,
 TYPE       = fc, 
 STREAM     = oper,
 EXPVER     = EXP,
 PARAM      = 35.212,
-LEVTYPE    = pl,
+LEVTYPE    = sfc,
 GRID       = 1.0/1.0,
 TIME       = 00:00:00,
-STEP       = 0/6/12/18,
-LEVELIST   = 10/20/30/50/70/100/150/200/250/300/400/500/700/850/925/1000,
+STEP       = 12,
 DATE       = SDATE/to/EDATE,
-target     = "MYPATH/MYFILE_C1_pl.grb"
+target     = "MYPATH/MYFILE_C1_sfc.grb"

qlc/config/nml/mars_C2_pl.nml

@@ -3,7 +3,7 @@
 # param_C2=("36.212" "35.212" "2.212" "1.212")
 # ncvar_C2=("var36"  "var35"  "var2"  "var1" )
 # myvar_C2=("NO3_as" "NH4_as" "SO4_as" "N_as")
-# model levels
+# pressure levels
 RETRIEVE,
 CLASS      = XCLASS,
 TYPE       = fc,