ana_grid.h from upwelling test case

Discussion of how to use ROMS on different regional and basin scale applications.

Moderators: arango, robertson

Post Reply
Message
Author
angelolemos
Posts: 34
Joined: Wed Mar 11, 2009 4:24 pm
Location: UFES

ana_grid.h from upwelling test case

#1 Unread post by angelolemos »

People,

Anyone can tell me what means the values 84.5, 66.526 (the complete equation too)and strings not declared JstrR, JendR, IstrR and IendR below from ana_grid.h?

#elif defined UPWELLING
DO j=JstrR,JendR
IF (j.le.Mm(ng)/2) THEN
val1=REAL(j,r8)
ELSE
val1=REAL(Mm(ng)+1-j,r8)
END IF
val2=MIN(depth,84.5_r8+66.526_r8*TANH((val1-10.0_r8)/7.0_r8))
DO i=IstrR,IendR
h(i,j)=val2
END DO
END DO


Thanks everybody.

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: ana_grid.h from upwelling test case

#2 Unread post by kate »

IstrR and friends are shown in a diagram here.

The shape of the bathymetry in UPWELLING is a flat interior, with hyperbolic tangent ramps up on either side. The parameters are chosen to make the tanh (almost) meet the flat interior, to have a certain width and minimum depth, etc. Note that the range of tanh is -1 to 1.

Post Reply