- Timestamp:
- 2010-07-29T17:05:35+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r2006_merge_TRA_TRC/NEMO/OPA_SRC/TRA/trazdf_imp.F90
r2024 r2034 43 43 # include "vectopt_loop_substitute.h90" 44 44 !!---------------------------------------------------------------------- 45 !! NEMO/OPA 3. 2 , LOCEAN-IPSL (2009)45 !! NEMO/OPA 3.3 , LOCEAN-IPSL (2010) 46 46 !! $Id$ 47 47 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) … … 49 49 CONTAINS 50 50 51 SUBROUTINE tra_zdf_imp( kt , cdtype, p2dt, & 52 & ptrab , ptraa , kjpt ) 51 SUBROUTINE tra_zdf_imp( kt, cdtype, p2dt, ptb, pta, kjpt ) 53 52 !!---------------------------------------------------------------------- 54 53 !! *** ROUTINE tra_zdf_imp *** … … 89 88 !! 90 89 !!--------------------------------------------------------------------- 91 !! * Modules used90 !! 92 91 USE oce , ONLY : zwd => ua ! ua used as workspace 93 92 USE oce , ONLY : zws => va ! va - - 94 !! * Arguments93 !! 95 94 INTEGER , INTENT(in ) :: kt ! ocean time-step index 96 95 CHARACTER(len=3), INTENT(in ) :: cdtype ! =TRA or TRC (tracer indicator) 97 96 INTEGER , INTENT(in ) :: kjpt ! number of tracers 98 97 REAL(wp) , INTENT(in ), DIMENSION(jpk) :: p2dt ! vertical profile of tracer time-step 99 REAL(wp) , INTENT(in ), DIMENSION(jpi,jpj,jpk,kjpt) :: pt rab ! before and now tracer fields100 REAL(wp) , INTENT(inout), DIMENSION(jpi,jpj,jpk,kjpt) :: pt raa ! tracer trend98 REAL(wp) , INTENT(in ), DIMENSION(jpi,jpj,jpk,kjpt) :: ptb ! before and now tracer fields 99 REAL(wp) , INTENT(inout), DIMENSION(jpi,jpj,jpk,kjpt) :: pta ! tracer trend 101 100 !! 102 101 INTEGER :: ji, jj, jk, jn ! dummy loop indices … … 265 264 ze3tb = ( 1. - znvvl ) + znvvl * fse3t_b(ji,jj,1) 266 265 ze3tn = ( 1. - znvvl ) + znvvl * fse3t(ji,jj,1) 267 pt raa(ji,jj,1,jn) = ze3tb * ptrab(ji,jj,1,jn) + p2dt(1) * ze3tn * ptraa(ji,jj,1,jn)266 pta(ji,jj,1,jn) = ze3tb * ptb(ji,jj,1,jn) + p2dt(1) * ze3tn * pta(ji,jj,1,jn) 268 267 END DO 269 268 END DO … … 273 272 ze3tb = ( 1. - znvvl ) + znvvl * fse3t_b(ji,jj,jk) 274 273 ze3tn = ( 1. - znvvl ) + znvvl * fse3t (ji,jj,jk) 275 zrhs = ze3tb * pt rab(ji,jj,jk,jn) + p2dt(jk) * ze3tn * ptraa(ji,jj,jk,jn) ! zrhs=right hand side276 pt raa(ji,jj,jk,jn) = zrhs - zwi(ji,jj,jk) / zwt(ji,jj,jk-1) * ptraa(ji,jj,jk-1,jn)274 zrhs = ze3tb * ptb(ji,jj,jk,jn) + p2dt(jk) * ze3tn * pta(ji,jj,jk,jn) ! zrhs=right hand side 275 pta(ji,jj,jk,jn) = zrhs - zwi(ji,jj,jk) / zwt(ji,jj,jk-1) * pta(ji,jj,jk-1,jn) 277 276 END DO 278 277 END DO … … 284 283 DO jj = 2, jpjm1 285 284 DO ji = fs_2, fs_jpim1 286 pt raa(ji,jj,jpkm1,jn) = ptraa(ji,jj,jpkm1,jn) / zwt(ji,jj,jpkm1) * tmask(ji,jj,jpkm1)285 pta(ji,jj,jpkm1,jn) = pta(ji,jj,jpkm1,jn) / zwt(ji,jj,jpkm1) * tmask(ji,jj,jpkm1) 287 286 END DO 288 287 END DO … … 290 289 DO jj = 2, jpjm1 291 290 DO ji = fs_2, fs_jpim1 292 pt raa(ji,jj,jk,jn) = ( ptraa(ji,jj,jk,jn) - zws(ji,jj,jk) * ptraa(ji,jj,jk+1,jn) ) &291 pta(ji,jj,jk,jn) = ( pta(ji,jj,jk,jn) - zws(ji,jj,jk) * pta(ji,jj,jk+1,jn) ) & 293 292 & / zwt(ji,jj,jk) * tmask(ji,jj,jk) 294 293 END DO
Note: See TracChangeset
for help on using the changeset viewer.