Changeset 7753 for trunk/NEMOGCM/NEMO/LIM_SRC_3/limitd_th.F90
- Timestamp:
- 2017-03-03T12:46:59+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/LIM_SRC_3/limitd_th.F90
r7698 r7753 106 106 CALL lim_column_sum (jpl, v_s, vt_s_init) 107 107 CALL lim_column_sum_energy (jpl, nlay_i, e_i, et_i_init) 108 DO jl = 1, jpl 109 !$OMP PARALLEL DO schedule(static) private(jj,ji) 110 DO jj = 1, jpj 111 DO ji = 1, jpi 112 dummy_es(ji,jj,jl) = e_s(ji,jj,1,jl) 113 END DO 114 END DO 115 END DO 108 dummy_es(:,:,:) = e_s(:,:,1,:) 116 109 CALL lim_column_sum (jpl, dummy_es(:,:,:) , et_s_init) 117 110 ENDIF … … 128 121 ENDIF 129 122 130 !$OMP PARALLEL 131 DO jl = 1, jpl 132 !$OMP DO schedule(static) private(jj,ji) 133 DO jj = 1, jpj 134 DO ji = 1, jpi 135 zdhice(ji,jj,jl) = 0._wp 136 END DO 137 END DO 138 END DO 123 zdhice(:,:,:) = 0._wp 139 124 DO jl = klbnd, kubnd 140 !$OMP DO schedule(static) private(jj,ji,rswitch)141 125 DO jj = 1, jpj 142 126 DO ji = 1, jpi … … 153 137 ! 2) Compute fractional ice area in each grid cell 154 138 !----------------------------------------------------------------------------------------------- 155 !$OMP DO schedule(static) private(jj,ji) 156 DO jj = 1, jpj 157 DO ji = 1, jpi 158 at_i(ji,jj) = 0._wp 159 END DO 160 END DO 139 at_i(:,:) = 0._wp 161 140 DO jl = klbnd, kubnd 162 !$OMP DO schedule(static) private(jj,ji) 163 DO jj = 1, jpj 164 DO ji = 1, jpi 165 at_i(ji,jj) = at_i(ji,jj) + a_i(ji,jj,jl) 166 END DO 167 END DO 168 END DO 169 !$OMP END PARALLEL 141 at_i(:,:) = at_i(:,:) + a_i(:,:,jl) 142 END DO 170 143 171 144 !----------------------------------------------------------------------------------------------- … … 190 163 !----------------------------------------------------------------------------------------------- 191 164 !- 4.1 Compute category boundaries 192 !$OMP PARALLEL 193 DO jl = 0, jpl 194 !$OMP DO schedule(static) private(jj,ji) 195 DO jj = 1, jpj 196 DO ji = 1, jpi 197 zhbnew(ji,jj,jl) = 0._wp 198 END DO 199 END DO 200 END DO 165 zhbnew(:,:,:) = 0._wp 201 166 202 167 DO jl = klbnd, kubnd - 1 203 !$OMP DO schedule(static) private(ji,ii,ij,zslope)204 168 DO ji = 1, nbrem 205 169 ii = nind_i(ji) … … 219 183 220 184 !- 4.2 Check that each zhbnew lies between adjacent values of ice thickness 221 !$OMP DO schedule(static) private(ji,ii,ij)222 185 DO ji = 1, nbrem 223 186 ii = nind_i(ji) … … 242 205 243 206 END DO 244 !$OMP END PARALLEL245 207 246 208 !----------------------------------------------------------------------------------------------- … … 261 223 ! 6) Fill arrays with lowermost / uppermost boundaries of 'new' categories 262 224 !----------------------------------------------------------------------------------------------- 263 !$OMP PARALLEL DO schedule(static) private(jj,ji)264 225 DO jj = 1, jpj 265 226 DO ji = 1, jpi … … 293 254 294 255 !- 7.2 Area lost due to melting of thin ice (first category, klbnd) 295 !$OMP PARALLEL DO schedule(static) private(ji,ii,ij,zdh0,zetamax,zx1,zx2,zda0,zdamax)296 256 DO ji = 1, nbrem 297 257 ii = nind_i(ji) … … 339 299 !----------------------------------------------------------------------------------------------- 340 300 341 !$OMP PARALLEL342 301 DO jl = klbnd, kubnd - 1 343 !$OMP DO schedule(static) private(jj,ji)344 302 DO jj = 1, jpj 345 303 DO ji = 1, jpi … … 350 308 END DO 351 309 352 !$OMP DO schedule(static) private(ji,ii,ij,zetamax,zetamin,zx1,zwk1,zwk2,zx2,zx3,nd)353 310 DO ji = 1, nbrem 354 311 ii = nind_i(ji) … … 385 342 END DO 386 343 END DO 387 !$OMP END PARALLEL388 344 389 345 !!---------------------------------------------------------------------------------------------- … … 396 352 !!---------------------------------------------------------------------------------------------- 397 353 398 !$OMP PARALLEL DO schedule(static) private(ji,ii,ij)399 354 DO ji = 1, nbrem 400 355 ii = nind_i(ji) … … 422 377 CALL lim_cons_check (vt_s_init, vt_s_final, 1.0e-6, fieldid) 423 378 424 DO jl = 1, jpl 425 !$OMP PARALLEL DO schedule(static) private(jj,ji) 426 DO jj = 1, jpj 427 DO ji = 1, jpi 428 dummy_es(ji,jj,jl) = e_s(ji,jj,1,jl) 429 END DO 430 END DO 431 END DO 379 dummy_es(:,:,:) = e_s(:,:,1,:) 432 380 CALL lim_column_sum (jpl, dummy_es(:,:,:) , et_s_final) 433 381 fieldid = ' e_s : limitd_th ' … … 473 421 !!------------------------------------------------------------------ 474 422 ! 475 !$OMP PARALLEL DO schedule(static) private(jj,ji,zh13,zh23,zdhr,zwk1,zwk2)476 423 DO jj = 1, jpj 477 424 DO ji = 1, jpi … … 553 500 554 501 DO jl = klbnd, kubnd 555 !$OMP PARALLEL DO schedule(static) private(jj,ji) 556 DO jj = 1, jpj 557 DO ji = 1, jpi 558 zaTsfn(ji,jj,jl) = a_i(ji,jj,jl) * t_su(ji,jj,jl) 559 END DO 560 END DO 502 zaTsfn(:,:,jl) = a_i(:,:,jl) * t_su(:,:,jl) 561 503 END DO 562 504 … … 577 519 END DO 578 520 579 !$OMP PARALLEL DO schedule(static) private(ji,ii,ij,jl1,jl2,rswitch,zdvsnow,zdesnow,zdo_aice,zdsm_vice,zdaTsf)580 521 DO ji = 1, nbrem 581 522 ii = nind_i(ji) … … 643 584 644 585 DO jk = 1, nlay_i 645 !$OMP PARALLEL DO schedule(static) private(ji,ii,ij,jl1,jl2,zdeice)646 586 DO ji = 1, nbrem 647 587 ii = nind_i(ji) … … 668 608 669 609 DO jl = klbnd, kubnd 670 !$OMP PARALLEL DO schedule(static) private(jj,ji)671 610 DO jj = 1, jpj 672 611 DO ji = 1, jpi … … 724 663 ! 1) Compute ice thickness. 725 664 !------------------------------------------------------------------------------ 726 !$OMP PARALLEL727 665 DO jl = klbnd, kubnd 728 !$OMP DO schedule(static) private(jj,ji,rswitch)729 666 DO jj = 1, jpj 730 667 DO ji = 1, jpi … … 743 680 !------------------------- 744 681 DO jl = klbnd, kubnd 745 !$OMP DO schedule(static) private(jj,ji) 746 DO jj = 1, jpj 747 DO ji = 1, jpi 748 zdonor(ji,jj,jl) = 0 749 zdaice(ji,jj,jl) = 0._wp 750 zdvice(ji,jj,jl) = 0._wp 751 END DO 752 END DO 753 END DO 754 !$OMP END PARALLEL 682 zdonor(:,:,jl) = 0 683 zdaice(:,:,jl) = 0._wp 684 zdvice(:,:,jl) = 0._wp 685 END DO 755 686 756 687 !------------------------- … … 765 696 zshiftflag = 0 766 697 767 !$OMP PARALLEL DO schedule(static) private(jj,ji) REDUCTION(MAX:zshiftflag)768 698 DO jj = 1, jpj 769 699 DO ji = 1, jpi … … 786 716 CALL lim_itd_shiftice( klbnd, kubnd, zdonor, zdaice, zdvice ) 787 717 ! Reset shift parameters 788 !$OMP PARALLEL DO schedule(static) private(jj,ji) 789 DO jj = 1, jpj 790 DO ji = 1, jpi 791 zdonor(ji,jj,jl) = 0 792 zdaice(ji,jj,jl) = 0._wp 793 zdvice(ji,jj,jl) = 0._wp 794 END DO 795 END DO 718 zdonor(:,:,jl) = 0 719 zdaice(:,:,jl) = 0._wp 720 zdvice(:,:,jl) = 0._wp 796 721 ENDIF 797 722 ! … … 809 734 zshiftflag = 0 810 735 811 !$OMP PARALLEL DO schedule(static) private(jj,ji) REDUCTION(MAX:zshiftflag)812 736 DO jj = 1, jpj 813 737 DO ji = 1, jpi 814 738 IF( a_i(ji,jj,jl+1) > epsi10 .AND. ht_i(ji,jj,jl+1) <= hi_max(jl) ) THEN 739 ! 815 740 zshiftflag = 1 816 741 zdonor(ji,jj,jl) = jl + 1 … … 826 751 CALL lim_itd_shiftice( klbnd, kubnd, zdonor, zdaice, zdvice ) 827 752 ! Reset shift parameters 828 !$OMP PARALLEL DO schedule(static) private(jj,ji) 829 DO jj = 1, jpj 830 DO ji = 1, jpi 831 zdonor(ji,jj,jl) = 0 832 zdaice(ji,jj,jl) = 0._wp 833 zdvice(ji,jj,jl) = 0._wp 834 END DO 835 END DO 753 zdonor(:,:,jl) = 0 754 zdaice(:,:,jl) = 0._wp 755 zdvice(:,:,jl) = 0._wp 836 756 ENDIF 837 757
Note: See TracChangeset
for help on using the changeset viewer.