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.
iscplini.F90 in branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/DOM – NEMO

source: branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/DOM/iscplini.F90 @ 9019

Last change on this file since 9019 was 9019, checked in by timgraham, 6 years ago

Merge of dev_CNRS_2017 into branch

File size: 4.3 KB
Line 
1MODULE iscplini
2   !!======================================================================
3   !!                       ***  MODULE  sbciscpl  ***
4   !! Ocean forcing:  river runoff
5   !!=====================================================================
6   !! History :  NEMO  ! 2015-01 P. Mathiot: original
7   !!----------------------------------------------------------------------
8
9   !!----------------------------------------------------------------------
10   !!   iscpl_init     : initialisation routine (namelist)
11   !!   iscpl_alloc    : allocation of correction variables
12   !!----------------------------------------------------------------------
13   USE oce             ! global tra/dyn variable
14   USE dom_oce         ! ocean space and time domain
15   !
16   USE lib_mpp         ! MPP library
17   USE lib_fortran     ! MPP library
18   USE in_out_manager  ! I/O manager
19
20   IMPLICIT NONE
21   PRIVATE
22   
23   PUBLIC   iscpl_init     
24   PUBLIC   iscpl_alloc 
25   !!                                                      !!* namsbc_iscpl namelist *
26   LOGICAL , PUBLIC                                        ::   ln_hsb
27   INTEGER , PUBLIC                                        ::   nn_fiscpl, nstp_iscpl
28   INTEGER , PUBLIC                                        ::   nn_drown
29   REAL(wp), PUBLIC                                        ::   rdt_iscpl
30   !!                                                      !!* namsbc_iscpl namelist *
31   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:  ) ::   hdiv_iscpl
32   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) ::   htsc_iscpl
33   !!----------------------------------------------------------------------
34   !! NEMO/OPA 3.3 , NEMO Consortium (2010)
35   !! $Id: sbcrnf.F90 4666 2014-06-11 12:52:23Z mathiot $
36   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
37   !!----------------------------------------------------------------------
38CONTAINS
39
40   INTEGER FUNCTION iscpl_alloc()
41      !!----------------------------------------------------------------------
42      !!                ***  ROUTINE sbc_iscpl_alloc  ***
43      !!----------------------------------------------------------------------
44      ALLOCATE( htsc_iscpl(jpi,jpj,jpk,jpts) , hdiv_iscpl(jpi,jpj,jpk) , STAT=iscpl_alloc )
45         !
46      IF( lk_mpp          )   CALL mpp_sum ( iscpl_alloc )
47      IF( iscpl_alloc > 0 )   CALL ctl_warn('iscpl_alloc: allocation of arrays failed')
48   END FUNCTION iscpl_alloc
49
50
51   SUBROUTINE iscpl_init()
52      !!----------------------------------------------------------------------
53      INTEGER ::   ios           ! Local integer output status for namelist read
54      NAMELIST/namsbc_iscpl/ nn_fiscpl, ln_hsb, nn_drown
55      !!----------------------------------------------------------------------
56      !
57      nn_fiscpl = 0
58      ln_hsb    = .FALSE.
59      REWIND( numnam_ref )              ! Namelist namsbc_iscpl in reference namelist : Ice sheet coupling
60      READ  ( numnam_ref, namsbc_iscpl, IOSTAT = ios, ERR = 901)
61901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_iscpl in reference namelist', lwp )
62
63      REWIND( numnam_cfg )              ! Namelist namsbc_iscpl in configuration namelist : Ice Sheet coupling
64      READ  ( numnam_cfg, namsbc_iscpl, IOSTAT = ios, ERR = 902 )
65902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_iscpl in configuration namelist', lwp )
66      IF(lwm) WRITE ( numond, namsbc_iscpl )
67      !
68      nstp_iscpl=MIN(nn_fiscpl, nitend-nit000+1) ! the coupling period have to be less or egal than the total number of time step
69      rdt_iscpl = nstp_iscpl * rn_rdt
70      !
71      IF (lwp) THEN
72         WRITE(numout,*) 'iscpl_rst:'
73         WRITE(numout,*) '~~~~~~~~~'
74         WRITE(numout,*) ' coupling     flag (ln_iscpl )            = ', ln_iscpl
75         WRITE(numout,*) ' conservation flag (ln_hsb   )            = ', ln_hsb
76         WRITE(numout,*) ' nb of stp for cons (rn_fiscpl)           = ', nstp_iscpl
77         IF (nstp_iscpl .NE. nn_fiscpl) WRITE(numout,*) 'W A R N I N G: nb of stp for cons has been modified &
78            &                                           (larger than run length)'
79         WRITE(numout,*) ' coupling time step                       = ', rdt_iscpl
80         WRITE(numout,*) ' number of call of the extrapolation loop = ', nn_drown
81      ENDIF
82      !
83   END SUBROUTINE iscpl_init
84
85   !!======================================================================
86END MODULE iscplini
Note: See TracBrowser for help on using the repository browser.