Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

GOTM plume model: (1) change two layer temperature below -2.0, what could happen? (2) Is the NN randomly chosen? #53

Open
YingpuRin opened this issue Oct 10, 2024 · 0 comments

Comments

@YingpuRin
Copy link

Hi Community,

I am doing a GOTM-Observation comparison. And my observation is from cold region, somewhere below -2.0 C. I am now using the GOTM plume model published in 2022. I noticed that when using the "two-layer" setting in temperature, you need to set the temperature between -2.0 and 40 C. May I know what would happen if I change/remove this threshold? Or what would happen if I modify the lower threshold from -2.0C into -3.0C?

Also, may I know if the "NN" for temperature and salinity is randomly chosen, or if there are references for these numbers? If there are references, could you please share with me?

My observation has temperatures below -2.0C. Also, now the stratification in the temperature and salinity setting ("NN") is not doing a good job, because the plume depth is too large when compared with observations. Therefore, I would like to modify the "NN" in temperature and salinity, but I do not know if there is a specific reason for you to pick NN: 1.00000000E-04 and NN: 2.56000000E-04 for the plume model. If there are any references related to these numbers, please let me know. Otherwise, I will consider modify these numbers.

Looking forward to your reply,

Best Regards,
Yingpu

temperature: # temperature profile used for initialization and optionally relaxation [Celsius]
method: two_layer # method [off, file=from file, constant, two_layer=two layers with linear gradient in between, buoyancy=from salinity and buoyancy frequency; default=off]
constant_value: 20.0 # value to use throughout the simulation [Celsius; default=0.0]
file: T_HWD2_input.dat # path to file with series of profiles [default=]
column: 1 # index of column to read from [default=1]
two_layer:
z_s: 5.0 # depth where upper layer ends [m; min=0.0; default=0.0]
t_s: -1.95 # upper layer temperature [Celsius; min=-2.0; max=40.0; default=0.0]
z_b: 5.0 # depth where lower layer begins [m; min=0.0; default=0.0]
t_b: -1.9 # lower layer temperature [Celsius; min=-2.0; max=40.0; default=0.0]
NN: 1.00000000E-04 # square of buoyancy frequency [s^-2; min=0.0; default=0.0]
relax: # relax model temperature to observed/prescribed value
tau: 1.00000000E+15 # time scale for interior layer [s; min=0.0; default=1.00000000E+15]
salinity: # salinity profile used for initialization and optionally relaxation [psu]
method: two_layer # method [off, file=from file, constant, two_layer=two layers with linear gradient in between, buoyancy=from temperature and buoyancy frequency; default=off]
constant_value: 20.0 # value to use throughout the simulation [psu; min=0.0; default=0.0]
file: S_HWD2_input.dat # path to file with series of profiles [default=]
column: 1 # index of column to read from [default=1]
two_layer:
z_s: 5.0 # depth where upper layer ends [m; min=0.0; default=0.0]
s_s: 32.0 # upper layer salinity [g/kg; min=0.0; max=40.0; default=0.0]
z_b: 5.0 # depth where lower layer begins [m; min=0.0; default=0.0]
s_b: 34.5 # lower layer salinity [g/kg; min=0.0; max=40.0; default=0.0]
NN: 2.56000000E-04 # square of buoyancy frequency [s^-2; min=0.0; default=0.0]
relax: # relax model salinity to observed/prescribed value
tau: 1.00000000E+15 # time scale for interior layer [s; min=0.0; default=1.00000000E+15]

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant