source: branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icestp.F90 @ 8515

Last change on this file since 8515 was 8515, checked in by clem, 3 years ago

changes in style - part5 - very nearly finished

File size: 22.5 KB
Line 
1MODULE icestp
2   !!======================================================================
3   !!                       ***  MODULE  icestp  ***
4   !! Surface module :  update the ocean surface boundary condition over ice
5   !!       &           covered area using LIM sea-ice model
6   !! Sea-Ice model  :  LIM-3 Sea ice model time-stepping
7   !!=====================================================================
8   !! History :  2.0  ! 2006-12  (M. Vancoppenolle) Original code
9   !!            3.0  ! 2008-02  (C. Talandier)  Surface module from icestp.F90
10   !!             -   ! 2008-04  (G. Madec)  sltyle and ice_ctl routine
11   !!            3.3  ! 2010-11  (G. Madec) ice-ocean stress always computed at each ocean time-step
12   !!            3.4  ! 2011-01  (A Porter)  dynamical allocation
13   !!             -   ! 2012-10  (C. Rousset)  add ice_dia
14   !!            3.6  ! 2014-07  (M. Vancoppenolle, G. Madec, O. Marti) revise coupled interface
15   !!            4.0  ! 2016-06  (L. Brodeau) new unified bulk routine (based on AeroBulk)
16   !!----------------------------------------------------------------------
17#if defined key_lim3
18   !!----------------------------------------------------------------------
19   !!   'key_lim3'                                    LIM 3.0 sea-ice model
20   !!----------------------------------------------------------------------
21   !!   ice_stp       : sea-ice model time-stepping and update ocean surf. boundary cond. over ice-covered area
22   !!   ice_init      :
23   !!   ice_run_init  :
24   !!----------------------------------------------------------------------
25   USE oce            ! ocean dynamics and tracers
26   USE dom_oce        ! ocean space and time domain
27   USE c1d            ! 1D vertical configuration
28   USE ice            ! sea-ice: variables
29   USE ice1D          ! sea-ice: thermodynamical 1D variables
30   !
31   USE sbc_oce        ! Surface boundary condition: ocean fields
32   USE sbc_ice        ! Surface boundary condition: ice   fields
33   USE iceforcing     ! sea-ice: Surface boundary condition       !!gm why not icesbc module name
34   !
35   USE phycst         ! Define parameters for the routines
36   USE eosbn2         ! equation of state
37   USE icerhg         ! sea-ice: rheology
38   USE iceadv         ! sea-ice: advection
39   USE icethd         ! sea-ice: thermodynamics
40   USE icerdgrft      ! sea-ice: ridging/rafting
41   USE iceupdate      ! sea-ice: sea surface boundary condition update
42   USE icedia         ! sea-ice: budget diagnostics
43   USE icewri         ! sea-ice: outputs
44   USE icerst         ! sea-ice: restarts
45   USE icecor         ! sea-ice: corrections
46   USE icevar         ! sea-ice: operations
47   USE icectl         ! sea-ice: control
48   ! MV MP 2016
49   USE limmp          ! sea-ice: melt ponds
50   ! END MV MP 2016
51   USE iceistate      ! sea-ice: initial state
52   USE icethd_sal     ! sea-ice: thermodynamics and salinity
53   USE iceitd         ! sea-ice: remapping thickness distribution
54   USE icealb         ! sea-ice: albedo
55   !
56   USE bdy_oce , ONLY : ln_bdy   ! flag for bdy
57   USE bdyice         ! unstructured open boundary data for sea-ice
58# if defined key_agrif
59   USE agrif_ice
60   USE agrif_lim3_update
61   USE agrif_lim3_interp
62# endif
63   !
64   USE in_out_manager ! I/O manager
65   USE iom            ! I/O manager library
66   USE prtctl         ! Print control
67   USE lib_fortran    !
68   USE lbclnk         ! lateral boundary condition - MPP link
69   USE lib_mpp        ! MPP library
70   USE timing         ! Timing
71
72   IMPLICIT NONE
73   PRIVATE
74
75   PUBLIC   ice_stp    ! called by sbcmod.F90
76   PUBLIC   ice_init   ! called by sbcmod.F90
77
78   !! * Substitutions
79#  include "vectopt_loop_substitute.h90"
80   !!----------------------------------------------------------------------
81   !! NEMO/ICE 4.0 , NEMO Consortium (2017)
82   !! $Id: icestp.F90 8319 2017-07-11 15:00:44Z clem $
83   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
84   !!----------------------------------------------------------------------
85CONTAINS
86
87   SUBROUTINE ice_stp( kt, ksbc )
88      !!---------------------------------------------------------------------
89      !!                  ***  ROUTINE ice_stp  ***
90      !!
91      !! ** Purpose :   sea-ice model time-stepping and update ocean surface
92      !!              boundary condition over ice-covered area
93      !!
94      !! ** Method  :   ice model time stepping
95      !!              - call the ice dynamics routine
96      !!              - call the ice advection/diffusion routine
97      !!              - call the ice thermodynamics routine
98      !!              - call the routine that computes mass and
99      !!                heat fluxes at the ice/ocean interface
100      !!              - save the outputs
101      !!              - save the outputs for restart when necessary
102      !!
103      !! ** Action  : - time evolution of the LIM sea-ice model
104      !!              - update all sbc variables below sea-ice:
105      !!                utau, vtau, taum, wndm, qns , qsr, emp , sfx
106      !!---------------------------------------------------------------------
107      INTEGER, INTENT(in) ::   kt      ! ocean time step
108      INTEGER, INTENT(in) ::   ksbc    ! flux formulation (user defined, bulk, or Pure Coupled)
109      !
110      INTEGER ::   jl   ! dummy loop index
111      !!----------------------------------------------------------------------
112      !
113      IF( nn_timing == 1 )   CALL timing_start('ice_stp')
114      !
115      !                                      !-----------------------!
116      IF( MOD( kt-1, nn_fsbc ) == 0 ) THEN   ! --- Ice time step --- !
117         !                                   !-----------------------!
118         !
119         kt_ice = kt                              ! -- Ice model time step
120         !
121# if defined key_agrif
122         IF( .NOT. Agrif_Root() )  lim_nbstep = MOD( lim_nbstep, Agrif_irhot() * Agrif_Parent(nn_fsbc) / nn_fsbc ) + 1
123# endif
124         u_oce(:,:) = ssu_m(:,:)                  ! -- mean surface ocean current
125         v_oce(:,:) = ssv_m(:,:)
126         !
127         CALL eos_fzp( sss_m(:,:) , t_bo(:,:) )   ! -- freezing temperature [Kelvin] (set to rt0 over land)
128         t_bo(:,:) = ( t_bo(:,:) + rt0 ) * tmask(:,:,1) + rt0 * ( 1._wp - tmask(:,:,1) )
129         !
130         CALL store_fields                        ! -- Store now ice values
131         !
132         !------------------------------------------------!
133         ! --- Dynamical coupling with the atmosphere --- !
134         !------------------------------------------------!
135         ! It provides the following fields used in sea ice model:
136         !    utau_ice, vtau_ice = surface ice stress [N/m2]
137         !------------------------------------------------!
138         CALL ice_forcing_tau( kt, ksbc, utau_ice, vtau_ice )
139                                     
140         !-------------------------------------!
141         ! --- ice dynamics and advection  --- !
142         !-------------------------------------!
143         CALL ice_diag0             ! set diag of mass, heat and salt fluxes to 0
144         CALL ice_rst_opn( kt )     ! Open Ice restart file (if necessary)
145         !
146         SELECT CASE( nice_dyn )   
147         CASE ( np_dynFULL )          !==  all dynamical processes  ==!
148            CALL ice_rhg   ( kt )          ! -- rheology 
149            CALL ice_adv   ( kt )          ! -- advection of ice
150            CALL ice_rdgrft( kt )          ! -- ridging/rafting
151            CALL ice_cor   ( kt , 1 )      ! -- Corrections
152         CASE ( np_dyn )              !==  pure dynamics only ==!   (no ridging/rafting)   (nono cat. case 2)
153            CALL ice_rhg   ( kt )          ! -- rheology 
154            CALL ice_adv   ( kt )          ! -- advection of ice
155            CALL ice_cor   ( kt , 1 )      ! -- Corrections
156         CASE ( np_dynPURE )          !==  pure dynamics only ==!   (nn_icedyn= 0 or 1 )
157            CALL ice_rhg   ( kt )          ! -- rheology 
158            CALL ice_adv   ( kt )          ! -- advection of ice
159         END SELECT
160         !
161         !                          !==  lateral boundary conditions  ==!
162#if defined key_agrif
163         IF( .NOT. Agrif_Root()     )   CALL agrif_interp_lim3('T')   ! -- AGRIF
164#endif
165         IF( ln_icethd .AND. ln_bdy )   CALL bdy_ice( kt )            ! -- bdy ice thermo
166         !
167         !
168         !                          !==  previous lead fraction and ice volume for flux calculations
169         !
170         CALL ice_var_glo2eqv            ! ht_i and ht_s for ice albedo calculation
171         CALL ice_var_agg(1)             ! at_i for coupling
172         CALL store_fields               ! Store now ice values
173
174         !------------------------------------------------------!
175         ! --- Thermodynamical coupling with the atmosphere --- !
176         !------------------------------------------------------!
177         ! It provides the following fields used in sea ice model:
178         !    fr1_i0  , fr2_i0     = 1sr & 2nd fraction of qsr penetration in ice  [%]
179         !    emp_oce , emp_ice    = E-P over ocean and sea ice                    [Kg/m2/s]
180         !    sprecip              = solid precipitation                           [Kg/m2/s]
181         !    evap_ice             = sublimation                                   [Kg/m2/s]
182         !    qsr_tot , qns_tot    = solar & non solar heat flux (total)           [W/m2]
183         !    qsr_ice , qns_ice    = solar & non solar heat flux over ice          [W/m2]
184         !    dqns_ice             = non solar  heat sensistivity                  [W/m2]
185         !    qemp_oce, qemp_ice,  = sensible heat (associated with evap & precip) [W/m2]
186         !    qprec_ice, qevap_ice
187         !------------------------------------------------------!
188         CALL ice_forcing_flx( kt, ksbc )
189
190         !----------------------------!
191         ! --- ice thermodynamics --- !
192         !----------------------------!
193         IF( ln_icethd )            CALL ice_thd( kt )          ! -- Ice thermodynamics     
194
195         ! MV MP 2016
196         IF ( ln_pnd )              CALL lim_mp( kt )           ! -- Melt ponds
197         ! END MV MP 2016
198
199         IF( ln_icethd )            CALL ice_cor( kt , 2 )      ! -- Corrections
200         ! ---
201# if defined key_agrif
202         IF( .NOT. Agrif_Root() )   CALL agrif_update_lim3( kt )
203# endif
204                                    CALL ice_var_glo2eqv        ! necessary calls (at least for coupling)
205                                    CALL ice_var_agg( 2 )       ! necessary calls (at least for coupling)
206                                    !
207!! clem: one should switch the calculation of the fluxes onto the parent grid but the following calls do not work
208!!       moreover it should only be called at the update frequency (as in agrif_lim3_update.F90)
209!!# if defined key_agrif
210!!         IF( .NOT. Agrif_Root() )   CALL Agrif_ChildGrid_To_ParentGrid()
211!!# endif
212                                    CALL ice_update_flx( kt )   ! -- Update ocean surface mass, heat and salt fluxes
213!!# if defined key_agrif
214!!         IF( .NOT. Agrif_Root() )   CALL Agrif_ParentGrid_To_ChildGrid()
215!!# endif
216         IF( ln_icediahsb )         CALL ice_dia( kt )        ! -- Diagnostics and outputs
217         !
218                                    CALL ice_wri( 1 )         ! -- Ice outputs
219         !
220         IF( kt == nit000 .AND. ln_rstart )   &                !!gm  I don't understand the ln_rstart, if needed, add a comment, please
221            &                       CALL iom_close( numrir )  ! close input ice restart file
222         !
223         IF( lrst_ice )             CALL ice_rst_write( kt )  ! -- Ice restart file
224         !
225         IF( ln_icectl )            CALL ice_ctl( kt )        ! alerts in case of model crash
226         !
227      ENDIF   ! End sea-ice time step only
228
229      !-------------------------!
230      ! --- Ocean time step --- !
231      !-------------------------!
232      ! Update surface ocean stresses (only in ice-dynamic case) otherwise the atm.-ocean stresses are used everywhere
233      !    using before instantaneous surf. currents
234      IF( ln_icedyn )               CALL ice_update_tau( kt, ub(:,:,1), vb(:,:,1) )
235!!gm   remark, the ocean-ice stress is not saved in ice diag call above .....  find a solution!!!
236      !
237      IF( nn_timing == 1 ) CALL timing_stop('ice_stp')
238      !
239   END SUBROUTINE ice_stp
240
241
242   SUBROUTINE ice_init
243      !!----------------------------------------------------------------------
244      !!                  ***  ROUTINE ice_init  ***
245      !!
246      !! ** purpose :   Allocate all the dynamic arrays of the LIM-3 modules
247      !!----------------------------------------------------------------------
248      INTEGER :: ji, jj, ierr
249      !!----------------------------------------------------------------------
250      IF(lwp) WRITE(numout,*)
251      IF(lwp) WRITE(numout,*) 'ice_init: Arrays allocation & Initialization off all routines & init state' 
252      IF(lwp) WRITE(numout,*) '~~~~~~~~'
253      !
254      !                                ! Open the reference and configuration namelist files and namelist output file
255      CALL ctl_opn( numnam_ice_ref, 'namelist_ice_ref',    'OLD',     'FORMATTED', 'SEQUENTIAL', -1, numout, lwp )
256      CALL ctl_opn( numnam_ice_cfg, 'namelist_ice_cfg',    'OLD',     'FORMATTED', 'SEQUENTIAL', -1, numout, lwp )
257      IF(lwm) CALL ctl_opn( numoni, 'output.namelist.ice', 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, 1 )
258      !
259      CALL ice_run_init                ! set some ice run parameters
260      !
261      !                                ! Allocate the ice arrays (sbc_ice already allocated in sbc_init)
262      ierr =        ice_alloc        ()      ! ice variables
263      ierr = ierr + sbc_ice_alloc    ()      ! surface forcing
264      ierr = ierr + ice1D_alloc      ()      ! thermodynamics
265      ierr = ierr + ice_rdgrft_alloc ()      ! ridging/rafting
266      !
267      IF( lk_mpp    )   CALL mpp_sum( ierr )
268      IF( ierr /= 0 )   CALL ctl_stop('STOP', 'ice_init : unable to allocate ice arrays')
269      !
270      CALL ice_itd_init                ! ice thickness distribution initialization
271      !
272      CALL ice_rdgrft_init             ! set ice ridging/rafting parameters
273      !
274      CALL ice_rhg_init                ! set ice rheology parameters
275      !
276      CALL ice_adv_init                ! set ice advection parameters
277      !
278      CALL ice_thd_init                ! set ice thermodynics parameters
279      !
280      CALL ice_thd_sal_init            ! set ice salinity parameters
281       
282      ! MV MP 2016
283      CALL lim_mp_init                 ! set melt ponds parameters
284      ! END MV MP 2016
285      !                                ! Initial sea-ice state
286      IF( .NOT. ln_rstart ) THEN              ! start from rest: sea-ice deduced from sst
287         CALL ice_istate_init
288         CALL ice_istate
289      ELSE                                    ! start from a restart file
290         CALL ice_rst_read
291      ENDIF
292      CALL ice_var_agg(2)
293      CALL ice_var_glo2eqv
294      !
295      CALL ice_update_init                 ! ice surface boundary condition
296      !
297      CALL ice_alb_init                    ! ice surface albedo
298      !
299      CALL ice_dia_init                    ! initialization for diags
300      !
301      fr_i  (:,:)   = at_i(:,:)         ! initialisation of sea-ice fraction
302      tn_ice(:,:,:) = t_su(:,:,:)       ! initialisation of surface temp for coupled simu
303      !
304      WHERE( gphit(:,:) > 0._wp )   ;   rn_amax_2d(:,:) = rn_amax_n  ! NH
305      ELSEWHERE                     ;   rn_amax_2d(:,:) = rn_amax_s  ! SH
306      END WHERE
307      !
308   END SUBROUTINE ice_init
309
310
311   SUBROUTINE ice_run_init
312      !!-------------------------------------------------------------------
313      !!                  ***  ROUTINE ice_run_init ***
314      !!
315      !! ** Purpose :   Definition some run parameter for ice model
316      !!
317      !! ** Method  :   Read the namice_run namelist and check the parameter
318      !!              values called at the first timestep (nit000)
319      !!
320      !! ** input   :   Namelist namice_run
321      !!-------------------------------------------------------------------
322      INTEGER  ::   ios                 ! Local integer output status for namelist read
323      NAMELIST/namice_run/ jpl, nlay_i, nlay_s, nn_monocat, rn_amax_n, rn_amax_s, cn_icerst_in, cn_icerst_indir,   &
324         &                cn_icerst_out, cn_icerst_outdir
325      !!-------------------------------------------------------------------
326      !
327      REWIND( numnam_ice_ref )      ! Namelist namice_run in reference namelist : Parameters for ice
328      READ  ( numnam_ice_ref, namice_run, IOSTAT = ios, ERR = 901)
329901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namice_run in reference namelist', lwp )
330
331      REWIND( numnam_ice_cfg )      ! Namelist namice_run in configuration namelist : Parameters for ice
332      READ  ( numnam_ice_cfg, namice_run, IOSTAT = ios, ERR = 902 )
333902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namice_run in configuration namelist', lwp )
334      IF(lwm) WRITE ( numoni, namice_run )
335      !
336      IF(lwp) THEN                  ! control print
337         WRITE(numout,*)
338         WRITE(numout,*) 'ice_run_init : ice share parameters for dynamics/advection/thermo of sea-ice'
339         WRITE(numout,*) ' ~~~~~~'
340         WRITE(numout,*) '   Namelist namice_run : '
341         WRITE(numout,*) '      number of ice  categories                              jpl    = ', jpl
342         WRITE(numout,*) '      number of ice  layers                                  nlay_i = ', nlay_i
343         WRITE(numout,*) '      number of snow layers                                  nlay_s = ', nlay_s
344         WRITE(numout,*) '      virtual ITD mono-category param (1-4) or not (0)   nn_monocat = ', nn_monocat
345         WRITE(numout,*) '      maximum ice concentration for NH                              = ', rn_amax_n 
346         WRITE(numout,*) '      maximum ice concentration for SH                              = ', rn_amax_s
347      ENDIF
348      !
349      !                                         !--- check consistency
350      IF ( jpl > 1 .AND. nn_monocat == 1 ) THEN
351         nn_monocat = 0
352         IF(lwp) WRITE(numout,*)
353         IF(lwp) WRITE(numout,*) '   nn_monocat forced to 0 as jpl>1, i.e. multi-category case is chosen'
354      ENDIF
355      IF ( jpl == 1 .AND. nn_monocat == 0 ) THEN
356         CALL ctl_stop( 'STOP', 'ice_run_init : if jpl=1 then nn_monocat should be between 1 and 4' )
357      ENDIF
358      !
359      IF( ln_bdy .AND. ln_icediachk )   CALL ctl_warn('ice_run_init: online conservation check does not work with BDY')
360      !
361      rdt_ice   = REAL(nn_fsbc) * rdt           !--- sea-ice timestep and inverse
362      r1_rdtice = 1._wp / rdt_ice
363      IF( lwp ) WRITE(numout,*) '   ice timestep rdt_ice  = ', rdt_ice
364      !
365      r1_nlay_i = 1._wp / REAL( nlay_i, wp )    !--- inverse of nlay_i and nlay_s
366      r1_nlay_s = 1._wp / REAL( nlay_s, wp )
367      !
368   END SUBROUTINE ice_run_init
369
370
371   SUBROUTINE store_fields
372      !!----------------------------------------------------------------------
373      !!                  ***  ROUTINE store_fields  ***
374      !!
375      !! ** purpose :  store ice variables at "before" time step
376      !!----------------------------------------------------------------------
377      INTEGER  ::   ji, jj, jl      ! dummy loop index
378      !!----------------------------------------------------------------------
379      !
380      a_i_b  (:,:,:)   = a_i  (:,:,:)     ! ice area
381      v_i_b  (:,:,:)   = v_i  (:,:,:)     ! ice volume
382      v_s_b  (:,:,:)   = v_s  (:,:,:)     ! snow volume
383      smv_i_b(:,:,:)   = smv_i(:,:,:)     ! salt content
384      oa_i_b (:,:,:)   = oa_i (:,:,:)     ! areal age content
385      e_s_b  (:,:,:,:) = e_s  (:,:,:,:)   ! snow thermal energy
386      e_i_b  (:,:,:,:) = e_i  (:,:,:,:)   ! ice thermal energy
387      WHERE( a_i_b(:,:,:) >= epsi20 )
388         ht_i_b(:,:,:) = v_i_b (:,:,:) / a_i_b(:,:,:)   ! ice thickness
389         ht_s_b(:,:,:) = v_s_b (:,:,:) / a_i_b(:,:,:)   ! snw thickness
390      ELSEWHERE
391         ht_i_b(:,:,:) = 0._wp
392         ht_s_b(:,:,:) = 0._wp
393      END WHERE
394     
395      ! ice velocities & total concentration
396      at_i_b(:,:)  = SUM( a_i_b(:,:,:), dim=3 )
397      u_ice_b(:,:) = u_ice(:,:)
398      v_ice_b(:,:) = v_ice(:,:)
399      !
400   END SUBROUTINE store_fields
401
402
403   SUBROUTINE ice_diag0
404      !!----------------------------------------------------------------------
405      !!                  ***  ROUTINE ice_diag0  ***
406      !!
407      !! ** purpose :  set ice-ocean and ice-atm. fluxes to zeros at the beggining
408      !!               of the time step
409      !!----------------------------------------------------------------------
410      INTEGER  ::   ji, jj      ! dummy loop index
411      !!----------------------------------------------------------------------
412      sfx    (:,:) = 0._wp   ;
413      sfx_bri(:,:) = 0._wp   ;   sfx_lam(:,:) = 0._wp
414      sfx_sni(:,:) = 0._wp   ;   sfx_opw(:,:) = 0._wp
415      sfx_bog(:,:) = 0._wp   ;   sfx_dyn(:,:) = 0._wp
416      sfx_bom(:,:) = 0._wp   ;   sfx_sum(:,:) = 0._wp
417      sfx_res(:,:) = 0._wp   ;   sfx_sub(:,:) = 0._wp
418      !
419      wfx_snw(:,:) = 0._wp   ;   wfx_ice(:,:) = 0._wp
420      wfx_sni(:,:) = 0._wp   ;   wfx_opw(:,:) = 0._wp
421      wfx_bog(:,:) = 0._wp   ;   wfx_dyn(:,:) = 0._wp
422      wfx_bom(:,:) = 0._wp   ;   wfx_sum(:,:) = 0._wp
423      wfx_res(:,:) = 0._wp   ;   wfx_sub(:,:) = 0._wp
424      wfx_spr(:,:) = 0._wp   ;   wfx_lam(:,:) = 0._wp 
425      wfx_snw_dyn(:,:) = 0._wp ; wfx_snw_sum(:,:) = 0._wp
426      wfx_snw_sub(:,:) = 0._wp ; wfx_ice_sub(:,:) = 0._wp
427      wfx_snw_sni(:,:) = 0._wp 
428      ! MV MP 2016
429      wfx_pnd(:,:) = 0._wp
430      ! END MV MP 2016
431
432      hfx_thd(:,:) = 0._wp   ;
433      hfx_snw(:,:) = 0._wp   ;   hfx_opw(:,:) = 0._wp
434      hfx_bog(:,:) = 0._wp   ;   hfx_dyn(:,:) = 0._wp
435      hfx_bom(:,:) = 0._wp   ;   hfx_sum(:,:) = 0._wp
436      hfx_res(:,:) = 0._wp   ;   hfx_sub(:,:) = 0._wp
437      hfx_spr(:,:) = 0._wp   ;   hfx_dif(:,:) = 0._wp
438      hfx_err(:,:) = 0._wp   ;   hfx_err_rem(:,:) = 0._wp
439      hfx_err_dif(:,:) = 0._wp
440      wfx_err_sub(:,:) = 0._wp
441      !
442      afx_tot(:,:) = 0._wp   ;
443      !
444      diag_heat(:,:) = 0._wp ;   diag_smvi(:,:) = 0._wp
445      diag_vice(:,:) = 0._wp ;   diag_vsnw(:,:) = 0._wp
446
447      ! SIMIP diagnostics
448      diag_fc_bo(:,:)    = 0._wp ; diag_fc_su(:,:)    = 0._wp
449
450      tau_icebfr(:,:) = 0._wp; ! landfast ice param only (clem: important to keep the init here)
451     
452   END SUBROUTINE ice_diag0
453
454#else
455   !!----------------------------------------------------------------------
456   !!   Default option           Dummy module      NO LIM 3.0 sea-ice model
457   !!----------------------------------------------------------------------
458CONTAINS
459   SUBROUTINE ice_stp ( kt, ksbc )     ! Dummy routine
460      INTEGER, INTENT(in) ::   kt, ksbc
461      WRITE(*,*) 'ice_stp: You should not have seen this print! error?', kt
462   END SUBROUTINE ice_stp
463   SUBROUTINE ice_init                 ! Dummy routine
464   END SUBROUTINE ice_init
465#endif
466
467   !!======================================================================
468END MODULE icestp
Note: See TracBrowser for help on using the repository browser.