Changeset 2819 for branches/2011/dev_r2787_LOCEAN3_TRA_TRP/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zmort.F90

Timestamp:

2011-08-09T10:29:53+02:00 (13 years ago)

Author:

cetlod

Message:

Improvment of branch dev_r2787_LOCEAN3_TRA_TRP

File:

• branches/2011/dev_r2787_LOCEAN3_TRA_TRP/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zmort.F90 (modified) (7 diffs)

branches/2011/dev_r2787_LOCEAN3_TRA_TRP/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zmort.F90

@@ -14 +14 @@
    !!   p4z_mort_init  :   Initialize the mortality params for phytoplankton
    !!----------------------------------------------------------------------
-   USE trc
-   USE oce_trc         !
-   USE trc         ! 
-   USE sms_pisces      ! 
-   USE p4zsink
-   USE prtctl_trc
+   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 p4zsink         !  vertical flux of particulate matter due to sinking
+   USE prtctl_trc      !  print control for debugging
 
    IMPLICIT NONE
@@ -26 +25 @@
    PUBLIC   p4z_mort    
    PUBLIC   p4z_mort_init    
-
+   PUBLIC   p4z_mort_alloc    
 
    !! * Shared module variables
@@ -81 +80 @@
       !!---------------------------------------------------------------------
 
-
-#if defined key_diatrc
-     prodcal(:,:,:) = 0.  !: Initialisation of calcite production variable
-#endif
-
+      prodcal(:,:,:) = 0.  !: calcite production variable set to zero
       DO jk = 1, jpkm1
          DO jj = 1, jpj
             DO ji = 1, jpi
-
                zcompaph = MAX( ( trn(ji,jj,jk,jpphy) - 1e-8 ), 0.e0 )
-
+               zstep   = xstep
 # if defined key_degrad
-               zstep =  xstep * facvol(ji,jj,jk)  
-# else
-               zstep =  xstep  
+               zstep = zstep * facvol(ji,jj,jk)
 # endif
                !     Squared mortality of Phyto similar to a sedimentation term during
@@ -117 +109 @@
                tra(ji,jj,jk,jpnfe) = tra(ji,jj,jk,jpnfe) - zmortp * zfactfe
                zprcaca = xfracal(ji,jj,jk) * zmortp
-#if defined key_diatrc
+
+               ! calcite production
                prodcal(ji,jj,jk) = prodcal(ji,jj,jk) + zprcaca  ! prodcal=prodcal(nanophy)+prodcal(microzoo)+prodcal(mesozoo)
-#endif
+               !
                zfracal = 0.5 * xfracal(ji,jj,jk)
                tra(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic) - zprcaca
@@ -177 +170 @@
                !    sticky and coagulate to sink quickly out of the euphotic zone
                !     ------------------------------------------------------------
-
+               zstep   = xstep
 # if defined key_degrad
-               zstep =  xstep * facvol(ji,jj,jk)  
-# else
-               zstep =  xstep  
+               zstep = zstep * facvol(ji,jj,jk)
 # endif
                !  Phytoplankton respiration 
@@ -243 +234 @@
       NAMELIST/nampismort/ wchl, wchld, mprat, mprat2, mpratm
 
-      REWIND( numnat )                     ! read numnat
-      READ  ( numnat, nampismort )
+      REWIND( numnatp )                     ! read numnatp
+      READ  ( numnatp, nampismort )
 
       IF(lwp) THEN                         ! control print
@@ -259 +250 @@
    END SUBROUTINE p4z_mort_init
 
+   INTEGER FUNCTION p4z_mort_alloc()
+      !!----------------------------------------------------------------------
+      !!                     ***  ROUTINE p4z_mort_alloc  ***
+      !!----------------------------------------------------------------------
+      ALLOCATE( prodcal(jpi,jpj,jpk), STAT=p4z_mort_alloc )
+      IF( p4z_mort_alloc /= 0 ) CALL ctl_warn('p4z_mort_alloc : failed to allocate arrays.')
+
+   END FUNCTION p4z_mort_alloc
+
 #else
    !!======================================================================