Changeset 2034 for branches/DEV_r2006_merge_TRA_TRC/NEMO/OPA_SRC/TRA/traldf_bilapg.F90

Timestamp:

2010-07-29T17:05:35+02:00 (14 years ago)

Author:

cetlod

Message:

cosmetic changes

File:

• branches/DEV_r2006_merge_TRA_TRC/NEMO/OPA_SRC/TRA/traldf_bilapg.F90 (modified) (10 diffs)

branches/DEV_r2006_merge_TRA_TRC/NEMO/OPA_SRC/TRA/traldf_bilapg.F90

@@ -36 +36 @@
 #  include "ldfeiv_substitute.h90"
    !!----------------------------------------------------------------------
-   !!   OPA 9.0 , LOCEAN-IPSL (2005) 
+   !! NEMO/OPA 3.3 , LOCEAN-IPSL (2010) 
    !! $Id$ 
    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
@@ -43 +43 @@
 CONTAINS
 
-   SUBROUTINE tra_ldf_bilapg( kt, cdtype, ptrab, ptraa, kjpt )
+   SUBROUTINE tra_ldf_bilapg( kt, cdtype, ptb, pta, kjpt )
       !!----------------------------------------------------------------------
       !!                 ***  ROUTINE tra_ldf_bilapg  ***
@@ -55 +55 @@
       !!      computed in routine inildf.
       !!         -1- compute the geopotential harmonic operator applied to
-      !!        ptrab and multiply it by the eddy diffusivity coefficient
+      !!        ptb and multiply it by the eddy diffusivity coefficient
       !!       (done by a call to ldfght routine, result in wk1 arrays).
       !!      Applied the domain lateral boundary conditions by call to lbc_lnk
@@ -62 +62 @@
       !!      arrays).
       !!         -3- Add this trend to the general trend 
-      !!            ptraa = ptraa + wk2
-      !!
-      !! ** Action : - Update ptraa arrays with the before geopotential 
+      !!            pta = pta + wk2
+      !!
+      !! ** Action : - Update pta arrays with the before geopotential 
       !!               biharmonic mixing trend.
       !!----------------------------------------------------------------------
@@ -71 +71 @@
       CHARACTER(len=3), INTENT(in   )                                ::   cdtype         ! =TRA or TRC (tracer indicator)
       INTEGER         , INTENT(in   )                                ::   kjpt            ! number of tracers
-      REAL(wp)        , INTENT(in   ), DIMENSION(jpi,jpj,jpk,kjpt)   ::   ptrab          ! before and now tracer fields
-      REAL(wp)        , INTENT(inout), DIMENSION(jpi,jpj,jpk,kjpt)   ::   ptraa          ! tracer trend 
+      REAL(wp)        , INTENT(in   ), DIMENSION(jpi,jpj,jpk,kjpt)   ::   ptb          ! before and now tracer fields
+      REAL(wp)        , INTENT(inout), DIMENSION(jpi,jpj,jpk,kjpt)   ::   pta          ! tracer trend 
       !! * Local declarations
       INTEGER ::   ji, jj, jk, jn                 ! dummy loop indices
-      REAL(wp), DIMENSION(jpi,jpj,jpk,SIZE(ptrab,4)) ::   &
+      REAL(wp), DIMENSION(jpi,jpj,jpk,SIZE(ptb,4)) ::   &
          wk1, wk2                ! work array used for rotated biharmonic
          !                       ! operator on tracers and/or momentum
@@ -88 +88 @@
       !
 
-      ! 1. Laplacian of ptrab * aht
+      ! 1. Laplacian of ptb * aht
       ! ----------------------------- 
-      ! rotated harmonic operator applied to ptrab and multiply by aht ; output in wk1 
-
-      CALL ldfght( kt, cdtype, ptrab, wk1, kjpt, 1 )
+      ! rotated harmonic operator applied to ptb and multiply by aht ; output in wk1 
+
+      CALL ldfght( kt, cdtype, ptb, wk1, kjpt, 1 )
 
       !
@@ -100 +100 @@
       END DO
 
-      ! 2. Bilaplacian of ptrab
+      ! 2. Bilaplacian of ptb
       ! -------------------------
       ! rotated harmonic operator applied to wk1 ; output in wk2
@@ -117 +117 @@
                DO ji = 2, jpim1
                   ! add it to the general tracer trends
-                  ptraa(ji,jj,jk,jn) = ptraa(ji,jj,jk,jn) + wk2(ji,jj,jk,jn)
+                  pta(ji,jj,jk,jn) = pta(ji,jj,jk,jn) + wk2(ji,jj,jk,jn)
                END DO
             END DO
@@ -166 +166 @@
       !!
       !!----------------------------------------------------------------------
-      !! * Modules used
+      !! 
       USE oce         , zftv => ua     ! use ua as workspace
-      !! * Arguments
+      !!
       INTEGER         , INTENT(in )                              ::  kt      ! ocean time-step index
       CHARACTER(len=3), INTENT(in )                              ::  cdtype  ! =TRA or TRC (tracer indicator) 
@@ -178 +178 @@
       INTEGER         , INTENT(in )                              ::  kaht    !: =1 multiply the laplacian by the eddy diffusivity coeff.
       !                                                             !: =2 no multiplication 
-      !! * Local declarations
+      !!
       INTEGER ::   ji, jj, jk,jn          ! dummy loop indices
       !                                   ! temporary scalars