Changeset 13218
- Timestamp:
- 2020-07-02T12:07:27+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/ABL/ablmod.F90
r13214 r13218 348 348 ! 349 349 zztmp1 = pcd_du(ji, jj) 350 zztmp2 = 0.5_wp * pcd_du(ji, jj) * ( pssu(ji-1, jj) + pssu(ji,jj) ) * rn_vfac350 zztmp2 = 0.5_wp * pcd_du(ji, jj) * ( pssu(ji-1, jj) + pssu(ji,jj) ) 351 351 #if defined key_si3 352 352 zztmp1 = zztmp1 * pfrac_oce(ji,jj) + (1._wp - pfrac_oce(ji,jj)) * pcd_du_ice(ji, jj) 353 zzice = 0.5_wp * ( pssu_ice(ji-1, jj) + pssu_ice(ji, jj) ) * rn_vfac353 zzice = 0.5_wp * ( pssu_ice(ji-1, jj) + pssu_ice(ji, jj) ) 354 354 zztmp2 = zztmp2 * pfrac_oce(ji,jj) + (1._wp - pfrac_oce(ji,jj)) * pcd_du_ice(ji, jj) * zzice 355 355 #endif … … 422 422 ! 423 423 zztmp1 = pcd_du(ji, jj) 424 zztmp2 = 0.5_wp * pcd_du(ji, jj) * ( pssv(ji, jj) + pssv(ji, jj-1) ) * rn_vfac424 zztmp2 = 0.5_wp * pcd_du(ji, jj) * ( pssv(ji, jj) + pssv(ji, jj-1) ) 425 425 #if defined key_si3 426 426 zztmp1 = zztmp1 * pfrac_oce(ji,jj) + (1._wp - pfrac_oce(ji,jj)) * pcd_du_ice(ji, jj) 427 zzice = 0.5_wp * ( pssv_ice(ji, jj) + pssv_ice(ji, jj-1) ) * rn_vfac427 zzice = 0.5_wp * ( pssv_ice(ji, jj) + pssv_ice(ji, jj-1) ) 428 428 zztmp2 = zztmp2 * pfrac_oce(ji,jj) + (1._wp - pfrac_oce(ji,jj)) * pcd_du_ice(ji, jj) * zzice 429 429 #endif … … 590 590 591 591 DO_2D_01_01 592 <<<<<<< .working593 592 zwnd_i(ji,jj) = u_abl(ji ,jj,2,nt_a) - 0.5_wp * ( pssu(ji ,jj) + pssu(ji-1,jj) ) 594 593 zwnd_j(ji,jj) = v_abl(ji,jj ,2,nt_a) - 0.5_wp * ( pssv(ji,jj ) + pssv(ji,jj-1) ) 595 ||||||| .merge-left.r13136596 zwnd_i(ji,jj) = u_abl(ji ,jj,2,nt_a) - 0.5_wp * rn_vfac * ( pssu(ji ,jj) + pssu(ji-1,jj) )597 zwnd_j(ji,jj) = v_abl(ji,jj ,2,nt_a) - 0.5_wp * rn_vfac * ( pssv(ji,jj ) + pssv(ji,jj-1) )598 =======599 zwnd_i(ji,jj) = u_abl(ji ,jj,2,nt_a) - 0.5_wp * rn_vfac * ( pssu(ji,jj) + pssu(ji-1,jj ) )600 zwnd_j(ji,jj) = v_abl(ji,jj ,2,nt_a) - 0.5_wp * rn_vfac * ( pssv(ji,jj) + pssv(ji ,jj-1) )601 >>>>>>> .merge-right.r13211602 594 END_2D 603 595 ! … … 638 630 639 631 #if defined key_si3 640 <<<<<<< .working641 632 ! ------------------------------------------------------------ ! 642 633 ! Wind stress relative to the moving ice ( U10m - U_ice ) ! … … 652 643 IF(sn_cfctl%l_prtctl) CALL prt_ctl( tab2d_1=ptaui_ice , clinfo1=' abl_stp: putaui : ' & 653 644 & , tab2d_2=ptauj_ice , clinfo2=' pvtaui : ' ) 654 ||||||| .merge-left.r13136655 ! ------------------------------------------------------------ !656 ! Wind stress relative to the moving ice ( U10m - U_ice ) !657 ! ------------------------------------------------------------ !658 DO_2D_00_00659 660 zztmp1 = 0.5_wp * ( u_abl(ji+1,jj,2,nt_a) + u_abl(ji,jj,2,nt_a) )661 zztmp2 = 0.5_wp * ( v_abl(ji,jj+1,2,nt_a) + v_abl(ji,jj,2,nt_a) )662 663 ptaui_ice(ji,jj) = 0.5_wp * ( rhoa(ji+1,jj) * pCd_du_ice(ji+1,jj) &664 & + rhoa(ji ,jj) * pCd_du_ice(ji ,jj) ) &665 & * ( zztmp1 - rn_vfac * pssu_ice(ji,jj) )666 ptauj_ice(ji,jj) = 0.5_wp * ( rhoa(ji,jj+1) * pCd_du_ice(ji,jj+1) &667 & + rhoa(ji,jj ) * pCd_du_ice(ji,jj ) ) &668 & * ( zztmp2 - rn_vfac * pssv_ice(ji,jj) )669 END_2D670 CALL lbc_lnk_multi( 'ablmod', ptaui_ice, 'U', -1., ptauj_ice, 'V', -1. )671 !672 IF(sn_cfctl%l_prtctl) CALL prt_ctl( tab2d_1=ptaui_ice , clinfo1=' abl_stp: putaui : ' &673 & , tab2d_2=ptauj_ice , clinfo2=' pvtaui : ' )674 =======675 645 ! ------------------------------------------------------------ ! 676 646 ! Wind stress relative to the moving ice ( U10m - U_ice ) ! … … 683 653 ptaui_ice(ji,jj) = 0.5_wp * ( rhoa(ji+1,jj) * pCd_du_ice(ji+1,jj) & 684 654 & + rhoa(ji ,jj) * pCd_du_ice(ji ,jj) ) & 685 & * ( zztmp1 - rn_vfac *pssu_ice(ji,jj) )655 & * ( zztmp1 - pssu_ice(ji,jj) ) 686 656 ptauj_ice(ji,jj) = 0.5_wp * ( rhoa(ji,jj+1) * pCd_du_ice(ji,jj+1) & 687 657 & + rhoa(ji,jj ) * pCd_du_ice(ji,jj ) ) & 688 & * ( zztmp2 - rn_vfac *pssv_ice(ji,jj) )658 & * ( zztmp2 - pssv_ice(ji,jj) ) 689 659 END_2D 690 660 CALL lbc_lnk_multi( 'ablmod', ptaui_ice, 'U', -1., ptauj_ice, 'V', -1. ) … … 694 664 & tab2d_2=ptauj_ice , clinfo2=' vtau_ice : ', mask2=vmask ) 695 665 END IF 696 >>>>>>> .merge-right.r13211697 666 #endif 698 667 ! !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
Note: See TracChangeset
for help on using the changeset viewer.