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 10115 for NEMO/branches/2018/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90 – NEMO

Ignore:
Timestamp:
2018-09-12T15:59:13+02:00 (6 years ago)
Author:
cbricaud
Message:

phase 3.6 coarsening branch with nemo_3.6_rev9192

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2018/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90

    r7256 r10115  
    3636   USE limvar          ! redistribute ice input into categories 
    3737#endif 
     38   USE sbc_oce 
    3839   USE sbcapr 
    3940 
     
    395396      ! end jchanut tschanges 
    396397 
    397       IF ( ln_apr_obc ) THEN 
    398          DO ib_bdy = 1, nb_bdy 
    399             IF (cn_tra(ib_bdy) /= 'runoff')THEN 
    400                igrd = 1                      ! meridional velocity 
    401                DO ib = 1, idx_bdy(ib_bdy)%nblenrim(igrd) 
    402                   ii   = idx_bdy(ib_bdy)%nbi(ib,igrd) 
    403                   ij   = idx_bdy(ib_bdy)%nbj(ib,igrd) 
    404                   dta_bdy(ib_bdy)%ssh(ib) = dta_bdy(ib_bdy)%ssh(ib) + ssh_ib(ii,ij) 
    405                ENDDO 
    406             ENDIF 
    407          ENDDO 
     398      IF( ln_apr_dyn )THEN 
     399         IF( ln_apr_obc ) THEN 
     400            DO ib_bdy = 1, nb_bdy 
     401               IF (cn_tra(ib_bdy) /= 'runoff')THEN 
     402                  igrd = 1                      ! meridional velocity 
     403                  DO ib = 1, idx_bdy(ib_bdy)%nblenrim(igrd) 
     404                     ii   = idx_bdy(ib_bdy)%nbi(ib,igrd) 
     405                     ij   = idx_bdy(ib_bdy)%nbj(ib,igrd) 
     406                     dta_bdy(ib_bdy)%ssh(ib) = dta_bdy(ib_bdy)%ssh(ib) + ssh_ib(ii,ij) 
     407                  ENDDO 
     408               ENDIF 
     409            ENDDO 
     410         ENDIF 
    408411      ENDIF 
    409412 
     
    440443      TYPE(OBC_DATA), POINTER                ::   dta           ! short cut 
    441444#if defined key_lim3 
    442       INTEGER, DIMENSION(3) ::   zdimsz   ! number of elements in each of the 4 dimensions (i.e. i,j,t,ice-cat) for an array 
    443445      INTEGER               ::   zndims   ! number of dimensions in an array (i.e. 3 = wo ice cat; 4 = w ice cat) 
    444446      INTEGER               ::   inum,id1 ! local integer 
     
    672674               ! Test for types of ice input (lim2 or lim3) 
    673675               ! Build file name to find dimensions  
    674                clname=TRIM(bn_a_i%clname) 
     676               clname=TRIM( cn_dir )//TRIM(bn_a_i%clname) 
    675677               IF( .NOT. bn_a_i%ln_clim ) THEN    
    676                                                   WRITE(clname, '(a,"_y",i4.4)' ) TRIM( bn_a_i%clname ), nyear    ! add year 
    677                   IF( bn_a_i%cltype /= 'yearly' ) WRITE(clname, '(a,"m" ,i2.2)' ) TRIM( clname        ), nmonth   ! add month 
     678                                                  WRITE(clname, '(a,"_y",i4.4)' ) TRIM( clname ), nyear    ! add year 
     679                  IF( bn_a_i%cltype /= 'yearly' ) WRITE(clname, '(a,"m" ,i2.2)' ) TRIM( clname ), nmonth   ! add month 
    678680               ELSE 
    679                   IF( bn_a_i%cltype /= 'yearly' ) WRITE(clname, '(a,"_m",i2.2)' ) TRIM( bn_a_i%clname ), nmonth   ! add month 
     681                  IF( bn_a_i%cltype /= 'yearly' ) WRITE(clname, '(a,"_m",i2.2)' ) TRIM( clname ), nmonth   ! add month 
    680682               ENDIF 
    681683               IF( bn_a_i%cltype == 'daily' .OR. bn_a_i%cltype(1:4) == 'week' ) & 
    682                &                                  WRITE(clname, '(a,"d" ,i2.2)' ) TRIM( clname        ), nday     ! add day 
     684               &                                  WRITE(clname, '(a,"d" ,i2.2)' ) TRIM( clname ), nday     ! add day 
    683685               ! 
    684686               CALL iom_open  ( clname, inum ) 
    685                id1 = iom_varid( inum, bn_a_i%clvar, kdimsz=zdimsz, kndims=zndims, ldstop = .FALSE. ) 
     687               id1 = iom_varid( inum, bn_a_i%clvar, kndims=zndims, ldstop = .FALSE. ) 
    686688               CALL iom_close ( inum ) 
    687689 
Note: See TracChangeset for help on using the changeset viewer.