Changeset 3446 for branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmicro.F90

Timestamp:

2012-08-10T13:13:55+02:00 (12 years ago)

Author:

cetlod

Message:

branch:2012/dev_r3438_LOCEAN15_PISLOB : 2nd step new PISCES updates from Olivier, see ticket #972

File:

• branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmicro.F90 (modified) (10 diffs)

branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmicro.F90

@@ -22 +22 @@
    USE p4zint          !  interpolation and computation of various fields
    USE p4zprod         !  production
+   USE iom             !  I/O manager
    USE prtctl_trc      !  print control for debugging
 
@@ -29 +30 @@
    PUBLIC   p4z_micro         ! called in p4zbio.F90
    PUBLIC   p4z_micro_init    ! called in trcsms_pisces.F90
-   PUBLIC   p4z_micro_alloc    ! called in trcsms_pisces.F90
 
    !! * Shared module variables
@@ -59 +59 @@
 CONTAINS
 
-   SUBROUTINE p4z_micro( kt )
+   SUBROUTINE p4z_micro( kt, jnt )
       !!---------------------------------------------------------------------
       !!                     ***  ROUTINE p4z_micro  ***
@@ -67 +67 @@
       !! ** Method  : - ???
       !!---------------------------------------------------------------------
-      INTEGER, INTENT(in) ::   kt ! ocean time step
+      INTEGER, INTENT(in) ::  kt  ! ocean time step
+      INTEGER, INTENT(in) ::  jnt 
+      !
       INTEGER  :: ji, jj, jk
       REAL(wp) :: zcompadi, zcompaz , zcompaph, zcompapoc
@@ -77 +79 @@
       REAL(wp) :: zgrazp, zgrazm, zgrazsd
       REAL(wp) :: zgrazmf, zgrazsf, zgrazpf
+      REAL(wp) :: zrfact2
+      REAL(wp), POINTER, DIMENSION(:,:,:) :: zgrazing
       CHARACTER (len=25) :: charout
       !!---------------------------------------------------------------------
@@ -82 +86 @@
       IF( nn_timing == 1 )  CALL timing_start('p4z_micro')
       !
-      grazing(:,:,:) = 0.  !: grazing set to zero
+      IF( ln_diatrc .AND. lk_iomput )  CALL wrk_alloc( jpi, jpj, jpk, zgrazing )
+      !
       DO jk = 1, jpkm1
          DO jj = 1, jpj
@@ -127 +132 @@
 
                ! Grazing by microzooplankton
-               grazing(ji,jj,jk) = grazing(ji,jj,jk) + zgraztot
+               IF( ln_diatrc .AND. lk_iomput )  zgrazing(ji,jj,jk) = zgraztot
 
                !    Various remineralization and excretion terms
@@ -185 +190 @@
       END DO
       !
+      IF( ln_diatrc .AND. lk_iomput .AND. jnt == nrdttrc ) THEN
+         zrfact2 = 1.e3 * rfact2r
+         CALL iom_put( "GRAZ1" , zgrazing(:,:,:) * zrfact2 * tmask(:,:,:) )  ! Total grazing of phyto by zooplankton
+      ENDIF
+      !
       IF(ln_ctl)   THEN  ! print mean trends (used for debugging)
          WRITE(charout, FMT="('micro')")
@@ -190 +200 @@
          CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm)
       ENDIF
+      !
+      IF( ln_diatrc .AND. lk_iomput )  CALL wrk_dealloc( jpi, jpj, jpk, zgrazing )
       !
       IF( nn_timing == 1 )  CALL timing_stop('p4z_micro')
@@ -240 +252 @@
    END SUBROUTINE p4z_micro_init
 
-   INTEGER FUNCTION p4z_micro_alloc()
-      !!----------------------------------------------------------------------
-      !!                     ***  ROUTINE p4z_micro_alloc  ***
-      !!----------------------------------------------------------------------
-      ALLOCATE( grazing(jpi,jpj,jpk), STAT=p4z_micro_alloc )
-      IF( p4z_micro_alloc /= 0 ) CALL ctl_warn('p4z_micro_alloc : failed to allocate arrays.')
-
-   END FUNCTION p4z_micro_alloc
-
 #else
    !!======================================================================