source: branches/2014/dev_r4621_NOC4_BDY_VERT_INTERP/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90 @ 6079

Last change on this file since 6079 was 6079, checked in by jamesharle, 6 years ago

merge to trunk@5936

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