Where to introduce obs error

[xavierC]

Hello,
Playing with some grep i realized that "obs_error" variable is never used in the code. Therefore i wonder where i should introduce observation error when known ?
Thanks

[arango]

Of course, you need the observation error covariance (diagonal matrix) in the ROMS 4D-Var algorithm! The NetCDF variable is obs_error, which is read into ROMS variable ObsErr. You are searching for the wrong variable. Check ROMS/Modules/mod_fourdvar.F for the description of all variables.

[xavierC]

thanks a lot for your fast answer
guess i got it on my way too:
CALL netcdf_get_fvar (ng, model, OBS(ng)%name, Vname(1,idOerr)