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.
iceini_2.F90 in branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/LIM_SRC_2 – NEMO

source: branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/LIM_SRC_2/iceini_2.F90 @ 4319

Last change on this file since 4319 was 4319, checked in by clevy, 10 years ago

change ctl_opn status of output.namelist files for compatibility with some compilers

  • Property svn:keywords set to Id
File size: 7.7 KB
RevLine 
[821]1MODULE iceini_2
[3]2   !!======================================================================
3   !!                       ***  MODULE iceini   ***
[821]4   !!   Sea-ice model : LIM 2.0 Sea ice model Initialization
[3]5   !!======================================================================
[2715]6   !! History :  1.0  ! 2002-08  (G. Madec)  F90: Free form and modules
7   !!            2.0  ! 2003-08  (C. Ethe)  add ice_run
8   !!            3.3  ! 2009-05  (G. Garric, C. Bricaud) addition of the lim2_evp case
9   !!            4.0  ! 2011-02  (G. Madec) dynamical allocation
[508]10   !!----------------------------------------------------------------------
[821]11#if defined key_lim2
[3]12   !!----------------------------------------------------------------------
[821]13   !!   'key_lim2' :                                  LIM 2.0 sea-ice model
[3]14   !!----------------------------------------------------------------------
[3625]15   !!   ice_init_2    : sea-ice model initialization
16   !!   ice_run_2     : Definition some run parameter for ice model
[3]17   !!----------------------------------------------------------------------
[3625]18   USE phycst         ! physical constants
19   USE dom_oce        ! ocean domain
20   USE sbc_oce        ! surface boundary condition: ocean
21   USE sbc_ice        ! LIM2 surface boundary condition
22   USE dom_ice_2      ! LIM2 ice domain
23   USE par_ice_2      ! LIM2 parameters
24   USE thd_ice_2      ! LIM2 thermodynamical variables
25   USE ice_2          ! LIM2 ice variable
26   USE limmsh_2       ! LIM2 mesh
27   USE limistate_2    ! LIM2 initial state
28   USE limrst_2       ! LIM2 restart
29   USE limsbc_2       ! LIM2 surface boundary condition
30   USE in_out_manager ! I/O manager
31   USE lib_mpp        ! MPP library
32   USE lib_fortran    ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined) 
[2715]33
[3]34   IMPLICIT NONE
35   PRIVATE
36
[1228]37   PUBLIC   ice_init_2               ! called by sbcice_lim_2.F90
[3]38
39   !!----------------------------------------------------------------------
[2715]40   !! NEMO/LIM2 4.0 , UCL - NEMO Consortium (2011)
[1156]41   !! $Id$
[2528]42   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
[3]43   !!----------------------------------------------------------------------
44CONTAINS
45
[821]46   SUBROUTINE ice_init_2
[3]47      !!----------------------------------------------------------------------
[821]48      !!                  ***  ROUTINE ice_init_2  ***
[3]49      !!
[2528]50      !! ** purpose :   initialisation of LIM-2 domain and variables 
[3]51      !!----------------------------------------------------------------------
[2715]52      INTEGER :: ierr
53      !!----------------------------------------------------------------------
[508]54      !
[2715]55      IF(lwp) THEN
56         WRITE(numout,*)
57         WRITE(numout,*) 'ice_init_2 : LIM-2 sea-ice - initialization'
58         WRITE(numout,*) '~~~~~~~~~~~   '
59      ENDIF
60      !                                ! Allocate the ice arrays
61      ierr =        ice_alloc_2    ()       ! ice variables
62      ierr = ierr + dom_ice_alloc_2()       ! domain
63      ierr = ierr + sbc_ice_alloc  ()       ! surface forcing
64      ierr = ierr + thd_ice_alloc_2()       ! thermodynamics
[3294]65
[2715]66      IF( lk_mpp    )   CALL mpp_sum( ierr )
67      IF( ierr /= 0 )   CALL ctl_stop( 'STOP', 'ice_init_2 : unable to allocate ice arrays' )
68
69      !                                ! adequation jpk versus ice/snow layers
70      IF( jpl > jpk  .OR.  jplayersp1 > jpk  )   CALL ctl_stop( 'STOP',           &
71         &     'ice_init: the 3rd dimension of workspace arrays is too small.',   &
72         &     'use more ocean levels or less ice layers/categories.' )
73
[4147]74      !                                ! Open the reference and configuration namelist files and namelist output file
75      CALL ctl_opn( numnam_ice_ref, 'namelist_ice_ref',    'OLD',     'FORMATTED', 'SEQUENTIAL', -1, numout, lwp ) 
76      CALL ctl_opn( numnam_ice_cfg, 'namelist_ice_cfg',    'OLD',     'FORMATTED', 'SEQUENTIAL', -1, numout, lwp )
[4319]77      CALL ctl_opn( numoni,         'output.namelist.ice', 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, 1 )
[4309]78
[4306]79      ! When Diurnal cycle, core bulk and LIM2  are activated,
80      ! a daily mean qsr is computed for tracer/biogeochemistery model                                 !
81      IF( ltrcdm2dc ) ltrcdm2dc_ice = .TRUE.
82
83      !                                ! Open the namelist file
[2528]84      !   
85      CALL ice_run_2                   ! read in namelist some run parameters
86      !         
87      rdt_ice = nn_fsbc * rdttra(1)    ! sea-ice time step
88      numit   = nit000 - 1
89      !
90      CALL lim_msh_2                   ! ice mesh initialization
91      !
92      !                                ! Initial sea-ice state
93      IF( .NOT.ln_rstart ) THEN   ;   CALL lim_istate_2     ! start from rest: sea-ice deduced from sst
94      ELSE                        ;   CALL lim_rst_read_2   ! start from a restart file
[3]95      ENDIF
[508]96      !
[2715]97      tn_ice(:,:,1) = sist(:,:)        ! ice temperature  known by the ocean
98      fr_i  (:,:)   = 1.0 - frld(:,:)  ! sea-ice fraction known by the ocean
[2528]99      !
[2715]100      CALL lim_sbc_init_2              ! ice surface boundary condition   
101      !
[2730]102      IF( lk_lim2_vp )   THEN   ;   IF(lwp) WRITE(numout,*) '                VP  rheology - B-grid case'
103      ELSE                      ;   IF(lwp) WRITE(numout,*) '                EVP rheology - C-grid case'
[2715]104      ENDIF
105      !
[821]106   END SUBROUTINE ice_init_2
[3]107
108
[821]109   SUBROUTINE ice_run_2
[3]110      !!-------------------------------------------------------------------
[821]111      !!                  ***  ROUTINE ice_run_2 ***
[3]112      !!                 
113      !! ** Purpose :   Definition some run parameter for ice model
114      !!
115      !! ** Method  :   Read the namicerun namelist and check the parameter
116      !!       values called at the first timestep (nit000)
117      !!
118      !! ** input   :   Namelist namicerun
119      !!-------------------------------------------------------------------
[1229]120      NAMELIST/namicerun/ cn_icerst_in, cn_icerst_out, ln_limdyn, ln_limdmp, acrit, hsndif, hicdif
[4147]121      INTEGER  ::   ios                 ! Local integer output status for namelist read
[3]122      !!-------------------------------------------------------------------
[508]123      !                   
[4147]124      REWIND( numnam_ice_ref )              ! Namelist namicerun in reference namelist : Parameters for ice
125      READ  ( numnam_ice_ref, namicerun, IOSTAT = ios, ERR = 901)
126901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicerun in reference namelist', lwp )
127
128      REWIND( numnam_ice_cfg )              ! Namelist namicerun in configuration namelist : Parameters for ice
129      READ  ( numnam_ice_cfg, namicerun, IOSTAT = ios, ERR = 902 )
130902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicerun in configuration namelist', lwp )
131      WRITE ( numoni, namicerun )
[2528]132      !
133      IF(lwp) THEN                              ! control print
[3]134         WRITE(numout,*)
135         WRITE(numout,*) 'ice_run : ice share parameters for dynamics/advection/thermo of sea-ice'
136         WRITE(numout,*) ' ~~~~~~'
[76]137         WRITE(numout,*) '   switch for ice dynamics (1) or not (0)      ln_limdyn   = ', ln_limdyn
[675]138         WRITE(numout,*) '   Ice damping                                 ln_limdmp   = ', ln_limdmp
[76]139         WRITE(numout,*) '   minimum fraction for leads in the NH (SH)  acrit(1/2)   = ', acrit(:)
[3]140         WRITE(numout,*) '   computation of temp. in snow (=0) or not (=9999) hsndif = ', hsndif
141         WRITE(numout,*) '   computation of temp. in ice  (=0) or not (=9999) hicdif = ', hicdif
142      ENDIF
[508]143      !
[821]144   END SUBROUTINE ice_run_2
[3]145
146#else
147   !!----------------------------------------------------------------------
[821]148   !!   Default option :        Empty module       NO LIM 2.0 sea-ice model
[3]149   !!----------------------------------------------------------------------
150CONTAINS
[821]151   SUBROUTINE ice_init_2      ! Empty routine
152   END SUBROUTINE ice_init_2
[3]153#endif
154
155   !!======================================================================
[821]156END MODULE iceini_2
Note: See TracBrowser for help on using the repository browser.