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.
sedsfc.F90 in NEMO/branches/2019/dev_r11351_fldread_with_XIOS/src/TOP/PISCES/SED – NEMO

source: NEMO/branches/2019/dev_r11351_fldread_with_XIOS/src/TOP/PISCES/SED/sedsfc.F90 @ 13463

Last change on this file since 13463 was 13463, checked in by andmirek, 4 years ago

Ticket #2195:update to trunk 13461

  • 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   !! * Substitutions
14#  include "do_loop_substitute.h90"
15   !! $Id$
16CONTAINS
17
18   SUBROUTINE sed_sfc( kt, Kbb )
19      !!---------------------------------------------------------------------
20      !!                  ***  ROUTINE sed_sfc ***
21      !!
22      !! ** Purpose :  Give data from sediment model to tracer model
23      !!
24      !!
25      !!   History :
26      !!        !  06-04 (C. Ethe)  Orginal code
27      !!----------------------------------------------------------------------
28      !!* Arguments
29      INTEGER, INTENT(in) ::  kt              ! time step
30      INTEGER, INTENT(in) ::  Kbb             ! time index
31
32      ! * local variables
33      INTEGER :: ji, jj, ikt     ! dummy loop indices
34
35      !------------------------------------------------------------------------
36      ! reading variables
37
38      IF( ln_timing )  CALL timing_start('sed_sfc')
39
40      CALL unpack_arr ( jpoce, trc_data(1:jpi,1:jpj,1), iarroce(1:jpoce), pwcp(1:jpoce,1,jwalk) )
41      CALL unpack_arr ( jpoce, trc_data(1:jpi,1:jpj,2), iarroce(1:jpoce), pwcp(1:jpoce,1,jwdic) )
42      CALL unpack_arr ( jpoce, trc_data(1:jpi,1:jpj,3), iarroce(1:jpoce), pwcp(1:jpoce,1,jwno3) )
43      CALL unpack_arr ( jpoce, trc_data(1:jpi,1:jpj,4), iarroce(1:jpoce), pwcp(1:jpoce,1,jwpo4) )
44      CALL unpack_arr ( jpoce, trc_data(1:jpi,1:jpj,5), iarroce(1:jpoce), pwcp(1:jpoce,1,jwoxy) )
45      CALL unpack_arr ( jpoce, trc_data(1:jpi,1:jpj,6), iarroce(1:jpoce), pwcp(1:jpoce,1,jwsil) )
46      CALL unpack_arr ( jpoce, trc_data(1:jpi,1:jpj,7), iarroce(1:jpoce), pwcp(1:jpoce,1,jwnh4) )
47      CALL unpack_arr ( jpoce, trc_data(1:jpi,1:jpj,8), iarroce(1:jpoce), pwcp(1:jpoce,1,jwfe2) )
48
49
50      DO_2D( 1, 1, 1, 1 )
51         ikt = mbkt(ji,jj)
52         IF ( tmask(ji,jj,ikt) == 1 ) THEN
53            tr(ji,jj,ikt,jptal,Kbb) = trc_data(ji,jj,1)
54            tr(ji,jj,ikt,jpdic,Kbb) = trc_data(ji,jj,2)
55            tr(ji,jj,ikt,jpno3,Kbb) = trc_data(ji,jj,3) * 7.625
56            tr(ji,jj,ikt,jppo4,Kbb) = trc_data(ji,jj,4) * 122.
57            tr(ji,jj,ikt,jpoxy,Kbb) = trc_data(ji,jj,5)
58            tr(ji,jj,ikt,jpsil,Kbb) = trc_data(ji,jj,6)
59            tr(ji,jj,ikt,jpnh4,Kbb) = trc_data(ji,jj,7) * 7.625
60            tr(ji,jj,ikt,jpfer,Kbb) = trc_data(ji,jj,8)
61         ENDIF
62      END_2D
63
64      IF( ln_timing )  CALL timing_stop('sed_sfc')
65
66   END SUBROUTINE sed_sfc
67
68END MODULE sedsfc
Note: See TracBrowser for help on using the repository browser.