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.
Changeset 2168 for branches/DEV_r1986_BDY_updates/NEMO/OPA_SRC/BDY/bdydta.F90 – NEMO

Ignore:
Timestamp:
2010-10-06T16:19:27+02:00 (13 years ago)
Author:
rblod
Message:

Cosmetic changes on BDY branch

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/DEV_r1986_BDY_updates/NEMO/OPA_SRC/BDY/bdydta.F90

    r2100 r2168  
    1616   !!---------------------------------------------------------------------- 
    1717   !!   bdy_dta    : read u, v, t, s data along open boundaries 
    18    !!   bdy_dta_bt : read depth-mean velocities and elevation along open 
    19    !!                boundaries         
     18   !!   bdy_dta_bt : read depth-mean velocities and elevation along open boundaries         
    2019   !!---------------------------------------------------------------------- 
    2120   USE oce             ! ocean dynamics and tracers 
     
    3736   PUBLIC   bdy_dta_bt  
    3837 
    39    INTEGER ::   numbdyt, numbdyu, numbdyv                      !: logical units for T-, U-, & V-points data file, resp. 
    40    INTEGER ::   ntimes_bdy                                     !: exact number of time dumps in data files 
    41    INTEGER ::   nbdy_b, nbdy_a                                 !: record of bdy data file for before and after model time step 
    42    INTEGER ::   numbdyt_bt, numbdyu_bt, numbdyv_bt             !: logical unit for T-, U- & V-points data file, resp. 
    43    INTEGER ::   ntimes_bdy_bt                                  !: exact number of time dumps in data files 
    44    INTEGER ::   nbdy_b_bt, nbdy_a_bt                           !: record of bdy data file for before and after model time step 
    45  
    46    INTEGER, DIMENSION (jpbtime) ::   istep, istep_bt           !: time array in seconds in each data file 
    47  
    48    REAL(wp) ::  zoffset                                        !: time offset between time origin in file & start time of model run 
    49  
    50    REAL(wp), DIMENSION(jpbdim,jpk,2) ::   tbdydta, sbdydta     !: time interpolated values of T and S bdy data    
    51    REAL(wp), DIMENSION(jpbdim,jpk,2) ::   ubdydta, vbdydta     !: time interpolated values of U and V bdy data  
    52    REAL(wp), DIMENSION(jpbdim,2)     ::   ubtbdydta, vbtbdydta !: Arrays used for time interpolation of bdy data    
    53    REAL(wp), DIMENSION(jpbdim,2)     ::   sshbdydta            !: bdy data of ssh 
     38   INTEGER ::   numbdyt, numbdyu, numbdyv                      ! logical units for T-, U-, & V-points data file, resp. 
     39   INTEGER ::   ntimes_bdy                                     ! exact number of time dumps in data files 
     40   INTEGER ::   nbdy_b, nbdy_a                                 ! record of bdy data file for before and after time step 
     41   INTEGER ::   numbdyt_bt, numbdyu_bt, numbdyv_bt             ! logical unit for T-, U- & V-points data file, resp. 
     42   INTEGER ::   ntimes_bdy_bt                                  ! exact number of time dumps in data files 
     43   INTEGER ::   nbdy_b_bt, nbdy_a_bt                           ! record of bdy data file for before and after time step 
     44 
     45   INTEGER, DIMENSION (jpbtime) ::   istep, istep_bt           ! time array in seconds in each data file 
     46 
     47   REAL(wp) ::  zoffset                                        ! time offset between time origin in file & start time of model run 
     48 
     49   REAL(wp), DIMENSION(jpbdim,jpk,2) ::   tbdydta, sbdydta     ! time interpolated values of T and S bdy data    
     50   REAL(wp), DIMENSION(jpbdim,jpk,2) ::   ubdydta, vbdydta     ! time interpolated values of U and V bdy data  
     51   REAL(wp), DIMENSION(jpbdim,2)     ::   ubtbdydta, vbtbdydta ! Arrays used for time interpolation of bdy data    
     52   REAL(wp), DIMENSION(jpbdim,2)     ::   sshbdydta            ! bdy data of ssh 
    5453 
    5554#if defined key_lim2 
    56    REAL(wp), DIMENSION(jpbdim,2)     ::   frld_bdydta          !: } 
    57    REAL(wp), DIMENSION(jpbdim,2)     ::   hicif_bdydta         !: } Arrays used for time interpolation of bdy data for ice variables 
    58    REAL(wp), DIMENSION(jpbdim,2)     ::   hsnif_bdydta         !: } 
     55   REAL(wp), DIMENSION(jpbdim,2)     ::   frld_bdydta          ! } 
     56   REAL(wp), DIMENSION(jpbdim,2)     ::   hicif_bdydta         ! } Arrays used for time interp. of ice bdy data  
     57   REAL(wp), DIMENSION(jpbdim,2)     ::   hsnif_bdydta         ! } 
    5958#endif 
    6059 
    6160   !!---------------------------------------------------------------------- 
    62    !! NEMO/OPA 3.0 , LOCEAN-IPSL (2008)  
     61   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    6362   !! $Id$  
    6463   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 
    6564   !!---------------------------------------------------------------------- 
    66  
    6765CONTAINS 
    6866 
     
    7876      !!                the file. If so read it in. Time interpolate. 
    7977      !!---------------------------------------------------------------------- 
    80       INTEGER, INTENT( in ) ::   kt                             ! ocean time-step index (for timesplitting option, otherwise zero) 
     78      INTEGER, INTENT( in ) ::   kt   ! ocean time-step index (for timesplitting option, otherwise zero) 
    8179      !! 
    8280      CHARACTER(LEN=80), DIMENSION(3) ::   clfile               ! names of input files 
     
    9088      INTEGER ::   itimer, totime 
    9189      INTEGER ::   ii, ij                                       ! array addresses 
    92       INTEGER ::   ipi, ipj, ipk, inum                          ! temporary integers (NetCDF read) 
     90      INTEGER ::   ipi, ipj, ipk, inum                          ! local integers (NetCDF read) 
    9391      INTEGER ::   iyear0, imonth0, iday0 
    9492      INTEGER ::   ihours0, iminutes0, isec0 
     
    102100      !!--------------------------------------------------------------------------- 
    103101 
    104       IF( ln_bdy_dyn_frs .OR. ln_bdy_tra_frs ) THEN  ! If these are both false then this routine  
    105                                                      ! does nothing. 
     102 
     103      IF( ln_bdy_dyn_frs .OR. ln_bdy_tra_frs    & 
     104         &               .OR. ln_bdy_ice_frs ) THEN  ! If these are both false then this routine does nothing 
    106105 
    107106      ! -------------------- ! 
     
    113112      ! Some time variables for monthly climatological forcing: 
    114113      ! ******************************************************* 
    115  !!gm  here  use directely daymod variables 
     114 
     115!!gm  here  use directely daymod calendar variables 
    116116  
    117117      iman = INT( raamo )      ! Number of months in a year 
     
    132132         !                                             !-------------------! 
    133133         istep(:) = 0 
    134          nbdy_b    = 0 
    135          nbdy_a    = 0 
     134         nbdy_b   = 0 
     135         nbdy_a   = 0 
    136136 
    137137         ! Get time information from bdy data file 
     
    162162            igrd_start = 1 
    163163            igrd_end   = 3 
    164             IF(.NOT. ln_bdy_tra_frs .AND. .NOT. ln_bdy_ice_frs) THEN 
    165                ! No T-grid file. 
     164            IF(.NOT. ln_bdy_tra_frs .AND. .NOT. ln_bdy_ice_frs) THEN       ! No T-grid file. 
    166165               igrd_start = 2 
    167             ELSEIF ( .NOT. ln_bdy_dyn_frs ) THEN 
    168                ! No U-grid or V-grid file. 
     166            ELSEIF ( .NOT. ln_bdy_dyn_frs ) THEN                           ! No U-grid or V-grid file. 
    169167               igrd_end   = 1          
    170168            ENDIF 
     
    176174 
    177175               SELECT CASE( igrd ) 
    178                   CASE (1)  
    179                      numbdyt = inum 
    180                   CASE (2)  
    181                      numbdyu = inum 
    182                   CASE (3)  
    183                      numbdyv = inum 
     176                  CASE (1)   ;   numbdyt = inum 
     177                  CASE (2)   ;   numbdyu = inum 
     178                  CASE (3)   ;   numbdyv = inum 
    184179               END SELECT 
    185180 
     
    216211 
    217212               ! Check that time array increases: 
    218          
    219213               it = 1 
    220                DO WHILE( zstepr(it+1) > zstepr(it) .AND. it /= ntimes_bdy - 1 ) 
    221                  it = it + 1 
     214               DO WHILE( zstepr(it+1) > zstepr(it) .AND. it /= ntimes_bdy - 1 )  
     215                  it = it + 1 
    222216               END DO 
    223  
    224                IF( it.NE.ntimes_bdy-1 .AND. ntimes_bdy > 1 ) THEN 
     217               ! 
     218               IF( it /= ntimes_bdy-1 .AND. ntimes_bdy > 1 ) THEN 
    225219                     WRITE(ctmp1,*) 'Check file: ', clfile(igrd) 
    226220                     CALL ctl_stop( 'Time array in unstructured boundary data files',   & 
     
    274268            ENDIF 
    275269 
    276             IF( igrd_start == 1 ) THEN 
    277                istep(:) = istept(:) 
    278             ELSE 
    279                istep(:) = istepu(:) 
     270            IF( igrd_start == 1 ) THEN   ;   istep(:) = istept(:) 
     271            ELSE                         ;   istep(:) = istepu(:) 
    280272            ENDIF 
    281273 
     
    302294            it = 1 
    303295            DO WHILE( istep(it+1) <= 0 .AND. it <= ntimes_bdy - 1 ) 
    304               it = it + 1 
     296               it = it + 1 
    305297            END DO 
    306298            nbdy_b = it 
     
    315307         ! ***************************************************************** 
    316308 
    317          IF( nbdy_dta == 0) THEN      ! boundary data arrays are filled with initial conditions 
     309         IF( nbdy_dta == 0 ) THEN      ! boundary data arrays are filled with initial conditions 
    318310            ! 
    319311            IF (ln_bdy_tra_frs) THEN 
    320               igrd = 1            ! T-points data  
    321               DO ib = 1, nblen(igrd) 
    322                 ii = nbi(ib,igrd) 
    323                 ij = nbj(ib,igrd) 
    324                 DO ik = 1, jpkm1 
    325                   tbdy(ib,ik) = tn(ii, ij, ik) 
    326                   sbdy(ib,ik) = sn(ii, ij, ik) 
    327                 ENDDO 
    328               END DO 
     312               igrd = 1            ! T-points data  
     313               DO ib = 1, nblen(igrd) 
     314                  ii = nbi(ib,igrd) 
     315                  ij = nbj(ib,igrd) 
     316                  DO ik = 1, jpkm1 
     317                     tbdy(ib,ik) = tn(ii,ij,ik) 
     318                     sbdy(ib,ik) = sn(ii,ij,ik) 
     319                  END DO 
     320               END DO 
    329321            ENDIF 
    330322 
    331323            IF(ln_bdy_dyn_frs) THEN 
    332               igrd = 2            ! U-points data  
    333               DO ib = 1, nblen(igrd) 
    334                 ii = nbi(ib,igrd) 
    335                 ij = nbj(ib,igrd) 
    336                 DO ik = 1, jpkm1 
    337                   ubdy(ib,ik) = un(ii, ij, ik) 
    338                 ENDDO 
    339               END DO 
    340  
    341               igrd = 3            ! V-points data  
    342               DO ib = 1, nblen(igrd)             
    343                 ii = nbi(ib,igrd) 
    344                 ij = nbj(ib,igrd) 
    345                 DO ik = 1, jpkm1 
    346                   vbdy(ib,ik) = vn(ii, ij, ik) 
    347                 ENDDO 
    348               END DO 
     324               igrd = 2            ! U-points data  
     325               DO ib = 1, nblen(igrd) 
     326                  ii = nbi(ib,igrd) 
     327                  ij = nbj(ib,igrd) 
     328                  DO ik = 1, jpkm1 
     329                     ubdy(ib,ik) = un(ii, ij, ik) 
     330                  END DO 
     331               END DO 
     332               ! 
     333               igrd = 3            ! V-points data  
     334               DO ib = 1, nblen(igrd)             
     335                  ii = nbi(ib,igrd) 
     336                  ij = nbj(ib,igrd) 
     337                  DO ik = 1, jpkm1 
     338                     vbdy(ib,ik) = vn(ii, ij, ik) 
     339                  END DO 
     340               END DO 
    349341            ENDIF 
    350342            ! 
    351343#if defined key_lim2 
    352             IF (ln_bdy_ice_frs) THEN 
     344            IF( ln_bdy_ice_frs ) THEN 
    353345               igrd = 1            ! T-points data 
    354346               DO ib = 1, nblen(igrd) 
    355                   frld_bdy(ib) =  frld(nbi(ib,igrd), nbj(ib,igrd)) 
     347                  frld_bdy (ib) =  frld(nbi(ib,igrd), nbj(ib,igrd)) 
    356348                  hicif_bdy(ib) = hicif(nbi(ib,igrd), nbj(ib,igrd)) 
    357349                  hsnif_bdy(ib) = hsnif(nbi(ib,igrd), nbj(ib,igrd)) 
     
    377369 
    378370            IF(ln_bdy_tra_frs) THEN 
     371               ! 
    379372               igrd = 1                                           ! Temperature 
    380373               IF( nblendta(igrd) <=  0 ) THEN  
     
    385378               ipi = nblendta(igrd) 
    386379               CALL iom_get ( numbdyt, jpdom_unknown, 'votemper', zdta(1:ipi,1:ipj,1:ipk), nbdy_a ) 
    387  
     380               ! 
    388381               DO ib = 1, nblen(igrd) 
    389382                  DO ik = 1, jpkm1 
     
    400393               ipi = nblendta(igrd) 
    401394               CALL iom_get ( numbdyt, jpdom_unknown, 'vosaline', zdta(1:ipi,1:ipj,1:ipk), nbdy_a ) 
    402  
     395               ! 
    403396               DO ib = 1, nblen(igrd) 
    404397                  DO ik = 1, jpkm1 
     
    408401            ENDIF  ! ln_bdy_tra_frs 
    409402  
    410             IF(ln_bdy_dyn_frs) THEN 
    411  
     403            IF( ln_bdy_dyn_frs ) THEN 
     404               ! 
    412405               igrd = 2                                           ! u-velocity 
    413406               IF ( nblendta(igrd) .le. 0 ) THEN  
     
    440433 
    441434#if defined key_lim2 
    442             IF(ln_bdy_ice_frs) THEN 
     435            IF( ln_bdy_ice_frs ) THEN 
    443436              ! 
    444437              igrd=1                                              ! leads fraction 
     
    468461#endif 
    469462 
    470             IF ((.NOT.ln_bdy_clim) .AND. (istep(1) > 0)) THEN 
    471                ! First data time is after start of run 
    472                ! Put first value in both time levels 
     463            IF( .NOT.ln_bdy_clim .AND. istep(1) > 0 ) THEN     ! First data time is after start of run 
     464               nbdy_b = nbdy_a                                 ! Put first value in both time levels 
     465               IF( ln_bdy_tra_frs ) THEN 
     466                 tbdydta(:,:,1) = tbdydta(:,:,2) 
     467                 sbdydta(:,:,1) = sbdydta(:,:,2) 
     468               ENDIF 
     469               IF( ln_bdy_dyn_frs ) THEN 
     470                 ubdydta(:,:,1) = ubdydta(:,:,2) 
     471                 vbdydta(:,:,1) = vbdydta(:,:,2) 
     472               ENDIF 
     473#if defined key_lim2 
     474               IF( ln_bdy_ice_frs ) THEN 
     475                  frld_bdydta (:,1) =  frld_bdydta(:,2) 
     476                  hicif_bdydta(:,1) = hicif_bdydta(:,2) 
     477                  hsnif_bdydta(:,1) = hsnif_bdydta(:,2) 
     478               ENDIF 
     479#endif 
     480            END IF 
     481            ! 
     482         END IF   ! nbdy_dta == 0/1 
     483  
     484         ! In the case of constant boundary forcing fill bdy arrays once for all 
     485         IF( ln_bdy_clim .AND. ntimes_bdy == 1 ) THEN 
     486            IF( ln_bdy_tra_frs ) THEN 
     487               tbdy  (:,:) = tbdydta  (:,:,2) 
     488               sbdy  (:,:) = sbdydta  (:,:,2) 
     489            ENDIF 
     490            IF( ln_bdy_dyn_frs) THEN 
     491               ubdy  (:,:) = ubdydta  (:,:,2) 
     492               vbdy  (:,:) = vbdydta  (:,:,2) 
     493            ENDIF 
     494#if defined key_lim2 
     495            IF( ln_bdy_ice_frs ) THEN 
     496               frld_bdy (:) = frld_bdydta (:,2) 
     497               hicif_bdy(:) = hicif_bdydta(:,2) 
     498               hsnif_bdy(:) = hsnif_bdydta(:,2) 
     499            ENDIF 
     500#endif 
     501 
     502            IF( ln_bdy_tra_frs .OR. ln_bdy_ice_frs) CALL iom_close( numbdyt ) 
     503            IF( ln_bdy_dyn_frs                    ) CALL iom_close( numbdyu ) 
     504            IF( ln_bdy_dyn_frs                    ) CALL iom_close( numbdyv ) 
     505         END IF 
     506         ! 
     507      ENDIF                                            ! End if nit000 
     508 
     509 
     510      !                                                !---------------------! 
     511      IF( nbdy_dta == 1 .AND. ntimes_bdy > 1 ) THEN    !  at each time step  ! 
     512         !                                             !---------------------! 
     513         ! Read one more record if necessary 
     514         !********************************** 
     515 
     516         IF( ln_bdy_clim .AND. imois /= nbdy_b ) THEN      ! remember that nbdy_b=0 for kt=nit000 
     517            nbdy_b = imois 
     518            nbdy_a = imois + 1 
     519            nbdy_b = MOD( nbdy_b, iman )   ;   IF( nbdy_b == 0 ) nbdy_b = iman 
     520            nbdy_a = MOD( nbdy_a, iman )   ;   IF( nbdy_a == 0 ) nbdy_a = iman 
     521            lect=.true. 
     522         ELSEIF( .NOT.ln_bdy_clim .AND. itimer >= istep(nbdy_a) ) THEN 
     523 
     524            IF( nbdy_a < ntimes_bdy ) THEN 
     525               nbdy_b = nbdy_a 
     526               nbdy_a = nbdy_a + 1 
     527               lect  =.true. 
     528            ELSE 
     529               ! We have reached the end of the file 
     530               ! put the last data time into both time levels 
    473531               nbdy_b = nbdy_a 
    474532               IF(ln_bdy_tra_frs) THEN 
    475                  tbdydta(:,:,1) = tbdydta(:,:,2) 
    476                  sbdydta(:,:,1) = sbdydta(:,:,2) 
     533                  tbdydta(:,:,1) = tbdydta(:,:,2) 
     534                  sbdydta(:,:,1) = sbdydta(:,:,2) 
    477535               ENDIF 
    478536               IF(ln_bdy_dyn_frs) THEN 
    479                  ubdydta(:,:,1) = ubdydta(:,:,2) 
    480                  vbdydta(:,:,1) = vbdydta(:,:,2) 
     537                  ubdydta(:,:,1) = ubdydta(:,:,2) 
     538                  vbdydta(:,:,1) = vbdydta(:,:,2) 
    481539               ENDIF 
    482540#if defined key_lim2 
    483              IF( ln_bdy_ice_frs ) THEN 
    484                frld_bdydta (:,1) =  frld_bdydta(:,2) 
    485                hicif_bdydta(:,1) = hicif_bdydta(:,2) 
    486                hsnif_bdydta(:,1) = hsnif_bdydta(:,2) 
    487              ENDIF 
    488 #endif 
    489             END IF 
    490  
    491          END IF ! nbdy_dta == 0/1 
    492   
    493          ! In the case of constant boundary forcing fill bdy arrays once for all 
    494          IF ((ln_bdy_clim).AND.(ntimes_bdy==1)) THEN 
    495             IF(ln_bdy_tra_frs) THEN 
    496               tbdy  (:,:) = tbdydta  (:,:,2) 
    497               sbdy  (:,:) = sbdydta  (:,:,2) 
    498             ENDIF 
    499             IF(ln_bdy_dyn_frs) THEN 
    500               ubdy  (:,:) = ubdydta  (:,:,2) 
    501               vbdy  (:,:) = vbdydta  (:,:,2) 
    502             ENDIF 
    503 #if defined key_lim2 
    504           IF(ln_bdy_ice_frs) THEN 
    505             frld_bdy (:) = frld_bdydta (:,2) 
    506             hicif_bdy(:) = hicif_bdydta(:,2) 
    507             hsnif_bdy(:) = hsnif_bdydta(:,2) 
    508           ENDIF 
    509 #endif 
    510  
    511             IF(ln_bdy_tra_frs .or. ln_bdy_ice_frs) CALL iom_close( numbdyt ) 
    512             IF(ln_bdy_dyn_frs) CALL iom_close( numbdyu ) 
    513             IF(ln_bdy_dyn_frs) CALL iom_close( numbdyv ) 
    514          END IF 
    515  
    516       ENDIF                                            ! End if nit000 
    517  
    518  
    519       !                                                !---------------------! 
    520       !                                                !  at each time step  ! 
    521       !                                                !---------------------! 
    522  
    523       IF( nbdy_dta == 1 .AND. ntimes_bdy > 1 ) THEN  
    524          ! 
    525          ! Read one more record if necessary 
    526          !********************************** 
    527  
    528         IF( ln_bdy_clim .AND. imois /= nbdy_b ) THEN      ! remember that nbdy_b=0 for kt=nit000 
    529            nbdy_b = imois 
    530            nbdy_a = imois + 1 
    531            nbdy_b = MOD( nbdy_b, iman )   ;   IF( nbdy_b == 0 ) nbdy_b = iman 
    532            nbdy_a = MOD( nbdy_a, iman )   ;   IF( nbdy_a == 0 ) nbdy_a = iman 
    533            lect=.true. 
    534         ELSEIF( .NOT.ln_bdy_clim .AND. itimer >= istep(nbdy_a) ) THEN 
    535  
    536            IF ( nbdy_a < ntimes_bdy ) THEN 
    537               nbdy_b = nbdy_a 
    538               nbdy_a = nbdy_a + 1 
    539               lect  =.true. 
    540            ELSE 
    541               ! We have reached the end of the file 
    542               ! put the last data time into both time levels 
    543               nbdy_b = nbdy_a 
    544               IF(ln_bdy_tra_frs) THEN 
    545                 tbdydta(:,:,1) =  tbdydta(:,:,2) 
    546                 sbdydta(:,:,1) =  sbdydta(:,:,2) 
    547               ENDIF 
    548               IF(ln_bdy_dyn_frs) THEN 
    549                 ubdydta(:,:,1) =  ubdydta(:,:,2) 
    550                 vbdydta(:,:,1) =  vbdydta(:,:,2) 
    551               ENDIF 
    552 #if defined key_lim2 
    553               IF(ln_bdy_ice_frs) THEN 
    554                 frld_bdydta (:,1) =  frld_bdydta (:,2) 
    555                 hicif_bdydta(:,1) =  hicif_bdydta(:,2) 
    556                 hsnif_bdydta(:,1) =  hsnif_bdydta(:,2) 
    557               ENDIF 
     541               IF(ln_bdy_ice_frs) THEN 
     542                  frld_bdydta (:,1) =  frld_bdydta (:,2) 
     543                  hicif_bdydta(:,1) =  hicif_bdydta(:,2) 
     544                  hsnif_bdydta(:,1) =  hsnif_bdydta(:,2) 
     545               ENDIF 
    558546#endif 
    559547            END IF ! nbdy_a < ntimes_bdy 
    560  
     548            ! 
    561549        END IF 
    562550          
    563         IF( lect ) THEN 
    564            ! Swap arrays 
    565            IF(ln_bdy_tra_frs) THEN 
     551        IF( lect ) THEN           ! Swap arrays 
     552           IF( ln_bdy_tra_frs ) THEN 
    566553             tbdydta(:,:,1) =  tbdydta(:,:,2) 
    567554             sbdydta(:,:,1) =  sbdydta(:,:,2) 
    568555           ENDIF 
    569            IF(ln_bdy_dyn_frs) THEN 
     556           IF( ln_bdy_dyn_frs ) THEN 
    570557             ubdydta(:,:,1) =  ubdydta(:,:,2) 
    571558             vbdydta(:,:,1) =  vbdydta(:,:,2) 
    572559           ENDIF 
    573560#if defined key_lim2 
    574            IF(ln_bdy_ice_frs) THEN 
     561           IF( ln_bdy_ice_frs ) THEN 
    575562             frld_bdydta (:,1) =  frld_bdydta (:,2) 
    576563             hicif_bdydta(:,1) =  hicif_bdydta(:,2) 
     
    582569           ipk  = jpk 
    583570 
    584            IF(ln_bdy_tra_frs) THEN 
     571           IF( ln_bdy_tra_frs ) THEN 
    585572              !  
    586573              igrd = 1                                   ! temperature 
     
    720707      ! 
    721708      ENDIF ! ln_bdy_dyn_frs .OR. ln_bdy_tra_frs 
    722  
     709      ! 
    723710   END SUBROUTINE bdy_dta 
    724711 
Note: See TracChangeset for help on using the changeset viewer.