source: trunk/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90 @ 7753

Last change on this file since 7753 was 7646, checked in by timgraham, 4 years ago

Merge of dev_merge_2016 into trunk. UPDATE TO ARCHFILES NEEDED for XIOS2.
LIM_SRC_s/limrhg.F90 to follow in next commit due to change of kind (I'm unable to do it in this commit).
Merged using the following steps:

1) svn merge —reintegrate svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/trunk .
2) Resolve minor conflicts in sette.sh and namelist_cfg for ORCA2LIM3 (due to a change in trunk after branch was created)
3) svn commit
4) svn switch svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/trunk
5) svn merge svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/branches/2016/dev_merge_2016 .
6) At this stage I checked out a clean copy of the branch to compare against what is about to be committed to the trunk.
6) svn commit #Commit code to the trunk

In this commit I have also reverted a change to Fcheck_archfile.sh which was causing problems on the Paris machine.

  • Property svn:keywords set to Id
File size: 40.5 KB
Line 
1MODULE bdydta
2   !!======================================================================
3   !!                       ***  MODULE bdydta  ***
4   !! Open boundary data : read the data for the unstructured open boundaries.
5   !!======================================================================
6   !! History :  1.0  !  2005-01  (J. Chanut, A. Sellar)  Original code
7   !!             -   !  2007-01  (D. Storkey) Update to use IOM module
8   !!             -   !  2007-07  (D. Storkey) add bdy_dta_fla
9   !!            3.0  !  2008-04  (NEMO team)  add in the reference version
10   !!            3.3  !  2010-09  (E.O'Dea) modifications for Shelf configurations
11   !!            3.3  !  2010-09  (D.Storkey) add ice boundary conditions
12   !!            3.4  !  2011     (D. Storkey) rewrite in preparation for OBC-BDY merge
13   !!            3.6  !  2012-01  (C. Rousset) add ice boundary conditions for lim3
14   !!----------------------------------------------------------------------
15   !!    bdy_dta        : read external data along open boundaries from file
16   !!    bdy_dta_init   : initialise arrays etc for reading of external data
17   !!----------------------------------------------------------------------
18   USE timing          ! Timing
19   USE oce             ! ocean dynamics and tracers
20   USE dom_oce         ! ocean space and time domain
21   USE phycst          ! physical constants
22   USE bdy_oce         ! ocean open boundary conditions 
23   USE bdytides        ! tidal forcing at boundaries
24   USE fldread         ! read input fields
25   USE iom             ! IOM library
26   USE in_out_manager  ! I/O logical units
27#if defined key_lim2
28   USE ice_2
29#elif defined key_lim3
30   USE ice
31   USE limvar          ! redistribute ice input into categories
32#endif
33   USE sbcapr
34   USE sbctide         ! Tidal forcing or not
35
36   IMPLICIT NONE
37   PRIVATE
38
39   PUBLIC   bdy_dta          ! routine called by step.F90 and dynspg_ts.F90
40   PUBLIC   bdy_dta_init     ! routine called by nemogcm.F90
41
42   INTEGER, ALLOCATABLE, DIMENSION(:)   ::   nb_bdy_fld        ! Number of fields to update for each boundary set.
43   INTEGER                              ::   nb_bdy_fld_sum    ! Total number of fields to update for all boundary sets.
44   LOGICAL,           DIMENSION(jp_bdy) ::   ln_full_vel_array ! =T => full velocities in 3D boundary conditions
45                                                               ! =F => baroclinic velocities in 3D boundary conditions
46!$AGRIF_DO_NOT_TREAT
47   TYPE(FLD), PUBLIC, ALLOCATABLE, DIMENSION(:), TARGET ::   bf        ! structure of input fields (file informations, fields read)
48!$AGRIF_END_DO_NOT_TREAT
49   TYPE(MAP_POINTER), ALLOCATABLE, DIMENSION(:) :: nbmap_ptr   ! array of pointers to nbmap
50
51#if defined key_lim3
52   LOGICAL :: ll_bdylim3                  ! determine whether ice input is lim2 (F) or lim3 (T) type
53   INTEGER :: jfld_hti, jfld_hts, jfld_ai ! indices of ice thickness, snow thickness and concentration in bf structure
54#endif
55
56   !!----------------------------------------------------------------------
57   !! NEMO/OPA 3.3 , NEMO Consortium (2010)
58   !! $Id$
59   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
60   !!----------------------------------------------------------------------
61CONTAINS
62
63      SUBROUTINE bdy_dta( kt, jit, time_offset )
64      !!----------------------------------------------------------------------
65      !!                   ***  SUBROUTINE bdy_dta  ***
66      !!                   
67      !! ** Purpose :   Update external data for open boundary conditions
68      !!
69      !! ** Method  :   Use fldread.F90
70      !!               
71      !!----------------------------------------------------------------------
72      INTEGER, INTENT(in)           ::   kt           ! ocean time-step index
73      INTEGER, INTENT(in), OPTIONAL ::   jit          ! subcycle time-step index (for timesplitting option)
74      INTEGER, INTENT(in), OPTIONAL ::   time_offset  ! time offset in units of timesteps. NB. if jit
75      !                                               ! is present then units = subcycle timesteps.
76      !                                               ! time_offset = 0 => get data at "now" time level
77      !                                               ! time_offset = -1 => get data at "before" time level
78      !                                               ! time_offset = +1 => get data at "after" time level
79      !                                               ! etc.
80      !
81      INTEGER ::  ib_bdy, jfld, jstart, jend, ib, ii, ij, ik, igrd, jl  ! local indices
82      INTEGER,          DIMENSION(jpbgrd) ::   ilen1 
83      INTEGER, POINTER, DIMENSION(:)      ::   nblen, nblenrim  ! short cuts
84      TYPE(OBC_DATA), POINTER             ::   dta              ! short cut
85      !!---------------------------------------------------------------------------
86      !
87      IF( nn_timing == 1 )   CALL timing_start('bdy_dta')
88      !
89      ! Initialise data arrays once for all from initial conditions where required
90      !---------------------------------------------------------------------------
91      IF( kt == nit000 .AND. .NOT.PRESENT(jit) ) THEN
92
93         ! Calculate depth-mean currents
94         !-----------------------------
95         
96         DO ib_bdy = 1, nb_bdy
97            !
98            nblen => idx_bdy(ib_bdy)%nblen
99            nblenrim => idx_bdy(ib_bdy)%nblenrim
100            dta => dta_bdy(ib_bdy)
101
102            IF( nn_dyn2d_dta(ib_bdy) == 0 ) THEN
103               ilen1(:) = nblen(:)
104               IF( dta%ll_ssh ) THEN
105                  igrd = 1
106                  DO ib = 1, ilen1(igrd)
107                     ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
108                     ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
109                     dta_bdy(ib_bdy)%ssh(ib) = sshn(ii,ij) * tmask(ii,ij,1)         
110                  END DO
111               END IF
112               IF( dta%ll_u2d ) THEN
113                  igrd = 2
114                  DO ib = 1, ilen1(igrd)
115                     ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
116                     ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
117                     dta_bdy(ib_bdy)%u2d(ib) = un_b(ii,ij) * umask(ii,ij,1)         
118                  END DO
119               END IF
120               IF( dta%ll_v2d ) THEN
121                  igrd = 3
122                  DO ib = 1, ilen1(igrd)
123                     ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
124                     ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
125                     dta_bdy(ib_bdy)%v2d(ib) = vn_b(ii,ij) * vmask(ii,ij,1)         
126                  END DO
127               END IF
128            ENDIF
129
130            IF( nn_dyn3d_dta(ib_bdy) == 0 ) THEN
131               ilen1(:) = nblen(:)
132               IF( dta%ll_u3d ) THEN
133                  igrd = 2 
134                  DO ib = 1, ilen1(igrd)
135                     DO ik = 1, jpkm1
136                        ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
137                        ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
138                        dta_bdy(ib_bdy)%u3d(ib,ik) =  ( un(ii,ij,ik) - un_b(ii,ij) ) * umask(ii,ij,ik)         
139                     END DO
140                  END DO
141               END IF
142               IF( dta%ll_v3d ) THEN
143                  igrd = 3 
144                  DO ib = 1, ilen1(igrd)
145                     DO ik = 1, jpkm1
146                        ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
147                        ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
148                        dta_bdy(ib_bdy)%v3d(ib,ik) =  ( vn(ii,ij,ik) - vn_b(ii,ij) ) * vmask(ii,ij,ik)         
149                        END DO
150                  END DO
151               END IF
152            ENDIF
153
154            IF( nn_tra_dta(ib_bdy) == 0 ) THEN
155               ilen1(:) = nblen(:)
156               IF( dta%ll_tem ) THEN
157                  igrd = 1 
158                  DO ib = 1, ilen1(igrd)
159                     DO ik = 1, jpkm1
160                        ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
161                        ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
162                        dta_bdy(ib_bdy)%tem(ib,ik) = tsn(ii,ij,ik,jp_tem) * tmask(ii,ij,ik)         
163                     END DO
164                  END DO
165               END IF
166               IF( dta%ll_sal ) THEN
167                  igrd = 1 
168                  DO ib = 1, ilen1(igrd)
169                     DO ik = 1, jpkm1
170                        ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
171                        ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
172                        dta_bdy(ib_bdy)%sal(ib,ik) = tsn(ii,ij,ik,jp_sal) * tmask(ii,ij,ik)         
173                     END DO
174                  END DO
175               END IF
176            ENDIF
177
178#if defined key_lim2
179            IF( nn_ice_lim_dta(ib_bdy) == 0 ) THEN
180               ilen1(:) = nblen(:)
181               IF( dta%ll_frld ) THEN
182                  igrd = 1 
183                  DO ib = 1, ilen1(igrd)
184                     ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
185                     ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
186                     dta_bdy(ib_bdy)%frld(ib) = frld(ii,ij) * tmask(ii,ij,1)         
187                  END DO
188               END IF
189               IF( dta%ll_hicif ) THEN
190                  igrd = 1 
191                  DO ib = 1, ilen1(igrd)
192                     ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
193                     ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
194                     dta_bdy(ib_bdy)%hicif(ib) = hicif(ii,ij) * tmask(ii,ij,1)         
195                  END DO
196               END IF
197               IF( dta%ll_hsnif ) THEN
198                  igrd = 1 
199                  DO ib = 1, ilen1(igrd)
200                     ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
201                     ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
202                     dta_bdy(ib_bdy)%hsnif(ib) = hsnif(ii,ij) * tmask(ii,ij,1)         
203                  END DO
204               END IF
205            ENDIF
206#elif defined key_lim3
207            IF( nn_ice_lim_dta(ib_bdy) == 0 ) THEN
208               ilen1(:) = nblen(:)
209               IF( dta%ll_a_i ) THEN
210                  igrd = 1   
211                  DO jl = 1, jpl
212                     DO ib = 1, ilen1(igrd)
213                        ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
214                        ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
215                        dta_bdy(ib_bdy)%a_i (ib,jl) =  a_i(ii,ij,jl) * tmask(ii,ij,1) 
216                     END DO
217                  END DO
218               ENDIF
219               IF( dta%ll_ht_i ) THEN
220                  igrd = 1   
221                  DO jl = 1, jpl
222                     DO ib = 1, ilen1(igrd)
223                        ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
224                        ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
225                        dta_bdy(ib_bdy)%ht_i (ib,jl) =  ht_i(ii,ij,jl) * tmask(ii,ij,1) 
226                     END DO
227                  END DO
228               ENDIF
229               IF( dta%ll_ht_s ) THEN
230                  igrd = 1   
231                  DO jl = 1, jpl
232                     DO ib = 1, ilen1(igrd)
233                        ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
234                        ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
235                        dta_bdy(ib_bdy)%ht_s (ib,jl) =  ht_s(ii,ij,jl) * tmask(ii,ij,1) 
236                     END DO
237                  END DO
238               ENDIF
239            ENDIF
240#endif
241         END DO ! ib_bdy
242         !
243      ENDIF ! kt == nit000
244
245      ! update external data from files
246      !--------------------------------
247     
248      jstart = 1
249      DO ib_bdy = 1, nb_bdy   
250         dta => dta_bdy(ib_bdy)
251         IF( nn_dta(ib_bdy) == 1 ) THEN ! skip this bit if no external data required
252     
253            IF( PRESENT(jit) ) THEN
254               ! Update barotropic boundary conditions only
255               ! jit is optional argument for fld_read and bdytide_update
256               IF( cn_dyn2d(ib_bdy) /= 'none' ) THEN
257                  IF( nn_dyn2d_dta(ib_bdy) == 2 ) THEN ! tidal harmonic forcing ONLY: initialise arrays
258                     IF( dta%ll_ssh ) dta%ssh(:) = 0._wp
259                     IF( dta%ll_u2d ) dta%u2d(:) = 0._wp
260                     IF( dta%ll_u3d ) dta%v2d(:) = 0._wp
261                  ENDIF
262                  IF (cn_tra(ib_bdy) /= 'runoff') THEN
263                     IF( nn_dyn2d_dta(ib_bdy) == 1 .OR. nn_dyn2d_dta(ib_bdy) == 3 ) THEN
264
265                        jend = jstart + dta%nread(2) - 1
266                        IF( ln_full_vel_array(ib_bdy) ) THEN
267                           CALL fld_read( kt=kt, kn_fsbc=1, sd=bf(jstart:jend), map=nbmap_ptr(jstart:jend),  &
268                                     & kit=jit, kt_offset=time_offset , jpk_bdy=nb_jpk_bdy, fvl=ln_full_vel_array(ib_bdy)  )
269                        ELSE
270                           CALL fld_read( kt=kt, kn_fsbc=1, sd=bf(jstart:jend), map=nbmap_ptr(jstart:jend),  &
271                                     & kit=jit, kt_offset=time_offset  )
272                        ENDIF
273
274                        ! If full velocities in boundary data then extract barotropic velocities from 3D fields
275                        IF( ln_full_vel_array(ib_bdy) .AND.                                             &
276                          &    ( nn_dyn2d_dta(ib_bdy) == 1 .OR. nn_dyn2d_dta(ib_bdy) == 3 .OR.  &
277                          &      nn_dyn3d_dta(ib_bdy) == 1 ) )THEN
278
279                           igrd = 2                      ! zonal velocity
280                           dta%u2d(:) = 0._wp
281                           DO ib = 1, idx_bdy(ib_bdy)%nblen(igrd)
282                              ii   = idx_bdy(ib_bdy)%nbi(ib,igrd)
283                              ij   = idx_bdy(ib_bdy)%nbj(ib,igrd)
284                              DO ik = 1, jpkm1
285                                 dta%u2d(ib) = dta%u2d(ib) &
286                       &                          + e3u_n(ii,ij,ik) * umask(ii,ij,ik) * dta%u3d(ib,ik)
287                              END DO
288                              dta%u2d(ib) =  dta%u2d(ib) * r1_hu_n(ii,ij)
289                           END DO
290                           igrd = 3                      ! meridional velocity
291                           dta%v2d(:) = 0._wp
292                           DO ib = 1, idx_bdy(ib_bdy)%nblen(igrd)
293                              ii   = idx_bdy(ib_bdy)%nbi(ib,igrd)
294                              ij   = idx_bdy(ib_bdy)%nbj(ib,igrd)
295                              DO ik = 1, jpkm1
296                                 dta%v2d(ib) = dta%v2d(ib) &
297                       &                       + e3v_n(ii,ij,ik) * vmask(ii,ij,ik) * dta%v3d(ib,ik)
298                              END DO
299                              dta%v2d(ib) =  dta%v2d(ib) * r1_hv_n(ii,ij)
300                           END DO
301                        ENDIF                   
302                     ENDIF
303                     IF( nn_dyn2d_dta(ib_bdy) .ge. 2 ) THEN ! update tidal harmonic forcing
304                        CALL bdytide_update( kt=kt, idx=idx_bdy(ib_bdy), dta=dta, td=tides(ib_bdy),   & 
305                          &                 jit=jit, time_offset=time_offset )
306                     ENDIF
307                  ENDIF
308               ENDIF
309            ELSE
310               IF (cn_tra(ib_bdy) == 'runoff') then      ! runoff condition
311                  jend = nb_bdy_fld(ib_bdy)
312                  CALL fld_read( kt=kt, kn_fsbc=1, sd=bf(jstart:jend),  &
313                               & map=nbmap_ptr(jstart:jend), kt_offset=time_offset )
314                  !
315                  igrd = 2                      ! zonal velocity
316                  DO ib = 1, idx_bdy(ib_bdy)%nblen(igrd)
317                     ii   = idx_bdy(ib_bdy)%nbi(ib,igrd)
318                     ij   = idx_bdy(ib_bdy)%nbj(ib,igrd)
319                     dta%u2d(ib) = dta%u2d(ib) / ( e2u(ii,ij) * hu_0(ii,ij) )
320                  END DO
321                  !
322                  igrd = 3                      ! meridional velocity
323                  DO ib = 1, idx_bdy(ib_bdy)%nblen(igrd)
324                     ii   = idx_bdy(ib_bdy)%nbi(ib,igrd)
325                     ij   = idx_bdy(ib_bdy)%nbj(ib,igrd)
326                     dta%v2d(ib) = dta%v2d(ib) / ( e1v(ii,ij) * hv_0(ii,ij) )
327                  END DO
328               ELSE
329                  IF( nn_dyn2d_dta(ib_bdy) == 2 ) THEN ! tidal harmonic forcing ONLY: initialise arrays
330                     IF( dta%ll_ssh ) dta%ssh(:) = 0._wp
331                     IF( dta%ll_u2d ) dta%u2d(:) = 0._wp
332                     IF( dta%ll_v2d ) dta%v2d(:) = 0._wp
333                  ENDIF
334                  IF( dta%nread(1) .gt. 0 ) THEN ! update external data
335                     jend = jstart + dta%nread(1) - 1
336                     CALL fld_read( kt=kt, kn_fsbc=1, sd=bf(jstart:jend), &
337                                  & map=nbmap_ptr(jstart:jend), kt_offset=time_offset, jpk_bdy=nb_jpk_bdy, fvl=ln_full_vel_array(ib_bdy) )
338                  ENDIF
339                  ! If full velocities in boundary data then split into barotropic and baroclinic data
340                  IF( ln_full_vel_array(ib_bdy) .and.                                             &
341                    & ( nn_dyn2d_dta(ib_bdy) == 1 .OR. nn_dyn2d_dta(ib_bdy) == 3 .OR. &
342                    &   nn_dyn3d_dta(ib_bdy) == 1 ) ) THEN
343                     igrd = 2                      ! zonal velocity
344                     dta%u2d(:) = 0._wp
345                     DO ib = 1, idx_bdy(ib_bdy)%nblen(igrd)
346                        ii   = idx_bdy(ib_bdy)%nbi(ib,igrd)
347                        ij   = idx_bdy(ib_bdy)%nbj(ib,igrd)
348                        DO ik = 1, jpkm1
349                           dta%u2d(ib) = dta%u2d(ib) &
350                 &                       + e3u_n(ii,ij,ik) * umask(ii,ij,ik) * dta%u3d(ib,ik)
351                        END DO
352                        dta%u2d(ib) =  dta%u2d(ib) * r1_hu_n(ii,ij)
353                        DO ik = 1, jpkm1
354                           dta%u3d(ib,ik) = dta%u3d(ib,ik) - dta%u2d(ib)
355                        END DO
356                     END DO
357                     igrd = 3                      ! meridional velocity
358                     dta%v2d(:) = 0._wp
359                     DO ib = 1, idx_bdy(ib_bdy)%nblen(igrd)
360                        ii   = idx_bdy(ib_bdy)%nbi(ib,igrd)
361                        ij   = idx_bdy(ib_bdy)%nbj(ib,igrd)
362                        DO ik = 1, jpkm1
363                           dta%v2d(ib) = dta%v2d(ib) &
364                 &                       + e3v_n(ii,ij,ik) * vmask(ii,ij,ik) * dta%v3d(ib,ik)
365                        END DO
366                        dta%v2d(ib) =  dta%v2d(ib) * r1_hv_n(ii,ij)
367                        DO ik = 1, jpkm1
368                           dta%v3d(ib,ik) = dta%v3d(ib,ik) - dta%v2d(ib)
369                        END DO
370                     END DO
371                  ENDIF
372
373               ENDIF
374#if defined key_lim3
375               IF( .NOT. ll_bdylim3 .AND. cn_ice_lim(ib_bdy) /= 'none' .AND. nn_ice_lim_dta(ib_bdy) == 1 ) THEN ! bdy ice input (case input is lim2 type)
376                CALL lim_var_itd ( bf(jfld_hti)%fnow(:,1,1), bf(jfld_hts)%fnow(:,1,1), bf(jfld_ai)%fnow(:,1,1), &
377                                  & dta_bdy(ib_bdy)%ht_i,     dta_bdy(ib_bdy)%ht_s,     dta_bdy(ib_bdy)%a_i     )
378               ENDIF
379#endif
380            ENDIF
381            jstart = jstart + dta%nread(1)
382         END IF ! nn_dta(ib_bdy) = 1
383      END DO  ! ib_bdy
384
385      IF ( ln_tide ) THEN
386         IF (ln_dynspg_ts) THEN      ! Fill temporary arrays with slow-varying bdy data                           
387            DO ib_bdy = 1, nb_bdy    ! Tidal component added in ts loop
388               IF ( nn_dyn2d_dta(ib_bdy) .ge. 2 ) THEN
389                  nblen => idx_bdy(ib_bdy)%nblen
390                  nblenrim => idx_bdy(ib_bdy)%nblenrim
391                  IF( cn_dyn2d(ib_bdy) == 'frs' ) THEN; ilen1(:)=nblen(:) ; ELSE ; ilen1(:)=nblenrim(:) ; ENDIF
392                  IF ( dta_bdy(ib_bdy)%ll_ssh ) dta_bdy_s(ib_bdy)%ssh(1:ilen1(1)) = dta_bdy(ib_bdy)%ssh(1:ilen1(1))
393                  IF ( dta_bdy(ib_bdy)%ll_u2d ) dta_bdy_s(ib_bdy)%u2d(1:ilen1(2)) = dta_bdy(ib_bdy)%u2d(1:ilen1(2))
394                  IF ( dta_bdy(ib_bdy)%ll_v2d ) dta_bdy_s(ib_bdy)%v2d(1:ilen1(3)) = dta_bdy(ib_bdy)%v2d(1:ilen1(3))
395               ENDIF
396            END DO
397         ELSE ! Add tides if not split-explicit free surface else this is done in ts loop
398            !
399            CALL bdy_dta_tides( kt=kt, time_offset=time_offset )
400         ENDIF
401      ENDIF
402
403      IF ( ln_apr_obc ) THEN
404         DO ib_bdy = 1, nb_bdy
405            IF (cn_tra(ib_bdy) /= 'runoff')THEN
406               igrd = 1                      ! meridional velocity
407               DO ib = 1, idx_bdy(ib_bdy)%nblenrim(igrd)
408                  ii   = idx_bdy(ib_bdy)%nbi(ib,igrd)
409                  ij   = idx_bdy(ib_bdy)%nbj(ib,igrd)
410                  dta_bdy(ib_bdy)%ssh(ib) = dta_bdy(ib_bdy)%ssh(ib) + ssh_ib(ii,ij)
411               END DO
412            ENDIF
413         END DO
414      ENDIF
415      !
416      IF( nn_timing == 1 ) CALL timing_stop('bdy_dta')
417      !
418   END SUBROUTINE bdy_dta
419
420
421   SUBROUTINE bdy_dta_init
422      !!----------------------------------------------------------------------
423      !!                   ***  SUBROUTINE bdy_dta_init  ***
424      !!                   
425      !! ** Purpose :   Initialise arrays for reading of external data
426      !!                for open boundary conditions
427      !!
428      !! ** Method  :   
429      !!               
430      !!----------------------------------------------------------------------
431      INTEGER ::   ib_bdy, jfld, jstart, jend, ierror, ios     ! Local integers
432      !
433      CHARACTER(len=100)                     ::   cn_dir        ! Root directory for location of data files
434      CHARACTER(len=100), DIMENSION(nb_bdy)  ::   cn_dir_array  ! Root directory for location of data files
435      CHARACTER(len = 256)::   clname                           ! temporary file name
436      LOGICAL                                ::   ln_full_vel   ! =T => full velocities in 3D boundary data
437                                                                ! =F => baroclinic velocities in 3D boundary data
438      INTEGER                                ::   ilen_global   ! Max length required for global bdy dta arrays
439      INTEGER, ALLOCATABLE, DIMENSION(:)     ::   ilen1, ilen3  ! size of 1st and 3rd dimensions of local arrays
440      INTEGER, ALLOCATABLE, DIMENSION(:)     ::   ibdy           ! bdy set for a particular jfld
441      INTEGER, ALLOCATABLE, DIMENSION(:)     ::   igrid         ! index for grid type (1,2,3 = T,U,V)
442      INTEGER, POINTER, DIMENSION(:)         ::   nblen, nblenrim  ! short cuts
443      TYPE(OBC_DATA), POINTER                ::   dta           ! short cut
444#if defined key_lim3
445      INTEGER, DIMENSION(3) ::   zdimsz   ! number of elements in each of the 4 dimensions (i.e. i,j,t,ice-cat) for an array
446      INTEGER               ::   zndims   ! number of dimensions in an array (i.e. 3 = wo ice cat; 4 = w ice cat)
447      INTEGER               ::   inum,id1 ! local integer
448#endif
449      TYPE(FLD_N), ALLOCATABLE, DIMENSION(:) ::   blf_i         !  array of namelist information structures
450      TYPE(FLD_N) ::   bn_tem, bn_sal, bn_u3d, bn_v3d   !
451      TYPE(FLD_N) ::   bn_ssh, bn_u2d, bn_v2d           ! informations about the fields to be read
452#if defined key_lim2
453      TYPE(FLD_N) ::   bn_frld, bn_hicif, bn_hsnif      !
454#elif defined key_lim3
455      TYPE(FLD_N) ::   bn_a_i, bn_ht_i, bn_ht_s     
456#endif
457      NAMELIST/nambdy_dta/ cn_dir, bn_tem, bn_sal, bn_u3d, bn_v3d, bn_ssh, bn_u2d, bn_v2d 
458#if defined key_lim2
459      NAMELIST/nambdy_dta/ bn_frld, bn_hicif, bn_hsnif
460#elif defined key_lim3
461      NAMELIST/nambdy_dta/ bn_a_i, bn_ht_i, bn_ht_s
462#endif
463      NAMELIST/nambdy_dta/ ln_full_vel, nb_jpk_bdy
464      !!---------------------------------------------------------------------------
465      !
466      IF( nn_timing == 1 )   CALL timing_start('bdy_dta_init')
467      !
468      IF(lwp) WRITE(numout,*)
469      IF(lwp) WRITE(numout,*) 'bdy_dta_ini : initialization of data at the open boundaries'
470      IF(lwp) WRITE(numout,*) '~~~~~~~~~~'
471      IF(lwp) WRITE(numout,*) ''
472
473      ! Set nn_dta
474      DO ib_bdy = 1, nb_bdy
475         nn_dta(ib_bdy) = MAX(  nn_dyn2d_dta(ib_bdy)       &
476                               ,nn_dyn3d_dta(ib_bdy)       &
477                               ,nn_tra_dta(ib_bdy)         &
478#if ( defined key_lim2 || defined key_lim3 )
479                              ,nn_ice_lim_dta(ib_bdy)    &
480#endif
481                              )
482         IF(nn_dta(ib_bdy) > 1) nn_dta(ib_bdy) = 1
483      END DO
484
485      ! Work out upper bound of how many fields there are to read in and allocate arrays
486      ! ---------------------------------------------------------------------------
487      ALLOCATE( nb_bdy_fld(nb_bdy) )
488      nb_bdy_fld(:) = 0
489      DO ib_bdy = 1, nb_bdy         
490         IF( cn_dyn2d(ib_bdy) /= 'none' .and. ( nn_dyn2d_dta(ib_bdy) == 1 .or. nn_dyn2d_dta(ib_bdy) == 3 ) ) THEN
491            nb_bdy_fld(ib_bdy) = nb_bdy_fld(ib_bdy) + 3
492         ENDIF
493         IF( cn_dyn3d(ib_bdy) /= 'none' .and. nn_dyn3d_dta(ib_bdy) == 1 ) THEN
494            nb_bdy_fld(ib_bdy) = nb_bdy_fld(ib_bdy) + 2
495         ENDIF
496         IF( cn_tra(ib_bdy) /= 'none' .and. nn_tra_dta(ib_bdy) == 1  ) THEN
497            nb_bdy_fld(ib_bdy) = nb_bdy_fld(ib_bdy) + 2
498         ENDIF
499#if ( defined key_lim2 || defined key_lim3 )
500         IF( cn_ice_lim(ib_bdy) /= 'none' .and. nn_ice_lim_dta(ib_bdy) == 1  ) THEN
501            nb_bdy_fld(ib_bdy) = nb_bdy_fld(ib_bdy) + 3
502         ENDIF
503#endif               
504         IF(lwp) WRITE(numout,*) 'Maximum number of files to open =',nb_bdy_fld(ib_bdy)
505      END DO           
506
507      nb_bdy_fld_sum = SUM( nb_bdy_fld )
508
509      ALLOCATE( bf(nb_bdy_fld_sum), STAT=ierror )
510      IF( ierror > 0 ) THEN   
511         CALL ctl_stop( 'bdy_dta: unable to allocate bf structure' )   ;   RETURN 
512      ENDIF
513      ALLOCATE( blf_i(nb_bdy_fld_sum), STAT=ierror )
514      IF( ierror > 0 ) THEN   
515         CALL ctl_stop( 'bdy_dta: unable to allocate blf_i structure' )   ;   RETURN 
516      ENDIF
517      ALLOCATE( nbmap_ptr(nb_bdy_fld_sum), STAT=ierror )
518      IF( ierror > 0 ) THEN   
519         CALL ctl_stop( 'bdy_dta: unable to allocate nbmap_ptr structure' )   ;   RETURN 
520      ENDIF
521      ALLOCATE( ilen1(nb_bdy_fld_sum), ilen3(nb_bdy_fld_sum) ) 
522      ALLOCATE( ibdy(nb_bdy_fld_sum) ) 
523      ALLOCATE( igrid(nb_bdy_fld_sum) ) 
524
525      ! Read namelists
526      ! --------------
527      REWIND(numnam_ref)
528      REWIND(numnam_cfg)
529      jfld = 0 
530      DO ib_bdy = 1, nb_bdy         
531         IF( nn_dta(ib_bdy) == 1 ) THEN
532            READ  ( numnam_ref, nambdy_dta, IOSTAT = ios, ERR = 901)
533901         IF( ios /= 0 )   CALL ctl_nam ( ios , 'nambdy_dta in reference namelist', lwp )
534
535            READ  ( numnam_cfg, nambdy_dta, IOSTAT = ios, ERR = 902 )
536902         IF( ios /= 0 )   CALL ctl_nam ( ios , 'nambdy_dta in configuration namelist', lwp )
537            IF(lwm) WRITE( numond, nambdy_dta )
538
539            cn_dir_array(ib_bdy) = cn_dir
540            ln_full_vel_array(ib_bdy) = ln_full_vel
541
542            nblen => idx_bdy(ib_bdy)%nblen
543            nblenrim => idx_bdy(ib_bdy)%nblenrim
544            dta => dta_bdy(ib_bdy)
545            dta%nread(2) = 0
546
547            ! Only read in necessary fields for this set.
548            ! Important that barotropic variables come first.
549            IF( nn_dyn2d_dta(ib_bdy) == 1 .or. nn_dyn2d_dta(ib_bdy) == 3 ) THEN
550
551               IF( dta%ll_ssh ) THEN
552                  if(lwp) write(numout,*) '++++++ reading in ssh field'
553                  jfld = jfld + 1
554                  blf_i(jfld) = bn_ssh
555                  ibdy(jfld) = ib_bdy
556                  igrid(jfld) = 1
557                  ilen1(jfld) = nblen(igrid(jfld))
558                  ilen3(jfld) = 1
559                  dta%nread(2) = dta%nread(2) + 1
560               ENDIF
561
562               IF( dta%ll_u2d .and. .not. ln_full_vel_array(ib_bdy) ) THEN
563                  if(lwp) write(numout,*) '++++++ reading in u2d field'
564                  jfld = jfld + 1
565                  blf_i(jfld) = bn_u2d
566                  ibdy(jfld) = ib_bdy
567                  igrid(jfld) = 2
568                  ilen1(jfld) = nblen(igrid(jfld))
569                  ilen3(jfld) = 1
570                  dta%nread(2) = dta%nread(2) + 1
571               ENDIF
572
573               IF( dta%ll_v2d .and. .not. ln_full_vel_array(ib_bdy) ) THEN
574                  if(lwp) write(numout,*) '++++++ reading in v2d field'
575                  jfld = jfld + 1
576                  blf_i(jfld) = bn_v2d
577                  ibdy(jfld) = ib_bdy
578                  igrid(jfld) = 3
579                  ilen1(jfld) = nblen(igrid(jfld))
580                  ilen3(jfld) = 1
581                  dta%nread(2) = dta%nread(2) + 1
582               ENDIF
583
584            ENDIF
585
586            ! read 3D velocities if baroclinic velocities require OR if
587            ! barotropic velocities required and ln_full_vel set to .true.
588            IF( nn_dyn3d_dta(ib_bdy) == 1 .OR. &
589           &  ( ln_full_vel_array(ib_bdy) .AND. ( nn_dyn2d_dta(ib_bdy) == 1 .or. nn_dyn2d_dta(ib_bdy) == 3 ) ) ) THEN
590
591               IF( dta%ll_u3d .OR. ( ln_full_vel_array(ib_bdy) .and. dta%ll_u2d ) ) THEN
592                  if(lwp) write(numout,*) '++++++ reading in u3d field'
593                  jfld = jfld + 1
594                  blf_i(jfld) = bn_u3d
595                  ibdy(jfld) = ib_bdy
596                  igrid(jfld) = 2
597                  ilen1(jfld) = nblen(igrid(jfld))
598                  ilen3(jfld) = jpk
599                  IF( ln_full_vel_array(ib_bdy) .and. dta%ll_u2d ) dta%nread(2) = dta%nread(2) + 1
600               ENDIF
601
602               IF( dta%ll_v3d .OR. ( ln_full_vel_array(ib_bdy) .and. dta%ll_v2d ) ) THEN
603                  if(lwp) write(numout,*) '++++++ reading in v3d field'
604                  jfld = jfld + 1
605                  blf_i(jfld) = bn_v3d
606                  ibdy(jfld) = ib_bdy
607                  igrid(jfld) = 3
608                  ilen1(jfld) = nblen(igrid(jfld))
609                  ilen3(jfld) = jpk
610                  IF( ln_full_vel_array(ib_bdy) .and. dta%ll_v2d ) dta%nread(2) = dta%nread(2) + 1
611               ENDIF
612
613            ENDIF
614
615            ! temperature and salinity
616            IF( nn_tra_dta(ib_bdy) == 1 ) THEN
617
618               IF( dta%ll_tem ) THEN
619                  if(lwp) write(numout,*) '++++++ reading in tem field'
620                  jfld = jfld + 1
621                  blf_i(jfld) = bn_tem
622                  ibdy(jfld) = ib_bdy
623                  igrid(jfld) = 1
624                  ilen1(jfld) = nblen(igrid(jfld))
625                  ilen3(jfld) = jpk
626               ENDIF
627
628               IF( dta%ll_sal ) THEN
629                  if(lwp) write(numout,*) '++++++ reading in sal field'
630                  jfld = jfld + 1
631                  blf_i(jfld) = bn_sal
632                  ibdy(jfld) = ib_bdy
633                  igrid(jfld) = 1
634                  ilen1(jfld) = nblen(igrid(jfld))
635                  ilen3(jfld) = jpk
636               ENDIF
637
638            ENDIF
639
640#if defined key_lim2
641            ! sea ice
642            IF( nn_ice_lim_dta(ib_bdy) == 1 ) THEN
643
644               IF( dta%ll_frld ) THEN
645                  jfld = jfld + 1
646                  blf_i(jfld) = bn_frld
647                  ibdy(jfld) = ib_bdy
648                  igrid(jfld) = 1
649                  ilen1(jfld) = nblen(igrid(jfld))
650                  ilen3(jfld) = 1
651               ENDIF
652
653               IF( dta%ll_hicif ) THEN
654                  jfld = jfld + 1
655                  blf_i(jfld) = bn_hicif
656                  ibdy(jfld) = ib_bdy
657                  igrid(jfld) = 1
658                  ilen1(jfld) = nblen(igrid(jfld))
659                  ilen3(jfld) = 1
660               ENDIF
661
662               IF( dta%ll_hsnif ) THEN
663                  jfld = jfld + 1
664                  blf_i(jfld) = bn_hsnif
665                  ibdy(jfld) = ib_bdy
666                  igrid(jfld) = 1
667                  ilen1(jfld) = nblen(igrid(jfld))
668                  ilen3(jfld) = 1
669               ENDIF
670
671            ENDIF
672#elif defined key_lim3
673            ! sea ice
674            IF( nn_ice_lim_dta(ib_bdy) == 1 ) THEN
675               ! Test for types of ice input (lim2 or lim3)
676               ! Build file name to find dimensions
677               clname=TRIM(bn_a_i%clname)
678               IF( .NOT. bn_a_i%ln_clim ) THEN   
679                                                  WRITE(clname, '(a,"_y",i4.4)' ) TRIM( bn_a_i%clname ), nyear    ! add year
680                  IF( bn_a_i%cltype /= 'yearly' ) WRITE(clname, '(a,"m" ,i2.2)' ) TRIM( clname        ), nmonth   ! add month
681               ELSE
682                  IF( bn_a_i%cltype /= 'yearly' ) WRITE(clname, '(a,"_m",i2.2)' ) TRIM( bn_a_i%clname ), nmonth   ! add month
683               ENDIF
684               IF( bn_a_i%cltype == 'daily' .OR. bn_a_i%cltype(1:4) == 'week' ) &
685               &                                  WRITE(clname, '(a,"d" ,i2.2)' ) TRIM( clname        ), nday     ! add day
686               !
687               CALL iom_open  ( clname, inum )
688               id1 = iom_varid( inum, bn_a_i%clvar, kdimsz=zdimsz, kndims=zndims, ldstop = .FALSE. )
689               CALL iom_close ( inum )
690
691                IF ( zndims == 4 ) THEN
692                 ll_bdylim3 = .TRUE.   ! lim3 input
693               ELSE
694                 ll_bdylim3 = .FALSE.  ! lim2 input     
695               ENDIF
696               ! End test
697
698               IF( dta%ll_a_i ) THEN
699                  jfld = jfld + 1
700                  blf_i(jfld) = bn_a_i
701                  ibdy(jfld) = ib_bdy
702                  igrid(jfld) = 1
703                  ilen1(jfld) = nblen(igrid(jfld))
704                  IF ( ll_bdylim3 ) THEN ; ilen3(jfld)=jpl ; ELSE ; ilen3(jfld)=1 ; ENDIF
705               ENDIF
706
707               IF( dta%ll_ht_i ) THEN
708                  jfld = jfld + 1
709                  blf_i(jfld) = bn_ht_i
710                  ibdy(jfld) = ib_bdy
711                  igrid(jfld) = 1
712                  ilen1(jfld) = nblen(igrid(jfld))
713                  IF ( ll_bdylim3 ) THEN ; ilen3(jfld)=jpl ; ELSE ; ilen3(jfld)=1 ; ENDIF
714               ENDIF
715
716               IF( dta%ll_ht_s ) THEN
717                  jfld = jfld + 1
718                   blf_i(jfld) = bn_ht_s
719                  ibdy(jfld) = ib_bdy
720                  igrid(jfld) = 1
721                  ilen1(jfld) = nblen(igrid(jfld))
722                  IF ( ll_bdylim3 ) THEN ; ilen3(jfld)=jpl ; ELSE ; ilen3(jfld)=1 ; ENDIF
723               ENDIF
724
725            ENDIF
726#endif
727            ! Recalculate field counts
728            !-------------------------
729            IF( ib_bdy == 1 ) THEN
730               nb_bdy_fld_sum = 0
731               nb_bdy_fld(ib_bdy) = jfld
732               nb_bdy_fld_sum     = jfld             
733            ELSE
734               nb_bdy_fld(ib_bdy) = jfld - nb_bdy_fld_sum
735               nb_bdy_fld_sum = nb_bdy_fld_sum + nb_bdy_fld(ib_bdy)
736            ENDIF
737
738            dta%nread(1) = nb_bdy_fld(ib_bdy)
739
740         ENDIF ! nn_dta == 1
741      ENDDO ! ib_bdy
742
743      DO jfld = 1, nb_bdy_fld_sum
744         ALLOCATE( bf(jfld)%fnow(ilen1(jfld),1,ilen3(jfld)) )
745         IF( blf_i(jfld)%ln_tint ) ALLOCATE( bf(jfld)%fdta(ilen1(jfld),1,ilen3(jfld),2) )
746         nbmap_ptr(jfld)%ptr => idx_bdy(ibdy(jfld))%nbmap(:,igrid(jfld))
747         nbmap_ptr(jfld)%ll_unstruc = ln_coords_file(ibdy(jfld))
748      ENDDO
749
750      ! fill bf with blf_i and control print
751      !-------------------------------------
752      jstart = 1
753      DO ib_bdy = 1, nb_bdy
754         jend = jstart - 1 + nb_bdy_fld(ib_bdy) 
755         CALL fld_fill( bf(jstart:jend), blf_i(jstart:jend), cn_dir_array(ib_bdy), 'bdy_dta',   &
756         &              'open boundary conditions', 'nambdy_dta' )
757         jstart = jend + 1
758      ENDDO
759
760      ! Initialise local boundary data arrays
761      ! nn_xxx_dta=0 : allocate space - will be filled from initial conditions later
762      ! nn_xxx_dta=1 : point to "fnow" arrays
763      !-------------------------------------
764
765      jfld = 0
766      DO ib_bdy=1, nb_bdy
767
768         nblen => idx_bdy(ib_bdy)%nblen
769         dta => dta_bdy(ib_bdy)
770
771         if(lwp) then
772            write(numout,*) '++++++ dta%ll_ssh = ',dta%ll_ssh
773            write(numout,*) '++++++ dta%ll_u2d = ',dta%ll_u2d
774            write(numout,*) '++++++ dta%ll_v2d = ',dta%ll_v2d
775            write(numout,*) '++++++ dta%ll_u3d = ',dta%ll_u3d
776            write(numout,*) '++++++ dta%ll_v3d = ',dta%ll_v3d
777            write(numout,*) '++++++ dta%ll_tem = ',dta%ll_tem
778            write(numout,*) '++++++ dta%ll_sal = ',dta%ll_sal
779         endif
780
781         IF ( nn_dyn2d_dta(ib_bdy) == 0 .or. nn_dyn2d_dta(ib_bdy) == 2 ) THEN
782            if(lwp) write(numout,*) '++++++ dta%ssh/u2d/u3d allocated space'
783            IF( dta%ll_ssh ) ALLOCATE( dta%ssh(nblen(1)) )
784            IF( dta%ll_u2d ) ALLOCATE( dta%u2d(nblen(2)) )
785            IF( dta%ll_v2d ) ALLOCATE( dta%v2d(nblen(3)) )
786         ENDIF
787         IF ( nn_dyn2d_dta(ib_bdy) == 1 .or. nn_dyn2d_dta(ib_bdy) == 3 ) THEN
788            IF( dta%ll_ssh ) THEN
789               if(lwp) write(numout,*) '++++++ dta%ssh pointing to fnow'
790               jfld = jfld + 1
791               dta%ssh => bf(jfld)%fnow(:,1,1)
792            ENDIF
793            IF ( dta%ll_u2d ) THEN
794               IF ( ln_full_vel_array(ib_bdy) ) THEN
795                  if(lwp) write(numout,*) '++++++ dta%u2d allocated space'
796                  ALLOCATE( dta%u2d(nblen(2)) )
797               ELSE
798                  if(lwp) write(numout,*) '++++++ dta%u2d pointing to fnow'
799                  jfld = jfld + 1
800                  dta%u2d => bf(jfld)%fnow(:,1,1)
801               ENDIF
802            ENDIF
803            IF ( dta%ll_v2d ) THEN
804               IF ( ln_full_vel_array(ib_bdy) ) THEN
805                  if(lwp) write(numout,*) '++++++ dta%v2d allocated space'
806                  ALLOCATE( dta%v2d(nblen(3)) )
807               ELSE
808                  if(lwp) write(numout,*) '++++++ dta%v2d pointing to fnow'
809                  jfld = jfld + 1
810                  dta%v2d => bf(jfld)%fnow(:,1,1)
811               ENDIF
812            ENDIF
813         ENDIF
814
815         IF ( nn_dyn3d_dta(ib_bdy) == 0 ) THEN
816            if(lwp) write(numout,*) '++++++ dta%u3d/v3d allocated space'
817            IF( dta%ll_u3d ) ALLOCATE( dta_bdy(ib_bdy)%u3d(nblen(2),jpk) )
818            IF( dta%ll_v3d ) ALLOCATE( dta_bdy(ib_bdy)%v3d(nblen(3),jpk) )
819         ENDIF
820         IF ( nn_dyn3d_dta(ib_bdy) == 1 .or. &
821           &  ( ln_full_vel_array(ib_bdy) .and. ( nn_dyn2d_dta(ib_bdy) == 1 .or. nn_dyn2d_dta(ib_bdy) == 3 ) ) ) THEN
822            IF ( dta%ll_u3d .or. ( ln_full_vel_array(ib_bdy) .and. dta%ll_u2d ) ) THEN
823               if(lwp) write(numout,*) '++++++ dta%u3d pointing to fnow'
824               jfld = jfld + 1
825               dta_bdy(ib_bdy)%u3d => bf(jfld)%fnow(:,1,:)
826            ENDIF
827            IF ( dta%ll_v3d .or. ( ln_full_vel_array(ib_bdy) .and. dta%ll_v2d ) ) THEN
828               if(lwp) write(numout,*) '++++++ dta%v3d pointing to fnow'
829               jfld = jfld + 1
830               dta_bdy(ib_bdy)%v3d => bf(jfld)%fnow(:,1,:)
831            ENDIF
832         ENDIF
833
834         IF( nn_tra_dta(ib_bdy) == 0 ) THEN
835            if(lwp) write(numout,*) '++++++ dta%tem/sal allocated space'
836            IF( dta%ll_tem ) ALLOCATE( dta_bdy(ib_bdy)%tem(nblen(1),jpk) )
837            IF( dta%ll_sal ) ALLOCATE( dta_bdy(ib_bdy)%sal(nblen(1),jpk) )
838         ELSE
839            IF( dta%ll_tem ) THEN
840               if(lwp) write(numout,*) '++++++ dta%tem pointing to fnow'
841               jfld = jfld + 1
842               dta_bdy(ib_bdy)%tem => bf(jfld)%fnow(:,1,:)
843            ENDIF
844            IF( dta%ll_sal ) THEN
845               if(lwp) write(numout,*) '++++++ dta%sal pointing to fnow'
846               jfld = jfld + 1
847               dta_bdy(ib_bdy)%sal => bf(jfld)%fnow(:,1,:)
848            ENDIF
849         ENDIF
850
851#if defined key_lim2
852         IF (cn_ice_lim(ib_bdy) /= 'none') THEN
853            IF( nn_ice_lim_dta(ib_bdy) == 0 ) THEN
854               ALLOCATE( dta_bdy(ib_bdy)%frld(nblen(1)) )
855               ALLOCATE( dta_bdy(ib_bdy)%hicif(nblen(1)) )
856               ALLOCATE( dta_bdy(ib_bdy)%hsnif(nblen(1)) )
857            ELSE
858               jfld = jfld + 1
859               dta_bdy(ib_bdy)%frld  => bf(jfld)%fnow(:,1,1)
860               jfld = jfld + 1
861               dta_bdy(ib_bdy)%hicif => bf(jfld)%fnow(:,1,1)
862               jfld = jfld + 1
863               dta_bdy(ib_bdy)%hsnif => bf(jfld)%fnow(:,1,1)
864            ENDIF
865         ENDIF
866#elif defined key_lim3
867         IF (cn_ice_lim(ib_bdy) /= 'none') THEN
868            IF( nn_ice_lim_dta(ib_bdy) == 0 ) THEN
869               ALLOCATE( dta_bdy(ib_bdy)%a_i (nblen(1),jpl) )
870               ALLOCATE( dta_bdy(ib_bdy)%ht_i(nblen(1),jpl) )
871               ALLOCATE( dta_bdy(ib_bdy)%ht_s(nblen(1),jpl) )
872            ELSE
873               IF ( ll_bdylim3 ) THEN ! case input is lim3 type
874                  jfld = jfld + 1
875                  dta_bdy(ib_bdy)%a_i  => bf(jfld)%fnow(:,1,:)
876                  jfld = jfld + 1
877                  dta_bdy(ib_bdy)%ht_i => bf(jfld)%fnow(:,1,:)
878                  jfld = jfld + 1
879                  dta_bdy(ib_bdy)%ht_s => bf(jfld)%fnow(:,1,:)
880               ELSE ! case input is lim2 type
881                  jfld_ai  = jfld + 1
882                  jfld_hti = jfld + 2
883                  jfld_hts = jfld + 3
884                  jfld     = jfld + 3
885                  ALLOCATE( dta_bdy(ib_bdy)%a_i (nblen(1),jpl) )
886                  ALLOCATE( dta_bdy(ib_bdy)%ht_i(nblen(1),jpl) )
887                  ALLOCATE( dta_bdy(ib_bdy)%ht_s(nblen(1),jpl) )
888                  dta_bdy(ib_bdy)%a_i (:,:) = 0._wp
889                  dta_bdy(ib_bdy)%ht_i(:,:) = 0._wp
890                  dta_bdy(ib_bdy)%ht_s(:,:) = 0._wp
891               ENDIF
892
893            ENDIF
894         ENDIF
895#endif
896         !
897      END DO ! ib_bdy
898      !
899      IF( nn_timing == 1 ) CALL timing_stop('bdy_dta_init')
900      !
901   END SUBROUTINE bdy_dta_init
902
903   !!==============================================================================
904END MODULE bdydta
Note: See TracBrowser for help on using the repository browser.