grids for defining OCEAN(ng)%ubar

Facts, news, and guidance about ROMS software

Moderators: arango, robertson

Post Reply
Message
Author
chenchai
Posts: 3
Joined: Wed Aug 03, 2016 6:29 pm
Location: SAIC/NRLMRY

grids for defining OCEAN(ng)%ubar

#1 Post by chenchai » Tue Aug 09, 2016 1:58 pm

Hi,

In esmf_roms.F, we have the following.

CASE ('Ubar') ! 2D u-momentum
found=.TRUE.
DO j=JstrR,JendR
DO i=Istr,IendR
romsExpData(ifield,ng)%field(i,j)= &
& OCEAN(ng)%ubar(i,j,Kexp(ng))*scale
Fmin=MIN(Fmin,romsExpData(ifield,ng)%field(i,j))
Fmax=MAX(Fmax,romsExpData(ifield,ng)%field(i,j))
END DO
END DO

Does this mean that grids for OCEAN(ng)%ubar are at the 'u' points, since the 'i' loop starts from Istr? I would think it should be at the 'rho' points, since all romsExpData are at the 'rho' points.

Thanks,

James

noah.eckman
Posts: 7
Joined: Tue Jun 28, 2016 7:00 pm

Re: grids for defining OCEAN(ng)%ubar

#2 Post by noah.eckman » Mon Aug 15, 2016 2:13 pm

Yes, indeed, ubar is defined at u-points, vbar at v-points, just as the 3d momentum are at u and v points. This is verified by the script c_initial.m, which defines the dimensions for each of those.

Post Reply