New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
bdydta.F90 in branches/UKMO/dev_5518_tide_analysis_restart/NEMOGCM/NEMO/OPA_SRC/BDY – NEMO

source: branches/UKMO/dev_5518_tide_analysis_restart/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90 @ 6061

Last change on this file since 6061 was 6061, checked in by deazer, 8 years ago

Removed SVN Keywords

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