Changeset 2443 for branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/TRA/trabbc.F90

Timestamp:

2010-11-29T08:52:36+01:00 (13 years ago)

Author:

gm

Message:

v3.3beta: #766 share the deepest ocean level indices (mbkt, mbku & mbkv)

File:

• branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/TRA/trabbc.F90 (modified) (9 diffs)

branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/TRA/trabbc.F90

@@ -8 +8 @@
    !!             -   ! 2002-11 (A. Bozec)  tra_bbc_init: original code
    !!            3.3  ! 2010-10 (G. Madec)  dynamical allocation + suppression of key_trabbc
+   !!             -   ! 2010-11 (G. Madec)  use mbkt array (deepest ocean t-level)
    !!----------------------------------------------------------------------
 
+   !!----------------------------------------------------------------------
    !!   tra_bbc      : update the tracer trend at ocean bottom 
    !!   tra_bbc_init : initialization of geothermal heat flux trend
    !!----------------------------------------------------------------------
-   USE oce             ! ocean dynamics and active tracers
-   USE dom_oce         ! ocean space and time domain
+   USE oce             ! ocean variables
+   USE dom_oce         ! domain: ocean
    USE phycst          ! physical constants
-   USE trdmod_oce      ! ocean trends 
-   USE trdtra      ! ocean trends 
+   USE trdmod_oce      ! trends: ocean variables 
+   USE trdtra          ! trends: active tracers 
    USE in_out_manager  ! I/O manager
    USE prtctl          ! Print control
@@ -32 +34 @@
    REAL(wp)        ::   rn_geoflx_cst = 86.4e-3_wp   !  Constant value of geothermal heat flux
 
-   INTEGER , DIMENSION(:,:), ALLOCATABLE ::   nbotlevt   ! ocean bottom level index at T-pt
    REAL(wp), PUBLIC, DIMENSION(:,:), ALLOCATABLE ::   qgh_trd0   ! geothermal heating trend
  
@@ -57 +58 @@
       !!       ocean bottom can be computed once and is added to the temperature
       !!       trend juste above the bottom at each time step:
-      !!            ta = ta + Qsf / (rau0 rcp e3T) for k= mbathy -1
+      !!            ta = ta + Qsf / (rau0 rcp e3T) for k= mbkt
       !!       Where Qsf is the geothermal heat flux.
       !!
@@ -66 +67 @@
       !!              Emile-Geay and Madec, 2009, Ocean Science.
       !!----------------------------------------------------------------------
-      !!
       INTEGER, INTENT( in ) ::   kt   ! ocean time-step index
       !!
       INTEGER  ::   ji, jj, ik    ! dummy loop indices
-      REAL(wp) ::   zqgh_trd  ! geothermal heat flux trend
+      REAL(wp) ::   zqgh_trd      ! geothermal heat flux trend
       REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  ztrdt
       !!----------------------------------------------------------------------
-
+      !
       IF( l_trdtra )   THEN         ! Save ta and sa trends
          ALLOCATE( ztrdt(jpi,jpj,jpk) )     ;   ztrdt(:,:,:) = tsa(:,:,:,jp_tem)
       ENDIF
-
+      !
       !                             !  Add the geothermal heat flux trend on temperature
 #if defined key_vectopt_loop
@@ -86 +86 @@
          DO ji = 2, jpim1
 #endif
-            ik = nbotlevt(ji,jj)
+            ik = mbkt(ji,jj)
             zqgh_trd = qgh_trd0(ji,jj) / fse3t(ji,jj,ik)
             tsa(ji,jj,ik,jp_tem) = tsa(ji,jj,ik,jp_tem) + zqgh_trd
          END DO
       END DO
-
+      !
       IF( l_trdtra ) THEN        ! Save the geothermal heat flux trend for diagnostics
          ztrdt(:,:,:) = tsa(:,:,:,jp_tem) - ztrdt(:,:,:)
@@ -117 +117 @@
       !!
       !! ** Action  : - read/fix the geothermal heat qgh_trd0
-      !!              - compute the bottom ocean level nbotlevt
       !!----------------------------------------------------------------------
       USE iom
@@ -127 +126 @@
       !!----------------------------------------------------------------------
 
-      REWIND ( numnam )              ! Read Namelist nambbc : bottom momentum boundary condition
-      READ   ( numnam, nambbc )
+      REWIND( numnam )                 ! Read Namelist nambbc : bottom momentum boundary condition
+      READ  ( numnam, nambbc )
 
       IF(lwp) THEN                     ! Control print
@@ -143 +142 @@
       IF( ln_trabbc ) THEN             !==  geothermal heating  ==!
          !
-         ALLOCATE( nbotlevt(jpi,jpj) )    ! allocation
-         ALLOCATE( qgh_trd0(jpi,jpj) )     
-         !             
-         DO jj = 1, jpj                   ! level of the ocean bottom at T-point
-            DO ji = 1, jpi
-               nbotlevt(ji,jj) = MAX( mbathy(ji,jj)-1, 1 )
-            END DO
-         END DO
+         ALLOCATE( qgh_trd0(jpi,jpj) )    ! allocation
          !
          SELECT CASE ( nn_geoflx )        ! geothermal heat flux / (rauO * Cp)
@@ -172 +164 @@
          !
       ELSE
-            IF(lwp) WRITE(numout,*) '      *** no geothermal heat flux'
+         IF(lwp) WRITE(numout,*) '      *** no geothermal heat flux'
       ENDIF
       !