source: trunk/NEMO/LIM_SRC/iceini.F90 @ 12

Last change on this file since 12 was 12, checked in by opalod, 17 years ago

CT : UPDATE001 : First major NEMO update

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 5.4 KB
Line 
1MODULE iceini
2   !!======================================================================
3   !!                       ***  MODULE iceini   ***
4   !!   Sea-ice model : LIM Sea ice model Initialization
5   !!======================================================================
6#if defined key_ice_lim
7   !!----------------------------------------------------------------------
8   !!   'key_ice_lim' :                                   LIM sea-ice model
9   !!----------------------------------------------------------------------
10   !!   ice_init       : sea-ice model initialization
11   !!----------------------------------------------------------------------
12   USE dom_oce
13   USE in_out_manager
14   USE ice_oce         ! ice variables
15   USE flx_oce
16   USE phycst          ! Define parameters for the routines
17   USE ocfzpt
18   USE ice
19   USE limmsh
20   USE limistate
21   USE limrst
22
23   IMPLICIT NONE
24   PRIVATE
25
26   !! * Routine accessibility
27   PUBLIC ice_init                 ! called by opa.F90
28
29   !! * Share Module variables
30   LOGICAL , PUBLIC  ::   & !!! ** init namelist (namicerun) **
31      ldyn   = .TRUE.        !: flag for ice dynamics (T) or not (F)
32   INTEGER , PUBLIC  ::   &  !:
33      nstart ,            &  !: iteration number of the begining of the run
34      nlast  ,            &  !: iteration number of the end of the run
35      nitrun ,            &  !: number of iteration
36      numit                  !: iteration number
37   REAL(wp), PUBLIC  ::   &  !:
38      hsndif = 0.0 ,      &  !: computation of temp. in snow (0) or not (9999)
39      hicdif = 0.0 ,      &  !: computation of temp. in ice (0) or not (9999)
40      tpstot                 !: time of the run in seconds
41   REAL(wp), PUBLIC, DIMENSION(2)  ::  &  !:
42      acrit  = (/ 1.0e-06 , 1.0e-06 /)    !: minimum fraction for leads in
43      !                                   !  north and south hemisphere
44   !!----------------------------------------------------------------------
45   !!   LIM 2.0   UCL-LODYC-IPSL (2003)
46   !!----------------------------------------------------------------------
47
48CONTAINS
49
50   SUBROUTINE ice_init
51      !!----------------------------------------------------------------------
52      !!                  ***  ROUTINE ice_init  ***
53      !!
54      !! ** purpose :   
55      !!
56      !! History :
57      !!   8.5  !  02-08  (G. Madec)  F90: Free form and modules
58      !!----------------------------------------------------------------------
59
60      ! Open the namelist file
61      OPEN( numnam_ice, FILE= 'namelist_ice', FORM='formatted', STATUS = 'old') 
62
63      CALL ice_run                    !  read in namelist some run parameters
64                 
65      ! Louvain la Neuve Ice model
66      IF( nacc == 1 ) THEN
67          dtsd2   = nfice * rdtmin * 0.5
68          rdt_ice = nfice * rdtmin
69      ELSE
70          dtsd2   = nfice * rdt * 0.5
71          rdt_ice = nfice * rdt
72      ENDIF
73
74      CALL lim_msh                    ! ice mesh initialization
75     
76      ! Initial sea-ice state
77      IF( .NOT.ln_rstart ) THEN
78         numit = 0
79         CALL lim_istate              ! start from rest: sea-ice deduced from sst
80      ELSE
81         CALL lim_rst_read( numit )   ! start from a restart file
82      ENDIF
83     
84      tn_ice(:,:) = sist(:,:)         ! initialisation of ice temperature   
85      freeze(:,:) = 1.0 - frld(:,:)   ! initialisation of sea/ice cover   
86# if defined key_coupled
87      alb_ice(:,:) = albege(:,:)      ! sea-ice albedo
88# endif
89     
90      nstart = numit  + nfice     
91      nitrun = nitend - nit000 + 1 
92      nlast  = numit  + nitrun 
93
94      IF( nstock == 0  )  nstock = nlast + 1
95
96   END SUBROUTINE ice_init
97
98
99   SUBROUTINE ice_run
100      !!-------------------------------------------------------------------
101      !!                  ***  ROUTINE ice_run ***
102      !!                 
103      !! ** Purpose :   Definition some run parameter for ice model
104      !!
105      !! ** Method  :   Read the namicerun namelist and check the parameter
106      !!       values called at the first timestep (nit000)
107      !!
108      !! ** input   :   Namelist namicerun
109      !!
110      !! history :
111      !!   2.0  !  03-08 (C. Ethe)  Original code
112      !!-------------------------------------------------------------------
113      NAMELIST/namicerun/ ldyn, acrit, hsndif, hicdif
114      !!-------------------------------------------------------------------
115
116      !                                           ! Read Namelist namicerun
117      REWIND ( numnam_ice )
118      READ   ( numnam_ice , namicerun )
119      IF(lwp) THEN
120         WRITE(numout,*)
121         WRITE(numout,*) 'ice_run : ice share parameters for dynamics/advection/thermo of sea-ice'
122         WRITE(numout,*) ' ~~~~~~'
123         WRITE(numout,*) '   switch for ice dynamics (1) or not (0)           ldyn   = ', ldyn
124         WRITE(numout,*) '   minimum fraction for leads in the NH (SH)  acrit(1/2)   = ', acrit
125         WRITE(numout,*) '   computation of temp. in snow (=0) or not (=9999) hsndif = ', hsndif
126         WRITE(numout,*) '   computation of temp. in ice  (=0) or not (=9999) hicdif = ', hicdif
127      ENDIF
128     
129   END SUBROUTINE ice_run
130
131#else
132   !!----------------------------------------------------------------------
133   !!   Default option :        Empty module           NO LIM sea-ice model
134   !!----------------------------------------------------------------------
135CONTAINS
136   SUBROUTINE ice_init        ! Empty routine
137   END SUBROUTINE ice_init
138#endif
139
140   !!======================================================================
141END MODULE iceini
Note: See TracBrowser for help on using the repository browser.