soil class extraction stops

[ShervanGharari]

@kasra-keshavarz faced some issues for FHIMP, the job failed for gistool soil class. I tried below to test the shapefile and the dataset. I am following:

./extract-gis.sh --dataset="soil_class" \
  --dataset-dir="/project/rpp-kshook/Model_Output/soil_classes/" \
  --variable="soil_classes" \
  --shape-file="/path/to/shapefile/domain_cat.shp" \
  --print-geotiff=true \
  --output-dir="/path/to/output" \
  --prefix="test_" \
  --stat="frac"

I get this:

Attaching package: ‘dplyr’

The following objects are masked from ‘package:stats’:

    filter, lag

The following objects are masked from ‘package:base’:

    intersect, setdiff, setequal, union

[1] "Transforming EPSG to 4326"
[1] "NOT including NA values"
^M  |                                                                            ^M  |          >
Error in `$<-.data.frame`(`*tmp*`, "base_value", value = NA_real_) :
  replacement has 1 row, data has 0
Calls: <Anonymous> ... tryCatchList -> .resultColumns -> $<- -> $<-.data.frame
Execution halted

the shapefile of the subbasins is domain.zip.

Thank you in advance for looking into this.

[ShervanGharari]

I believe this is resolved after the replacement of the global geotiff?

[kasra-keshavarz]

Yes, correct. Sorry I missed this.