...
Ensure that the initial fix is in your branch (here: commit with fix ).
If your simulation uses an initial condition file for land (finidat), replace the NaNs in TWS_MONTH_BEGIN with the fill value of 1.e+36. Below are methods to perform the conversion :(note that they overwrite the input file, so make a copy).
Using NCO functions:
Code Block language none ncatted -a _FillValue,TWS_MONTH_BEGIN,o,f,NaN ${infile.nc} ncatted -a _FillValue,TWS_MONTH_BEGIN,m,f,1.0e36 ${infile.nc}
Using a Python script:
Code Block language py from netCDF4 import Dataset import numpy as np ofile = Dataset('infile.nc','r+') var_array = f.variables['TWS_MONTH_BEGIN'] var_array[:][np.isnan(var_array[:])] = 1.e+36 ofile.close()
...