Changeset 503 for trunk/NEMO/OPA_SRC/TRA/traldf_iso.F90

Timestamp:

2006-09-27T10:52:29+02:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_064 : CT : general trends update including the addition of mean windows analysis possibility in the mixed layer

File:

• trunk/NEMO/OPA_SRC/TRA/traldf_iso.F90 (modified) (10 diffs)

trunk/NEMO/OPA_SRC/TRA/traldf_iso.F90

@@ -1 +1 @@
 MODULE traldf_iso
-   !!==============================================================================
+   !!======================================================================
    !!                   ***  MODULE  traldf_iso  ***
    !! Ocean active tracers:  horizontal component of the lateral tracer mixing trend
-   !!==============================================================================
+   !!======================================================================
+   !! History :        !  94-08  (G. Madec, M. Imbard)
+   !!                  !  97-05  (G. Madec)  split into traldf and trazdf
+   !!             8.5  !  02-08  (G. Madec)  Free form, F90
+   !!             9.0  !  05-11  (G. Madec)  merge traldf and trazdf :-)
+   !!----------------------------------------------------------------------
 #if   defined key_ldfslp   ||   defined key_esopa
    !!----------------------------------------------------------------------
    !!   'key_ldfslp'               slope of the lateral diffusive direction
+   !!----------------------------------------------------------------------
    !!----------------------------------------------------------------------
    !!   tra_ldf_iso  : update the tracer trend with the horizontal 
@@ -14 +20 @@
    !!                  vector optimization, use k-j-i loops.
    !!----------------------------------------------------------------------
-   !! * Modules used
    USE oce             ! ocean dynamics and active tracers
    USE dom_oce         ! ocean space and time domain
@@ -26 +31 @@
    USE prtctl          ! Print control
 
-
    IMPLICIT NONE
    PRIVATE
 
-   !! * Accessibility
-   PUBLIC tra_ldf_iso   ! routine called by step.F90
+   PUBLIC   tra_ldf_iso   ! routine called by step.F90
 
    !! * Substitutions
@@ -39 +42 @@
    !!----------------------------------------------------------------------
    !!   OPA 9.0 , LOCEAN-IPSL (2005) 
+   !! $Header$ 
+   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
 
@@ -81 +86 @@
       !!         ta = ta + difft
       !!
-      !! ** Action :
-      !!         Update (ta,sa) arrays with the before rotated diffusion trend
-      !!            (except the dk[ dk[.] ] term)
-      !!
-      !! History :
-      !!        !  94-08  (G. Madec, M. Imbard)
-      !!        !  97-05  (G. Madec)  split into traldf and trazdf
-      !!   8.5  !  02-08  (G. Madec)  Free form, F90
-      !!   9.0  !  04-08  (C. Talandier) New trends organization
-      !!        !  05-11  (G. Madec)  merge traldf and trazdf :-)
-      !!----------------------------------------------------------------------
-      !! * Modules used
-      USE oce           , zftv => ua,   &  ! use ua as workspace
-         &                zfsv => va       ! use va as workspace
-
-      !! * Arguments
-      INTEGER, INTENT( in ) ::   kt        ! ocean time-step index
-
-      !! * Local declarations
-      INTEGER ::   ji, jj, jk              ! dummy loop indices
-      INTEGER ::   iku, ikv                ! temporary integer
-      REAL(wp) ::   &
-         zabe1, zabe2, zcof1, zcof2,    &  ! temporary scalars
-         zmsku, zta,                    &  !    "           "
-         zmskv, zsa, zbtr,              &  !    "           "
-         zcoef0, zcoef3, zcoef4            !    "         "
-      REAL(wp), DIMENSION(jpi,jpj) ::   &  ! temporary workspace
-         zdkt , zdk1t, zftu,            &  !    "           "
-         zdks , zdk1s, zfsu                !    "           "
-      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   & 
-         zdit, zdjt, ztfw,              &  ! temporary workspace
-         zdis, zdjs, zsfw                  !    "           "
+      !! ** Action :   Update (ta,sa) arrays with the before rotated diffusion
+      !!            trend (except the dk[ dk[.] ] term)
+      !!----------------------------------------------------------------------
+      USE oce           , zftv => ua   ! use ua as workspace
+      USE oce           , zfsv => va   ! use va as workspace
+      !!
+      INTEGER, INTENT( in ) ::   kt    ! ocean time-step index
+      !!
+      INTEGER  ::   ji, jj, jk   ! dummy loop indices
+      INTEGER  ::   iku, ikv     ! temporary integer
+      REAL(wp) ::   zmsku, zabe1, zcof1, zcoef3, zta   ! temporary scalars
+      REAL(wp) ::   zmskv, zabe2, zcof2, zcoef4, zsa   !    "         "
+      REAL(wp) ::   zcoef0, zbtr                       !    "         "
+      REAL(wp), DIMENSION(jpi,jpj)     ::   zdkt , zdk1t, zftu   ! 2D workspace
+      REAL(wp), DIMENSION(jpi,jpj)     ::   zdks , zdk1s, zfsu   !    "           "
+      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zdit, zdjt, ztfw     ! 3D workspace
+      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zdis, zdjs, zsfw     !  "      "
       !!----------------------------------------------------------------------
 
@@ -130 +120 @@
       zdjs (1,:,:) = 0.e0     ;     zdjs (jpi,:,:) = 0.e0
 !!end
-
 
       ! Horizontal temperature and salinity gradient 
@@ -217 +206 @@
          ! II.4 Second derivative (divergence) and add to the general trend
          ! ----------------------------------------------------------------
-
          DO jj = 2 , jpjm1
             DO ji = fs_2, fs_jpim1   ! vector opt.
@@ -231 +219 @@
       !                                             ! ===============
 
-      IF( ln_diaptr .AND. ( MOD( kt, nf_ptr ) == 0 ) ) THEN
-         ! "zonal" mean lateral diffusive heat and salt transports
+      IF( ln_diaptr .AND. ( MOD( kt, nf_ptr ) == 0 ) ) THEN   ! Poleward diffusive heat and salt transports
          pht_ldf(:) = ptr_vj( zftv(:,:,:) )
          pst_ldf(:) = ptr_vj( zfsv(:,:,:) )
@@ -297 +284 @@
          END DO
       END DO
-
+      !
    END SUBROUTINE tra_ldf_iso
 
@@ -306 +293 @@
 CONTAINS
    SUBROUTINE tra_ldf_iso( kt )               ! Empty routine
-      INTEGER :: kt
       WRITE(*,*) 'tra_ldf_iso: You should not have seen this print! error?', kt
    END SUBROUTINE tra_ldf_iso