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

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

Last change on this file since 5930 was 5930, checked in by jchanut, 8 years ago

#1620 Merge free surface simplification into trunk

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