- Timestamp:
- 2016-10-18T15:32:04+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90
r6748 r7037 325 325 ! need to be like this to compute the pressure gradient with ISF. If not, level beneath the ISF are not aligned (sum(e3t) /= depth) 326 326 ! define e3t_0 and e3w_0 as the differences between gdept and gdepw respectively 327 !$OMP PARALLEL DO schedule(static) private(jk)328 327 DO jk = 1, jpkm1 329 328 e3t_1d(jk) = gdepw_1d(jk+1)-gdepw_1d(jk) 330 329 END DO 331 330 e3t_1d(jpk) = e3t_1d(jpk-1) ! we don't care because this level is masked in NEMO 332 !$OMP PARALLEL DO schedule(static) private(jk)333 331 DO jk = 2, jpk 334 332 e3w_1d(jk) = gdept_1d(jk) - gdept_1d(jk-1) … … 429 427 !$OMP END PARALLEL WORKSHARE 430 428 ELSE ! z-coordinate (zco or zps): step-like topography 431 !$OMP PARALLEL WORKSHARE 429 !$OMP PARALLEL 430 !$OMP WORKSHARE 432 431 idta(:,:) = jpkm1 433 !$OMP END PARALLELWORKSHARE432 !$OMP END WORKSHARE 434 433 DO jk = 1, jpkm1 434 !$OMP WORKSHARE 435 435 WHERE( gdept_1d(jk) < zdta(:,:) .AND. zdta(:,:) <= gdept_1d(jk+1) ) idta(:,:) = jk 436 !$OMP END WORKSHARE 436 437 END DO 438 !$OMP END PARALLEL 437 439 ENDIF 438 440 ELSE … … 506 508 507 509 ! ! local domain level and meter bathymetries (mbathy,bathy) 508 !$OMP PARALLEL WORKSHARE 510 !$OMP PARALLEL 511 !$OMP WORKSHARE 509 512 mbathy(:,:) = 0 ! set to zero extra halo points 510 513 bathy (:,:) = 0._wp ! (require for mpp case) 511 !$OMP END PARALLELWORKSHARE512 !$OMP PARALLELDO schedule(static) private(jj, ji)514 !$OMP END WORKSHARE 515 !$OMP DO schedule(static) private(jj, ji) 513 516 DO jj = 1, nlcj ! interior values 514 517 DO ji = 1, nlci … … 517 520 END DO 518 521 END DO 519 !$OMP PARALLEL WORKSHARE 522 !$OMP END DO NOWAIT 523 !$OMP WORKSHARE 520 524 risfdep(:,:)=0.e0 521 525 misfdep(:,:)=1 522 !$OMP END PARALLEL WORKSHARE 526 !$OMP END WORKSHARE NOWAIT 527 !$OMP END PARALLEL 523 528 ! 524 529 DEALLOCATE( idta, zdta ) … … 532 537 CALL iom_get ( inum, jpdom_data, 'Bathy_level', bathy ) 533 538 CALL iom_close( inum ) 534 !$OMP PARALLEL WORKSHARE535 539 mbathy(:,:) = INT( bathy(:,:) ) 536 540 ! initialisation isf variables 541 !$OMP PARALLEL WORKSHARE 537 542 risfdep(:,:)=0._wp ; misfdep(:,:)=1 538 543 !$OMP END PARALLEL WORKSHARE … … 1196 1201 END DO 1197 1202 ELSE ! no cavity 1203 !$OMP PARALLEL 1204 !$OMP WORKSHARE 1198 1205 gde3w_0(:,:,1) = 0.5_wp * e3w_0(:,:,1) 1206 !$OMP END WORKSHARE 1199 1207 DO jk = 2, jpk 1200 !$OMP PARALLELDO schedule(static) private(jj, ji)1208 !$OMP DO schedule(static) private(jj, ji) 1201 1209 DO jj =1, jpj 1202 1210 DO ji=1, jpi … … 1205 1213 END DO 1206 1214 END DO 1215 !$OMP END PARALLEL 1207 1216 END IF 1208 1217 !
Note: See TracChangeset
for help on using the changeset viewer.