Changeset 3192 for branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/PISCES/trcsms_pisces.F90

Timestamp:

2011-12-05T14:55:12+01:00 (12 years ago)

Author:

cetlod

Message:

Perform the initialisation phase of PISCES at nit000 rather than nittrc000

File:

• branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/PISCES/trcsms_pisces.F90 (modified) (4 diffs)

branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/PISCES/trcsms_pisces.F90

@@ -18 +18 @@
    USE p4zbio          !  Biological model
    USE p4zche          !  Chemical model
-   USE p4zsink         !  vertical flux of particulate matter due to sinking
-   USE p4zopt          !  optical model
-   USE p4zlim          !  Co-limitations of differents nutrients
-   USE p4zprod         !  Growth rate of the 2 phyto groups
-   USE p4zmort         !  Mortality terms for phytoplankton
-   USE p4zmicro        !  Sources and sinks of microzooplankton
-   USE p4zmeso         !  Sources and sinks of mesozooplankton
-   USE p4zrem          !  Remineralisation of organic matter
    USE p4zlys          !  Calcite saturation
    USE p4zflx          !  Gas exchange
@@ -74 +66 @@
       IF( nn_timing == 1 )  CALL timing_start('trc_sms_pisces')
       !
-      IF( kt == nittrc000 )                                        CALL trc_sms_pisces_init       ! Initialization (first time-step only)
       IF( ln_pisdmp .AND. MOD( kt - nn_dttrc, nn_pisdmp ) == 0 )   CALL trc_sms_pisces_dmp( kt )  ! Relaxation of some tracers
-                                                                   CALL trc_sms_pisces_mass_conserv( kt )
+                                                                   CALL trc_sms_pisces_mass_conserv( kt ) ! Mass conservation checking
 
       IF( ndayflxtr /= nday_year ) THEN      ! New days
@@ -176 +167 @@
 
    END SUBROUTINE trc_sms_pisces_dmp
-
-   SUBROUTINE trc_sms_pisces_init
-      !!----------------------------------------------------------------------
-      !!                  ***  ROUTINE trc_sms_pisces_init  ***
-      !!
-      !! ** Purpose :   Initialization of PH variable
-      !!
-      !!----------------------------------------------------------------------
-      INTEGER  ::  ji, jj, jk
-      REAL(wp) ::  zcaralk, zbicarb, zco3
-      REAL(wp) ::  ztmas, 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 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
 
    SUBROUTINE trc_sms_pisces_mass_conserv ( kt )
@@ -277 +220 @@
          ENDIF
        ENDIF
- 9500  FORMAT(i6,e18.10)     
+ 9500  FORMAT(i10,e18.10)     
        !
    END SUBROUTINE trc_sms_pisces_mass_conserv