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 14169 – NEMO

Changeset 14169


Ignore:
Timestamp:
2020-12-14T19:32:36+01:00 (3 years ago)
Author:
jchanut
Message:

#2592, ensure all processors read 2d tidal data

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/BDY/bdytides.F90

    r13286 r14169  
    162162               ! 
    163163               ! SSH fields 
    164                IF( ASSOCIATED(dta%ssh) ) THEN   ! we use bdy ssh on this mpi subdomain 
    165                   clfile = TRIM(filtide)//'_grid_T.nc' 
    166                   CALL iom_open( clfile , inum )  
    167                   igrd = 1                       ! Everything is at T-points here 
    168                   DO itide = 1, nb_harmo 
    169                      CALL iom_get( inum, jpdom_auto, TRIM(tide_harmonics(itide)%cname_tide)//'_z1', ztr(:,:) ) 
    170                      CALL iom_get( inum, jpdom_auto, TRIM(tide_harmonics(itide)%cname_tide)//'_z2', zti(:,:) )  
     164               clfile = TRIM(filtide)//'_grid_T.nc' 
     165               CALL iom_open( clfile , inum )  
     166               igrd = 1                       ! Everything is at T-points here 
     167               DO itide = 1, nb_harmo 
     168                  CALL iom_get( inum, jpdom_auto, TRIM(tide_harmonics(itide)%cname_tide)//'_z1', ztr(:,:) ) 
     169                  CALL iom_get( inum, jpdom_auto, TRIM(tide_harmonics(itide)%cname_tide)//'_z2', zti(:,:) )  
     170                  IF( ASSOCIATED(dta%ssh) ) THEN   ! we use bdy ssh on this mpi subdomain 
    171171                     DO ib = 1, SIZE(dta%ssh) 
    172172                        ii = idx_bdy(ib_bdy)%nbi(ib,igrd) 
     
    175175                        td%ssh0(ib,itide,2) = zti(ii,ij) 
    176176                     END DO 
    177                   END DO 
    178                   CALL iom_close( inum ) 
    179                ENDIF 
     177                  ENDIF 
     178               END DO 
     179               CALL iom_close( inum ) 
    180180               ! 
    181181               ! U fields 
    182                IF( ASSOCIATED(dta%u2d) ) THEN   ! we use bdy u2d on this mpi subdomain 
    183                   clfile = TRIM(filtide)//'_grid_U.nc' 
    184                   CALL iom_open( clfile , inum )  
    185                   igrd = 2                       ! Everything is at U-points here 
    186                   DO itide = 1, nb_harmo 
    187                      CALL iom_get(inum, jpdom_auto, TRIM(tide_harmonics(itide)%cname_tide)//'_u1', ztr(:,:),cd_type='U',psgn=-1._wp) 
    188                      CALL iom_get(inum, jpdom_auto, TRIM(tide_harmonics(itide)%cname_tide)//'_u2', zti(:,:),cd_type='U',psgn=-1._wp) 
     182               clfile = TRIM(filtide)//'_grid_U.nc' 
     183               CALL iom_open( clfile , inum )  
     184               igrd = 2                       ! Everything is at U-points here 
     185               DO itide = 1, nb_harmo 
     186                  CALL iom_get(inum, jpdom_auto, TRIM(tide_harmonics(itide)%cname_tide)//'_u1', ztr(:,:),cd_type='U',psgn=-1._wp) 
     187                  CALL iom_get(inum, jpdom_auto, TRIM(tide_harmonics(itide)%cname_tide)//'_u2', zti(:,:),cd_type='U',psgn=-1._wp) 
     188                  IF( ASSOCIATED(dta%u2d) ) THEN   ! we use bdy u2d on this mpi subdomain 
    189189                     DO ib = 1, SIZE(dta%u2d) 
    190190                        ii = idx_bdy(ib_bdy)%nbi(ib,igrd) 
     
    193193                        td%u0(ib,itide,2) = zti(ii,ij) 
    194194                     END DO 
    195                   END DO 
    196                   CALL iom_close( inum ) 
    197                ENDIF 
     195                  ENDIF  
     196               END DO  
     197               CALL iom_close( inum ) 
    198198               ! 
    199199               ! V fields 
    200                IF( ASSOCIATED(dta%v2d) ) THEN   ! we use bdy v2d on this mpi subdomain 
    201                   clfile = TRIM(filtide)//'_grid_V.nc' 
    202                   CALL iom_open( clfile , inum )  
    203                   igrd = 3                       ! Everything is at V-points here 
    204                   DO itide = 1, nb_harmo 
    205                      CALL iom_get(inum, jpdom_auto, TRIM(tide_harmonics(itide)%cname_tide)//'_v1', ztr(:,:),cd_type='V',psgn=-1._wp) 
    206                      CALL iom_get(inum, jpdom_auto, TRIM(tide_harmonics(itide)%cname_tide)//'_v2', zti(:,:),cd_type='V',psgn=-1._wp) 
     200               clfile = TRIM(filtide)//'_grid_V.nc' 
     201               CALL iom_open( clfile , inum )  
     202               igrd = 3                       ! Everything is at V-points here 
     203               DO itide = 1, nb_harmo 
     204                  CALL iom_get(inum, jpdom_auto, TRIM(tide_harmonics(itide)%cname_tide)//'_v1', ztr(:,:),cd_type='V',psgn=-1._wp) 
     205                  CALL iom_get(inum, jpdom_auto, TRIM(tide_harmonics(itide)%cname_tide)//'_v2', zti(:,:),cd_type='V',psgn=-1._wp) 
     206                  IF( ASSOCIATED(dta%v2d) ) THEN   ! we use bdy v2d on this mpi subdomain 
    207207                     DO ib = 1, SIZE(dta%v2d) 
    208208                        ii = idx_bdy(ib_bdy)%nbi(ib,igrd) 
     
    211211                        td%v0(ib,itide,2) = zti(ii,ij) 
    212212                     END DO 
    213                   END DO 
    214                   CALL iom_close( inum ) 
    215                ENDIF 
     213                  ENDIF  
     214               END DO 
     215               CALL iom_close( inum ) 
    216216               ! 
    217217               DEALLOCATE( ztr, zti )  
Note: See TracChangeset for help on using the changeset viewer.