# Ocean Modeling Discussion

ROMS/TOMS

Search for:
 It is currently Sat Jan 19, 2019 11:03 pm

 Page 1 of 1 [ 3 posts ]
 All times are UTC
Author Message
 Post subject: Keep net volume as a constant without VolconsPosted: Thu Dec 20, 2018 4:36 am

Joined: Mon Sep 24, 2018 7:40 pm
Posts: 11
Location: CSIRO - Hobart Site
Dear All,

My purpose:
Keep the net volume as a constant without turning on the VolCons

My design:
Code:
LBC(isFsur) ==   Clo     Cha     Cha     Cha            ! free-surface
LBC(isUbar) ==   Clo     Cla     Cla     Cla            ! 2D U-momentum
LBC(isVbar) ==   Clo     Cla     Cla     Cla            ! 2D V-momentum
LBC(isMtke) ==   Clo     Gra     Gra     Gra            ! mixing TKE

LBC(isTvar) ==   Clo     Cla     Cla     Cla \          ! temperature
Clo     Cla     Cla     Cla            ! salinity

I Clamped the 2D momentum, which is balanced at every input time step (the frequency is 5 days) in the BOUNDARY.NC so that I can control the net flux from the boundaries as zero all the time, and gave the freedom for 3D momentum to adjust.

My problem:
The net volume increases at the speed of 1cm per year approximately, and I think the key point is that I still did not provide 'balanced flux' in the boundary file because of i,j indices offset.

I calculated the boundary flux in this way (the ubar/vbar points used are indicated with black circle, the h/pn/pm/mask points used are indicated with black square in the attached figure):

EAST_Flux = ubar_east(1:M-1) * (h(L-1,1:M-1)+h(L,1:M-1)) * 0.5 * (2/(pn(L-1,1:M-1)+pn(L,1:M-1))) * mask_u(L,1:M-1)
NORTH_FLUX = vbar_north(1:L-1) * (h(1:L-1,M-1)+h(1:L-1,M)) * 0.5 * (2/(pm(1:L-1,M-1)+pm(1:L-1,M))) * mask_v(1:L-1,M)
SOUTH_FLUX = vbar_south(1:L-1) * (h(1:L-1,0)+h(1:L-1,1)) * 0.5 * (2/(pm(1:L-1,0)+pm(1:L-1,1))) * mask_v(1:L-1,1)

Did my design make sense for my purpose?
Why I cannot get a constant net volume?

Thank you very much!

 Attachments: File comment: c_grid staggered_grid_rho_cells_indicated.png [ 129.47 KiB | Viewed 172 times ]
Top

 Post subject: Re: Keep net volume as a constant without VolconsPosted: Sun Jan 06, 2019 4:39 am

Joined: Mon Sep 24, 2018 7:40 pm
Posts: 11
Location: CSIRO - Hobart Site
Could anyone help me please, I am still confused.

Thank you very much!

Top

 Post subject: Re: Keep net volume as a constant without VolconsPosted: Mon Jan 14, 2019 3:26 am

Joined: Mon Sep 24, 2018 7:40 pm
Posts: 11
Location: CSIRO - Hobart Site
The net volume is changed when I just used different SST/SSS for QCORRECTION and SCORRECTION.

So temperature and salinity will influence the net volume? How the stflx makes difference to zeta? ROMS is a Boussinesq model.

Top

 Display posts from previous: All posts1 day7 days2 weeks1 month3 months6 months1 year Sort by AuthorPost timeSubject AscendingDescending
 Page 1 of 1 [ 3 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