Skip to content

Add the possibility to read a time-dependant reference pressure in a file

The contribution of the atmospheric pressure into the sea surface height is estimated through the inverse barometer relation :

ssh_ib = - 1 / (rho * g) * (Patm - Pr)

The computation depends on the reference pressure Pr. Currently, there are 2 ways to estimate Pr :

  • considered it as a constant
  • compute it directly as the average of the atmospheric pressure forcing file over ocean

The second is certainly the better option but makes sense only in global modelling.

Here we propose a third option for regional modelling which consists in reading directly an atmospheric pressure average calculated upstream and stored in a 1D file.

In order to make NEMO able to read a such file, a preliminary change in fldread.F90 has already been done as explained here #552 (closed).

This option changes slightly the namelist :

ln_ref_apr = .false. ! ref. pressure: global mean Patm (T) or a constant (F)

is replaced by

nn_ref_apr = 0 ! ref. pressure: 0: constant, 1: global mean or 2: read in a file

Edited by Guillaume Samson
To upload designs, you'll need to enable LFS and have an admin enable hashed storage. More information