New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 13218 for NEMO/trunk/src/ABL – NEMO

Changeset 13218 for NEMO/trunk/src/ABL


Ignore:
Timestamp:
2020-07-02T12:07:27+02:00 (4 years ago)
Author:
smasson
Message:

trunk: fix merge issues in ablmod

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/ABL/ablmod.F90

    r13214 r13218  
    348348            ! 
    349349            zztmp1  = pcd_du(ji, jj) 
    350             zztmp2  = 0.5_wp * pcd_du(ji, jj) * ( pssu(ji-1, jj) + pssu(ji,jj) ) * rn_vfac 
     350            zztmp2  = 0.5_wp * pcd_du(ji, jj) * ( pssu(ji-1, jj) + pssu(ji,jj) ) 
    351351#if defined key_si3 
    352352            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_vfac 
     353            zzice  = 0.5_wp * ( pssu_ice(ji-1, jj) + pssu_ice(ji, jj) ) 
    354354            zztmp2 = zztmp2 * pfrac_oce(ji,jj) + (1._wp - pfrac_oce(ji,jj)) * pcd_du_ice(ji, jj) * zzice 
    355355#endif 
     
    422422            ! 
    423423            zztmp1 = pcd_du(ji, jj) 
    424             zztmp2 = 0.5_wp * pcd_du(ji, jj) * ( pssv(ji, jj) + pssv(ji, jj-1) ) * rn_vfac 
     424            zztmp2 = 0.5_wp * pcd_du(ji, jj) * ( pssv(ji, jj) + pssv(ji, jj-1) ) 
    425425#if defined key_si3 
    426426            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_vfac 
     427            zzice  = 0.5_wp * ( pssv_ice(ji, jj) + pssv_ice(ji, jj-1) ) 
    428428            zztmp2 = zztmp2 * pfrac_oce(ji,jj) + (1._wp - pfrac_oce(ji,jj)) * pcd_du_ice(ji, jj) * zzice 
    429429#endif 
     
    590590 
    591591      DO_2D_01_01 
    592 <<<<<<< .working 
    593592         zwnd_i(ji,jj) = u_abl(ji  ,jj,2,nt_a) - 0.5_wp * ( pssu(ji  ,jj) + pssu(ji-1,jj) )   
    594593         zwnd_j(ji,jj) = v_abl(ji,jj  ,2,nt_a) - 0.5_wp * ( pssv(ji,jj  ) + pssv(ji,jj-1) )  
    595 ||||||| .merge-left.r13136 
    596          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.r13211 
    602594      END_2D 
    603595      ! 
     
    638630 
    639631#if defined key_si3 
    640 <<<<<<< .working 
    641632      ! ------------------------------------------------------------ ! 
    642633      !    Wind stress relative to the moving ice ( U10m - U_ice )   ! 
     
    652643      IF(sn_cfctl%l_prtctl)   CALL prt_ctl( tab2d_1=ptaui_ice  , clinfo1=' abl_stp: putaui : '   & 
    653644         &                                , tab2d_2=ptauj_ice  , clinfo2='          pvtaui : ' ) 
    654 ||||||| .merge-left.r13136 
    655          ! ------------------------------------------------------------ ! 
    656          !    Wind stress relative to the moving ice ( U10m - U_ice )   ! 
    657          ! ------------------------------------------------------------ ! 
    658          DO_2D_00_00 
    659              
    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_2D 
    670          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 ======= 
    675645      ! ------------------------------------------------------------ ! 
    676646      !    Wind stress relative to the moving ice ( U10m - U_ice )   ! 
     
    683653         ptaui_ice(ji,jj) = 0.5_wp * (  rhoa(ji+1,jj) * pCd_du_ice(ji+1,jj)             & 
    684654            &                      +    rhoa(ji  ,jj) * pCd_du_ice(ji  ,jj)  )          & 
    685             &         * ( zztmp1 - rn_vfac * pssu_ice(ji,jj) ) 
     655            &         * ( zztmp1 - pssu_ice(ji,jj) ) 
    686656         ptauj_ice(ji,jj) = 0.5_wp * (  rhoa(ji,jj+1) * pCd_du_ice(ji,jj+1)             & 
    687657            &                      +    rhoa(ji,jj  ) * pCd_du_ice(ji,jj  )  )          & 
    688             &         * ( zztmp2 - rn_vfac * pssv_ice(ji,jj) ) 
     658            &         * ( zztmp2 - pssv_ice(ji,jj) ) 
    689659      END_2D 
    690660      CALL lbc_lnk_multi( 'ablmod', ptaui_ice, 'U', -1., ptauj_ice, 'V', -1. ) 
     
    694664            &          tab2d_2=ptauj_ice , clinfo2='          vtau_ice : ', mask2=vmask ) 
    695665      END IF 
    696 >>>>>>> .merge-right.r13211 
    697666#endif 
    698667      !                            !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
Note: See TracChangeset for help on using the changeset viewer.