Fortran runtime error: Bad integer while reading precipitation data

[klekarkar]

Hi mHM community,
I am setting up an mHM model for my simulation domain. The model reads the morphological inputs without any problems, but it runs into a runtime error while reading the precipitation file. My meteorological data consists of grids that lie within irregular boundaries of a country, so I have allocated all the cells outside of this boundary -9999 values to match the extents of the morphological data (which contains values that extend outside of the country boundaries). The meteorological product is not available outside these boundaries. I am not sure if the problem comes from this mismatch in the data.
My model files and mhm.nml are here:
https://drive.google.com/drive/folders/1aVNG6u_Qoy_W9RVcfMPnxA8ua4BsD__V?usp=sharing

What are your suggestions to fix this problem?

The images are attached below.

mHM error

precipitation map and DEM in the background

As always, thanks for your feedback!

[pallavhydro]

Hi @klekarkar

Thanks for sharing the issue.

I have two suggestions:

1. Just in case you have some pixel of pre without data inside the DEM, I would suggest you to enter 0 instead of -9999 while preparing your pre.nc. That will ensure the missing pre is not an issue.

2. I also have a guess that the error message originates from the structure of your pre.nc. Could you please share the output of "ncdump -h pre.nc" here, so that we can assess the structure further.

3. Also please directly upload the mhm.nml file here. It's much easier to access and view in that way. Thanks.

Pallav

[klekarkar]

Thanks a lot @pallavhydro for the suggestions. Here is the output of ncdump -h pre.nc and my mhm.nml file

mhm_nml.zip.

In the meantime, I will follow your first suggestion and I will post the result here.

Thanks for your feedback

[klekarkar]

I changed the values of pre.nc and pet.nc outside the DEM from -9999 to 0 but the error still remains.

[pallavhydro]

Hi @klekarkar

Thanks for following my suggestions.
I checked your setup, mainly your meteo files and have the following suggestions:

1. The error that you are getting in this case is because mHM is getting confused while reading the reference time from your netcdf file. I am referring to the time:units attribute. In your case it is : days since 1994-01-01T23:00:00.000000000. Please try to convert this to days since 1994-01-01 23:00:00. E.g., see the screenshot (test_domain data)

2. Please ensure the _FillValue and missing_value attributes for tavg and pet follow what you have done for pre.nc. Due to the missing attributes/ correct attributes, your pet and temp data currently have -9999 in the value pool. (Note: You can cdo infon tavg.nc | less to confirm this).

Let me know if you could progress with these hints.

Regards
Pallav

p.s. Nice model setup by the way.

[klekarkar]

Hi @pallavhydro ,
Thanks a lot for your feedback. The error was indeed coming from the time attribute of my meteo files. (I also just realized the error message contained a clue to the source of the error since the time attribute is the only attribute with integer format in the files.)
I fixed them and unified the other attributes as you have suggested, and it worked like magic. The model now runs!

Thanks for your helpful suggestions

[pallavhydro]

Hi @klekarkar

Congratulations on your model run.

I too recall the first time I was able to run mHM for my setup. I confirm, it was magical ;)

Cheers
Pallav