Hi Matt,
Thanks very much for your warm welcome. I finally solved my issue - "num_depths" was originally defined as "1" for the initial condition setting. I was thinking of modeling the inital water depth from ~0m. It looks like that will not be simulated by GLM. That is all fixed now.
I am cross-checking simulated water levels by GLM and another program used by the project. It looks like GLM results are about 10m lower when the system water level approaches stable. I am looking at the lake file ("csv_lake_fname"), which displays a list of parameters with no units. Where can I find more information (e.g. data source and unit) if I need details on those parameters? One parameter I am paying extra attention to at this stage is "Evaporation". There is a long list of parameters from ECMWF website (https://apps.ecmwf.int/codes/grib/param-db/). Would you mind giving me some advices on which one to look at?
Thank you very much.
Kind regards,
Jade