Changes between Version 11 and Version 12 of 2011WP/2011Stream2/Open Boundaries


Ignore:
Timestamp:
2011-03-29T13:47:08+02:00 (10 years ago)
Author:
davestorkey
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • 2011WP/2011Stream2/Open Boundaries

    v11 v12  
    223223   END IF 
    224224 
    225  
     225   jstart = 1 
     226   jfld = 0 
    226227   DO ib_set = 1, nb_set    
    227     
    228       IF( PRESENT(jit) ) THEN 
    229          ! Update barotropic boundary conditions only 
    230          ! jit is optional argument for fld_read 
    231          IF( nn_barotropic(ib_set) .gt. 0 ) THEN 
    232             CALL fld_read( kt, nn_fobc, bf(ib_set, 1:3), jit=jit ) 
    233          ENDIF 
    234       ELSE 
    235          IF ( nn_barotropic(ib_set) .gt. 0 .OR. & 
    236         &     nn_tra(ib_set)        .gt. 0 .OR. & 
    237         &     nn_dyn(ib_set)        .gt. 0 )      THEN 
     228      IF( nn_dtactl(ib_set) = 1 ) THEN 
     229       
     230         IF( PRESENT(jit) ) THEN 
     231            ! Update barotropic boundary conditions only 
     232            ! jit is optional argument for fld_read 
     233            IF( nn_dyn2d(ib_set) .gt. 0 ) THEN 
     234               jend = jstart + 2 
     235               CALL fld_read( kt, nn_fobc, bf(jstart:jend), jit=jit ) 
     236            ENDIF 
     237         ELSE 
     238            jend = jstart + nb_obc(ib_set) - 1 
    238239            CALL fld_read( kt, nn_fobc, bf(ib_set, 1:nb_dta(ib_set) ) 
    239240         ENDIF 
    240       ENDIF 
    241  
    242       IF( nn_barotropic .gt. 0 ) THEN 
    243          sshbdy(ib_set,:) = bf(ib_set, jp_ssh)%fnow 
    244          ubar(ib_set,:) = bf(ib_set, jp_ubar)%fnow 
    245          ... 
    246       ENDIF 
     241 
     242         IF( nn_dyn2d .gt. 0 ) THEN 
     243            jpfld = jpfld + 1 
     244            sshbdy(ib_set,:) = bf(jpfld)%fnow 
     245            jpfld = jpfld + 1 
     246            ubar(ib_set,:) = bf(jpfld)%fnow 
     247            ... 
     248         ENDIF 
    247249       
    248       IF( .NOT. PRESENT(jit) ) THEN 
    249          IF( nn_tra(ib_set) .gt. 0 ) THEN 
    250          .... 
     250         IF( .NOT. PRESENT(jit) ) THEN 
     251            IF( nn_tra(ib_set) .gt. 0 ) THEN 
     252               jpfld = jpfld + 1 
     253               sshbdy(ib_set,:) = bf(jpfld)%fnow 
     254               .... 
    251255          
    252256 
     257      END IF ! nn_dtactl(ib_set) = 1 
    253258   END DO  ! ib_set 
    254259