Changeset 71 for trunk/NEMO/OPA_SRC/LDF/ldfdyn.F90

Timestamp:

2004-04-22T14:21:33+02:00 (20 years ago)

Author:

opalod

Message:

CT : BUGFIX045 : Addition of a generic interface ldf_zpf_1d_3d

File:

• trunk/NEMO/OPA_SRC/LDF/ldfdyn.F90 (modified) (2 diffs)

trunk/NEMO/OPA_SRC/LDF/ldfdyn.F90

@@ -28 +28 @@
 
   INTERFACE ldf_zpf
-     MODULE PROCEDURE ldf_zpf_1d, ldf_zpf_3d
+     MODULE PROCEDURE ldf_zpf_1d, ldf_zpf_1d_3d, ldf_zpf_3d
   END INTERFACE
 
@@ -297 +297 @@
 
 
+   SUBROUTINE ldf_zpf_1d_3d( ld_print, pdam, pwam, pbot, pdep, pah )
+      !!----------------------------------------------------------------------
+      !!                  ***  ROUTINE ldf_zpf  ***
+      !!
+      !! ** Purpose :   vertical adimensional profile for eddy coefficient
+      !!
+      !! ** Method  :   1D eddy viscosity coefficients ( depth )
+      !!
+      !!----------------------------------------------------------------------
+      !! * Arguments
+      LOGICAL , INTENT (in   ) :: ld_print   ! If true, output arrays on numout
+      REAL(wp), INTENT (in   ) ::   &
+          pdam,     &  ! depth of the inflection point
+          pwam,     &  ! width of inflection
+          pbot         ! battom value (0<pbot<= 1)
+      REAL(wp), INTENT (in   ), DIMENSION(jpk) ::   &
+          pdep         ! depth of the gridpoint (T, U, V, F)
+      REAL(wp), INTENT (inout), DIMENSION(jpi,jpj,jpk) ::   &
+          pah          ! adimensional vertical profile
+
+      !! * Local variables
+      INTEGER  ::   jk           ! dummy loop indices
+      REAL(wp) ::   zm00, zm01, zmhb, zmhs, zcf  ! temporary scalars
+      !!----------------------------------------------------------------------
+
+      zm00 = TANH( ( pdam - gdept(1    ) ) / pwam )
+      zm01 = TANH( ( pdam - gdept(jpkm1) ) / pwam )
+      zmhs = zm00 / zm01
+      zmhb = ( 1.e0 - pbot ) / ( 1.e0 - zmhs ) / zm01
+
+      DO jk = 1, jpk
+         zcf = 1.e0 + zmhb * ( zm00 - TANH( ( pdam - pdep(jk) ) / pwam )  )
+         pah(:,:,jk) = zcf
+      END DO
+
+      ! Control print
+      IF(lwp .AND. ld_print ) THEN
+         WRITE(numout,*)
+         WRITE(numout,*) '         ahm profile : '
+         WRITE(numout,*)
+         WRITE(numout,'("  jk      ahm       ","  depth t-level " )')
+         DO jk = 1, jpk
+            WRITE(numout,'(i6,2f12.4,3x,2f12.4)') jk, pah(:,:,jk), pdep(jk)
+         END DO
+      ENDIF
+
+   END SUBROUTINE ldf_zpf_1d_3d
+
+
    SUBROUTINE ldf_zpf_3d( ld_print, pdam, pwam, pbot, pdep, pah )
       !!----------------------------------------------------------------------