- Timestamp:
- 2021-04-29T16:25:38+02: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
r14680 r14765 647 647 ! ! Compute lateral diffusive coefficient at T-point 648 648 IF( ln_traldf_triad ) THEN 649 ! [comm_cleanup] ! DO_3D( 0, 0, 0, 0, 1, jpk )650 649 DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 1, jpk ) 651 650 ! Take the max of N^2 and zero then take the vertical sum … … 662 661 END_3D 663 662 ELSE 664 ! [comm_cleanup] ! DO_3D( 0, 0, 0, 0, 1, jpk )665 663 DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 1, jpk ) 666 664 ! Take the max of N^2 and zero then take the vertical sum … … 679 677 ENDIF 680 678 681 ! [comm_cleanup] ! DO_2D( 0, 0, 0, 0 )682 679 DO_2D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1 ) 683 680 zfw = MAX( ABS( 2. * omega * SIN( rad * gphit(ji,jj) ) ) , 1.e-10 ) … … 690 687 ! !== Bound on eiv coeff. ==! 691 688 z1_f20 = 1._wp / ( 2._wp * omega * sin( rad * 20._wp ) ) 692 ! [comm_cleanup] ! DO_2D( 0, 0, 0, 0 )693 689 DO_2D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1 ) 694 690 zzaei = MIN( 1._wp, ABS( ff_t(ji,jj) * z1_f20 ) ) * zaeiw(ji,jj) ! tropical decrease … … 697 693 CALL lbc_lnk( 'ldftra', zaeiw(:,:), 'W', 1.0_wp ) ! lateral boundary condition 698 694 ! 699 ! [comm_cleanup] ! DO_2D( 0, 0, 0, 0 )700 695 DO_2D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1 ) 701 696 paeiu(ji,jj,1) = 0.5_wp * ( zaeiw(ji,jj) + zaeiw(ji+1,jj ) ) * umask(ji,jj,1) … … 757 752 zpsi_uw(:,:,jpk) = 0._wp ; zpsi_vw(:,:,jpk) = 0._wp 758 753 ! 759 ! [comm_cleanup] ! DO_3D( 1, 0, 1, 0, 2, jpkm1 )760 754 DO_3D( nn_hls, nn_hls-1, nn_hls, nn_hls-1, 2, jpkm1 ) 761 755 zpsi_uw(ji,jj,jk) = - r1_4 * e2u(ji,jj) * ( wslpi(ji,jj,jk ) + wslpi(ji+1,jj,jk) ) & … … 765 759 END_3D 766 760 ! 767 ! [comm_cleanup] ! DO_3D( 1, 0, 1, 0, 1, jpkm1 )768 761 DO_3D_OVR( nn_hls, nn_hls-1, nn_hls, nn_hls-1, 1, jpkm1 ) 769 762 pu(ji,jj,jk) = pu(ji,jj,jk) - ( zpsi_uw(ji,jj,jk) - zpsi_uw(ji,jj,jk+1) ) 770 763 pv(ji,jj,jk) = pv(ji,jj,jk) - ( zpsi_vw(ji,jj,jk) - zpsi_vw(ji,jj,jk+1) ) 771 764 END_3D 772 ! [comm_cleanup] ! DO_3D( 0, 0, 0, 0, 1, jpkm1 )773 765 DO_3D_OVR( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 1, jpkm1 ) 774 766 pw(ji,jj,jk) = pw(ji,jj,jk) + ( zpsi_uw(ji,jj,jk) - zpsi_uw(ji-1,jj ,jk) &
Note: See TracChangeset
for help on using the changeset viewer.