Changeset 13208 for NEMO/trunk/src/ABL


Ignore:
Timestamp:
2020-07-02T10:54:35+02:00 (7 months ago)
Author:
smasson
Message:

trunk: Mid-year merge, merge back dev_r12472_ASINTER-05_Masson_CurrentFeedback

Location:
NEMO/trunk/src/ABL
Files:
2 edited

Legend:

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

    r12814 r13208  
    1818   USE dom_oce, ONLY  : tmask   
    1919   USE sbc_oce, ONLY  : ght_abl, ghw_abl, e3t_abl, e3w_abl, jpka, jpkam1, rhoa 
    20    USE sbcblk         ! use rn_?fac 
     20   USE sbcblk, ONLY   : rn_efac 
    2121   USE sbcblk_phy     ! use some physical constants for flux computation 
    2222   ! 
     
    530530       
    531531      DO_2D_01_01 
    532          zwnd_i(ji,jj) = u_abl(ji  ,jj,2,nt_a) - 0.5_wp * rn_vfac * ( pssu(ji  ,jj) + pssu(ji-1,jj) )   
    533          zwnd_j(ji,jj) = v_abl(ji,jj  ,2,nt_a) - 0.5_wp * rn_vfac * ( pssv(ji,jj  ) + pssv(ji,jj-1) )  
     532         zwnd_i(ji,jj) = u_abl(ji  ,jj,2,nt_a) - 0.5_wp * ( pssu(ji  ,jj) + pssu(ji-1,jj) )   
     533         zwnd_j(ji,jj) = v_abl(ji,jj  ,2,nt_a) - 0.5_wp * ( pssv(ji,jj  ) + pssv(ji,jj-1) )  
    534534      END_2D 
    535535      !  
     
    570570 
    571571#if defined key_si3 
    572          ! ------------------------------------------------------------ ! 
    573          !    Wind stress relative to the moving ice ( U10m - U_ice )   ! 
    574          ! ------------------------------------------------------------ ! 
    575          DO_2D_00_00 
    576              
    577             zztmp1 = 0.5_wp * ( u_abl(ji+1,jj,2,nt_a) + u_abl(ji,jj,2,nt_a) ) 
    578             zztmp2 = 0.5_wp * ( v_abl(ji,jj+1,2,nt_a) + v_abl(ji,jj,2,nt_a) ) 
    579     
    580             ptaui_ice(ji,jj) = 0.5_wp * (  rhoa(ji+1,jj) * pCd_du_ice(ji+1,jj)             & 
    581                &                      +    rhoa(ji  ,jj) * pCd_du_ice(ji  ,jj)  )          & 
    582                &         * ( zztmp1 - rn_vfac * pssu_ice(ji,jj) ) 
    583             ptauj_ice(ji,jj) = 0.5_wp * (  rhoa(ji,jj+1) * pCd_du_ice(ji,jj+1)             & 
    584                &                      +    rhoa(ji,jj  ) * pCd_du_ice(ji,jj  )  )          & 
    585                &         * ( zztmp2 - rn_vfac * pssv_ice(ji,jj) ) 
    586          END_2D 
    587          CALL lbc_lnk_multi( 'ablmod', ptaui_ice, 'U', -1., ptauj_ice, 'V', -1. ) 
    588          ! 
    589          IF(sn_cfctl%l_prtctl)   CALL prt_ctl( tab2d_1=ptaui_ice  , clinfo1=' abl_stp: putaui : '   & 
    590             &                                , tab2d_2=ptauj_ice  , clinfo2='          pvtaui : ' ) 
     572      ! ------------------------------------------------------------ ! 
     573      !    Wind stress relative to the moving ice ( U10m - U_ice )   ! 
     574      ! ------------------------------------------------------------ ! 
     575      DO_2D_00_00            
     576         ptaui_ice(ji,jj) = 0.5_wp * ( rhoa(ji+1,jj) * pCd_du_ice(ji+1,jj) + rhoa(ji,jj) * pCd_du_ice(ji,jj)      )   & 
     577            &                      * ( 0.5_wp * ( u_abl(ji+1,jj,2,nt_a) + u_abl(ji,jj,2,nt_a) ) - pssu_ice(ji,jj) ) 
     578         ptauj_ice(ji,jj) = 0.5_wp * ( rhoa(ji,jj+1) * pCd_du_ice(ji,jj+1) + rhoa(ji,jj) * pCd_du_ice(ji,jj)      )   & 
     579            &                      * ( 0.5_wp * ( v_abl(ji,jj+1,2,nt_a) + v_abl(ji,jj,2,nt_a) ) - pssv_ice(ji,jj) ) 
     580      END_2D 
     581      CALL lbc_lnk_multi( 'ablmod', ptaui_ice, 'U', -1., ptauj_ice, 'V', -1. ) 
     582      ! 
     583      IF(sn_cfctl%l_prtctl)   CALL prt_ctl( tab2d_1=ptaui_ice  , clinfo1=' abl_stp: putaui : '   & 
     584         &                                , tab2d_2=ptauj_ice  , clinfo2='          pvtaui : ' ) 
    591585#endif 
    592586      !                            !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
     
    599593   END SUBROUTINE abl_stp 
    600594!=================================================================================================== 
    601  
    602  
    603  
    604  
    605  
    606  
    607  
    608  
    609  
    610  
    611  
    612  
    613  
    614595 
    615596 
  • NEMO/trunk/src/ABL/sbcabl.F90

    r12814 r13208  
    335335            &                tq_abl(:,:,2,nt_n,jp_ta), tq_abl(:,:,2,nt_n,jp_qa),   &   !   <<= in 
    336336            &                sf(jp_slp )%fnow(:,:,1) , sst_m, ssu_m, ssv_m     ,   &   !   <<= in 
     337            &                sf(jp_uoatm)%fnow(:,:,1), sf(jp_voatm)%fnow(:,:,1),   &   !   <<= in 
    337338            &                sf(jp_qsr )%fnow(:,:,1) , sf(jp_qlw )%fnow(:,:,1) ,   &   !   <<= in 
    338339            &                tsk_m, zssq, zcd_du, zsen, zevp                       )   !   =>> out 
Note: See TracChangeset for help on using the changeset viewer.