- Timestamp:
- 2020-01-27T15:31:53+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11943_MERGE_2019/src/ICE/iceitd.F90
r11960 r12340 48 48 REAL(wp), DIMENSION(0:100) :: rn_catbnd ! ice categories bounds 49 49 ! 50 !! * Substitutions 51 # include "do_loop_substitute.h90" 50 52 !!---------------------------------------------------------------------- 51 53 !! NEMO/ICE 4.0 , NEMO Consortium (2018) … … 96 98 ! 97 99 npti = 0 ; nptidx(:) = 0 98 DO jj = 1, jpj 99 DO ji = 1, jpi 100 IF ( at_i(ji,jj) > epsi10 ) THEN 101 npti = npti + 1 102 nptidx( npti ) = (jj - 1) * jpi + ji 103 ENDIF 104 END DO 105 END DO 100 DO_2D_11_11 101 IF ( at_i(ji,jj) > epsi10 ) THEN 102 npti = npti + 1 103 nptidx( npti ) = (jj - 1) * jpi + ji 104 ENDIF 105 END_2D 106 106 107 107 !----------------------------------------------------------------------------------------------- … … 597 597 ! !--------------------------------------- 598 598 npti = 0 ; nptidx(:) = 0 599 DO jj = 1, jpj 600 DO ji = 1, jpi 601 IF( a_i(ji,jj,jl) > 0._wp .AND. v_i(ji,jj,jl) > (a_i(ji,jj,jl) * hi_max(jl)) ) THEN 602 npti = npti + 1 603 nptidx( npti ) = (jj - 1) * jpi + ji 604 ENDIF 605 END DO 606 END DO 599 DO_2D_11_11 600 IF( a_i(ji,jj,jl) > 0._wp .AND. v_i(ji,jj,jl) > (a_i(ji,jj,jl) * hi_max(jl)) ) THEN 601 npti = npti + 1 602 nptidx( npti ) = (jj - 1) * jpi + ji 603 ENDIF 604 END_2D 607 605 ! 608 606 !!clem CALL tab_2d_1d( npti, nptidx(1:npti), h_i_1d(1:npti), h_i(:,:,jl) ) … … 638 636 ! !----------------------------------------- 639 637 npti = 0 ; nptidx(:) = 0 640 DO jj = 1, jpj 641 DO ji = 1, jpi 642 IF( a_i(ji,jj,jl+1) > 0._wp .AND. v_i(ji,jj,jl+1) <= (a_i(ji,jj,jl+1) * hi_max(jl)) ) THEN 643 npti = npti + 1 644 nptidx( npti ) = (jj - 1) * jpi + ji 645 ENDIF 646 END DO 647 END DO 638 DO_2D_11_11 639 IF( a_i(ji,jj,jl+1) > 0._wp .AND. v_i(ji,jj,jl+1) <= (a_i(ji,jj,jl+1) * hi_max(jl)) ) THEN 640 npti = npti + 1 641 nptidx( npti ) = (jj - 1) * jpi + ji 642 ENDIF 643 END_2D 648 644 ! 649 645 CALL tab_2d_1d( npti, nptidx(1:npti), a_i_1d(1:npti), a_i(:,:,jl+1) ) ! jl+1 is ok
Note: See TracChangeset
for help on using the changeset viewer.