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.
Changeset 2814 for branches/2011/UKMO_MERCATOR_obc_bdy_merge/NEMOGCM/NEMO/OPA_SRC/OBC/obcdyn2d.F90 – NEMO

Ignore:
Timestamp:
2011-07-27T14:41:28+02:00 (13 years ago)
Author:
davestorkey
Message:
  1. Implement tidal harmonics forcing (UKMO version) in new structure.
  2. Other bug fixes and updates.
File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2011/UKMO_MERCATOR_obc_bdy_merge/NEMOGCM/NEMO/OPA_SRC/OBC/obcdyn2d.F90

    r2800 r2814  
    2424   USE phycst          ! physical constants 
    2525   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
    26    USE obctides        ! for tidal harmonic forcing at boundary 
    2726   USE in_out_manager  ! 
    2827 
     
    3029   PRIVATE 
    3130 
    32    PUBLIC   obc_dyn2d     ! routine called in dynspg_ts (time splitting case ONLY) 
     31   PUBLIC   obc_dyn2d     ! routine called in dynspg_ts and dyn_nxt and dynspg_flt 
    3332 
    3433   !!---------------------------------------------------------------------- 
     
    7776      !!               topography. Tellus, 365-382. 
    7877      !!---------------------------------------------------------------------- 
    79       INTEGER,         INTENT(in) :: kt 
     78      INTEGER,         INTENT(in) ::   kt 
    8079      TYPE(OBC_INDEX), INTENT(in) ::   idx  ! OBC indices 
    8180      TYPE(OBC_DATA),  INTENT(in) ::   dta  ! OBC external data 
     
    114113      !!              
    115114      !!              - Apply Flather boundary conditions on normal barotropic velocities  
    116       !!                (ln_dyn_fla=.true. or ln_tides=.true.) 
    117115      !! 
    118116      !! ** WARNINGS about FLATHER implementation: 
     
    143141      ! ---------------------------------!  
    144142      
    145 !!!!!!!!!!!! SOME WORK TO DO ON THE TIDES HERE !!!!!!!!!!!!! 
    146  
    147143!!! REPLACE spgu with nemo_wrk work space 
    148144 
     
    154150         ij = idx%nbj(jb,igrd) 
    155151         spgu(ii, ij) = dta%ssh(jb) 
    156 !!$         IF( ln_tides )   spgu(ii, ij) = spgu(ii, ij) + sshtide(jb) 
    157152      END DO 
    158153      ! 
     
    167162         ! 
    168163         zcorr = - idx%flagu(jb) * SQRT( grav * phur(ii, ij) ) * ( pssh(iim1, ij) - spgu(iip1,ij) ) 
    169 !!$         zforc = dta%u2d(jb) + utide(jb) 
    170164         zforc = dta%u2d(jb) 
    171165         pu2d(ii,ij) = zforc + zcorr * umask(ii,ij,1)  
     
    181175         ! 
    182176         zcorr = - idx%flagv(jb) * SQRT( grav * phvr(ii, ij) ) * ( pssh(ii, ijm1) - spgu(ii,ijp1) ) 
    183 !!$         zforc = dta%v2d(jb) + vtide(jb) 
    184177         zforc = dta%v2d(jb) 
    185178         pv2d(ii,ij) = zforc + zcorr * vmask(ii,ij,1) 
Note: See TracChangeset for help on using the changeset viewer.