New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
trcsms_pisces.F90 in tags/nemo_v3_2/nemo_v3_2/NEMO/TOP_SRC/PISCES – NEMO

source: tags/nemo_v3_2/nemo_v3_2/NEMO/TOP_SRC/PISCES/trcsms_pisces.F90 @ 1878

Last change on this file since 1878 was 1878, checked in by flavoni, 14 years ago

initial test for nemogcm

File size: 5.5 KB
Line 
1MODULE trcsms_pisces
2   !!======================================================================
3   !!                         ***  MODULE trcsms_pisces  ***
4   !! TOP :   PISCES Source Minus Sink manager
5   !!======================================================================
6   !! History :   1.0  !  2004-03 (O. Aumont) Original code
7   !!             2.0  !  2007-12  (C. Ethe, G. Madec)  F90
8   !!----------------------------------------------------------------------
9#if defined key_pisces
10   !!----------------------------------------------------------------------
11   !!   'key_pisces'                                       PISCES bio-model
12   !!----------------------------------------------------------------------
13   !!   trcsms_pisces        :  Time loop of passive tracers sms
14   !!----------------------------------------------------------------------
15   USE oce_trc         !
16   USE trc
17   USE sms_pisces
18   USE lbclnk
19   USE lib_mpp
20   
21   USE p4zint          !
22   USE p4zche          !
23   USE p4zbio          !
24   USE p4zsed          !
25   USE p4zlys          !
26   USE p4zflx          !
27
28   USE trdmld_trc_oce
29   USE trdmld_trc
30
31   USE sedmodel
32
33   IMPLICIT NONE
34   PRIVATE
35
36   PUBLIC   trc_sms_pisces    ! called in trcsms.F90
37
38   !!----------------------------------------------------------------------
39   !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007)
40   !! $Id: trcsms_pisces.F90 1753 2009-11-25 12:35:09Z cetlod $
41   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
42   !!----------------------------------------------------------------------
43
44CONTAINS
45
46   SUBROUTINE trc_sms_pisces( kt )
47      !!---------------------------------------------------------------------
48      !!                     ***  ROUTINE trc_sms_pisces  ***
49      !!
50      !! ** Purpose :   Managment of the call to Biological sources and sinks
51      !!              routines of PISCES bio-model
52      !!
53      !! ** Method  : - at each new day ...
54      !!              - several calls of bio and sed ???
55      !!              - ...
56      !!---------------------------------------------------------------------
57      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index     
58      !!
59      INTEGER ::   jnt, jn
60      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   ztrpis   ! used for pisces sms trends
61      !!---------------------------------------------------------------------
62
63      IF( kt == nittrc000  .AND. .NOT. ln_rsttr )   CALL trc_sms_pisces_init    ! Initialization (first time-step only)
64
65      IF( ndayflxtr /= nday ) THEN      ! New days
66         !
67         ndayflxtr = nday
68
69         CALL p4z_che          ! computation of chemical constants
70         CALL p4z_int          ! computation of various rates for biogeochemistry
71         !
72      ENDIF
73
74
75      DO jnt = 1, nrdttrc          ! Potential time splitting if requested
76         !
77         CALL p4z_bio (kt, jnt)    ! Compute soft tissue production (POC)
78         CALL p4z_sed (kt, jnt)    ! compute soft tissue remineralisation
79         !
80         trb(:,:,:,:) = trn(:,:,:,:)
81         !
82      END DO
83
84      CALL p4z_lys( kt )             ! Compute CaCO3 saturation
85      CALL p4z_flx( kt )             ! Compute surface fluxes
86
87      DO jn = jp_pcs0, jp_pcs1
88        CALL lbc_lnk( trn(:,:,:,jn), 'T', 1. )
89        CALL lbc_lnk( trb(:,:,:,jn), 'T', 1. )
90        CALL lbc_lnk( tra(:,:,:,jn), 'T', 1. )
91      END DO
92
93      IF( l_trdtrc ) THEN
94          DO jn = jp_pcs0, jp_pcs1
95            ztrpis(:,:,:) = tra(:,:,:,jn)
96            CALL trd_mod_trc( ztrpis, jn, jptrc_trd_sms, kt )   ! save trends
97          END DO
98      END IF
99
100      IF( lk_sed ) THEN 
101         !
102         CALL sed_model( kt )     !  Main program of Sediment model
103         !
104         DO jn = jp_pcs0, jp_pcs1
105           CALL lbc_lnk( trn(:,:,:,jn), 'T', 1. )
106         END DO
107         !
108      ENDIF
109
110   END SUBROUTINE trc_sms_pisces
111
112   SUBROUTINE trc_sms_pisces_init
113      !!----------------------------------------------------------------------
114      !!                  ***  ROUTINE trc_sms_pisces_init  ***
115      !!
116      !! ** Purpose :   Initialization of PH variable
117      !!
118      !!----------------------------------------------------------------------
119      INTEGER  ::  ji, jj, jk
120      REAL(wp) ::  zcaralk, zbicarb, zco3
121      REAL(wp) ::  ztmas, ztmas1
122
123      ! Initialization of chemical variables of the carbon cycle
124      ! --------------------------------------------------------
125      DO jk = 1, jpk
126         DO jj = 1, jpj
127            DO ji = 1, jpi
128               ztmas   = tmask(ji,jj,jk)
129               ztmas1  = 1. - tmask(ji,jj,jk)
130               zcaralk = trn(ji,jj,jk,jptal) - borat(ji,jj,jk) / (  1. + 1.E-8 / ( rtrn + akb3(ji,jj,jk) )  )
131               zco3    = ( zcaralk - trn(ji,jj,jk,jpdic) ) * ztmas + 0.5e-3 * ztmas1
132               zbicarb = ( 2. * trn(ji,jj,jk,jpdic) - zcaralk )
133               hi(ji,jj,jk) = ( ak23(ji,jj,jk) * zbicarb / zco3 ) * ztmas + 1.e-9 * ztmas1
134            END DO
135         END DO
136      END DO
137
138   END SUBROUTINE trc_sms_pisces_init
139
140#else
141   !!======================================================================
142   !!  Dummy module :                                   No PISCES bio-model
143   !!======================================================================
144CONTAINS
145   SUBROUTINE trc_sms_pisces( kt )                   ! Empty routine
146      INTEGER, INTENT( in ) ::   kt
147      WRITE(*,*) 'trc_sms_pisces: You should not have seen this print! error?', kt
148   END SUBROUTINE trc_sms_pisces
149#endif 
150
151   !!======================================================================
152END MODULE trcsms_pisces 
Note: See TracBrowser for help on using the repository browser.