- Timestamp:
- 2021-03-03T17:04:07+01:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2021/dev_r14273_HPC-02_Daley_Tiling/src/OCE/LDF/ldftra.F90
r14574 r14576 647 647 ! ! Compute lateral diffusive coefficient at T-point 648 648 IF( ln_traldf_triad ) THEN 649 DO_3D( 0, 0, 0, 0, 1, jpk ) 649 ! [comm_cleanup] ! DO_3D( 0, 0, 0, 0, 1, jpk ) 650 DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 1, jpk ) 650 651 ! Take the max of N^2 and zero then take the vertical sum 651 652 ! of the square root of the resulting N^2 ( required to compute … … 661 662 END_3D 662 663 ELSE 663 DO_3D( 0, 0, 0, 0, 1, jpk ) 664 ! [comm_cleanup] ! DO_3D( 0, 0, 0, 0, 1, jpk ) 665 DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 1, jpk ) 664 666 ! Take the max of N^2 and zero then take the vertical sum 665 667 ! of the square root of the resulting N^2 ( required to compute … … 677 679 ENDIF 678 680 679 DO_2D( 0, 0, 0, 0 ) 681 ! [comm_cleanup] ! DO_2D( 0, 0, 0, 0 ) 682 DO_2D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1 ) 680 683 zfw = MAX( ABS( 2. * omega * SIN( rad * gphit(ji,jj) ) ) , 1.e-10 ) 681 684 ! Rossby radius at w-point taken betwenn 2 km and 40km … … 687 690 ! !== Bound on eiv coeff. ==! 688 691 z1_f20 = 1._wp / ( 2._wp * omega * sin( rad * 20._wp ) ) 689 DO_2D( 0, 0, 0, 0 ) 692 ! [comm_cleanup] ! DO_2D( 0, 0, 0, 0 ) 693 DO_2D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1 ) 690 694 zzaei = MIN( 1._wp, ABS( ff_t(ji,jj) * z1_f20 ) ) * zaeiw(ji,jj) ! tropical decrease 691 695 zaeiw(ji,jj) = MIN( zzaei , paei0 ) ! Max value = paei0 … … 693 697 CALL lbc_lnk( 'ldftra', zaeiw(:,:), 'W', 1.0_wp ) ! lateral boundary condition 694 698 ! 695 DO_2D( 0, 0, 0, 0 ) !== aei at u- and v-points ==! 699 ! [comm_cleanup] ! DO_2D( 0, 0, 0, 0 ) 700 DO_2D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1 ) 696 701 paeiu(ji,jj,1) = 0.5_wp * ( zaeiw(ji,jj) + zaeiw(ji+1,jj ) ) * umask(ji,jj,1) 697 702 paeiv(ji,jj,1) = 0.5_wp * ( zaeiw(ji,jj) + zaeiw(ji ,jj+1) ) * vmask(ji,jj,1) … … 751 756 zpsi_uw(:,:,jpk) = 0._wp ; zpsi_vw(:,:,jpk) = 0._wp 752 757 ! 753 DO_3D( 1, 0, 1, 0, 2, jpkm1 ) 758 ! [comm_cleanup] ! DO_3D( 1, 0, 1, 0, 2, jpkm1 ) 759 DO_3D( nn_hls, nn_hls-1, nn_hls, nn_hls-1, 2, jpkm1 ) 754 760 zpsi_uw(ji,jj,jk) = - r1_4 * e2u(ji,jj) * ( wslpi(ji,jj,jk ) + wslpi(ji+1,jj,jk) ) & 755 761 & * ( aeiu (ji,jj,jk-1) + aeiu (ji ,jj,jk) ) * wumask(ji,jj,jk) … … 758 764 END_3D 759 765 ! 760 DO_3D( 1, 0, 1, 0, 1, jpkm1 ) 766 ! [comm_cleanup] ! DO_3D( 1, 0, 1, 0, 1, jpkm1 ) 767 DO_3D( nn_hls, nn_hls-1, nn_hls, nn_hls-1, 1, jpkm1 ) 761 768 pu(ji,jj,jk) = pu(ji,jj,jk) - ( zpsi_uw(ji,jj,jk) - zpsi_uw(ji,jj,jk+1) ) 762 769 pv(ji,jj,jk) = pv(ji,jj,jk) - ( zpsi_vw(ji,jj,jk) - zpsi_vw(ji,jj,jk+1) ) 763 770 END_3D 764 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 771 ! [comm_cleanup] ! DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 772 DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 1, jpkm1 ) 765 773 pw(ji,jj,jk) = pw(ji,jj,jk) + ( zpsi_uw(ji,jj,jk) - zpsi_uw(ji-1,jj ,jk) & 766 774 & + zpsi_vw(ji,jj,jk) - zpsi_vw(ji ,jj-1,jk) )
Note: See TracChangeset
for help on using the changeset viewer.