source: branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zsed.F90 @ 7041

Last change on this file since 7041 was 7041, checked in by cetlod, 4 years ago

ROBUST5_CNRS : implementation of part I of new TOP interface - 1st step -, see ticket #1782

  • Property svn:keywords set to Id
File size: 7.3 KB
Line 
1MODULE p2zsed
2   !!======================================================================
3   !!                         ***  MODULE p2zsed  ***
4   !! TOP :   PISCES Compute loss of organic matter in the sediments
5   !!======================================================================
6   !! History :    -   !  1995-06 (M. Levy)  original code
7   !!              -   !  2000-12 (E. Kestenare)  clean up
8   !!             2.0  !  2007-12  (C. Deltel, G. Madec)  F90 + simplifications
9   !!----------------------------------------------------------------------
10#if defined key_pisces
11   !!----------------------------------------------------------------------
12   !!   'key_pisces'                                     LOBSTER bio-model
13   !!----------------------------------------------------------------------
14   !!   p2z_sed        :  Compute loss of organic matter in the sediments
15   !!----------------------------------------------------------------------
16   USE oce_trc         !
17   USE trc
18   USE sms_pisces
19   USE lbclnk
20   USE trd_oce
21   USE trdtrc
22   USE iom
23   USE prtctl_trc      ! Print control for debbuging
24
25   IMPLICIT NONE
26   PRIVATE
27
28   PUBLIC   p2z_sed         ! called in ???
29   PUBLIC   p2z_sed_init    ! called in ???
30
31   REAL(wp), PUBLIC ::   sedlam      !: time coefficient of POC remineralization in sediments
32   REAL(wp), PUBLIC ::   sedlostpoc  ! mass of POC lost in sediments
33   REAL(wp), PUBLIC ::   vsed        ! detritus sedimentation speed [m/s]
34   REAL(wp), PUBLIC ::   xhr         ! coeff for martin''s remineralisation profile
35
36   !!----------------------------------------------------------------------
37   !! NEMO/TOP 3.3 , NEMO Consortium (2010)
38   !! $Id$
39   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
40   !!----------------------------------------------------------------------
41
42CONTAINS
43
44   SUBROUTINE p2z_sed( kt )
45      !!---------------------------------------------------------------------
46      !!                     ***  ROUTINE p2z_sed  ***
47      !!
48      !! ** Purpose :   compute the now trend due to the vertical sedimentation of
49      !!              detritus and add it to the general trend of detritus equations
50      !!
51      !! ** Method  :   this ROUTINE compute not exactly the advection but the
52      !!              transport term, i.e.  dz(wt) and dz(ws)., dz(wtr)
53      !!              using an upstream scheme
54      !!              the now vertical advection of tracers is given by:
55      !!                      dz(trn wn) = 1/bt dk+1( e1t e2t vsed (trn) )
56      !!              add this trend now to the general trend of tracer (ta,sa,tra):
57      !!                             tra = tra + dz(trn wn)
58      !!       
59      !!              IF 'key_diabio' is defined, the now vertical advection
60      !!              trend of passive tracers is saved for futher diagnostics.
61      !!---------------------------------------------------------------------
62      !!
63      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index     
64      !!
65      INTEGER  ::   ji, jj, jk, jl, ierr
66      CHARACTER (len=25) :: charout
67      REAL(wp), POINTER, DIMENSION(:,:  ) :: zw2d
68      REAL(wp), POINTER, DIMENSION(:,:,:) :: zwork, ztra
69      !!---------------------------------------------------------------------
70      !
71      IF( nn_timing == 1 )  CALL timing_start('p2z_sed')
72      !
73      IF( kt == nittrc000 ) THEN
74         IF(lwp) WRITE(numout,*)
75         IF(lwp) WRITE(numout,*) ' p2z_sed: LOBSTER sedimentation'
76         IF(lwp) WRITE(numout,*) ' ~~~~~~~'
77      ENDIF
78
79      ! Allocate temporary workspace
80      CALL wrk_alloc( jpi, jpj, jpk, zwork, ztra )
81
82      ! sedimentation of detritus  : upstream scheme
83      ! --------------------------------------------
84
85      ! for detritus sedimentation only - jpdet
86      zwork(:,:,1  ) = 0.e0      ! surface value set to zero
87      zwork(:,:,jpk) = 0.e0      ! bottom value  set to zero
88
89      ! tracer flux at w-point: we use -vsed (downward flux)  with simplification : no e1*e2
90      DO jk = 2, jpkm1
91         zwork(:,:,jk) = -vsed * trn(:,:,jk-1,jpdet)
92      END DO
93
94      ! tracer flux divergence at t-point added to the general trend
95      DO jk = 1, jpkm1
96         DO jj = 1, jpj
97            DO ji = 1, jpi
98               ztra(ji,jj,jk)  = - ( zwork(ji,jj,jk) - zwork(ji,jj,jk+1) ) / e3t_n(ji,jj,jk)
99               tra(ji,jj,jk,jpdet) = tra(ji,jj,jk,jpdet) + ztra(ji,jj,jk) 
100            END DO
101         END DO
102      END DO
103
104      IF( lk_iomput )  THEN
105         IF( iom_use( "TDETSED" ) ) THEN
106            CALL wrk_alloc( jpi, jpj, zw2d )
107            zw2d(:,:) =  ztra(:,:,1) * e3t_n(:,:,1) * 86400._wp
108            DO jk = 2, jpkm1
109               zw2d(:,:) = zw2d(:,:) + ztra(:,:,jk) * e3t_n(:,:,jk) * 86400._wp
110            END DO
111            CALL iom_put( "TDETSED", zw2d )
112            CALL wrk_dealloc( jpi, jpj, zw2d )
113         ENDIF
114      ENDIF
115      !
116      CALL wrk_dealloc( jpi, jpj, jpk, zwork, ztra )
117      !
118
119      IF(ln_ctl)   THEN  ! print mean trends (used for debugging)
120         WRITE(charout, FMT="('sed')")
121         CALL prt_ctl_trc_info(charout)
122         CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm)
123      ENDIF
124      !
125      IF( nn_timing == 1 )  CALL timing_stop('p2z_sed')
126      !
127   END SUBROUTINE p2z_sed
128
129   SUBROUTINE p2z_sed_init
130      !!----------------------------------------------------------------------
131      !!                  ***  ROUTINE p2z_sed_init  ***
132      !!
133      !! ** Purpose :   Parameters from aphotic layers to sediment
134      !!
135      !! ** Method  :   Read the namlobsed namelist and check the parameters
136      !!
137      !!----------------------------------------------------------------------
138      NAMELIST/namlobsed/ sedlam, sedlostpoc, vsed, xhr
139      INTEGER :: ios                 ! Local integer output status for namelist read
140
141      REWIND( numnatp_ref )              ! Namelist namlobsed in reference namelist : Lobster sediments
142      READ  ( numnatp_ref, namlobsed, IOSTAT = ios, ERR = 901)
143901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlosed in reference namelist', lwp )
144
145      REWIND( numnatp_cfg )              ! Namelist namlobsed in configuration namelist : Lobster sediments
146      READ  ( numnatp_cfg, namlobsed, IOSTAT = ios, ERR = 902 )
147902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobsed in configuration namelist', lwp )
148      IF(lwm) WRITE ( numonp, namlobsed )
149
150      IF(lwp) THEN
151          WRITE(numout,*) ' Namelist namlobsed'
152          WRITE(numout,*) '    time coeff of POC in sediments                       sedlam    =', sedlam
153          WRITE(numout,*) '    Sediment geol loss for POC                           sedlostpoc=', sedlostpoc
154          WRITE(numout,*) '    detritus sedimentation speed                         vsed      =', 86400 * vsed  , ' d'
155          WRITE(numout,*) '    coeff for martin''s remineralistion                  xhr       =', xhr
156          WRITE(numout,*) ' '
157      ENDIF
158      !
159   END SUBROUTINE p2z_sed_init
160
161#else
162   !!======================================================================
163   !!  Dummy module :                                   No PISCES bio-model
164   !!======================================================================
165CONTAINS
166   SUBROUTINE p2z_sed( kt )                   ! Empty routine
167      INTEGER, INTENT( in ) ::   kt
168      WRITE(*,*) 'p2z_sed: You should not have seen this print! error?', kt
169   END SUBROUTINE p2z_sed
170#endif 
171
172   !!======================================================================
173END MODULE  p2zsed
Note: See TracBrowser for help on using the repository browser.