- Timestamp:
- 2017-03-17T08:46:30+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/TOP_SRC/TRP/trcadv.F90
r5602 r7806 89 89 r2dt(:) = 2. * rdttrc(:) ! = 2 rdttrc (leapfrog) 90 90 ENDIF 91 ! ! effective transport 92 DO jk = 1, jpkm1 93 ! ! eulerian transport only 94 zun(:,:,jk) = e2u (:,:) * fse3u(:,:,jk) * un(:,:,jk) 95 zvn(:,:,jk) = e1v (:,:) * fse3v(:,:,jk) * vn(:,:,jk) 96 zwn(:,:,jk) = e1e2t(:,:) * wn(:,:,jk) 97 ! 98 END DO 99 ! 100 IF( ln_vvl_ztilde .OR. ln_vvl_layer ) THEN 101 zun(:,:,:) = zun(:,:,:) + un_td(:,:,:) 102 zvn(:,:,:) = zvn(:,:,:) + vn_td(:,:,:) 91 ! 92 IF( lk_offline ) THEN 93 zun(:,:,:) = un(:,:,:) ! effective transport already in un/vn/wn 94 zvn(:,:,:) = vn(:,:,:) 95 zwn(:,:,:) = wn(:,:,:) 96 ELSE 97 ! ! effective transport 98 DO jk = 1, jpkm1 99 ! ! eulerian transport only 100 zun(:,:,jk) = e2u (:,:) * fse3u(:,:,jk) * un(:,:,jk) 101 zvn(:,:,jk) = e1v (:,:) * fse3v(:,:,jk) * vn(:,:,jk) 102 zwn(:,:,jk) = e1e2t(:,:) * wn(:,:,jk) 103 ! 104 END DO 105 ! 106 IF( ln_vvl_ztilde .OR. ln_vvl_layer ) THEN 107 zun(:,:,:) = zun(:,:,:) + un_td(:,:,:) 108 zvn(:,:,:) = zvn(:,:,:) + vn_td(:,:,:) 109 ENDIF 110 ! 111 zun(:,:,jpk) = 0._wp ! no transport trough the bottom 112 zvn(:,:,jpk) = 0._wp ! no transport trough the bottom 113 zwn(:,:,jpk) = 0._wp ! no transport trough the bottom 114 ! 115 116 IF( lk_traldf_eiv .AND. .NOT. ln_traldf_grif ) & ! add the eiv transport (if necessary) 117 & CALL tra_adv_eiv( kt, nittrc000, zun, zvn, zwn, 'TRC' ) 118 ! 119 IF( ln_mle ) CALL tra_adv_mle( kt, nittrc000, zun, zvn, zwn, 'TRC' ) ! add the mle transport (if necessary) 120 ! 103 121 ENDIF 104 !105 zun(:,:,jpk) = 0._wp ! no transport trough the bottom106 zvn(:,:,jpk) = 0._wp ! no transport trough the bottom107 zwn(:,:,jpk) = 0._wp ! no transport trough the bottom108 109 IF( lk_traldf_eiv .AND. .NOT. ln_traldf_grif ) & ! add the eiv transport (if necessary)110 & CALL tra_adv_eiv( kt, nittrc000, zun, zvn, zwn, 'TRC' )111 !112 IF( ln_mle ) CALL tra_adv_mle( kt, nittrc000, zun, zvn, zwn, 'TRC' ) ! add the mle transport (if necessary)113 122 ! 114 123 SELECT CASE ( nadv ) !== compute advection trend and add it to general trend ==!
Note: See TracChangeset
for help on using the changeset viewer.