New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
2011WP/2011Stream2/OpenBoundaries (diff) – NEMO

Changes between Version 15 and Version 16 of 2011WP/2011Stream2/OpenBoundaries


Ignore:
Timestamp:
2011-03-29T17:03:24+02:00 (13 years ago)
Author:
davestorkey
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • 2011WP/2011Stream2/OpenBoundaries

    v15 v16  
    6666   ln_vol         = .false.,.false.,        !  total volume correction (see volbdy parameter) 
    6767   ln_tides       = .false.,.false.,        !  Apply tidal harmonic forcing with Flather condition 
    68    nn_rimwidth    =  9, 1,                  !  width of the relaxation zone 
     68   nn_rimwidth    =  9, 1,                  !  width of the relaxation zone for Flow Relaxation 
     69   nn_dmp2d_in    =                         !  damping timescale (days) for 2D variables for inward radiation velocity or Flow Relaxation 
     70   nn_dmp2d_out   =                         !  damping timescale (days) for 2D variables outward radiation velocity 
     71   nn_dmp3d_in    =                         !  damping timescale (days) for 3D variables inward radiation velocity or Flow Relaxation 
     72   nn_dmp3d_out   =                         !  damping timescale (days) for 3D variables outward radiation velocity 
    6973   nn_dtactl      =  1, 1,                  !  = 0, bdy data are equal to the initial state 
    7074                                            !  = 1, bdy data are read in 'bdydata   .nc' files 
     
    280284== Boundary geometry and initialisation == 
    281285The internal code of the new module will use 1D unstructured arrays to specify the boundary as is currently done in BDY. For the T-points along the boundary these are '''nbit''', '''nbjt''' and '''nbrt''', which specify the (i,j) coordinates of each point and the discrete distance from the boundary. Each boundary set can be initialised using straight-line segments in '''obc_par.F90''' (ln_read=.false.) or read in from a '''coordinates_bdy.nc''' file (ln_read=.true.). If the boundary coordinates are read in from '''obc_par.F90''' then the '''nbi''', '''nbj''', '''nbr''' arrays will be generated internally (and if necessary a '''coordinates_bdy.nc''' file could be written out).  
     286 
     287For the Flow Relaxation Scheme we require relaxation to external data within the model, not just at the boundary. The BDY module includes these points in the specification of the boundary arrays and uses '''nbr''' to tell the discrete distance from the model boundary. The obc_par.F90 specification only specifies the points along the edge of the model domain. If '''nn_rimwidth > 1''' the internal points will be added to the BDY-style arrays internally. 
    282288 
    283289The '''obc_par.F90''' file will be the easiest way to define rectangular boundaries. More than one segment can be specified for each of the north, east, south and west boundaries as follows. Where there are more than one boundary set, the '''jpieob'''-type arrays are specified as 1D arrays and the '''nobcsege'''-type arrays are used to split the numbers between the different boundary sets.