Ocean Modeling Discussion


Search for:
It is currently Tue Sep 17, 2019 5:07 am

Post new topic Reply to topic  [ 3 posts ] 

All times are UTC

Author Message
PostPosted: Tue Jul 11, 2017 11:06 am 

Joined: Wed Jul 02, 2014 2:50 pm
Posts: 1
Location: First Institute of Oceanography
something wrong happened ,when i make 'LuvSponge == T ;LtracerSponge == T T'...

GET_GRID - error while reading variable: spherical
in grid NetCDF file: ./input/GRID.nc

I use the code "roms_kata_svn". When i close the 'LuvSponge' and 'LtracerSponge', it works well , is thers some bug in the code,or i use the wrong option?

Reply with quote  
PostPosted: Tue Jul 11, 2017 5:40 pm 
User avatar

Joined: Wed Jul 02, 2003 5:29 pm
Posts: 3665
Location: IMS/UAF, USA
It looks like my version of get_grid.F got a bit scrambled. Try again? I've never used the sponge options, so that code hasn't been tested.

Reply with quote  
PostPosted: Wed Jul 12, 2017 7:16 am 
User avatar

Joined: Mon Apr 28, 2003 5:44 pm
Posts: 502
Location: Rutgers University
Hernan Arango's add_sponge.m function in the myroms.org Matlab tools will add (or overwrite) the sponge to an existing grid file:

add_sponge:  Adds sponge data to a ROMS Grid NetCDF file
  add_sponge(ncfile, visc_factor, diff_factor)
  Adds enhanced viscosity and diffusion scaling variables (visc_factor
  and diff_factor) to an existing ROMS Grid NetCDF file.  These scales
  are used in an application to set sponge areas with larger horizontal
  mixing coefficients for the damping of high frequency noise coming
  from open boundary conditions or nesting.  In ROMS, these scales are
  used as follows:
     visc2_r(i,j) = visc2_r(i,j) * visc_factor(i,j)
     visc4_r(i,j) = visc4_r(i,j) * visc_factor(i,j)
     diff2(i,j,itrc) = diff2(i,j,itrc) * diff_factor(i,j)
     diff4(i,j,itrc) = diff4(i,j,itrc) * diff_factor(i,j)
  where the variables 'visc_factor' and 'diff_factor' are defined at
  RHO-points.  Usually, sponges are linearly tapered over several grid
  points adjacent to the open boundaries.  Its positive values linearly
  increase from the inner to outer edges of the sponge. At the interior
  of the grid we can values of zero (no mixing) or one (regular mixing).
  NOTICE that it is more advantageous to specify these scaling factors
  in the ROMS Grid that coding it inside ROMS.  We can plot and adjust
  their values in an easy way in Matlab.
  On Input:
     ncfile        GRID NetCDF file name (string)
     visc_factor   Viscosity factor (2D array; nondimensional, positive)
     diff_factor   Diffusivity factor (2D array; nondimensional, positive)

John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu

Reply with quote  
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 3 posts ] 

All times are UTC

Who is online

Users browsing this forum: No registered users and 3 guests

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