Hi everyone, I am conducting research on Streamflow Nudging Data Assimilation. After completing all the preparatory work and running the program, there were no output files related to Nudging, and no error messages were reported either. I suspect that the assimilation was not successful. The terminal output information after running the program was "Ndg: Start initializing reach gages (netcdf)
Ndg: Starting read_reach_gage_collocation
Ndg: Finish read_reach_gage_collocation
Ndg: nGagesDomain: 3
Ndg: nLinksL 2406
Ndg: size(nodeGageStr%nodeId): 3
Ndg: nodeGageStr%usgsId((size(nodeGageStr%nodeId)-nGagesDomain+1):(size(nodeGageStr%nodeId))): 13 12 11
Ndg: Finish initializing reach gages (netcdf)
getting dimension from file: './domain/nudgingParams.nc'
Ndg: nParamGages: 3
Ndg: nParamMonth 0
Ndg: nParamThresh 0
Ndg: nParamThreshCat 0
Ndg: temporalPersistence: F
Ndg: start read_nudging_param_file
Ndg: paramFile: './domain/nudgingParams.nc'
Ndg: finish read_nudging_param_file
Ndg: nudgingParamsStr%usgsId 11
Ndg: nudgingParamsStr%R 0.25
Ndg: nudgingParamsStr%G 1.00
Ndg: nudgingParamsStr%tau 15.0 Ndg: obsResolution:15
Ndg: maxTau: 15.000000 Ndg: nObsTimes: 4
Ndg: start obs_static_to_struct
Ndg: nudgeSpatial = F
Finnish obs_static_to_struct
Ndg: Finish init_stream_nudging
Finnish obs_static_to_struct
Finnish obs_static_to_struct
Finnish obs_static_to_struct
output hydrology at time : 2020-06-27_08:00:00 0
before call output_chrt
before call output_chrt
before call output_chrt
before call output_chrt
end calling output functions
end calling output functions
end calling output functions"
I noticed that nParamMonth, nParamThresh, and nParamThreshCat seem to have issues; they are all set to 0. I will attach nudgingParams.nc. The final diag file only outputs "Ndg: start finish_stream_nudgingNdg: end of finish_stream_nudging
The model finished successfully......." at the end. I would like to know where the problem lies and how I can correct it. I have attached the namelist.hrldas and hydro.namelist files. I look forward to your reply. Thank you very much.