Changeset 2800 for branches/2011/UKMO_MERCATOR_obc_bdy_merge/NEMOGCM/NEMO/OPA_SRC/OBC/obcdyn2d.F90

Timestamp:

2011-07-13T15:31:05+02:00 (13 years ago)

Author:

davestorkey

Message:

1. Application of boundary conditions to barotropic and baroclinic velocities clearly separated.
2. Option to input full velocities in boundary data (default expects barotropic and baroclinic velocities separately).
3. Option to use initial conditions as boundary conditions coded.

File:

• branches/2011/UKMO_MERCATOR_obc_bdy_merge/NEMOGCM/NEMO/OPA_SRC/OBC/obcdyn2d.F90 (modified) (7 diffs)

branches/2011/UKMO_MERCATOR_obc_bdy_merge/NEMOGCM/NEMO/OPA_SRC/OBC/obcdyn2d.F90

@@ -39 +39 @@
 CONTAINS
 
-   SUBROUTINE obc_dyn2d( kt, pssh )
+   SUBROUTINE obc_dyn2d( kt )
       !!----------------------------------------------------------------------
       !!                  ***  SUBROUTINE obc_dyn2d  ***
@@ -47 +47 @@
       !!----------------------------------------------------------------------
       INTEGER,                      INTENT(in) ::   kt   ! Main time step counter
-      REAL(wp), DIMENSION(jpi,jpj), INTENT(in) ::   pssh
       !!
       INTEGER                                  ::   ib_obc ! Loop counter
@@ -56 +55 @@
          CASE(jp_none)
             CYCLE
+         CASE(jp_frs)
+            CALL obc_dyn2d_frs( idx_obc(ib_obc), dta_obc(ib_obc), kt )
          CASE(jp_flather)
-            CALL obc_dyn2d_fla( idx_obc(ib_obc), dta_obc(ib_obc), pssh )
+            CALL obc_dyn2d_fla( idx_obc(ib_obc), dta_obc(ib_obc) )
          CASE DEFAULT
             CALL ctl_stop( 'obc_dyn3d : unrecognised option for open boundaries for barotropic variables' )
@@ -65 +66 @@
    END SUBROUTINE obc_dyn2d
 
-   SUBROUTINE obc_dyn2d_fla( idx, dta, pssh )
-      !!----------------------------------------------------------------------
-      !!                 ***  SUBROUTINE obc_dyn_fla  ***
+   SUBROUTINE obc_dyn2d_frs( idx, dta, kt )
+      !!----------------------------------------------------------------------
+      !!                  ***  SUBROUTINE obc_dyn2d_frs  ***
+      !!
+      !! ** Purpose : - Apply the Flow Relaxation Scheme for barotropic velocities
+      !!                at open boundaries.
+      !!
+      !! References :- Engedahl H., 1995: Use of the flow relaxation scheme in 
+      !!               a three-dimensional baroclinic ocean model with realistic
+      !!               topography. Tellus, 365-382.
+      !!----------------------------------------------------------------------
+      INTEGER,         INTENT(in) :: kt
+      TYPE(OBC_INDEX), INTENT(in) ::   idx  ! OBC indices
+      TYPE(OBC_DATA),  INTENT(in) ::   dta  ! OBC external data
+      !!
+      INTEGER  ::   jb, jk         ! dummy loop indices
+      INTEGER  ::   ii, ij, igrd   ! local integers
+      REAL(wp) ::   zwgt           ! boundary weight
+      !!----------------------------------------------------------------------
+      !
+      !
+      igrd = 2                      ! Relaxation of zonal velocity
+      DO jb = 1, idx%nblen(igrd)
+         ii   = idx%nbi(jb,igrd)
+         ij   = idx%nbj(jb,igrd)
+         zwgt = idx%nbw(jb,igrd)
+         pu2d(ii,ij) = ( pu2d(ii,ij) + zwgt * ( dta%u2d(jb) - pu2d(ii,ij) ) ) * umask(ii,ij,1)
+      END DO
+      !
+      igrd = 3                      ! Relaxation of meridional velocity
+      DO jb = 1, idx%nblen(igrd)
+         ii   = idx%nbi(jb,igrd)
+         ij   = idx%nbj(jb,igrd)
+         zwgt = idx%nbw(jb,igrd)
+         pv2d(ii,ij) = ( pv2d(ii,ij) + zwgt * ( dta%v2d(jb) - pv2d(ii,ij) ) ) * vmask(ii,ij,1)
+      END DO 
+      CALL lbc_lnk( pu2d, 'U', -1. ) 
+      CALL lbc_lnk( pv2d, 'V', -1. )   ! Boundary points should be updated
+      !
+
+   END SUBROUTINE obc_dyn2d_frs
+
+
+   SUBROUTINE obc_dyn2d_fla( idx, dta )
+      !!----------------------------------------------------------------------
+      !!                 ***  SUBROUTINE obc_dyn2d_fla  ***
       !!             
       !!              - Apply Flather boundary conditions on normal barotropic velocities 
@@ -100 +144 @@
      
 !!!!!!!!!!!! SOME WORK TO DO ON THE TIDES HERE !!!!!!!!!!!!!
+
+!!! REPLACE spgu with nemo_wrk work space
 
       ! Fill temporary array with ssh data (here spgu):
@@ -120 +166 @@
          iip1 = ii - MIN( 0, INT( idx%flagu(jb) ) )   ! T pts i-indice outside the boundary 
          !
-         zcorr = - idx%flagu(jb) * SQRT( grav * hur_e(ii, ij) ) * ( pssh(iim1, ij) - spgu(iip1,ij) )
+         zcorr = - idx%flagu(jb) * SQRT( grav * phur(ii, ij) ) * ( pssh(iim1, ij) - spgu(iip1,ij) )
 !!$         zforc = dta%u2d(jb) + utide(jb)
          zforc = dta%u2d(jb)
-         ua_e(ii,ij) = zforc + zcorr * umask(ii,ij,1) 
+         pu2d(ii,ij) = zforc + zcorr * umask(ii,ij,1) 
       END DO
       !
@@ -134 +180 @@
          ijp1 = ij - MIN( 0, INT( idx%flagv(jb) ) )   ! T pts j-indice outside the boundary 
          !
-         zcorr = - idx%flagv(jb) * SQRT( grav * hvr_e(ii, ij) ) * ( pssh(ii, ijm1) - spgu(ii,ijp1) )
+         zcorr = - idx%flagv(jb) * SQRT( grav * phvr(ii, ij) ) * ( pssh(ii, ijm1) - spgu(ii,ijp1) )
 !!$         zforc = dta%v2d(jb) + vtide(jb)
          zforc = dta%v2d(jb)
-         va_e(ii,ij) = zforc + zcorr * vmask(ii,ij,1)
-      END DO
-      CALL lbc_lnk( ua_e, 'U', -1. )   ! Boundary points should be updated
-      CALL lbc_lnk( va_e, 'V', -1. )   !
+         pv2d(ii,ij) = zforc + zcorr * vmask(ii,ij,1)
+      END DO
+      CALL lbc_lnk( pu2d, 'U', -1. )   ! Boundary points should be updated
+      CALL lbc_lnk( pv2d, 'V', -1. )   !
       !
       !