- Timestamp:
- 2011-07-27T14:41:28+02:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/UKMO_MERCATOR_obc_bdy_merge/NEMOGCM/NEMO/OPA_SRC/OBC/obcdyn2d.F90
r2800 r2814 24 24 USE phycst ! physical constants 25 25 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 26 USE obctides ! for tidal harmonic forcing at boundary27 26 USE in_out_manager ! 28 27 … … 30 29 PRIVATE 31 30 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 33 32 34 33 !!---------------------------------------------------------------------- … … 77 76 !! topography. Tellus, 365-382. 78 77 !!---------------------------------------------------------------------- 79 INTEGER, INTENT(in) :: kt78 INTEGER, INTENT(in) :: kt 80 79 TYPE(OBC_INDEX), INTENT(in) :: idx ! OBC indices 81 80 TYPE(OBC_DATA), INTENT(in) :: dta ! OBC external data … … 114 113 !! 115 114 !! - Apply Flather boundary conditions on normal barotropic velocities 116 !! (ln_dyn_fla=.true. or ln_tides=.true.)117 115 !! 118 116 !! ** WARNINGS about FLATHER implementation: … … 143 141 ! ---------------------------------! 144 142 145 !!!!!!!!!!!! SOME WORK TO DO ON THE TIDES HERE !!!!!!!!!!!!!146 147 143 !!! REPLACE spgu with nemo_wrk work space 148 144 … … 154 150 ij = idx%nbj(jb,igrd) 155 151 spgu(ii, ij) = dta%ssh(jb) 156 !!$ IF( ln_tides ) spgu(ii, ij) = spgu(ii, ij) + sshtide(jb)157 152 END DO 158 153 ! … … 167 162 ! 168 163 zcorr = - idx%flagu(jb) * SQRT( grav * phur(ii, ij) ) * ( pssh(iim1, ij) - spgu(iip1,ij) ) 169 !!$ zforc = dta%u2d(jb) + utide(jb)170 164 zforc = dta%u2d(jb) 171 165 pu2d(ii,ij) = zforc + zcorr * umask(ii,ij,1) … … 181 175 ! 182 176 zcorr = - idx%flagv(jb) * SQRT( grav * phvr(ii, ij) ) * ( pssh(ii, ijm1) - spgu(ii,ijp1) ) 183 !!$ zforc = dta%v2d(jb) + vtide(jb)184 177 zforc = dta%v2d(jb) 185 178 pv2d(ii,ij) = zforc + zcorr * vmask(ii,ij,1)
Note: See TracChangeset
for help on using the changeset viewer.