source: branches/UKMO/r6232_tracer_advection/NEMOGCM/NEMO/TOP_SRC/PISCES/SED/sedsfc.F90 @ 9295

Last change on this file since 9295 was 9295, checked in by jcastill, 3 years ago

Remove svn keywords

File size: 2.7 KB
Line 
1MODULE sedsfc
2   !!======================================================================
3   !!              ***  MODULE  sedsfc  ***
4   !!    Sediment : Data at sediment surface
5   !!=====================================================================
6#if defined key_sed && ! defined key_sed_off
7   !! * Modules used
8   USE sed     ! sediment global variable
9   USE sedarr
10   USE seddta
11
12   PUBLIC sed_sfc
13
14   !! $Id$
15CONTAINS
16
17   SUBROUTINE sed_sfc( kt )
18      !!---------------------------------------------------------------------
19      !!                  ***  ROUTINE sed_sfc ***
20      !!
21      !! ** Purpose :  Give data from sediment model to tracer model
22      !!
23      !!
24      !!   History :
25      !!        !  06-04 (C. Ethe)  Orginal code
26      !!----------------------------------------------------------------------
27      !!* Arguments
28      INTEGER, INTENT(in) ::  kt              ! time step
29
30      ! * local variables
31      INTEGER :: ji, jj, ikt     ! dummy loop indices
32
33      !------------------------------------------------------------------------
34
35      IF( kt == nitsed000 ) THEN
36         WRITE(numsed,*) ' sed_sfc : Give data from sediment model to tracer model  '
37         WRITE(numsed,*) ' '
38      ENDIF
39
40      ! reading variables
41
42      CALL unpack_arr ( jpoce, trc_data(1:jpi,1:jpj,1), iarroce(1:jpoce), pwcp(1:jpoce,1,jwalk) )
43      CALL unpack_arr ( jpoce, trc_data(1:jpi,1:jpj,2), iarroce(1:jpoce), pwcp(1:jpoce,1,jwdic) )
44      CALL unpack_arr ( jpoce, trc_data(1:jpi,1:jpj,3), iarroce(1:jpoce), pwcp(1:jpoce,1,jwno3) )
45      CALL unpack_arr ( jpoce, trc_data(1:jpi,1:jpj,4), iarroce(1:jpoce), pwcp(1:jpoce,1,jwpo4) )
46      CALL unpack_arr ( jpoce, trc_data(1:jpi,1:jpj,5), iarroce(1:jpoce), pwcp(1:jpoce,1,jwoxy) )
47      CALL unpack_arr ( jpoce, trc_data(1:jpi,1:jpj,6), iarroce(1:jpoce), pwcp(1:jpoce,1,jwsil) )
48
49
50      DO jj = 1,jpj
51         DO ji = 1, jpi
52            ikt = mbkt(ji,jj)
53            IF ( tmask(ji,jj,ikt) == 1 ) THEN
54               trn(ji,jj,ikt,jptal) = trc_data(ji,jj,1)
55               trn(ji,jj,ikt,jpdic) = trc_data(ji,jj,2)
56               trn(ji,jj,ikt,jpno3) = trc_data(ji,jj,3) * 7.6
57               trn(ji,jj,ikt,jppo4) = trc_data(ji,jj,4) * 122.
58               trn(ji,jj,ikt,jpoxy) = trc_data(ji,jj,5)
59               trn(ji,jj,ikt,jpsil) = trc_data(ji,jj,6)
60            ENDIF
61         ENDDO
62      ENDDO
63
64   END SUBROUTINE sed_sfc
65
66#else
67   !!======================================================================
68   !! MODULE sedsfc  :   Dummy module
69   !!======================================================================
70   !! $Id$
71CONTAINS
72   SUBROUTINE sed_sfc ( kt )
73      INTEGER, INTENT(in) :: kt
74      WRITE(*,*) 'sed_sfc: You should not have seen this print! error?', kt
75   END SUBROUTINE sed_sfc
76#endif
77
78END MODULE sedsfc
Note: See TracBrowser for help on using the repository browser.