Hi Robert,
I check the model output. All model variables are there. If you look at the the error that it gives, it specifically has a problem in creating the “timestamp” variable.
However, I was suspecting if this has anything to do with skipping the first time stamp, and sure enough, when I set SKIPIND to “1 1” instead of to “0 0”, then the script works. What is this index for? As far as I understand, first 1 skips the first model timestamp when it reads all timestamps available in the first file and the second 1 is for the other model output files. If I use WRF model output, then I might want to skip certain timestamps at the beginning of a file since I prefer to dump model results to a single output file for WRF results. But, for MPAS, I don’t write all model results into a single output file since the output file size gets very large. Therefore, I want to read all hours that I generate and I don’t want to skip any timestamp and set it to “0 0”. To reiterate the issue that I am having, “matching_surface.csh” that uses MET_model.read.R works fine with SKINPIND set to “0 0”, but not “matching_raob.R” that reads exactly the same MET_model.read.R when SKIPIND set to “0 0”. I checked how MET_matching_surface.R and MET_matching_raob.R uses SKIPIND but couldn’t see any difference and understand how one works but not the other.
In short, the matching_raob.csh script works if I set to SKIPIND “1 1”, but not “0 0”. I would appreciate if you could kindly tell me if the way that I understand and explained above is the right way that that variable works. Thank you Robert.
Kemal.