Where to introduce obs error

Discussion about tangent linear and adjoint models, variational data assimilation, and other related issues.

Moderators: arango, robertson

Post Reply
Message
Author
xavierC
Posts: 24
Joined: Thu Jul 16, 2009 5:17 pm
Location: Madeira University

Where to introduce obs error

#1 Unread post by 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

User avatar
arango
Site Admin
Posts: 1347
Joined: Wed Feb 26, 2003 4:41 pm
Location: DMCS, Rutgers University
Contact:

Re: Where to introduce obs error

#2 Unread post by 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
Posts: 24
Joined: Thu Jul 16, 2009 5:17 pm
Location: Madeira University

Re: Where to introduce obs error

#3 Unread post by 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)

Post Reply