Changeset 1287 for trunk/NEMO/TOP_SRC/PISCES/trcsms_pisces.F90

Timestamp:

2009-02-03T14:57:43+01:00 (15 years ago)

Author:

cetlod

Message:

Initialization of PH in a better way in PISCEs model, see ticket:316

File:

• trunk/NEMO/TOP_SRC/PISCES/trcsms_pisces.F90 (modified) (3 diffs)

trunk/NEMO/TOP_SRC/PISCES/trcsms_pisces.F90

@@ -61 +61 @@
       !!---------------------------------------------------------------------
 
+      IF( kt == nittrc000  .AND. .NOT. lrsttr )   CALL trc_sms_pisces_init    ! Initialization (first time-step only)
+
       IF( ndayflxtr /= nday ) THEN      ! New days
          !
@@ -83 +85 @@
       CALL p4z_flx( kt )             ! Compute surface fluxes
 
-
       DO jn = jp_pcs0, jp_pcs1
         CALL lbc_lnk( trn(:,:,:,jn), 'T', 1. )
@@ -97 +98 @@
       END IF
 
-#if defined key_sed
+      IF( lk_sed ) THEN 
+         !
+         CALL sed_model( kt )     !  Main program of Sediment model
+         !
+         DO jn = jp_pcs0, jp_pcs1
+           CALL lbc_lnk( trn(:,:,:,jn), 'T', 1. )
+         END DO
+         !
+      ENDIF
 
-      CALL sed_model( kt )     !  Main program of Sediment model
+   END SUBROUTINE trc_sms_pisces
 
-      DO jn = jp_pcs0, jp_pcs1
-        CALL lbc_lnk( trn(:,:,:,jn), 'T', 1. )
+   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
+
+      ! 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
 
-#endif
-
-      !
-   END SUBROUTINE trc_sms_pisces
+   END SUBROUTINE trc_sms_pisces_init
 
 #else