Changeset 3318 for branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/TOP_SRC/PISCES/trcsms_pisces.F90

Timestamp:

2012-02-25T16:50:01+01:00 (12 years ago)

Author:

gm

Message:

Ediag branche: #927 split TRA/DYN trd computation

File:

• branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/TOP_SRC/PISCES/trcsms_pisces.F90 (modified) (11 diffs)

branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/TOP_SRC/PISCES/trcsms_pisces.F90

@@ -4 +4 @@
    !! TOP :   PISCES Source Minus Sink manager
    !!======================================================================
-   !! History :   1.0  !  2004-03 (O. Aumont) Original code
-   !!             2.0  !  2007-12  (C. Ethe, G. Madec)  F90
+   !! History :  1.0  !  2004-03 (O. Aumont) Original code
+   !!            2.0  !  2007-12  (C. Ethe, G. Madec)  F90
    !!----------------------------------------------------------------------
 #if defined key_pisces
@@ -11 +11 @@
    !!   'key_pisces'                                       PISCES bio-model
    !!----------------------------------------------------------------------
-   !!   trcsms_pisces        :  Time loop of passive tracers sms
-   !!----------------------------------------------------------------------
-   USE oce_trc         !  shared variables between ocean and passive tracers
-   USE trc             !  passive tracers common variables 
-   USE sms_pisces      !  PISCES Source Minus Sink variables
-   USE p4zbio          !  Biological model
-   USE p4zche          !  Chemical model
-   USE p4zlys          !  Calcite saturation
-   USE p4zflx          !  Gas exchange
-   USE p4zsed          !  Sedimentation
-   USE p4zint          !  time interpolation
-   USE trdmod_oce      !  Ocean trends variables
-   USE trdmod_trc      !  TOP trends variables
-   USE sedmodel        !  Sediment model
-   USE prtctl_trc      !  print control for debugging
+   !!   trcsms_pisces :  Time loop of passive tracers sms
+   !!----------------------------------------------------------------------
+   USE oce_trc        !  shared variables between ocean and passive tracers
+   USE trc            !  passive tracers common variables 
+   USE sms_pisces     !  PISCES Source Minus Sink variables
+   USE p4zbio         !  Biological model
+   USE p4zche         !  Chemical model
+   USE p4zlys         !  Calcite saturation
+   USE p4zflx         !  Gas exchange
+   USE p4zsed         !  Sedimentation
+   USE p4zint         !  time interpolation
+   USE trd_oce        ! trends: ocean variables
+   USE trdmod_trc     !  TOP trends variables
+   USE sedmodel       !  Sediment model
+   USE prtctl_trc     !  print control for debugging
 
    IMPLICIT NONE
    PRIVATE
 
-   PUBLIC   trc_sms_pisces    ! called in trcsms.F90
-
-   LOGICAL ::  ln_check_mass = .false.       !: Flag to check mass conservation 
-
-   INTEGER ::  numno3  !: logical unit for NO3 budget
-   INTEGER ::  numalk  !: logical unit for talk budget
-   INTEGER ::  numsil  !: logical unit for Si budget
+   PUBLIC   trc_sms_pisces   ! called in trcsms.F90
+
+   LOGICAL ::   ln_check_mass = .false.   ! Flag to check mass conservation 
+
+   INTEGER ::   numno3   ! logical unit for NO3 budget
+   INTEGER ::   numalk   ! logical unit for talk budget
+   INTEGER ::   numsil   ! logical unit for Si budget
 
    !!----------------------------------------------------------------------
@@ -43 +43 @@
    !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
-
 CONTAINS
 
@@ -57 +56 @@
       !!              - ...
       !!---------------------------------------------------------------------
-      !
-      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index      
-      !!
-      INTEGER ::   jnt, jn
+      INTEGER, INTENT( in ) ::   kt   ! ocean time-step index      
+      !
+      INTEGER ::   jnt, jn   ! dummy loop indices
       CHARACTER (len=25) :: charout
       !!---------------------------------------------------------------------
       !
-      IF( nn_timing == 1 )  CALL timing_start('trc_sms_pisces')
+      IF( nn_timing == 1 )   CALL timing_start('trc_sms_pisces')
       !
       IF( ln_pisdmp .AND. MOD( kt - nn_dttrc, nn_pisdmp ) == 0 )   CALL trc_sms_pisces_dmp( kt )  ! Relaxation of some tracers
@@ -102 +100 @@
       IF( l_trdtrc ) THEN
           DO jn = jp_pcs0, jp_pcs1
-            CALL trd_mod_trc( tra(:,:,:,jn), jn, jptra_trd_sms, kt )   ! save trends
+            CALL trd_mod_trc( tra(:,:,:,jn), jn, jptra_sms, kt )   ! save trends
           END DO
       END IF
@@ -120 +118 @@
    END SUBROUTINE trc_sms_pisces
 
+
    SUBROUTINE trc_sms_pisces_dmp( kt )
       !!----------------------------------------------------------------------
@@ -126 +125 @@
       !! ** purpose  : Relaxation of some tracers
       !!----------------------------------------------------------------------
-      !
       INTEGER, INTENT( in )  ::     kt ! time step
       !
-      REAL(wp) ::  alkmean = 2426.     ! mean value of alkalinity ( Glodap ; for Goyet 2391. )
-      REAL(wp) ::  po4mean = 2.165     ! mean value of phosphates
-      REAL(wp) ::  no3mean = 30.90     ! mean value of nitrate
-      REAL(wp) ::  silmean = 91.51     ! mean value of silicate
-      !
-      REAL(wp) :: zarea, zalksum, zpo4sum, zno3sum, zsilsum
-      !!---------------------------------------------------------------------
-
+      REAL(wp) ::   alkmean = 2426.   ! mean value of alkalinity ( Glodap ; for Goyet 2391. )
+      REAL(wp) ::   po4mean = 2.165   ! mean value of phosphates
+      REAL(wp) ::   no3mean = 30.90   ! mean value of nitrate
+      REAL(wp) ::   silmean = 91.51   ! mean value of silicate
+      !
+      REAL(wp) ::   zarea, zalksum, zpo4sum, zno3sum, zsilsum   ! local scalars
+      !!---------------------------------------------------------------------
 
       IF(lwp)  WRITE(numout,*)
@@ -165 +162 @@
          !
       ENDIF
-
+      !
    END SUBROUTINE trc_sms_pisces_dmp
+
 
    SUBROUTINE trc_sms_pisces_mass_conserv ( kt )
@@ -173 +171 @@
       !!
       !! ** Purpose :  Mass conservation check 
-      !!
-      !!---------------------------------------------------------------------
-      !
+      !!---------------------------------------------------------------------
       INTEGER, INTENT( in ) ::   kt      ! ocean time-step index      
-      !!
+      !
       REAL(wp) :: zalkbudget, zno3budget, zsilbudget
       !
       NAMELIST/nampismass/ ln_check_mass
       !!---------------------------------------------------------------------
-
+      !
       IF( kt == nittrc000 ) THEN 
          REWIND( numnatp )       
@@ -199 +195 @@
          ENDIF
       ENDIF
-
+      !
       IF( ln_check_mass ) THEN      !   Compute the budget of NO3, ALK, Si
          zno3budget = glob_sum( (   trn(:,:,:,jpno3) + trn(:,:,:,jpnh4)  &
@@ -206 +202 @@
             &                     + trn(:,:,:,jppoc) + trn(:,:,:,jpgoc)  &
             &                     + trn(:,:,:,jpdoc)                     ) * cvol(:,:,:)  ) 
-         ! 
+            ! 
          zsilbudget = glob_sum( (   trn(:,:,:,jpsil) + trn(:,:,:,jpgsi)  &
             &                     + trn(:,:,:,jpdsi)                     ) * cvol(:,:,:)  )
-         ! 
+            ! 
          zalkbudget = glob_sum( (   trn(:,:,:,jpno3) * rno3              &
             &                     + trn(:,:,:,jptal)                     &
             &                     + trn(:,:,:,jpcal) * 2.                ) * cvol(:,:,:)  )
-
+            !
          IF( lwp ) THEN
             WRITE(numno3,9500) kt,  zno3budget / areatot