source: trunk/NEMO/LIM_SRC_2/iceini_2.F90 @ 900

Last change on this file since 900 was 900, checked in by rblod, 13 years ago

Update 1D configuration according to SBC and LIM3, see ticket #117

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 5.5 KB
Line 
1MODULE iceini_2
2   !!======================================================================
3   !!                       ***  MODULE iceini   ***
4   !!   Sea-ice model : LIM 2.0 Sea ice model Initialization
5   !!======================================================================
6   !! History :   1.0  !  02-08  (G. Madec)  F90: Free form and modules
7   !!             2.0  !  03-08  (C. Ethe)  add ice_run
8   !!----------------------------------------------------------------------
9#if defined key_lim2
10   !!----------------------------------------------------------------------
11   !!   'key_lim2' :                                  LIM 2.0 sea-ice model
12   !!----------------------------------------------------------------------
13   !!----------------------------------------------------------------------
14   !!   ice_init_2       : sea-ice model initialization
15   !!   ice_run_2        : Definition some run parameter for ice model
16   !!----------------------------------------------------------------------
17   USE dom_oce
18   USE dom_ice_2
19   USE ice_oce         ! ice variables
20   USE sbc_oce         ! surface boundary condition: ocean
21   USE sbc_ice         ! surface boundary condition: ice
22   USE phycst          ! Define parameters for the routines
23   USE ocfzpt
24   USE ice_2
25   USE limmsh_2
26   USE limistate_2
27   USE limrst_2   
28   USE in_out_manager
29     
30   IMPLICIT NONE
31   PRIVATE
32
33   PUBLIC   ice_init_2               ! called by opa.F90
34
35   LOGICAL , PUBLIC               ::   ln_limdyn = .TRUE.   !: flag for ice dynamics (T) or not (F)
36   REAL(wp), PUBLIC               ::   hsndif = 0.e0        !: computation of temp. in snow (0) or not (9999)
37   REAL(wp), PUBLIC               ::   hicdif = 0.e0        !: computation of temp. in ice (0) or not (9999)
38   REAL(wp), PUBLIC, DIMENSION(2) ::   acrit  = (/ 1.e-06 , 1.e-06 /)    !: minimum fraction for leads in
39      !                                                                  !  north and south hemisphere
40   !!----------------------------------------------------------------------
41   !!   LIM 2.0,  UCL-LOCEAN-IPSL (2005)
42   !! $ Id: $
43   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt
44   !!----------------------------------------------------------------------
45
46CONTAINS
47
48   SUBROUTINE ice_init_2
49      !!----------------------------------------------------------------------
50      !!                  ***  ROUTINE ice_init_2  ***
51      !!
52      !! ** purpose :   
53      !!----------------------------------------------------------------------
54      CHARACTER(len=80) :: namelist_icename
55      !!----------------------------------------------------------------------
56      !
57      ! Open the namelist file
58      namelist_icename = 'namelist_ice'
59      CALL ctlopn(numnam_ice,namelist_icename,'OLD', 'FORMATTED', 'SEQUENTIAL',   &
60                     1,numout,.FALSE.,1)     
61      CALL ice_run_2                    !  read in namelist some run parameters
62                 
63      ! Louvain la Neuve Ice model
64      dtsd2   = nn_fsbc * rdttra(1) * 0.5
65      rdt_ice = nn_fsbc * rdttra(1)
66
67      CALL lim_msh_2                  ! ice mesh initialization
68     
69      ! Initial sea-ice state
70      IF( .NOT.ln_rstart ) THEN
71         CALL lim_istate_2            ! start from rest: sea-ice deduced from sst
72      ELSE
73         CALL lim_rst_read_2          ! start from a restart file
74      ENDIF
75     
76      tn_ice(:,:) = sist(:,:)         ! initialisation of ice temperature   
77      freeze(:,:) = 1.0 - frld(:,:)   ! initialisation of sea/ice cover   
78# if defined key_coupled
79      alb_ice(:,:) = albege(:,:)      ! sea-ice albedo
80# endif
81      !
82   END SUBROUTINE ice_init_2
83
84
85   SUBROUTINE ice_run_2
86      !!-------------------------------------------------------------------
87      !!                  ***  ROUTINE ice_run_2 ***
88      !!                 
89      !! ** Purpose :   Definition some run parameter for ice model
90      !!
91      !! ** Method  :   Read the namicerun namelist and check the parameter
92      !!       values called at the first timestep (nit000)
93      !!
94      !! ** input   :   Namelist namicerun
95      !!-------------------------------------------------------------------
96      NAMELIST/namicerun/ ln_limdyn, ln_limdmp, acrit, hsndif, hicdif
97      !!-------------------------------------------------------------------
98      !                   
99      REWIND ( numnam_ice )                       ! Read Namelist namicerun
100      READ   ( numnam_ice , namicerun )
101
102      IF(lwp) THEN
103         WRITE(numout,*)
104         WRITE(numout,*) 'ice_run : ice share parameters for dynamics/advection/thermo of sea-ice'
105         WRITE(numout,*) ' ~~~~~~'
106         WRITE(numout,*) '   switch for ice dynamics (1) or not (0)      ln_limdyn   = ', ln_limdyn
107         WRITE(numout,*) '   Ice damping                                 ln_limdmp   = ', ln_limdmp
108         WRITE(numout,*) '   minimum fraction for leads in the NH (SH)  acrit(1/2)   = ', acrit(:)
109         WRITE(numout,*) '   computation of temp. in snow (=0) or not (=9999) hsndif = ', hsndif
110         WRITE(numout,*) '   computation of temp. in ice  (=0) or not (=9999) hicdif = ', hicdif
111      ENDIF
112      !
113   END SUBROUTINE ice_run_2
114
115#else
116   !!----------------------------------------------------------------------
117   !!   Default option :        Empty module       NO LIM 2.0 sea-ice model
118   !!----------------------------------------------------------------------
119CONTAINS
120   SUBROUTINE ice_init_2      ! Empty routine
121   END SUBROUTINE ice_init_2
122#endif
123
124   !!======================================================================
125END MODULE iceini_2
Note: See TracBrowser for help on using the repository browser.