- Timestamp:
- 2015-05-18T17:19:50+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DYN/divcur.F90
r4990 r5282 17 17 !! 3.3 ! 2010-09 (D.Storkey and E.O'Dea) bug fixes for BDY module 18 18 !! - ! 2010-10 (R. Furner, G. Madec) runoff and cla added directly here 19 !! 3.6 ! 2014-15 cla removed 19 20 !!---------------------------------------------------------------------- 20 21 … … 28 29 USE sbcrnf ! river runoff 29 30 USE sbcisf ! ice shelf 30 USE cla ! cross land advection (cla_div routine)31 31 USE in_out_manager ! I/O manager 32 32 USE lbclnk ! ocean lateral boundary conditions (or mpp link) … … 68 68 !! hdivn = 1/(e1t*e2t*e3t) ( di[e2u*e3u un] + dj[e1v*e3v vn] ) 69 69 !! correct hdiv with runoff inflow (div_rnf), ice shelf melting (div_isf) 70 !! and cross land flow (div_cla)70 !! [REMOVED cross land flow (div_cla)] 71 71 !! II. vorticity : 72 72 !! - save the curl computed at the previous time-step … … 229 229 IF( ln_rnf ) CALL sbc_rnf_div( hdivn ) ! runoffs (update hdivn field) 230 230 IF( ln_divisf .AND. (nn_isf /= 0) ) CALL sbc_isf_div( hdivn ) ! ice shelf (update hdivn field) 231 IF( nn_cla == 1 ) CALL cla_div ( kt ) ! Cross Land Advection (Update Hor. divergence) 232 231 233 232 ! 4. Lateral boundary conditions on hdivn and rotn 234 233 ! ---------------------------------=======---====== … … 259 258 !! - compute the now divergence given by : 260 259 !! hdivn = 1/(e1t*e2t*e3t) ( di[e2u*e3u un] + dj[e1v*e3v vn] ) 261 !! correct hdiv with runoff inflow (div_rnf) and cross land flow (div_cla) 260 !! correct hdiv with runoff inflow (div_rnf) 261 !! [REMOVED cross land flow (div_cla)] 262 262 !! - Relavtive Vorticity : 263 263 !! - save the curl computed at the previous time-step (rotb = rotn) … … 328 328 IF( ln_rnf ) CALL sbc_rnf_div( hdivn ) ! runoffs (update hdivn field) 329 329 IF( ln_divisf .AND. (nn_isf .GT. 0) ) CALL sbc_isf_div( hdivn ) ! ice shelf (update hdivn field) 330 IF( nn_cla == 1 ) CALL cla_div ( kt ) ! Cross Land Advection (update hdivn field)331 330 ! 332 331 CALL lbc_lnk( hdivn, 'T', 1. ) ; CALL lbc_lnk( rotn , 'F', 1. ) ! lateral boundary cond. (no sign change)
Note: See TracChangeset
for help on using the changeset viewer.