Issue reading ISMN data into LVT #1685

alexamarco · 2025-02-13T14:56:21Z

alexamarco
Feb 13, 2025

Hello All! I'm running into some trouble reading ISMN data into LVT 7.5.4. I am running LVT with LIS output as the first datastream and ISMN as the second, and trying to compare soil moisture between the two. Based on my output data, it looks like LVT is reading in the LIS Output data, but not the ISMN data. I have anom, mean, stdev, and variance output for the LIS Output (datastream 1), but null for ISMN (datastream 2). The acorr, armse, bias, mae, rmse, and ubrmse are all missing values because they require both datastream 1 and 2 to be calculated. However, in the logfile it looks like each ISMN file is being read in.

Has anyone experienced a similar issue before? I have included my lvt.config information below and attached logfile and an example output file for reference.

LVT running mode: "Data intercomparison"
Map projection of the LVT analysis: "latlon"
LVT output format: "netcdf"
LVT output methodology: "2d gridspace"
Analysis data sources: "LIS output" "ISMN"

Start mode: coldstart
LVT output restart file: 1
LVT restart output interval: "1mo"
LVT restart filename: LVT_RST
Starting year: 2015
Starting month: 3
Starting day: 31
Starting hour: 0
Starting minute: 0
Starting second: 0
Ending year: 2015
Ending month: 12
Ending day: 31
Ending hour: 21
Ending minute: 00
Ending second: 0
LVT clock timestep: "3hr"
Undefined value: -9999
LVT diagnostic file: ./LOGDIR/LVT/lvtlog_ismn

#LIS domain
Run domain lower left lat: 25.00
Run domain lower left lon: -124.00
Run domain upper right lat: 50.00
Run domain upper right lon: -67.00
Run domain resolution (dx): 0.1
Run domain resolution (dy): 0.1

Input domain and mask data file: ./LDTOUT/lis_input.10km_noahmp.nc

LVT datastream attributes table::
SoilMoist 1 1 m3/m3 - 0 4 SoilMoist 1 1 m3/m3 - 0 1
::

LVT surface soil layer thickness: 0.10
LVT root zone soil layer thickness: 1.0

#stats specifications (only applicable to retrospective mode)
Apply external mask: 0
External mask directory: none #OBSMASK.bin
Compute information theory metrics: 0
Compute ensemble metrics: 0
Metrics attributes file: ./LVT/METRICS.TBL
Observation count threshold: 0
Metrics computation frequency: "3hr"
Spatial averaging mode: "pixel-by-pixel"
Starting month if a shifted year definition is used in temporal averaging: 1
Metrics output directory: ./LVT/STATS.ISMN #STATS.SMAP.ARS
Metrics output frequency: "1da"
Time series location file: ./LVT/TS_LOCATIONS_ISMN_CONUS_2015_limited.TXT
Variable-based stratification: 0
Compute LSM diagnostics: 0
Confidence interval (%): 95
External data-based stratification: 0
Stratification attributes file: none
Seasonal cycle minimum count threshold: 0
Seasonal cycle interval type: "monthly"
Average diurnal cycle minimum count threshold: 0
Apply temporal smoothing to obs: 0
Averaging window for computing mean values in anomaly calculations: 'monthly'

ISMN observation directory: ../ISMN_limited

LIS output analysis data class: "LSM"
LIS output number of surface model types: 2
LIS output surface model types: "LSM" "Openwater"
LIS output model name: "Noah-MP.4.0.1"
LIS output domain and parameter file: ./LDTOUT/lis_input.10km_noahmp.nc
LIS output directory: ./OUTDIR/NoahMP/NoahAnalysis1
LIS output naming style: "3 level hierarchy"
LIS output methodology: "2d gridspace"
LIS output format: "netcdf"
LIS output interval: "3hr"
LIS output attributes file: ./configs/MODEL_OUTPUT_LIST.TBL
LIS output maximum number of surface type tiles per grid: 1
LIS output minimum cutoff percentage (surface type tiles): 0.05
LIS output maximum number of soil texture tiles per grid: 1
LIS output minimum cutoff percentage (soil texture tiles): 0.05
LIS output maximum number of soil fraction tiles per grid: 1
LIS output minimum cutoff percentage (soil fraction tiles): 0.05
LIS output maximum number of elevation bands per grid: 1
LIS output minimum cutoff percentage (elevation bands): 0.05
LIS output maximum number of slope bands per grid: 1
LIS output minimum cutoff percentage (slope bands): 0.05
LIS output maximum number of aspect bands per grid: 1
LIS output minimum cutoff percentage (aspect bands): 0.05
LIS output number of ensembles per tile: 12
LIS output nest index: 1
LIS output elevation data source: "LDT"
LIS output slope data source: "LDT"
LIS output aspect data source: "LDT"
LIS output soil texture data source:"LDT"
LIS output soil fraction data source: none
LIS output number of soil moisture layers: 4
LIS output number of soil temperature layers: 4
LIS output soil moisture layer thickness: 0.1 0.3 0.6 1.0
LIS output soil temperature layer thickness: 0.1 0.3 0.6 1.0