Different output grids when interpolating to r360x180 and subsetting region
Hello, I found a potential bug when interpolating data: I read two different datasets and interpolate them to 'r360x180' grid, and also request to subset the region. Using the same parameters in both datasets I get two different output grids.
tmp1 <- Load("g500", exp = list(CMCC=list(name="CMCC",path="/esarchive/exp/cmcc/system3_m1-c3s/monthly_mean/g500_f12h/$VAR_NAME$_$START_DATE$.nc")), obs = NULL, sdates= "19931101", latmin=20, latmax=80, lonmin=270, lonmax=60, output="lonlat", grid="r360x180", nprocs=4)
...
* - Experimental data: ( 1 x 40 x 1 x 6 x 60 x 150 ) x 8 bytes = 17280000 bytes.
tmp2 <- Load("g500", exp = list(ECMWF=list(name="ECMWF",path="/esarchive/exp/ecmwf/system5c3s/monthly_mean/g500_f12h/$VAR_NAME$_$START_DATE$.nc")), obs = NULL, sdates= "19931101", latmin=20, latmax=80, lonmin=270, lonmax=60, output="lonlat", grid="r360x180", nprocs=4)
...
* - Experimental data: ( 1 x 25 x 1 x 6 x 60 x 151 ) x 8 bytes = 10872000 bytes.
I have a workaround, which is loading both datasets at the same time, then it returns the 60x151 grid for both datasets, no matter the dataset order.
tmp3 <- Load("g500", exp = list(CMCC=list(name="CMCC",path="/esarchive/exp/cmcc/system3_m1-c3s/monthly_mean/g500_f12h/$VAR_NAME$_$START_DATE$.nc"),ECMWF=list(name="ECMWF",path="/esarchive/exp/ecmwf/system5c3s/monthly_mean/g500_f12h/$VAR_NAME$_$START_DATE$.nc")), obs = NULL, sdates= "19931101", latmin=20, latmax=80, lonmin=270, lonmax=60, output="lonlat", grid="r360x180", nprocs=4)
tmp4 <- Load("g500", exp = list(ECMWF=list(name="ECMWF",path="/esarchive/exp/ecmwf/system5c3s/monthly_mean/g500_f12h/$VAR_NAME$_$START_DATE$.nc"),CMCC=list(name="CMCC",path="/esarchive/exp/cmcc/system3_m1-c3s/monthly_mean/g500_f12h/$VAR_NAME$_$START_DATE$.nc")), obs = NULL, sdates= "19931101", latmin=20, latmax=80, lonmin=270, lonmax=60, output="lonlat", grid="r360x180", nprocs=4)