Changeset 14168 for NEMO/releases/r4.0/r4.0-HEAD/src/OCE/BDY/bdytides.F90
- Timestamp:
- 2020-12-14T19:32:04+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/releases/r4.0/r4.0-HEAD/src/OCE/BDY/bdytides.F90
r12910 r14168 147 147 ! 148 148 ! SSH fields 149 IF( ASSOCIATED(dta%ssh) ) THEN ! we use bdy ssh on this mpi subdomain150 clfile = TRIM(filtide)//'_grid_T.nc'151 CALL iom_open( clfile , inum )152 igrd = 1 ! Everything is at T-points here153 DO itide = 1, nb_harmo154 CALL iom_get( inum, jpdom_autoglo, TRIM(Wave(ntide(itide))%cname_tide)//'_z1', ztr(:,:) )155 CALL iom_get( inum, jpdom_autoglo, TRIM(Wave(ntide(itide))%cname_tide)//'_z2', zti(:,:) )149 clfile = TRIM(filtide)//'_grid_T.nc' 150 CALL iom_open( clfile , inum ) 151 igrd = 1 ! Everything is at T-points here 152 DO itide = 1, nb_harmo 153 CALL iom_get( inum, jpdom_autoglo, TRIM(Wave(ntide(itide))%cname_tide)//'_z1', ztr(:,:) ) 154 CALL iom_get( inum, jpdom_autoglo, TRIM(Wave(ntide(itide))%cname_tide)//'_z2', zti(:,:) ) 155 IF( ASSOCIATED(dta%ssh) ) THEN ! we use bdy ssh on this mpi subdomain 156 156 DO ib = 1, SIZE(dta%ssh) 157 157 ii = idx_bdy(ib_bdy)%nbi(ib,igrd) … … 160 160 td%ssh0(ib,itide,2) = zti(ii,ij) 161 161 END DO 162 END DO163 CALL iom_close( inum )164 END IF162 ENDIF 163 END DO 164 CALL iom_close( inum ) 165 165 ! 166 166 ! U fields 167 IF( ASSOCIATED(dta%u2d) ) THEN ! we use bdy u2d on this mpi subdomain168 clfile = TRIM(filtide)//'_grid_U.nc'169 CALL iom_open( clfile , inum )170 igrd = 2 ! Everything is at U-points here171 DO itide = 1, nb_harmo172 CALL iom_get ( inum, jpdom_autoglo, TRIM(Wave(ntide(itide))%cname_tide)//'_u1', ztr(:,:) )173 CALL iom_get ( inum, jpdom_autoglo, TRIM(Wave(ntide(itide))%cname_tide)//'_u2', zti(:,:) )167 clfile = TRIM(filtide)//'_grid_U.nc' 168 CALL iom_open( clfile , inum ) 169 igrd = 2 ! Everything is at U-points here 170 DO itide = 1, nb_harmo 171 CALL iom_get ( inum, jpdom_autoglo, TRIM(Wave(ntide(itide))%cname_tide)//'_u1', ztr(:,:) ) 172 CALL iom_get ( inum, jpdom_autoglo, TRIM(Wave(ntide(itide))%cname_tide)//'_u2', zti(:,:) ) 173 IF( ASSOCIATED(dta%u2d) ) THEN ! we use bdy u2d on this mpi subdomain 174 174 DO ib = 1, SIZE(dta%u2d) 175 175 ii = idx_bdy(ib_bdy)%nbi(ib,igrd) … … 178 178 td%u0(ib,itide,2) = zti(ii,ij) 179 179 END DO 180 END DO181 CALL iom_close( inum )182 END IF180 ENDIF 181 END DO 182 CALL iom_close( inum ) 183 183 ! 184 184 ! V fields 185 IF( ASSOCIATED(dta%v2d) ) THEN ! we use bdy v2d on this mpi subdomain186 clfile = TRIM(filtide)//'_grid_V.nc'187 CALL iom_open( clfile , inum )188 igrd = 3 ! Everything is at V-points here189 DO itide = 1, nb_harmo190 CALL iom_get ( inum, jpdom_autoglo, TRIM(Wave(ntide(itide))%cname_tide)//'_v1', ztr(:,:) )191 CALL iom_get ( inum, jpdom_autoglo, TRIM(Wave(ntide(itide))%cname_tide)//'_v2', zti(:,:) )185 clfile = TRIM(filtide)//'_grid_V.nc' 186 CALL iom_open( clfile , inum ) 187 igrd = 3 ! Everything is at V-points here 188 DO itide = 1, nb_harmo 189 CALL iom_get ( inum, jpdom_autoglo, TRIM(Wave(ntide(itide))%cname_tide)//'_v1', ztr(:,:) ) 190 CALL iom_get ( inum, jpdom_autoglo, TRIM(Wave(ntide(itide))%cname_tide)//'_v2', zti(:,:) ) 191 IF( ASSOCIATED(dta%v2d) ) THEN ! we use bdy v2d on this mpi subdomain 192 192 DO ib = 1, SIZE(dta%v2d) 193 193 ii = idx_bdy(ib_bdy)%nbi(ib,igrd) … … 196 196 td%v0(ib,itide,2) = zti(ii,ij) 197 197 END DO 198 END DO199 CALL iom_close( inum )200 END IF198 ENDIF 199 END DO 200 CALL iom_close( inum ) 201 201 ! 202 202 DEALLOCATE( ztr, zti )
Note: See TracChangeset
for help on using the changeset viewer.