Add one check to ensure the grid point has all NAs or no NA along time dim. Refine the documentation.

Add '-L' option in cdo command. It prevents potential segmentation fault in the underlying hdf5 library.