[MITgcm-support] generate a stable current field

[guopu at scsio.ac.cn]

Hello Sir,

I feel the OBCS package is so complicate, so I start from generating a simple current field to understand this package. Now I have 3 questions:

1. I build a 2D region with uniform depth 1000m, and set OB_Iwest=1, OB_Ieast=-1. Now I am wondering how to generate a stable barotropic current in this region. The only effective way I found is to use files containing uniform velocity for parameters: "uVelInitFile", "OBWuFile", and "OBEuFile". Is it right? or is there another "standard procedure" ?

2. I want to generate above current field by only one prescribe file for "OBWuFile", so I make this (1*Nr*50)-dimension file with gradually increased velocity, from 0 (or 0.01) to 0.1m/s in 50 "time levels". However, the model results show that U at west boundary keep to be zero (or 0.01m/s). What could be the reason for this? and what is the meaning of "time levels"?

3. I add a Gaussian shape obstacle on the bottom and hope to see what will happen in this region. Should I turn on one or some of these flags: useOBCSsponge, useOrlanskiEast/West, useOBCS_balance? I thought the "OB_Iwest=1, OB_Ieast=-1" would be enough to let water/waves flow/propagate out freely.

Studying this package has cost me several weeks, I am hungry for help, thanks!

Pu






-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20150706/9ad34478/attachment.htm>