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/AMM15_v3_6_STABLE_package/NEMOGCM/NEMO/OPA_SRC/BDY – NEMO

source: branches/UKMO/AMM15_v3_6_STABLE_package/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90 @ 8790

Last change on this file since 8790 was 8790, checked in by deazer, 6 years ago

Add optional ssh switch to del with Baltic in CO7
Also added lbclnk in dynkeg for reproducible runs.

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