opendrift icon indicating copy to clipboard operation
opendrift copied to clipboard

How to implement "pressure gradient force" in Openberg module?

Open dayexingyang opened this issue 10 months ago • 4 comments

Hi all,

I have noticed that in the Openberg module, the pressure gradient force calculation needs the sea surface slope data. However, I didn't find any public forcing data related to the sea surface slope. Do you have any recommendation databases for this data? Or is the sea surface slope calculated manually from the sea surface height? If it is obtained by taking the derivative of the sea surface height with respect to position, how is this achieved in Openberg?

Looking forward to your help.

Best wishes, Daye

dayexingyang avatar Jan 16 '25 10:01 dayexingyang

Hi, perhaps you can answer this question @AchrefAO ?

knutfrode avatar Jan 16 '25 10:01 knutfrode

Hi Daye, Yes, indeed, the sea surface slope can be calculated manually from the sea surface height and read into OpenBerg from an external file. However, its contribution is very small compared to other forces and even smaller than the Coriolis force. Nevertheless, you can calculate and include it if needed. Regards, Achref

AchrefAO avatar Jan 16 '25 10:01 AchrefAO

Hi Daye, Yes, indeed, the sea surface slope can be calculated manually from the sea surface height and read into OpenBerg from an external file. However, its contribution is very small compared to other forces and even smaller than the Coriolis force. Nevertheless, you can calculate and include it if needed. Regards, Achref

Hi Achref,

Thanks a lot for your quick reply!

I have a follow-up question. Does Openberg provide any functions to calculate the variable's derivative with respect to position? This would help me to calculate the sea surface slope manually from the sea surface height.

Thanks again.

Best wishes, Daye

dayexingyang avatar Jan 16 '25 10:01 dayexingyang

No unfortunately but it can be done manually using a simple formula derived from the concept of a gradient. Achref

AchrefAO avatar Jan 16 '25 10:01 AchrefAO