Changeset 5972 for branches/2014/dev_r4650_UKMO14.5_SST_BIAS_CORRECTION/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv_cen2.F90
- Timestamp:
- 2015-12-02T09:52:20+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO14.5_SST_BIAS_CORRECTION/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv_cen2.F90
r5967 r5972 90 90 DO jj = 2, jpjm1 ! divergence of horizontal momentum fluxes 91 91 DO ji = fs_2, fs_jpim1 ! vector opt. 92 zbu = e1 u(ji,jj) *e2u(ji,jj) * fse3u(ji,jj,jk)93 zbv = e1 v(ji,jj) *e2v(ji,jj) * fse3v(ji,jj,jk)92 zbu = e1e2u(ji,jj) * fse3u(ji,jj,jk) 93 zbv = e1e2v(ji,jj) * fse3v(ji,jj,jk) 94 94 ! 95 95 ua(ji,jj,jk) = ua(ji,jj,jk) - ( zfu_t(ji+1,jj ,jk) - zfu_t(ji ,jj ,jk) & … … 114 114 DO jk = 1, jpkm1 ! ==================== ! 115 115 ! ! Vertical volume fluxesÊ 116 zfw(:,:,jk) = 0.25 * e1 t(:,:) *e2t(:,:) * wn(:,:,jk)116 zfw(:,:,jk) = 0.25 * e1e2t(:,:) * wn(:,:,jk) 117 117 ! 118 118 IF( jk == 1 ) THEN ! surface/bottom advective fluxes … … 144 144 DO ji = fs_2, fs_jpim1 ! vector opt. 145 145 ua(ji,jj,jk) = ua(ji,jj,jk) - ( zfu_uw(ji,jj,jk) - zfu_uw(ji,jj,jk+1) ) & 146 & / ( e1 u(ji,jj) *e2u(ji,jj) * fse3u(ji,jj,jk) )146 & / ( e1e2u(ji,jj) * fse3u(ji,jj,jk) ) 147 147 va(ji,jj,jk) = va(ji,jj,jk) - ( zfv_vw(ji,jj,jk) - zfv_vw(ji,jj,jk+1) ) & 148 & / ( e1 v(ji,jj) *e2v(ji,jj) * fse3v(ji,jj,jk) )148 & / ( e1e2v(ji,jj) * fse3v(ji,jj,jk) ) 149 149 END DO 150 150 END DO
Note: See TracChangeset
for help on using the changeset viewer.