Ocean Modeling Discussion

ROMS/TOMS

Search for:
It is currently Tue Sep 17, 2019 4:58 am




Post new topic Reply to topic  [ 9 posts ] 

All times are UTC

Author Message
PostPosted: Mon Jan 21, 2019 6:26 pm 
Offline

Joined: Thu Feb 27, 2014 7:51 pm
Posts: 5
Location: Instituto Mexicano de Tecnología del Agua
Dear ROMS users,

I had an issue with all the variables, like temperature, that has in all levels in the vertical a granular aspect (attached figures). I still can figured out how to fix this.

I run ROMS 3.7 for 15 years simulation. Am using HYCOM initial conditions and climatology for each month. Forcing from COADS (heat, wind stress, sss, sst and fresh water). My open boundaries Che for free surface, Shc for 2D, RadNud for 3D and tracers, Rad for TKE.

Thanks

Jorge Ibarra


Attachments:
File comment: CPP options
pacific_tides.h [1.51 KiB]
Downloaded 69 times
File comment: Bottom temperature
pacificMX_temperature_bottom.png
pacificMX_temperature_bottom.png [ 369.31 KiB | Viewed 1234 times ]
File comment: Surface temperature
pacificMX_temperature_surface.png
pacificMX_temperature_surface.png [ 329.55 KiB | Viewed 1234 times ]
Top
 Profile  
Reply with quote  
PostPosted: Wed Jan 30, 2019 9:14 pm 
Offline

Joined: Wed Dec 31, 2003 6:16 pm
Posts: 795
Location: USGS, USA
does it start out like that or evolve to that? check the fields in your grid file like dmde dndx. check angle.
are you using N2S2_HORAVG?

-j


Top
 Profile  
Reply with quote  
PostPosted: Wed Jan 30, 2019 10:22 pm 
Offline

Joined: Thu Feb 27, 2014 7:51 pm
Posts: 5
Location: Instituto Mexicano de Tecnología del Agua
Hello,

Thank you for the reply. The first time steps looks smooth, and it evolve to the figure i attached.
I not used N2S2_HORAVG. Do you recomend to use it?

I'm doing some test and founded that using advection schemes for tracers TS_U3HADVECTION with TS_C4VADVECTION and for momentum UV_U3HADVECTION and UV_C4VADVECTION, i got better results.

Cheers

Jorge Ibarra


Top
 Profile  
Reply with quote  
PostPosted: Thu Jan 31, 2019 1:10 pm 
Offline

Joined: Wed Dec 31, 2003 6:16 pm
Posts: 795
Location: USGS, USA
yes and try
#define N2S2_HORAVG


Top
 Profile  
Reply with quote  
PostPosted: Thu Jan 31, 2019 6:40 pm 
Offline

Joined: Thu Feb 27, 2014 7:51 pm
Posts: 5
Location: Instituto Mexicano de Tecnología del Agua
Hello,

With respect to the variables dmde, dndx and the angle, i attached the figures.

I will activate N2S2_HORAVG and will report it as soon as possible.

Cheers

Jorge Ibarra


Attachments:
File comment: DNDX
pacificMX_DNDX.png
pacificMX_DNDX.png [ 96.03 KiB | Viewed 1097 times ]
File comment: DMDE
pacificMX_DMDE.png
pacificMX_DMDE.png [ 104.11 KiB | Viewed 1097 times ]
File comment: Angle
pacificMX_angle.png
pacificMX_angle.png [ 136.02 KiB | Viewed 1097 times ]
Top
 Profile  
Reply with quote  
PostPosted: Thu Jan 31, 2019 9:09 pm 
Offline
Site Admin
User avatar

Joined: Wed Feb 26, 2003 4:41 pm
Posts: 1081
Location: IMCS, Rutgers University
I guess that the more fundamental question is what type of horizontal advection are you using? Are you specifying explicit diffusion and viscosity? What do you specify for vertical diffusion/viscosity parameterization closure? You should compute your Reynolds Number (U dx/nu), where U is the maximum velocity magnitude for your application. dx is grid spacing, and nu is diffusion/viscosity coefficient.


Top
 Profile  
Reply with quote  
PostPosted: Fri Feb 01, 2019 6:00 pm 
Offline

Joined: Thu Feb 27, 2014 7:51 pm
Posts: 5
Location: Instituto Mexicano de Tecnología del Agua
Hello,

I used 4th order centered advection for tracers and momentum. For explicit viscosity TS_DIF4=2 and UV_VIS2=15. For the vertical parametrization closure: Mellor & Yamada, with Kantha and Clayson stability function, also to mention I added RI_SPLINES.

About the Reynolds Number from the log file, where U = 1.9 m/s, dx = 7000 m and nu = 15 m2/s; Re = 886.

Cheers.

Jorge Ibarra


Top
 Profile  
Reply with quote  
PostPosted: Fri Feb 01, 2019 6:51 pm 
Offline
Site Admin
User avatar

Joined: Wed Feb 26, 2003 4:41 pm
Posts: 1081
Location: IMCS, Rutgers University
:shock: Well, that's your problem. You are using biharmonic diffusivity (m4/s), which scales as 1/12 * ABS(U) * (dx ** 3). A value for TS_DIF4 of 2 m4/s is too small :!: You need values of the order of 2.0d+8 m4/s. That is, you are using a value that is eight orders of magnitude too small :!: No wonder why your solution is so granular. I hope that you are not using MIX_S_TS, you will need, at least, mixing along geopotentials (MIX_GEO_TS) to avoid spurious diapycnal diffusion.


Top
 Profile  
Reply with quote  
PostPosted: Fri Feb 01, 2019 11:44 pm 
Offline

Joined: Thu Feb 27, 2014 7:51 pm
Posts: 5
Location: Instituto Mexicano de Tecnología del Agua
Hello,

Thank you for the reply. Will do my new tests with that values.

Yes, I'm using MIX_GEO_TS.

Cheers

Jorge Ibarra


Top
 Profile  
Reply with quote  
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 9 posts ] 

All times are UTC


Who is online

Users browsing this forum: No registered users and 0 guests


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
Powered by phpBB® Forum Software © phpBB Group