Changeset 11637 for NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk.F90
- Timestamp:
- 2019-10-02T12:48:31+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk.F90
r11623 r11637 153 153 !! *** ROUTINE sbc_blk_cswl_alloc *** 154 154 !!------------------------------------------------------------------- 155 ! PRINT *,'*** LB: allocating tsk!'155 !WRITE(numout,*) '*** LB: allocating tsk!' 156 156 ALLOCATE( tsk(jpi,jpj), STAT=sbc_blk_cswl_alloc ) 157 ! PRINT *,'*** LB: done!'157 !WRITE(numout,*) '*** LB: done!' 158 158 CALL mpp_sum ( 'sbcblk', sbc_blk_cswl_alloc ) 159 159 IF( sbc_blk_cswl_alloc /= 0 ) CALL ctl_stop( 'STOP', 'sbc_blk_cswl_alloc: failed to allocate arrays' ) … … 510 510 zqair(:,:) = sf(jp_humi)%fnow(:,:,1) ! what we read in file is already a spec. humidity! 511 511 CASE( np_humi_dpt ) 512 IF (lwp) PRINT *,' *** LB(sbcblk.F90) => computing q_air out of d_air and slp !'512 IF (lwp) WRITE(numout,*) ' *** LB(sbcblk.F90) => computing q_air out of d_air and slp !' 513 513 zqair(:,:) = q_sat( sf(jp_humi)%fnow(:,:,1), sf(jp_slp)%fnow(:,:,1) ) 514 514 CASE( np_humi_rlh ) 515 IF (lwp) PRINT *,' *** LB(sbcblk.F90) => computing q_air out of RH, t_air and slp !'515 IF (lwp) WRITE(numout,*) ' *** LB(sbcblk.F90) => computing q_air out of RH, t_air and slp !' 516 516 zqair(:,:) = q_air_rh( 0.01_wp*sf(jp_humi)%fnow(:,:,1), sf(jp_tair)%fnow(:,:,1), sf(jp_slp)%fnow(:,:,1) ) 517 517 END SELECT … … 529 529 530 530 CASE( np_COARE_3p0 ) 531 IF (lwp) PRINT *, ' *** LB(sbcblk.F90) => calling "turb_coare3p0" WITH CSWL options!!!, nsec_day=', nsec_day531 IF (lwp) WRITE(numout,*) ' *** LB(sbcblk.F90) => calling "turb_coare3p0" WITH CSWL options!!!, nsec_day, gdept_1d(1)=', nsec_day, gdept_1d(1) 532 532 CALL turb_coare3p0 ( kt, rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl,& ! COARE v3.0 533 533 & Cd_atm, Ch_atm, Ce_atm, t_zu, q_zu, zU_zu, cdn_oce, chn_oce, cen_oce, & … … 536 536 537 537 CASE( np_COARE_3p6 ) 538 IF (lwp) PRINT *,' *** LB(sbcblk.F90) => calling "turb_coare3p6" WITH CSWL options!!!, nsec_day, gdept_1d(1)=', nsec_day, gdept_1d(1)538 IF (lwp) WRITE(numout,*) ' *** LB(sbcblk.F90) => calling "turb_coare3p6" WITH CSWL options!!!, nsec_day, gdept_1d(1)=', nsec_day, gdept_1d(1) 539 539 CALL turb_coare3p6 ( kt, rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl,& ! COARE v3.6 540 540 & Cd_atm, Ch_atm, Ce_atm, t_zu, q_zu, zU_zu, cdn_oce, chn_oce, cen_oce, & … … 543 543 544 544 CASE( np_ECMWF ) 545 IF (lwp) PRINT *, ' *** LB(sbcblk.F90) => calling "turb_ecmwf" WITH CSWL options!!!, nsec_day=', nsec_day545 IF (lwp) WRITE(numout,*) ' *** LB(sbcblk.F90) => calling "turb_ecmwf" WITH CSWL options!!!, nsec_day, gdept_1d(1)=', nsec_day, gdept_1d(1) 546 546 CALL turb_ecmwf ( rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl, & ! ECMWF 547 547 & Cd_atm, Ch_atm, Ce_atm, t_zu, q_zu, zU_zu, cdn_oce, chn_oce, cen_oce, & … … 576 576 577 577 CASE( np_COARE_3p0 ) 578 IF (lwp) PRINT *,' *** LB(sbcblk.F90) => calling "turb_coare3p0" WITHOUT CSWL optional arrays!!!'578 IF (lwp) WRITE(numout,*) ' *** LB(sbcblk.F90) => calling "turb_coare3p0" WITHOUT CSWL optional arrays!!!' 579 579 CALL turb_coare3p0 ( kt, rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl,& ! COARE v3.0 580 580 & Cd_atm, Ch_atm, Ce_atm, t_zu, q_zu, zU_zu, cdn_oce, chn_oce, cen_oce ) 581 581 582 582 CASE( np_COARE_3p6 ) 583 IF (lwp) PRINT *,' *** LB(sbcblk.F90) => calling "turb_coare3p6" WITHOUT CSWL optional arrays!!!'583 IF (lwp) WRITE(numout,*) ' *** LB(sbcblk.F90) => calling "turb_coare3p6" WITHOUT CSWL optional arrays!!!' 584 584 CALL turb_coare3p6 ( kt, rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl,& ! COARE v3.6 585 585 & Cd_atm, Ch_atm, Ce_atm, t_zu, q_zu, zU_zu, cdn_oce, chn_oce, cen_oce ) 586 586 587 587 CASE( np_ECMWF ) 588 IF (lwp) PRINT *,' *** LB(sbcblk.F90) => calling "turb_ecmwf" WITHOUT CSWL optional arrays!!!'588 IF (lwp) WRITE(numout,*) ' *** LB(sbcblk.F90) => calling "turb_ecmwf" WITHOUT CSWL optional arrays!!!' 589 589 CALL turb_ecmwf ( rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl, & ! ECMWF 590 590 & Cd_atm, Ch_atm, Ce_atm, t_zu, q_zu, zU_zu, cdn_oce, chn_oce, cen_oce ) … … 854 854 zqair(:,:) = sf(jp_humi)%fnow(:,:,1) ! what we read in file is already a spec. humidity! 855 855 CASE( np_humi_dpt ) 856 IF (lwp) PRINT *,' *** LB(sbcblk.F90) => ICE !!! computing q_air out of d_air and slp !'856 IF (lwp) WRITE(numout,*) ' *** LB(sbcblk.F90) => ICE !!! computing q_air out of d_air and slp !' 857 857 zqair(:,:) = q_sat( sf(jp_humi)%fnow(:,:,1), sf(jp_slp)%fnow(:,:,1) ) 858 858 CASE( np_humi_rlh ) 859 IF (lwp) PRINT *,' *** LB(sbcblk.F90) => ICE !!! computing q_air out of RH, t_air and slp !'859 IF (lwp) WRITE(numout,*) ' *** LB(sbcblk.F90) => ICE !!! computing q_air out of RH, t_air and slp !' 860 860 zqair(:,:) = q_air_rh( 0.01_wp*sf(jp_humi)%fnow(:,:,1), sf(jp_tair)%fnow(:,:,1), sf(jp_slp)%fnow(:,:,1) ) 861 861 END SELECT
Note: See TracChangeset
for help on using the changeset viewer.