Changeset 2528 for trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/trcsms_pisces.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/PISCES/trcsms_pisces.F90 (modified) (6 diffs, 1 prop)

trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/trcsms_pisces.F90

@@ -22 +22 @@
    USE p4zche          ! 
    USE p4zbio          ! 
+   USE p4zsink         ! 
+   USE p4zopt          ! 
+   USE p4zlim          ! 
+   USE p4zprod         !
+   USE p4zmort         !
+   USE p4zmicro        ! 
+   USE p4zmeso         ! 
+   USE p4zrem          ! 
    USE p4zsed          ! 
    USE p4zlys          ! 
    USE p4zflx          ! 
 
-   USE trdmld_trc_oce
-   USE trdmld_trc
+   USE prtctl_trc
+
+   USE trdmod_oce
+   USE trdmod_trc
 
    USE sedmodel
@@ -37 +47 @@
 
    !!----------------------------------------------------------------------
-   !! 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)
    !!----------------------------------------------------------------------
 
@@ -59 +69 @@
       INTEGER ::   jnt, jn
       REAL(wp), DIMENSION(jpi,jpj,jpk) ::   ztrpis   ! used for pisces sms trends
+      CHARACTER (len=25) :: charout
       !!---------------------------------------------------------------------
 
-      IF( kt == nittrc000  .AND. .NOT. ln_rsttr )   CALL trc_sms_pisces_init    ! Initialization (first time-step only)
+      IF( kt == nit000 )   CALL trc_sms_pisces_init    ! Initialization (first time-step only)
 
-      IF( ndayflxtr /= nday ) THEN      ! New days
+      IF( ndayflxtr /= nday_year ) THEN      ! New days
          !
-         ndayflxtr = nday
+         ndayflxtr = nday_year
+
+         IF(lwp) write(numout,*)
+         IF(lwp) write(numout,*) ' New chemical constants and various rates for biogeochemistry at new day : ', nday_year
+         IF(lwp) write(numout,*) '~~~~~~'
 
          CALL p4z_che          ! computation of chemical constants
@@ -71 +86 @@
          !
       ENDIF
-
 
       DO jnt = 1, nrdttrc          ! Potential time splitting if requested
@@ -91 +105 @@
       END DO
 
+
       IF( l_trdtrc ) THEN
           DO jn = jp_pcs0, jp_pcs1
             ztrpis(:,:,:) = tra(:,:,:,jn)
-            CALL trd_mod_trc( ztrpis, jn, jptrc_trd_sms, kt )   ! save trends
+            CALL trd_mod_trc( ztrpis, jn, jptra_trd_sms, kt )   ! save trends
           END DO
       END IF
@@ -121 +136 @@
       REAL(wp) ::  ztmas, ztmas1
 
-      ! Initialization of chemical variables of the carbon cycle
-      ! --------------------------------------------------------
-      DO jk = 1, jpk
-         DO jj = 1, jpj
-            DO ji = 1, jpi
-               ztmas   = tmask(ji,jj,jk)
-               ztmas1  = 1. - tmask(ji,jj,jk)
-               zcaralk = trn(ji,jj,jk,jptal) - borat(ji,jj,jk) / (  1. + 1.E-8 / ( rtrn + akb3(ji,jj,jk) )  )
-               zco3    = ( zcaralk - trn(ji,jj,jk,jpdic) ) * ztmas + 0.5e-3 * ztmas1
-               zbicarb = ( 2. * trn(ji,jj,jk,jpdic) - zcaralk )
-               hi(ji,jj,jk) = ( ak23(ji,jj,jk) * zbicarb / zco3 ) * ztmas + 1.e-9 * ztmas1
+      IF( .NOT. ln_rsttr ) THEN
+         ! Initialization of chemical variables of the carbon cycle
+         ! --------------------------------------------------------
+         DO jk = 1, jpk
+            DO jj = 1, jpj
+               DO ji = 1, jpi
+                  ztmas   = tmask(ji,jj,jk)
+                  ztmas1  = 1. - tmask(ji,jj,jk)
+                  zcaralk = trn(ji,jj,jk,jptal) - borat(ji,jj,jk) / (  1. + 1.E-8 / ( rtrn + akb3(ji,jj,jk) )  )
+                  zco3    = ( zcaralk - trn(ji,jj,jk,jpdic) ) * ztmas + 0.5e-3 * ztmas1
+                  zbicarb = ( 2. * trn(ji,jj,jk,jpdic) - zcaralk )
+                  hi(ji,jj,jk) = ( ak23(ji,jj,jk) * zbicarb / zco3 ) * ztmas + 1.e-9 * ztmas1
+               END DO
             END DO
          END DO
-      END DO
+         !
+      END IF
+
+      ! Time step duration for biology
+      xstep = rfact2 / rday
+
+      CALL p4z_sink_init      ! vertical flux of particulate organic matter
+      CALL p4z_opt_init       ! Optic: PAR in the water column
+      CALL p4z_lim_init       ! co-limitations by the various nutrients
+      CALL p4z_prod_init      ! phytoplankton growth rate over the global ocean. 
+      CALL p4z_rem_init       ! remineralisation
+      CALL p4z_mort_init      ! phytoplankton mortality
+      CALL p4z_micro_init     ! microzooplankton
+      CALL p4z_meso_init      ! mesozooplankton
+      CALL p4z_sed_init       ! sedimentation
+      CALL p4z_lys_init       ! calcite saturation
+      CALL p4z_flx_init       ! gas exchange
+
+      ndayflxtr = 0
 
    END SUBROUTINE trc_sms_pisces_init