tide velocity problem

Discussion on computers, ROMS installation and compiling

Moderators: arango, robertson

Post Reply
Message
Author
Rocky_Wang
Posts: 45
Joined: Tue Mar 01, 2016 1:38 am
Location: the Ocean University of China

tide velocity problem

#1 Post by Rocky_Wang » Thu Jan 19, 2017 5:58 am

Hello, I run the roms model with only tide forcing. We have put an ADCP into the modeled region to record the tide current velocity. When I compare the computed current velocity with observed data,I find the simulated result is smaller. which method or parameter I should adjust to increse the simulated currtrn velocity? Attachments are my H file and some paremeters about bottom boundary layer.
Attachments
bhseatide.h
(1.76 KiB) Downloaded 113 times
捕获.PNG
捕获.PNG (6.21 KiB) Viewed 2187 times

lfparejar
Posts: 12
Joined: Sat Oct 05, 2013 5:21 am
Location: Rutgers University

Re: tide velocity problem

#2 Post by lfparejar » Wed Apr 05, 2017 1:02 am

One way could be to decrease the value of Zob in ocean.in. Zos is the surface roughness length, and you wouldn't need to tweak that here.

Another would be to decrease rdrg2, also in the ocean.in you show.

And another way could be to prescribe a constant, lower, bottom drag coefficient to increase the velocity. To do that (and if you want to keep using a quadratic drag law with the #define UV_QDRAG you have in the header file), you need to add these lines:

#define ANA_DRAG
#define UV_DRAG_GRID

Then you will need to modify and use ana_drag.h and add something like this:

#if defined YOUR_HEADER_NAME
# elif defined UV_QDRAG
DO j=JstrT,JendT ! based on Chezy coefficient (g/c^2)
DO i=IstrT,IendT
rdrag2(i,j)= *your new drag value*
END DO
END DO
# endif

I'd expect quite some testing if you want to tweak rdrag2 to match your ADCP data.

Rocky_Wang
Posts: 45
Joined: Tue Mar 01, 2016 1:38 am
Location: the Ocean University of China

Re: tide velocity problem

#3 Post by Rocky_Wang » Fri May 05, 2017 11:21 am

lfparejar wrote:One way could be to decrease the value of Zob in ocean.in. Zos is the surface roughness length, and you wouldn't need to tweak that here.

Another would be to decrease rdrg2, also in the ocean.in you show.

And another way could be to prescribe a constant, lower, bottom drag coefficient to increase the velocity. To do that (and if you want to keep using a quadratic drag law with the #define UV_QDRAG you have in the header file), you need to add these lines:

#define ANA_DRAG
#define UV_DRAG_GRID

Then you will need to modify and use ana_drag.h and add something like this:

#if defined YOUR_HEADER_NAME
# elif defined UV_QDRAG
DO j=JstrT,JendT ! based on Chezy coefficient (g/c^2)
DO i=IstrT,IendT
rdrag2(i,j)= *your new drag value*
END DO
END DO
# endif

I'd expect quite some testing if you want to tweak rdrag2 to match your ADCP data.


Thanks very much, it is really helpful!
Now I want to set a constant, spatially varying bottom roughness length Z0B (m) associated with grain roughness z0N, I have the median grain size distribution data of the region that I modeled.
I know that
Z0N = 2.5 * D50 / 30
So how to modify and use ana_drag.h to achieve the goal.


#if defined MY_HEADER_NAME
# elif defined UV_LOGDRAG
DO j=JstrT,JendT ! based on Chezy coefficient (g/c^2)
DO i=IstrT,IendT
rdrag2(i,j)= *your new drag value* !Sorry I do not know how to modify there, can you give me more detaild reminder?
END DO
END DO
# endif

Chang
Posts: 20
Joined: Sun Dec 04, 2016 1:03 am
Location: ocean university of china

Re: tide velocity problem

#4 Post by Chang » Sun May 21, 2017 2:57 pm

hello Rocky_Wang
I meet the same problem with you.so if you have figured it out,could you please show me how to do that?
Thanks in advance.
my qq is 2454518664.

Post Reply