Changeset 2715 for trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_lap.F90

Timestamp:

2011-03-30T17:58:35+02:00 (13 years ago)

Author:

rblod

Message:

First attempt to put dynamic allocation on the trunk

File:

• trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_lap.F90 (modified) (8 diffs)

trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_lap.F90

@@ -4 +4 @@
    !! Ocean dynamics:  lateral viscosity trend
    !!======================================================================
+   !! History :  OPA  ! 1990-09 (G. Madec) Original code
+   !!            4.0  ! 1991-11 (G. Madec)
+   !!            6.0  ! 1996-01 (G. Madec) statement function for e3 and ahm
+   !!   NEMO     1.0  ! 2002-06 (G. Madec)  F90: Free form and module
+   !!             -   ! 2004-08 (C. Talandier) New trends organization
+   !!----------------------------------------------------------------------
 
    !!----------------------------------------------------------------------
@@ -9 +15 @@
    !!                  using an iso-level harmonic operator
    !!----------------------------------------------------------------------
-   !! * Modules used
    USE oce             ! ocean dynamics and tracers
    USE dom_oce         ! ocean space and time domain
@@ -22 +27 @@
    PRIVATE
 
-   !! * Routine accessibility
    PUBLIC dyn_ldf_lap  ! called by step.F90
 
@@ -32 +36 @@
    !! NEMO/OPA 3.3 , NEMO Consortium (2010)
    !! $Id$ 
-   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
+   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
-
 CONTAINS
 
@@ -58 +61 @@
       !! ** Action : - Update (ua,va) with the before iso-level harmonic 
       !!               mixing trend.
-      !!
-      !! History :
-      !!        !  90-09 (G. Madec) Original code
-      !!        !  91-11 (G. Madec)
-      !!        !  96-01 (G. Madec) statement function for e3 and ahm
-      !!   8.5  !  02-06 (G. Madec)  F90: Free form and module
-      !!   9.0  !  04-08 (C. Talandier) New trends organization
       !!----------------------------------------------------------------------
-      !! * Arguments
-      INTEGER, INTENT( in ) ::   kt       ! ocean time-step index
-
-      !! * Local declarations
-      INTEGER  ::   ji, jj, jk            ! dummy loop indices
-      REAL(wp) ::   &
-         zua, zva, ze2u, ze1v             ! temporary scalars
+      INTEGER, INTENT( in ) ::   kt   ! ocean time-step index
+      !
+      INTEGER  ::   ji, jj, jk             ! dummy loop indices
+      REAL(wp) ::   zua, zva, ze2u, ze1v   ! local scalars
       !!----------------------------------------------------------------------
-
+      !
       IF( kt == nit000 ) THEN
          IF(lwp) WRITE(numout,*)
@@ -80 +73 @@
          IF(lwp) WRITE(numout,*) '~~~~~~~ '
       ENDIF
-
       !                                                ! ===============
       DO jk = 1, jpkm1                                 ! Horizontal slab
@@ -86 +78 @@
          DO jj = 2, jpjm1
             DO ji = fs_2, fs_jpim1   ! vector opt.
-               ze2u = rotb (ji,jj,jk)*fsahmf(ji,jj,jk)*fse3f(ji,jj,jk)
-               ze1v = hdivb(ji,jj,jk)*fsahmt(ji,jj,jk)
+               ze2u = rotb (ji,jj,jk) * fsahmf(ji,jj,jk) * fse3f(ji,jj,jk)
+               ze1v = hdivb(ji,jj,jk) * fsahmt(ji,jj,jk)
                ! horizontal diffusive trends
                zua = - ( ze2u - rotb (ji,jj-1,jk)*fsahmf(ji,jj-1,jk)*fse3f(ji,jj-1,jk) ) / ( e2u(ji,jj) * fse3u(ji,jj,jk) )   &
@@ -103 +95 @@
       END DO                                           !   End of slab
       !                                                ! ===============
-
    END SUBROUTINE dyn_ldf_lap