New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 2034 for branches/DEV_r2006_merge_TRA_TRC/NEMO/OPA_SRC/TRA/trazdf_imp.F90 – NEMO

Ignore:
Timestamp:
2010-07-29T17:05:35+02:00 (14 years ago)
Author:
cetlod
Message:

cosmetic changes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/DEV_r2006_merge_TRA_TRC/NEMO/OPA_SRC/TRA/trazdf_imp.F90

    r2024 r2034  
    4343#  include "vectopt_loop_substitute.h90" 
    4444   !!---------------------------------------------------------------------- 
    45    !! NEMO/OPA 3.2 , LOCEAN-IPSL (2009)  
     45   !! NEMO/OPA 3.3 , LOCEAN-IPSL (2010)  
    4646   !! $Id$ 
    4747   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 
     
    4949CONTAINS 
    5050  
    51    SUBROUTINE tra_zdf_imp( kt    , cdtype, p2dt,    & 
    52       &                    ptrab , ptraa , kjpt     ) 
     51   SUBROUTINE tra_zdf_imp( kt, cdtype, p2dt, ptb, pta, kjpt )  
    5352      !!---------------------------------------------------------------------- 
    5453      !!                  ***  ROUTINE tra_zdf_imp  *** 
     
    8988      !! 
    9089      !!--------------------------------------------------------------------- 
    91       !! * Modules used 
     90      !!  
    9291      USE oce    , ONLY :   zwd   => ua   ! ua used as workspace 
    9392      USE oce    , ONLY :   zws   => va   ! va  -          - 
    94       !! * Arguments 
     93      !!  
    9594      INTEGER         , INTENT(in   )                                ::   kt             ! ocean time-step index 
    9695      CHARACTER(len=3), INTENT(in   )                                ::   cdtype         ! =TRA or TRC (tracer indicator) 
    9796      INTEGER         , INTENT(in   )                                ::   kjpt            ! number of tracers 
    9897      REAL(wp)        , INTENT(in   ), DIMENSION(jpk)                ::   p2dt        ! vertical profile of tracer time-step 
    99       REAL(wp)        , INTENT(in   ), DIMENSION(jpi,jpj,jpk,kjpt)   ::   ptrab          ! before and now tracer fields 
    100       REAL(wp)        , INTENT(inout), DIMENSION(jpi,jpj,jpk,kjpt)   ::   ptraa          ! tracer trend  
     98      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  
    101100      !! 
    102101      INTEGER  ::  ji, jj, jk, jn        ! dummy loop indices 
     
    265264               ze3tb = ( 1. - znvvl ) + znvvl * fse3t_b(ji,jj,1) 
    266265               ze3tn = ( 1. - znvvl ) + znvvl * fse3t(ji,jj,1) 
    267                ptraa(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) 
    268267            END DO 
    269268         END DO 
     
    273272                  ze3tb = ( 1. - znvvl ) + znvvl * fse3t_b(ji,jj,jk) 
    274273                  ze3tn = ( 1. - znvvl ) + znvvl * fse3t  (ji,jj,jk) 
    275                   zrhs = ze3tb * ptrab(ji,jj,jk,jn) + p2dt(jk) * ze3tn * ptraa(ji,jj,jk,jn)   ! zrhs=right hand side  
    276                   ptraa(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) 
    277276               END DO 
    278277            END DO 
     
    284283         DO jj = 2, jpjm1 
    285284            DO ji = fs_2, fs_jpim1 
    286                ptraa(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) 
    287286            END DO 
    288287         END DO 
     
    290289            DO jj = 2, jpjm1 
    291290               DO ji = fs_2, fs_jpim1 
    292                   ptraa(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) ) & 
    293292                  &                    / zwt(ji,jj,jk) * tmask(ji,jj,jk) 
    294293               END DO 
Note: See TracChangeset for help on using the changeset viewer.