Changeset 13947
- Timestamp:
- 2020-12-01T18:13:02+01:00 (4 years ago)
- Location:
- NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/OCE
- Files:
-
- 9 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/OCE/DYN/dynspg.F90
r13710 r13947 6 6 !! History : 1.0 ! 2005-12 (C. Talandier, G. Madec, V. Garnier) Original code 7 7 !! 3.2 ! 2009-07 (R. Benshila) Suppression of rigid-lid option 8 !! 4.2 ! 2020-12 (G. Madec, E. Clementi) add Bernoulli Head for 9 !! wave coupling 8 10 !!---------------------------------------------------------------------- 9 11 -
NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/OCE/DYN/dynvor.F90
r13710 r13947 21 21 !! - ! 2018-03 (G. Madec) add two new schemes (ln_dynvor_enT and ln_dynvor_eet) 22 22 !! - ! 2018-04 (G. Madec) add pre-computed gradient for metric term calculation 23 !! 4.2 ! 2020- 05(G. Madec, E. Clementi) add vortex force trends (ln_vortex_force=T)23 !! 4.2 ! 2020-12 (G. Madec, E. Clementi) add vortex force trends (ln_vortex_force=T) 24 24 !!---------------------------------------------------------------------- 25 25 -
NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/OCE/SBC/sbc_oce.F90
r13710 r13947 12 12 !! 4.0 ! 2016-06 (L. Brodeau) new unified bulk routine (based on AeroBulk) 13 13 !! 4.0 ! 2019-03 (F. Lemarié, G. Samson) add compatibility with ABL mode 14 !! 4.2 ! 2020- 05 9Madec, E. Clementi) modified wave parameters in namelist14 !! 4.2 ! 2020-12 (G. Madec, E. Clementi) modified wave parameters in namelist 15 15 !!---------------------------------------------------------------------- 16 16 -
NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/OCE/SBC/sbcblk_algo_ecmwf.F90
r13710 r13947 17 17 !!---------------------------------------------------------------------- 18 18 !! History : 4.0 ! 2016-02 (L.Brodeau) Original code 19 !! 4.2 ! 2020- 05 (G. MAdec, E. Clementi) Charnock coeff from wave model19 !! 4.2 ! 2020-12 (G. Madec, E. Clementi) Charnock coeff from wave model 20 20 !!---------------------------------------------------------------------- 21 21 … … 32 32 USE in_out_manager ! I/O manager 33 33 USE prtctl ! Print control 34 USE sbcwave, ONLY : charn , ln_charn! wave module34 USE sbcwave, ONLY : charn ! wave module 35 35 #if defined key_si3 || defined key_cice 36 36 USE sbc_ice ! Surface boundary condition: ice fields … … 235 235 236 236 IF (ln_charn) THEN ! Charnock value if wave coupling 237 z0 = charn*u_star*u_star/grav + 0.11 *znu_a/u_star237 z0 = charn*u_star*u_star/grav + 0.11_wp*znu_a/u_star 238 238 ELSE 239 z0 = charn0*u_star*u_star/grav + 0.11 *znu_a/u_star239 z0 = charn0*u_star*u_star/grav + 0.11_wp*znu_a/u_star 240 240 ENDIF 241 241 -
NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/OCE/SBC/sbccpl.F90
r13759 r13947 8 8 !! 3.1 ! 2009_02 (G. Madec, S. Masson, E. Maisonave, A. Caubel) generic coupled interface 9 9 !! 3.4 ! 2011_11 (C. Harris) more flexibility + multi-category fields 10 !! 4.2 ! 2020- 05(G. Madec, E. Clementi) wave coupling updates10 !! 4.2 ! 2020-12 (G. Madec, E. Clementi) wave coupling updates 11 11 !!---------------------------------------------------------------------- 12 12 … … 670 670 cpl_wdrag = .TRUE. 671 671 ENDIF 672 srcv(jpr_charn)%clname = 'O_Charn' ! neutral surface dragcoefficient672 srcv(jpr_charn)%clname = 'O_Charn' ! Chranock coefficient 673 673 IF( TRIM(sn_rcv_charn%cldes ) == 'coupled' ) THEN 674 674 srcv(jpr_charn)%laction = .TRUE. 675 675 cpl_charn = .TRUE. 676 676 ENDIF 677 srcv(jpr_bhd)%clname = 'O_Bhd' ! neutral surface drag coefficient677 srcv(jpr_bhd)%clname = 'O_Bhd' ! Bernoulli head. waves' induced surface pressure 678 678 IF( TRIM(sn_rcv_bhd%cldes ) == 'coupled' ) THEN 679 679 srcv(jpr_bhd)%laction = .TRUE. 680 680 cpl_bhd = .TRUE. 681 681 ENDIF 682 srcv(jpr_tusd)%clname = 'O_Tusd' ! neutral surface drag coefficient682 srcv(jpr_tusd)%clname = 'O_Tusd' ! zonal stokes transport 683 683 IF( TRIM(sn_rcv_tusd%cldes ) == 'coupled' ) THEN 684 684 srcv(jpr_tusd)%laction = .TRUE. 685 685 cpl_tusd = .TRUE. 686 686 ENDIF 687 srcv(jpr_tvsd)%clname = 'O_Tvsd' ! neutral surface drag coefficient687 srcv(jpr_tvsd)%clname = 'O_Tvsd' ! meridional stokes tranmport 688 688 IF( TRIM(sn_rcv_tvsd%cldes ) == 'coupled' ) THEN 689 689 srcv(jpr_tvsd)%laction = .TRUE. … … 1389 1389 ! 1390 1390 ! ! ========================= ! 1391 ! ! Wave drag coefficient!1391 ! ! Chranock coefficient ! 1392 1392 ! ! ========================= ! 1393 1393 IF( srcv(jpr_charn)%laction .AND. ln_charn ) charn(:,:) = frcv(jpr_charn)%z3(:,:,1) 1394 1394 ! 1395 ! 1395 ! ! ========================= ! 1396 ! ! net wave-supported stress ! 1397 ! ! ========================= ! 1396 1398 IF( srcv(jpr_tawx)%laction .AND. ln_taw ) tawx(:,:) = frcv(jpr_tawx)%z3(:,:,1) 1397 1399 IF( srcv(jpr_tawy)%laction .AND. ln_taw ) tawy(:,:) = frcv(jpr_tawy)%z3(:,:,1) 1400 ! 1401 ! ! ========================= ! 1402 ! !wave to ocean momentum flux! 1403 ! ! ========================= ! 1398 1404 IF( srcv(jpr_twox)%laction .AND. ln_taw ) twox(:,:) = frcv(jpr_twox)%z3(:,:,1) 1399 1405 IF( srcv(jpr_twoy)%laction .AND. ln_taw ) twoy(:,:) = frcv(jpr_twoy)%z3(:,:,1) -
NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/OCE/SBC/sbcmod.F90
r13852 r13947 16 16 !! 4.0 ! 2016-06 (L. Brodeau) new general bulk formulation 17 17 !! 4.0 ! 2019-03 (F. Lemarié & G. Samson) add ABL compatibility (ln_abl=TRUE) 18 !! 4.2 ! 2020- 05 (G. Madec, E. Clementi) modified wave forcing and coipling18 !! 4.2 ! 2020-12 (G. Madec, E. Clementi) modified wave forcing and coupling 19 19 !!---------------------------------------------------------------------- 20 20 -
NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/OCE/SBC/sbcwave.F90
r13759 r13947 9 9 !! - ! 2016-12 (G. Madec, E. Clementi) update Stoke drift computation 10 10 !! + add sbc_wave_ini routine 11 !! 4.2 ! 2020- 06 (G. Madec, E. Clement) updates, new Stoke drift computation12 !! according to Co velard et al.,201911 !! 4.2 ! 2020-12 (G. Madec, E. Clementi) updates, new Stoke drift computation 12 !! according to Couvelard et al.,2019 13 13 !!---------------------------------------------------------------------- 14 14 -
NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/OCE/ZDF/zdfsh2.F90
r13710 r13947 6 6 !! History : - ! 2014-10 (A. Barthelemy, G. Madec) original code 7 7 !! NEMO 4.0 ! 2017-04 (G. Madec) remove u-,v-pts avm 8 !! NEMO 4.2 ! 2020-12 (G. Madec, E. Clementi) add Stokes Drift Shear 9 ! ! for wave coupling 8 10 !!---------------------------------------------------------------------- 9 11 -
NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/OCE/ZDF/zdftke.F90
r13873 r13947 29 29 !! 4.0 ! 2017-04 (G. Madec) remove CPP ddm key & avm at t-point only 30 30 !! - ! 2017-05 (G. Madec) add top/bottom friction as boundary condition 31 !! 4.2 ! 2020-12 (G. Madec, E. Clementi) add wave coupling 32 ! ! following Couvelard et al., 2019 31 33 !!---------------------------------------------------------------------- 32 34 … … 397 399 ! !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 398 400 ! ! Surface boundary condition on tke if 399 ! ! coupl eing with waves401 ! ! coupling with waves 400 402 ! !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 401 403 ! … … 425 427 ! 426 428 ! !* Matrix inversion from level 2 (tke prescribed at level 1) 427 ! DO_3D( 0, 0, 0, 0, 3, jpkm1 ) ! First recurrence : Dk = Dk - Lk * Uk-1 / Dk-1428 429 DO_3D( 0, 0, 0, 0, 2, jpkm1 ) ! First recurrence : Dk = Dk - Lk * Uk-1 / Dk-1 429 430 zdiag(ji,jj,jk) = zdiag(ji,jj,jk) - zd_lw(ji,jj,jk) * zd_up(ji,jj,jk-1) / zdiag(ji,jj,jk-1) … … 433 434 ! zd_lw(ji,jj,2) = en(ji,jj,2) - zd_lw(ji,jj,2) * en(ji,jj,1) ! Surface boudary conditions on tke 434 435 ! END_2D 435 ! DO_3D( 0, 0, 0, 0, 3, jpkm1 )436 436 DO_3D( 0, 0, 0, 0, 2, jpkm1 ) 437 437 zd_lw(ji,jj,jk) = en(ji,jj,jk) - zd_lw(ji,jj,jk) / zdiag(ji,jj,jk-1) *zd_lw(ji,jj,jk-1)
Note: See TracChangeset
for help on using the changeset viewer.