Changeset 11997 for NEMO/branches
- Timestamp:
- 2019-11-28T12:49:14+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk.F90
r11996 r11997 207 207 208 208 209 !LB: 209 210 210 ! !** initialization of the cool-skin / warm-layer parametrization 211 IF( ln_NCAR .AND. (ln_skin_cs .OR. ln_skin_wl) ) & 212 & CALL ctl_stop( 'sbc_blk_init: Cool-skin/warm-layer param. not compatible with NCAR algorithm!' ) 213 ! 211 IF( ln_skin_cs .OR. ln_skin_wl ) THEN 212 !! Some namelist sanity tests: 213 IF( ln_NCAR ) & 214 & CALL ctl_stop( 'sbc_blk_init: Cool-skin/warm-layer param. not compatible with NCAR algorithm' ) 215 IF( nn_fsbc /= 1 ) & 216 & CALL ctl_stop( 'sbc_blk_init: Please set "nn_fsbc" to 1 when using cool-skin/warm-layer param.') 217 END IF 218 214 219 ioptio = 0 215 220 IF( ln_humi_sph ) THEN … … 532 537 CALL ctl_stop( 'STOP', 'sbc_oce: non-existing bulk formula selected' ) 533 538 END SELECT 534 539 535 540 IF( ln_skin_cs .OR. ln_skin_wl ) THEN 536 541 !! In the presence of sea-ice we forget about the cool-skin/warm-layer update of zst and zsq: … … 546 551 tsk(:,:) = zst(:,:) !#LB: Update of tsk, the "official" array for skin temperature 547 552 END IF 548 549 !ELSE !IF( ln_skin_cs .OR. ln_skin_wl )550 !551 ! SELECT CASE( nblk ) !== transfer coefficients ==! Cd, Ch, Ce at T-point552 !553 !CASE( np_COARE_3p0 )554 ! CALL turb_coare3p0 ( kt, rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl, &555 ! & Cd_atm, Ch_atm, Ce_atm, t_zu, q_zu, zU_zu, cdn_oce, chn_oce, cen_oce )556 !557 ! CASE( np_COARE_3p6 )558 !CALL turb_coare3p6 ( kt, rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl, &559 ! & Cd_atm, Ch_atm, Ce_atm, t_zu, q_zu, zU_zu, cdn_oce, chn_oce, cen_oce )560 !561 ! CASE( np_ECMWF )562 ! CALL turb_ecmwf ( kt, rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl, &563 ! & Cd_atm, Ch_atm, Ce_atm, t_zu, q_zu, zU_zu, cdn_oce, chn_oce, cen_oce )564 ! !565 ! CASE DEFAULT566 ! CALL ctl_stop( 'STOP', 'sbc_oce: non-existing bulk formula selected' )567 ! END SELECT!568 !569 !ENDIF ! IF( ln_skin_cs .OR. ln_skin_wl )570 553 571 554 !! CALL iom_put( "Cd_oce", Cd_atm) ! output value of pure ocean-atm. transfer coef.
Note: See TracChangeset
for help on using the changeset viewer.