Changeset 2528 for trunk/NEMOGCM/NEMO/TOP_SRC/LOBSTER/trcsed.F90

Timestamp:

2010-12-27T18:33:53+01:00 (13 years ago)

Author:

rblod

Message:

Update NEMOGCM from branch nemo_v3_3_beta

File:

• trunk/NEMOGCM/NEMO/TOP_SRC/LOBSTER/trcsed.F90 (modified) (9 diffs)

trunk/NEMOGCM/NEMO/TOP_SRC/LOBSTER/trcsed.F90

@@ -18 +18 @@
    USE sms_lobster
    USE lbclnk
-   USE trdmld_trc
-   USE trdmld_trc_oce
+   USE trdmod_oce
+   USE trdmod_trc
    USE iom
    USE prtctl_trc      ! Print control for debbuging
@@ -31 +31 @@
 #  include "top_substitute.h90"
    !!----------------------------------------------------------------------
-   !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007) 
+   !! NEMO/TOP 3.3 , NEMO Consortium (2010)
    !! $Id$ 
-   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
+   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
 
@@ -53 +53 @@
       !!                             tra = tra + dz(trn wn)
       !!        
-      !!              IF 'key_trc_diabio' is defined, the now vertical advection
+      !!              IF 'key_diabio' is defined, the now vertical advection
       !!              trend of passive tracers is saved for futher diagnostics.
       !!---------------------------------------------------------------------
@@ -61 +61 @@
       REAL(wp) ::   ztra
       REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zwork
-#if defined key_trc_diaadd && defined key_iomput
+#if defined key_diatrc && defined key_iomput
       REAL(wp), DIMENSION(jpi,jpj) ::  zw2d
 #endif
@@ -77 +77 @@
       ! --------------------------------------------
 
-      ! for detritus sedimentation only - jpdet
+      ! for detritus sedimentation only - jp_lob_det
       zwork(:,:,1  ) = 0.e0      ! surface value set to zero
       zwork(:,:,jpk) = 0.e0      ! bottom value  set to zero
 
-#if defined key_trc_diaadd && defined key_iomput
+#if defined key_diatrc && defined key_iomput
       zw2d(:,:) = 0.
 # endif
@@ -87 +87 @@
       IF( l_trdtrc )THEN
          ALLOCATE( ztrbio(jpi,jpj,jpk) )
-         ztrbio(:,:,:) = tra(:,:,:,jpdet)
+         ztrbio(:,:,:) = tra(:,:,:,jp_lob_det)
       ENDIF
 
       ! tracer flux at w-point: we use -vsed (downward flux)  with simplification : no e1*e2
       DO jk = 2, jpkm1
-         zwork(:,:,jk) = -vsed * trn(:,:,jk-1,jpdet)
+         zwork(:,:,jk) = -vsed * trn(:,:,jk-1,jp_lob_det)
       END DO
 
@@ -100 +100 @@
             DO ji = 1,jpi
                ztra  = - ( zwork(ji,jj,jk) - zwork(ji,jj,jk+1) ) / fse3t(ji,jj,jk)
-               tra(ji,jj,jk,jpdet) = tra(ji,jj,jk,jpdet) + ztra
-#if defined key_trc_diabio
+               tra(ji,jj,jk,jp_lob_det) = tra(ji,jj,jk,jp_lob_det) + ztra
+#if defined key_diabio
                trbio(ji,jj,jk,jp_lob0_trd + 7) = ztra
 #endif
-#if defined key_trc_diaadd
+#if defined key_diatrc
 # if ! defined key_iomput
                trc2d(ji,jj,jp_lob0_2d + 7) = trc2d(ji,jj,jp_lob0_2d + 7) + ztra * fse3t(ji,jj,jk) * 86400.
@@ -115 +115 @@
       END DO
 
-#if defined key_trc_diabio
+#if defined key_diabio
       jl = jp_lob0_trd + 7
       CALL lbc_lnk (trbio(:,:,1,jl), 'T', 1. )    ! Lateral boundary conditions on trcbio
 #endif
-#if defined key_trc_diaadd
+#if defined key_diatrc
 # if ! defined key_iomput
       jl = jp_lob0_2d + 7
@@ -131 +131 @@
 
       IF( l_trdtrc ) THEN
-         ztrbio(:,:,:) = tra(:,:,:,jpdet) - ztrbio(:,:,:)
+         ztrbio(:,:,:) = tra(:,:,:,jp_lob_det) - ztrbio(:,:,:)
          jl = jp_lob0_trd + 7
          CALL trd_mod_trc( ztrbio, jl, kt )   ! handle the trend