Changeset 14200 for NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/OCE/BDY/bdytides.F90
- Timestamp:
- 2020-12-17T15:36:44+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/OCE/BDY/bdytides.F90
r13286 r14200 162 162 ! 163 163 ! SSH fields 164 IF( ASSOCIATED(dta%ssh) ) THEN ! we use bdy ssh on this mpi subdomain165 clfile = TRIM(filtide)//'_grid_T.nc'166 CALL iom_open( clfile , inum )167 igrd = 1 ! Everything is at T-points here168 DO itide = 1, nb_harmo169 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 171 171 DO ib = 1, SIZE(dta%ssh) 172 172 ii = idx_bdy(ib_bdy)%nbi(ib,igrd) … … 175 175 td%ssh0(ib,itide,2) = zti(ii,ij) 176 176 END DO 177 END DO178 CALL iom_close( inum )179 ENDIF177 ENDIF 178 END DO 179 CALL iom_close( inum ) 180 180 ! 181 181 ! U fields 182 IF( ASSOCIATED(dta%u2d) ) THEN ! we use bdy u2d on this mpi subdomain183 clfile = TRIM(filtide)//'_grid_U.nc'184 CALL iom_open( clfile , inum )185 igrd = 2 ! Everything is at U-points here186 DO itide = 1, nb_harmo187 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 189 189 DO ib = 1, SIZE(dta%u2d) 190 190 ii = idx_bdy(ib_bdy)%nbi(ib,igrd) … … 193 193 td%u0(ib,itide,2) = zti(ii,ij) 194 194 END DO 195 END DO196 CALL iom_close( inum )197 ENDIF195 ENDIF 196 END DO 197 CALL iom_close( inum ) 198 198 ! 199 199 ! V fields 200 IF( ASSOCIATED(dta%v2d) ) THEN ! we use bdy v2d on this mpi subdomain201 clfile = TRIM(filtide)//'_grid_V.nc'202 CALL iom_open( clfile , inum )203 igrd = 3 ! Everything is at V-points here204 DO itide = 1, nb_harmo205 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 207 207 DO ib = 1, SIZE(dta%v2d) 208 208 ii = idx_bdy(ib_bdy)%nbi(ib,igrd) … … 211 211 td%v0(ib,itide,2) = zti(ii,ij) 212 212 END DO 213 END DO214 CALL iom_close( inum )215 ENDIF213 ENDIF 214 END DO 215 CALL iom_close( inum ) 216 216 ! 217 217 DEALLOCATE( ztr, zti )
Note: See TracChangeset
for help on using the changeset viewer.