Assuring the spatial domain is fully covered (#50)

In this commit, a value of 0.1 degrees is added/subtracted to/from the
minimum and maximum values of latitude and longitude extents. This is to
assure that the spatial domain of the requested area is fully covered in
the subsetted dataset files. Previously, the extents were not accurately
covered.

The value of 0.1 degress reflect the resolution of RDRSv2.1.

Signed-off-by: Kasra Keshavarz <[email protected]>

Author: kasra-keshavarz

scripts/eccc-rdrs/rdrs.sh

@@ -206,14 +206,22 @@ maxLon=$(echo $lonLims | cut -d ',' -f 2)
 # unload and load necessary modules
 unload_core_modules
 load_ncl_module
+
 # choose a sample file as all files share the same grid
 domainFile="$(find "${datasetDir}/" -type f -name "*.nc" | head -n 1)"
+
 # parse the upper and lower bounds of a given spatial limit
 minLat=$(echo $latLims | cut -d ',' -f 1)
 maxLat=$(echo $latLims | cut -d ',' -f 2)
 minLon=$(echo $lonLims | cut -d ',' -f 1)
 maxLon=$(echo $lonLims | cut -d ',' -f 2)
 
+# adding/subtracting 0.1 degree to/from max/min values
+minLat=$(bc <<< "$minLat - 0.1")
+maxLat=$(bc <<< "$maxLat + 0.1")
+minLon=$(bc <<< "$minLon - 0.1")
+maxLon=$(bc <<< "$maxLon + 0.1")
+
 # extract the associated indices corresponding to $latLims and $lonLims
 coordIdx="$(ncl -nQ 'coord_file='\"$domainFile\" 'minlat='"$minLat" 'maxlat='"$maxLat" 'minlon='"$minLon" 'maxlon='"$maxLon" "$coordIdxScript")"
 

scripts/ouranos-espo-g6-r2/espo-g6-r2.sh

@@ -229,14 +229,22 @@ maxLon=$(echo $lonLims | cut -d ',' -f 2)
 # unload and load necessary modules
 unload_core_modules
 load_ncl_module
+
 # choose a sample file as all files share the same grid
 domainFile=$(find ${datasetDir} -type f -name "*.nc" | head -n 1)
+
 # parse the upper and lower bounds of a given spatial limit
 minLat=$(echo $latLims | cut -d ',' -f 1)
 maxLat=$(echo $latLims | cut -d ',' -f 2)
 minLon=$(echo $lonLims | cut -d ',' -f 1)
 maxLon=$(echo $lonLims | cut -d ',' -f 2)
 
+# adding/subtracting 0.1 degree to/from max/min values
+minLat=$(bc <<< "$minLat - 0.1")
+maxLat=$(bc <<< "$maxLat + 0.1")
+minLon=$(bc <<< "$minLon - 0.1")
+maxLon=$(bc <<< "$maxLon + 0.1")
+
 # extract the associated indices corresponding to $latLims and $lonLims
 coordIdx="$(ncl -nQ 'coord_file='\"$domainFile\" 'minlat='"$minLat" 'maxlat='"$maxLat" 'minlon='"$minLon" 'maxlon='"$maxLon" "$coordIdxScript")"