SMAPDA_OUTPUT/DAOBS Not Generating Files in LIS Despite Configuration Settings #1689
Unanswered
guesswhoam
asked this question in
Q&A
Replies: 0 comments
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
-
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