I am currently conducting research on predicting storm surge with NEMO-surge. A sensitivity test was conducted by changing the reference pressure in the NAMELIST, but there was no change in the sea surface height (SSH). It was confirmed that the model calculated the equivalent inverse barometer sea surface height(ssh_ib), but there was no change in the final SSH result. I cannot find the part where the reference pressure directly affects the sea surface height (not ocean momentum equation) in the source codes except sbcssm.F90. Is there any other process for calculating SSH?
Thanks in advance,
Further precisions are needed:
I guess you have in your namelist :
ln_apr_dyn = .true. ! Patm gradient added in ocean & ice Eqs., right ?
Is your model global (or do you have open boundaries making use of the Flather scheme) ?
When you say “no change in the final SSH result”, are you referring to the total SSH (hence including pressure effects if any) or corrected by the inverse barometer estimate ?
You should use
ln_apr_obc=.true. in order to add the inverse barometer component to your elevation data used by the Flather scheme.
Hence, It’s normal that you don’t see anything by changing the reference pressure with
ln_apr_obc=.false.. (since this is the gradient of the inverse barometer that is actually used in momentum equations).
I got it. I will check this part.
Thanks for your help.