SMAPDA_OUTPUT/DAOBS Not Generating Files in LIS Despite Configuration Settings #1688
guesswhoam
started this conversation in
General
Replies: 0 comments
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Uh oh!
There was an error while loading. Please reload this page.
-
Hello LISF team,
I am working through the LISF Public Testcase Walkthrough (Step 6: LSM Data Assimilation Experiment) using the support/lisf-public-7.4 branch, but I am encountering an issue where the SMAPDA_OUTPUT/DAOBS/ directory is not generating any files, despite having configured the data assimilation diagnostics output settings in lis.config_noah36_smapda.
Environment
LISF Version: Compiled from support/lisf-public-7.4 branch.
System: Linux (CentOS/Ubuntu, [please specify your system if known]).
Compilers: [e.g., Intel Fortran Compiler 2021 or GNU gfortran, specify if known].
Dependencies: NetCDF, HDF5, ESMF installed with default settings.
Configuration Details
I have followed the tutorial and set up the configuration as follows:
Number of data assimilation instances: 1
Data assimilation algorithm: "EnKF"
Data assimilation set: "SMAP(NASA) soil moisture"
Data assimilation output interval for diagnostics: "1mn" # Originally "1da", modified to test
Data assimilation output ensemble spread: 1
Data assimilation output processed observations: 1
Data assimilation output innovations: 1
SMAP(NASA) soil moisture data directory: ./RS_DATA/SMAP/SPL3SMP.005
SMAP(NASA) soil moisture apply SMAP QC flags: 1 # Also tested with 0
Simulation period: 2017-04-30 06:00 to 2017-05-01 06:00.
Restart file: ./DA_ensrst/LIS_EnRST_NOAH36_201801010000.d01.nc.
Observed Behavior
The LIS run completes successfully (LIS Run completed in lislog_smapda.0000).
Outputs are generated in SMAPDA_OUTPUT/SURFACEMODEL/ (e.g., LIS_RST_NOAH36_201705010600.d01.nc) and SMAPDA_OUTPUT/DAPERT/ (e.g., LIS_DAPERT_201705010600.d01.bin).
SMAP data is read successfully (e.g., [INFO] Reading descending pass ./RS_DATA/SMAP/SPL3SMP.005/2017.04.30/SMAP_L3_SM_P_20170430_R16010_001.h5).
However, the DAOBS directory remains empty, and no diagnostics files (innovations, processed observations) are generated.
Logs
Here are relevant excerpts from lislog_smapda.0000:
[INFO] Successfully created observations state
[INFO] Created the States to hold the SMAP NASA observations data
[INFO] Reading descending pass ./RS_DATA/SMAP/SPL3SMP.005/2017.04.30/SMAP_L3_SM_P_20170430_R16010_001.h5
[INFO] Writing Perturbations restart ./SMAPDA_OUTPUT/DAPERT/201705/LIS_DAPERT_201705010600.d01.bin
[INFO] Noah-3.6 archive restart written: ./SMAPDA_OUTPUT/SURFACEMODEL/201705/LIS_RST_NOAH36_201705010600.d01.nc
[INFO] LIS Run completed.
Questions
Why is the DAOBS folder not being populated despite the diagnostics output settings being enabled?
Is there a specific condition (e.g., time alignment, spatial matching) that might prevent EnKF from performing updates and generating DAOBS files?
Could this be related to a compilation issue or version mismatch?
Additional Information
Model domain: 34.375°N to 39.625°N, -102.875°W to -96.125°W (from lis_input.nldas.noah36.d01.nc).
I can provide the full lislog_smapda.0000 or other files if needed.
Any insights or suggestions would be greatly appreciated! Thank you for your support.
Best regards,
Beta Was this translation helpful? Give feedback.
All reactions