Changeset 9923 for NEMO/branches/2018/dev_r9838_ENHANCE04_MLF/src/OCE/BDY/bdylib.F90

Timestamp:

2018-07-11T10:24:17+02:00 (6 years ago)

Author:

gm

Message:

#1911 (ENHANCE-04): step I.2: dev_r9838_ENHANCE04_MLF

File:

• NEMO/branches/2018/dev_r9838_ENHANCE04_MLF/src/OCE/BDY/bdylib.F90 (modified) (7 diffs)

NEMO/branches/2018/dev_r9838_ENHANCE04_MLF/src/OCE/BDY/bdylib.F90

@@ -4 +4 @@
    !! Unstructured Open Boundary Cond. :  Library module of generic boundary algorithms.
    !!======================================================================
-   !! History :  3.6  !  2013     (D. Storkey) original code
-   !!            4.0  !  2014     (T. Lovato) Generalize OBC structure
+   !! History :  3.6  !  2013     (D. Storkey)  original code
+   !!            4.0  !  2014     (T. Lovato)  Generalize OBC structure
    !!----------------------------------------------------------------------
+   
    !!----------------------------------------------------------------------
-   !!   bdy_orlanski_2d
-   !!   bdy_orlanski_3d
+   !!  bdy_frs        : Apply the Flow Relaxation Scheme (tracers)
+   !!  bdy_spe        : Apply a specified value (tracers)
+   !!  bdy_orl        : Apply Orlanski radiation (tracers)
+   !!  bdy_orlanski_2d:   2D      -        -         -
+   !!  bdy_orlanski_3d:   3D      -        -         -
+   !!  bdy_nmn        : Duplicate the value at open boundaries (zero gradient)
    !!----------------------------------------------------------------------
    USE oce            ! ocean dynamics and tracers 
@@ -22 +27 @@
    PRIVATE
 
-   PUBLIC   bdy_frs, bdy_spe, bdy_nmn, bdy_orl
-   PUBLIC   bdy_orlanski_2d
-   PUBLIC   bdy_orlanski_3d
+   PUBLIC   bdy_frs, bdy_spe, bdy_nmn
+   PUBLIC   bdy_orl, bdy_orlanski_2d, bdy_orlanski_3d
 
    !!----------------------------------------------------------------------
@@ -230 +234 @@
          ! Note no rdt factor in expression for zdt because it cancels in the expressions for 
          ! zrx and zry.
-         zdt = phia(iibm1,ijbm1) - phib(iibm1,ijbm1)
+         zdt =     phia(iibm1,ijbm1) - phib(iibm1,ijbm1)
          zdx = ( ( phia(iibm1,ijbm1) - phia(iibm2,ijbm2) ) / zex2 ) * zmask_x 
          zdy_1 = ( ( phib(iibm1   ,ijbm1   ) - phib(iibm1jm1,ijbm1jm1) ) / zey1 ) * zmask_y1    
@@ -247 +251 @@
          zout = sign( 1., zrx )
          zout = 0.5*( zout + abs(zout) )
-         zwgt = 2.*rdt*( (1.-zout) * idx%nbd(jb,igrd) + zout * idx%nbdout(jb,igrd) )
+         zwgt = rDt*( (1.-zout) * idx%nbd(jb,igrd) + zout * idx%nbdout(jb,igrd) )
          ! only apply radiation on outflow points 
          if( ll_npo ) then     !! NPO version !!
@@ -385 +389 @@
             ! Centred derivative is calculated as average of "left" and "right" derivatives for 
             ! this reason. 
-            zdt = phia(iibm1,ijbm1,jk) - phib(iibm1,ijbm1,jk)
+            zdt =     phia(iibm1,ijbm1,jk) - phib(iibm1,ijbm1,jk)
             zdx = ( ( phia(iibm1,ijbm1,jk) - phia(iibm2,ijbm2,jk) ) / zex2 ) * zmask_x                  
             zdy_1 = ( ( phib(iibm1   ,ijbm1   ,jk) - phib(iibm1jm1,ijbm1jm1,jk) ) / zey1 ) * zmask_y1  
@@ -403 +407 @@
             zout = sign( 1., zrx )
             zout = 0.5*( zout + abs(zout) )
-            zwgt = 2.*rdt*( (1.-zout) * idx%nbd(jb,igrd) + zout * idx%nbdout(jb,igrd) )
+            zwgt = rDt*( (1.-zout) * idx%nbd(jb,igrd) + zout * idx%nbdout(jb,igrd) )
             ! only apply radiation on outflow points 
             if( ll_npo ) then     !! NPO version !!
@@ -426 +430 @@
       !
    END SUBROUTINE bdy_orlanski_3d
+
 
    SUBROUTINE bdy_nmn( idx, igrd, phia )