Changeset 11637 for NEMO/branches
- Timestamp:
- 2019-10-02T12:48:31+02:00 (5 years ago)
- Location:
- NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk
- Files:
-
- 11 added
- 2 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 -
NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk_skin_ecmwf.F90
r11623 r11637 151 151 !!--------------------------------------------------------------------- 152 152 153 IF (lwp) WRITE(numout,*) 'LOLO: WL_ECMWF@sbcblk_skin_ecmwf.F90: depth for SST bulk takent at', gdept_1d(1), ' m'154 155 153 DO jj = 1, jpj 156 154 DO ji = 1, jpi … … 164 162 flg = 0.5_wp + SIGN( 0.5_wp , gdept_1d(1)-zdz ) ! => 1 when gdept_1d(1)>zdz (pdT(ji,jj) = dT_wl) | 0 when z_s$ 165 163 dT_wl = pdT(ji,jj) / ( flg + (1._wp-flg)*gdept_1d(1)/zdz ) 166 !PRINT *, 'LOLO/mod_wl_ecmwf.f90: dT_wl2=', dT_wl167 !PRINT *, ''168 164 169 165 zalpha_w = alpha_sw( pSST(ji,jj) ) ! thermal expansion coefficient of sea-water (SST accurate enough!)
Note: See TracChangeset
for help on using the changeset viewer.