Changeset 6140 for trunk/NEMOGCM/NEMO/OPA_SRC/DYN/divhor.F90
- Timestamp:
- 2015-12-21T12:35:23+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/DYN/divhor.F90
r5836 r6140 20 20 USE oce ! ocean dynamics and tracers 21 21 USE dom_oce ! ocean space and time domain 22 USE sbc_oce, ONLY : ln_rnf, nn_isf ! surface boundary condition: ocean22 USE sbc_oce, ONLY : ln_rnf, ln_isf ! surface boundary condition: ocean 23 23 USE sbcrnf ! river runoff 24 24 USE sbcisf ! ice shelf 25 USE iscplhsb ! ice sheet / ocean coupling 26 USE iscplini ! ice sheet / ocean coupling 25 27 ! 26 28 USE in_out_manager ! I/O manager … … 36 38 37 39 !! * Substitutions 38 # include "domzgr_substitute.h90"39 40 # include "vectopt_loop_substitute.h90" 40 41 !!---------------------------------------------------------------------- … … 74 75 DO jj = 2, jpjm1 75 76 DO ji = fs_2, fs_jpim1 ! vector opt. 76 hdivn(ji,jj,jk) = ( e2u(ji ,jj) * fse3u_n(ji ,jj,jk) * un(ji ,jj,jk) &77 & - e2u(ji-1,jj) * fse3u_n(ji-1,jj,jk) * un(ji-1,jj,jk) &78 & + e1v(ji,jj ) * fse3v_n(ji,jj ,jk) * vn(ji,jj ,jk) &79 & - e1v(ji,jj-1) * fse3v_n(ji,jj-1,jk) * vn(ji,jj-1,jk) ) &80 & / ( e1e2t(ji,jj) * fse3t_n(ji,jj,jk) )77 hdivn(ji,jj,jk) = ( e2u(ji ,jj) * e3u_n(ji ,jj,jk) * un(ji ,jj,jk) & 78 & - e2u(ji-1,jj) * e3u_n(ji-1,jj,jk) * un(ji-1,jj,jk) & 79 & + e1v(ji,jj ) * e3v_n(ji,jj ,jk) * vn(ji,jj ,jk) & 80 & - e1v(ji,jj-1) * e3v_n(ji,jj-1,jk) * vn(ji,jj-1,jk) ) & 81 & / ( e1e2t(ji,jj) * e3t_n(ji,jj,jk) ) 81 82 END DO 82 83 END DO … … 89 90 END DO 90 91 ! 91 IF( ln_rnf 92 IF( ln_rnf ) CALL sbc_rnf_div( hdivn ) !== runoffs ==! (update hdivn field) 92 93 ! 93 IF( ln_ divisf .AND. nn_isf > 0) CALL sbc_isf_div( hdivn ) !== ice shelf ==! (update hdivn field)94 IF( ln_isf ) CALL sbc_isf_div( hdivn ) !== ice shelf ==! (update hdivn field) 94 95 ! 95 CALL lbc_lnk( hdivn, 'T', 1. ) !== lateral boundary cond. ==! (no sign change) 96 IF( ln_iscpl .AND. ln_hsb ) CALL iscpl_div( hdivn ) !== ice sheet ==! (update hdivn field) 97 ! 98 CALL lbc_lnk( hdivn, 'T', 1. ) !== lateral boundary cond. ==! (no sign change) 96 99 ! 97 100 IF( nn_timing == 1 ) CALL timing_stop('div_hor')
Note: See TracChangeset
for help on using the changeset viewer.