# Ocean Modeling Discussion

ROMS/TOMS

Search for:
 It is currently Wed Jun 26, 2019 8:06 am

 Page 1 of 1 [ 5 posts ]
 All times are UTC
Author Message
 Post subject: seamountPosted: Thu Oct 09, 2014 8:02 am

Joined: Tue Jun 18, 2013 1:10 pm
Posts: 11
Location: Tsinghua university
I run the SEAMOUNT case, but the seamount shape could not be found when the flow field (temperature, u, v) was plot. can someone know how to obtain the seamount shape in the seamount.h or ocean_seamount.in.
thankyou very much!

Top

 Post subject: Re: seamountPosted: Thu Oct 09, 2014 8:45 pm

Joined: Wed Jul 02, 2003 5:29 pm
Posts: 3614
Location: IMS/UAF, USA
How are you plotting it? Your plotting tools need to know about the z coordinate stretching.

The depth is set in ana_grid.h:
Code:
#elif defined SEAMOUNT
DO j=JstrT,JendT
DO i=IstrT,IendT
val1=(xr(i,j)-0.5_r8*Xsize)/40000.0_r8
val2=(yr(i,j)-0.5_r8*Esize)/40000.0_r8
h(i,j)=depth-4500.0_r8*EXP(-(val1*val1+val2*val2))
END DO
END DO

Top

 Post subject: Re: seamountPosted: Fri Oct 10, 2014 1:55 pm

Joined: Tue Jun 18, 2013 1:10 pm
Posts: 11
Location: Tsinghua university
Thankyou for your help! i want to create a flow from west to east with a constant eastward flow velocity of 0.5m/s for this seamount case. so i set a Cla boundary condition on western side and Rad boundary condition on eastern side (LBC(isUvel)== Cla Clo Rad Clo !3D u-momentum). The flow velocity was given through ANA_M3OBC (the velocity was 0.5). the details are as follows,

#if defined SEAMOUNT
IF (LBC(iwest,isUvel,ng)%acquire.and. &
& LBC(iwest,isVvel,ng)%acquire.and. &
& DOMAIN(ng)%Western_Edge(tile)) THEN
DO k=1,N(ng)
DO j=JstrT,JendT
BOUNDARY(ng)%u_west(j,k)=0.5_r8
END DO
DO j=JstrP,JendT
BOUNDARY(ng)%v_west(j,k)=0.0_r8
END DO
END DO
END IF
i want to know what i have set is right or wrong to realize the constant eastward flow? the results shows the velocity is very small and the flow is also irregular, i hope someone can help me make the right set. thankyou very much!

Top

 Post subject: Re: seamountPosted: Fri Oct 10, 2014 9:34 pm

Joined: Wed Jul 02, 2003 5:29 pm
Posts: 3614
Location: IMS/UAF, USA
What are your ubar boundary conditions? You need to set up a consistent ubar for ANA_M2OBC.

Top

 Post subject: Re: seamountPosted: Sat Oct 11, 2014 3:00 am

Joined: Tue Jun 18, 2013 1:10 pm
Posts: 11
Location: Tsinghua university
kate, thankyou i got the expected flow field.

Top

 Display posts from previous: All posts1 day7 days2 weeks1 month3 months6 months1 year Sort by AuthorPost timeSubject AscendingDescending
 Page 1 of 1 [ 5 posts ]

 All times are UTC

#### Who is online

Users browsing this forum: No registered users and 1 guest

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

Search for:
 Jump to:  Select a forum ------------------ News, Events & Job Opportunities    Meetings/Workshops    Job Opportunities    Ocean News ROMS/TOMS    ROMS Adjoint    ROMS Benchmarks    ROMS Bugs    ROMS Discussion    ROMS Documentation    ROMS Ecosystem    ROMS FAQ    ROMS Ice    ROMS Information    ROMS Installation    ROMS Known Problems    ROMS Messages    ROMS Problems    ROMS Releases    ROMS Results    ROMS Sediment    ROMS Source    ROMS Tips    ROMS Tools and Techniques    ROMS Trivia    ROMS Usage    ROMS Webinar    ROMS Wish List ROMS/TOMS Applications    User Applications    Adriatic Sea Ocean Modeling    Ocean Modeling FAQ