- Timestamp:
- 2015-09-13T09:42:41+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/NEMO/OPA_SRC/TRD/trdken.F90
r5656 r5737 96 96 nkstp = kt 97 97 DO jk = 1, jpkm1 98 bu (:,:,jk) = e1u(:,:) *e2u(:,:) * fse3u_n(:,:,jk)99 bv (:,:,jk) = e1v(:,:) *e2v(:,:) * fse3v_n(:,:,jk)98 bu (:,:,jk) = e1e2u(:,:) * fse3u_n(:,:,jk) 99 bv (:,:,jk) = e1e2v(:,:) * fse3v_n(:,:,jk) 100 100 r1_bt(:,:,jk) = 1._wp / ( e1e2t(:,:) * fse3t_n(:,:,jk) ) * tmask(:,:,jk) 101 101 END DO … … 263 263 ENDIF 264 264 ! ! allocate box volume arrays 265 IF 265 IF( trd_ken_alloc() /= 0 ) CALL ctl_stop('trd_ken_alloc: failed to allocate arrays') 266 266 ! 267 267 !!gm IF( .NOT. (ln_hpg_zco.OR.ln_hpg_zps) ) & 268 268 !!gm & CALL ctl_stop('trd_ken_init : only full and partial cells are coded for conversion rate') 269 269 ! 270 IF ( .NOT.lk_vvl ) THEN! constant volume: bu, bv, 1/bt computed one for all270 IF( .NOT.lk_vvl ) THEN ! constant volume: bu, bv, 1/bt computed one for all 271 271 DO jk = 1, jpkm1 272 bu (:,:,jk) = e1u(:,:) *e2u(:,:) * fse3u_n(:,:,jk)273 bv (:,:,jk) = e1v(:,:) *e2v(:,:) * fse3v_n(:,:,jk)272 bu (:,:,jk) = e1e2u(:,:) * fse3u_n(:,:,jk) 273 bv (:,:,jk) = e1e2v(:,:) * fse3v_n(:,:,jk) 274 274 r1_bt(:,:,jk) = 1._wp / ( e1e2t(:,:) * fse3t_n(:,:,jk) ) 275 275 END DO
Note: See TracChangeset
for help on using the changeset viewer.