Ocean Modeling Discussion

ROMS/TOMS

Search for:
It is currently Mon Oct 23, 2017 12:34 am




Post new topic Reply to topic  [ 5 posts ] 

All times are UTC

Author Message
 Post subject: seamount
PostPosted: Thu Oct 09, 2014 8:02 am 
Offline

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
 Profile  
Reply with quote  
 Post subject: Re: seamount
PostPosted: Thu Oct 09, 2014 8:45 pm 
Offline
User avatar

Joined: Wed Jul 02, 2003 5:29 pm
Posts: 3236
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
 Profile  
Reply with quote  
 Post subject: Re: seamount
PostPosted: Fri Oct 10, 2014 1:55 pm 
Offline

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
 Profile  
Reply with quote  
 Post subject: Re: seamount
PostPosted: Fri Oct 10, 2014 9:34 pm 
Offline
User avatar

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


Top
 Profile  
Reply with quote  
 Post subject: Re: seamount
PostPosted: Sat Oct 11, 2014 3:00 am 
Offline

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


Top
 Profile  
Reply with quote  
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 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 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