- 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/LIM_SRC_2/limthd_2.F90
r6140 r7037 121 121 122 122 !!gm needed? yes at least for some of these arrays 123 !$OMP PARALLEL 124 !$OMP WORKSHARE 123 125 zdvosif(:,:) = 0.e0 ! variation of ice volume at surface 124 126 zdvobif(:,:) = 0.e0 ! variation of ice volume at bottom … … 137 139 rdq_ice(:,:) = 0.e0 ! heat content associated with rdm_ice 138 140 zmsk (:,:,:) = 0.e0 141 !$OMP END WORKSHARE NOWAIT 139 142 140 143 ! set to zero snow thickness smaller than epsi04 144 !$OMP DO schedule(static) private(jj,ji) 141 145 DO jj = 1, jpj 142 146 DO ji = 1, jpi … … 144 148 END DO 145 149 END DO 150 !$OMP END DO NOWAIT 151 !$OMP END PARALLEL 146 152 !!gm better coded (do not use SIGN...) 147 153 ! WHERE( hsnif(:,:) < epsi04 ) hsnif(:,:) = 0.e0 … … 154 160 !-----------------------------------! 155 161 162 !$OMP PARALLEL DO schedule(static) private(jj,ji,zindg,zthsnice,zindb,za,zh) 156 163 DO jj = 1, jpj 157 164 DO ji = 1, jpi … … 196 203 !-------------------------------------------------------------------------- 197 204 205 !$OMP PARALLEL DO schedule(static) private(jj,ji,zthsnice,zindb,zinda,zfricp,zfric_u,zfntlat,zpareff) 198 206 DO jj = 1, jpj 199 207 DO ji = 1, jpi … … 337 345 ! Up-date sea ice thickness 338 346 !-------------------------- 347 !$OMP PARALLEL DO schedule(static) private(jj,ji) 339 348 DO jj = 1, jpj 340 349 DO ji = 1, jpi … … 348 357 !-------------------------------------------------------- 349 358 IF( fcor(1,1) < 0.e0 ) THEN 359 !$OMP PARALLEL DO schedule(static) private(jj,ji) 350 360 DO jj = 1, njeqm1 351 361 DO ji = 1, jpi … … 418 428 ! Update daily thermodynamic ice production. 419 429 !------------------------------------------------------------------------------ 430 !$OMP PARALLEL DO schedule(static) private(jj,ji) 420 431 DO jj = 1, jpj 421 432 DO ji = 1, jpi … … 453 464 IF( iom_use('u_imasstr') ) THEN 454 465 zztmp = 0.25 * rhoic 466 !$OMP PARALLEL DO schedule(static) private(jj,ji,zuice_m,zhice_u) 455 467 DO jj = 1, jpjm1 456 468 DO ji = 1, jpim1 ! NO vector opt. … … 466 478 IF( iom_use('v_imasstr') ) THEN 467 479 zztmp = 0.25 * rhoic 480 !$OMP PARALLEL DO schedule(static) private(jj,ji,zvice_m,zhice_v) 468 481 DO jj = 1, jpjm1 469 482 DO ji = 1, jpim1 ! NO vector opt. … … 480 493 !! Fram Strait sea-ice transport (sea-ice + snow) (in ORCA2 = 5 points) 481 494 IF( iom_use('fram_trans') .and. cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN ! ORCA R2 configuration 495 !$OMP PARALLEL DO schedule(static) private(jj,ji,zrhoij,zrhoijm1,ztr_fram) 482 496 DO jj = mj0(137), mj1(137) ! B grid 483 497 IF( mj0(jj-1) >= nldj ) THEN … … 497 511 !! ce A big warning because the model crashes on IDRIS/IBM SP6 with xlf 13.1.0.3, see ticket #761 498 512 !! ce We Unroll the loop and everything works fine 513 !$OMP PARALLEL DO schedule(static) private(jj,ji) 499 514 DO jj = 1, jpj 500 515 DO ji = 1, jpi
Note: See TracChangeset
for help on using the changeset viewer.