source: NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP/PISCES/SED/sedsfc.F90 @ 10345

Last change on this file since 10345 was 10345, checked in by smasson, 2 years ago

dev_r10164_HPC09_ESIWACE_PREP_MERGE: merge with trunk@10344, see #2133

  • Property svn:keywords set to Id
File size: 2.6 KB
Line 
1MODULE sedsfc
2   !!======================================================================
3   !!              ***  MODULE  sedsfc  ***
4   !!    Sediment : Data at sediment surface
5   !!=====================================================================
6   !! * Modules used
7   USE sed     ! sediment global variable
8   USE sedarr
9   USE seddta
10
11   PUBLIC sed_sfc
12
13   !! $Id$
14CONTAINS
15
16   SUBROUTINE sed_sfc( kt )
17      !!---------------------------------------------------------------------
18      !!                  ***  ROUTINE sed_sfc ***
19      !!
20      !! ** Purpose :  Give data from sediment model to tracer model
21      !!
22      !!
23      !!   History :
24      !!        !  06-04 (C. Ethe)  Orginal code
25      !!----------------------------------------------------------------------
26      !!* Arguments
27      INTEGER, INTENT(in) ::  kt              ! time step
28
29      ! * local variables
30      INTEGER :: ji, jj, ikt     ! dummy loop indices
31
32      !------------------------------------------------------------------------
33      ! reading variables
34
35      IF( ln_timing )  CALL timing_start('sed_sfc')
36
37      CALL unpack_arr ( jpoce, trc_data(1:jpi,1:jpj,1), iarroce(1:jpoce), pwcp(1:jpoce,1,jwalk) )
38      CALL unpack_arr ( jpoce, trc_data(1:jpi,1:jpj,2), iarroce(1:jpoce), pwcp(1:jpoce,1,jwdic) )
39      CALL unpack_arr ( jpoce, trc_data(1:jpi,1:jpj,3), iarroce(1:jpoce), pwcp(1:jpoce,1,jwno3) )
40      CALL unpack_arr ( jpoce, trc_data(1:jpi,1:jpj,4), iarroce(1:jpoce), pwcp(1:jpoce,1,jwpo4) )
41      CALL unpack_arr ( jpoce, trc_data(1:jpi,1:jpj,5), iarroce(1:jpoce), pwcp(1:jpoce,1,jwoxy) )
42      CALL unpack_arr ( jpoce, trc_data(1:jpi,1:jpj,6), iarroce(1:jpoce), pwcp(1:jpoce,1,jwsil) )
43      CALL unpack_arr ( jpoce, trc_data(1:jpi,1:jpj,7), iarroce(1:jpoce), pwcp(1:jpoce,1,jwnh4) )
44      CALL unpack_arr ( jpoce, trc_data(1:jpi,1:jpj,8), iarroce(1:jpoce), pwcp(1:jpoce,1,jwfe2) )
45
46
47      DO jj = 1,jpj
48         DO ji = 1, jpi
49            ikt = mbkt(ji,jj)
50            IF ( tmask(ji,jj,ikt) == 1 ) THEN
51               trb(ji,jj,ikt,jptal) = trc_data(ji,jj,1)
52               trb(ji,jj,ikt,jpdic) = trc_data(ji,jj,2)
53               trb(ji,jj,ikt,jpno3) = trc_data(ji,jj,3) * 7.625
54               trb(ji,jj,ikt,jppo4) = trc_data(ji,jj,4) * 122.
55               trb(ji,jj,ikt,jpoxy) = trc_data(ji,jj,5)
56               trb(ji,jj,ikt,jpsil) = trc_data(ji,jj,6)
57               trb(ji,jj,ikt,jpnh4) = trc_data(ji,jj,7) * 7.625
58               trb(ji,jj,ikt,jpfer) = trc_data(ji,jj,8)
59            ENDIF
60         ENDDO
61      ENDDO
62
63      IF( ln_timing )  CALL timing_stop('sed_sfc')
64
65   END SUBROUTINE sed_sfc
66
67END MODULE sedsfc
Note: See TracBrowser for help on using the repository browser.