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

Last change on this file since 6140 was 6140, checked in by timgraham, 5 years ago

Merge of branches/2015/dev_merge_2015 back into trunk. Merge excludes NEMOGCM/TOOLS/OBSTOOLS/ for now due to issues with the change of file type. Will sort these manually with further commits.

Branch merged as follows:
In the working copy of branch ran:
svn merge svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/trunk@HEAD
Small conflicts due to bug fixes applied to trunk since the dev_merge_2015 was copied. Bug fixes were applied to the branch as well so these were easy to resolve.
Branch committed at this stage

In working copy run:
svn switch svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/trunk
to switch working copy

Run:
svn merge —reintegrate svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/branches/2015/dev_merge_2015
to merge the branch into the trunk and then commit - no conflicts at this stage.

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