Changeset 11962
- Timestamp:
- 2019-11-25T23:31:07+01:00 (5 years ago)
- Location:
- NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk
- Files:
-
- 7 edited
-
src/OCE/SBC/sbcblk.F90 (modified) (14 diffs)
-
src/OCE/SBC/sbcblk_algo_coare3p0.F90 (modified) (5 diffs)
-
src/OCE/SBC/sbcblk_algo_coare3p6.F90 (modified) (5 diffs)
-
src/OCE/SBC/sbcblk_algo_ecmwf.F90 (modified) (6 diffs)
-
src/OCE/SBC/sbcblk_phy.F90 (modified) (3 diffs)
-
src/OCE/SBC/sbcblk_skin_coare.F90 (modified) (1 diff)
-
tests/STATION_ASF/EXPREF/plot_station_asf.py (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk.F90
r11841 r11962 142 142 !!------------------------------------------------------------------- 143 143 ALLOCATE( Cd_atm (jpi,jpj), Ch_atm (jpi,jpj), Ce_atm (jpi,jpj), t_zu(jpi,jpj), q_zu(jpi,jpj), & 144 & cdn_oce(jpi,jpj), chn_oce(jpi,jpj), cen_oce(jpi,jpj), STAT=sbc_blk_alloc )144 & cdn_oce(jpi,jpj), chn_oce(jpi,jpj), cen_oce(jpi,jpj), tsk(jpi,jpj), STAT=sbc_blk_alloc ) 145 145 ! 146 146 CALL mpp_sum ( 'sbcblk', sbc_blk_alloc ) 147 147 IF( sbc_blk_alloc /= 0 ) CALL ctl_stop( 'STOP', 'sbc_blk_alloc: failed to allocate arrays' ) 148 148 END FUNCTION sbc_blk_alloc 149 150 !LB:151 INTEGER FUNCTION sbc_blk_cswl_alloc()152 !!-------------------------------------------------------------------153 !! *** ROUTINE sbc_blk_cswl_alloc ***154 !!-------------------------------------------------------------------155 !WRITE(numout,*) '*** LB: allocating tsk!'156 ALLOCATE( tsk(jpi,jpj), STAT=sbc_blk_cswl_alloc )157 !WRITE(numout,*) '*** LB: done!'158 CALL mpp_sum ( 'sbcblk', sbc_blk_cswl_alloc )159 IF( sbc_blk_cswl_alloc /= 0 ) CALL ctl_stop( 'STOP', 'sbc_blk_cswl_alloc: failed to allocate arrays' )160 END FUNCTION sbc_blk_cswl_alloc161 !LB.162 149 163 150 … … 222 209 !LB: 223 210 ! !** initialization of the cool-skin / warm-layer parametrization 224 IF( ln_skin_cs .OR. ln_skin_wl ) THEN 225 IF ( ln_NCAR ) CALL ctl_stop( 'sbc_blk_init: Cool-skin/warm-layer param. not compatible with NCAR algorithm!' ) 226 ! ! allocate array(s) for cool-skin/warm-layer param. 227 IF( sbc_blk_cswl_alloc() /= 0 ) CALL ctl_stop( 'STOP', 'sbc_blk : unable to allocate standard arrays' ) 228 END IF 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!' ) 229 213 ! 230 214 ioptio = 0 … … 322 306 ! 323 307 WRITE(numout,*) 324 WRITE(numout,*) ' use cool-skin parameterization (SSST) ln_skin_cs = ', ln_skin_cs !LB325 WRITE(numout,*) ' use warm-layer parameterization (SSST) ln_skin_wl = ', ln_skin_wl !LB308 WRITE(numout,*) ' use cool-skin parameterization (SSST) ln_skin_cs = ', ln_skin_cs 309 WRITE(numout,*) ' use warm-layer parameterization (SSST) ln_skin_wl = ', ln_skin_wl 326 310 ! 327 !LB:328 311 WRITE(numout,*) 329 312 SELECT CASE( nhumi ) !* Print the choice of air humidity … … 332 315 CASE( np_humi_rlh ) ; WRITE(numout,*) ' ==>>> air humidity is RELATIVE HUMIDITY [%]' 333 316 END SELECT 334 !LB.335 317 ! 336 318 ENDIF … … 390 372 ENDIF 391 373 tatm_ice(:,:) = sf(jp_tair)%fnow(:,:,1) 392 !LB: 374 393 375 SELECT CASE( nhumi ) 394 376 CASE( np_humi_sph ) … … 399 381 qatm_ice(:,:) = q_air_rh( 0.01_wp*sf(jp_humi)%fnow(:,:,1), sf(jp_tair)%fnow(:,:,1), sf(jp_slp)%fnow(:,:,1)) !LB: 0.01 => RH is % percent in file 400 382 END SELECT 401 !LB. 383 402 384 tprecip(:,:) = sf(jp_prec)%fnow(:,:,1) * rn_pfac 403 385 sprecip(:,:) = sf(jp_snow)%fnow(:,:,1) * rn_pfac … … 445 427 REAL(wp), DIMENSION(jpi,jpj) :: zU_zu ! bulk wind speed at height zu [m/s] 446 428 REAL(wp), DIMENSION(jpi,jpj) :: ztpot ! potential temperature of air at z=rn_zqt [K] 447 REAL(wp), DIMENSION(jpi,jpj) :: zqair ! specific humidity of air at z=rn_zqt [kg/kg] !LB429 REAL(wp), DIMENSION(jpi,jpj) :: zqair ! specific humidity of air at z=rn_zqt [kg/kg] 448 430 !!--------------------------------------------------------------------- 449 431 ! … … 510 492 zqair(:,:) = sf(jp_humi)%fnow(:,:,1) ! what we read in file is already a spec. humidity! 511 493 CASE( np_humi_dpt ) 512 IF (lwp) WRITE(numout,*) ' *** blk_oce => computing q_air out of d_air and slp !' !LBrm494 !IF (lwp) WRITE(numout,*) ' *** blk_oce => computing q_air out of d_air and slp !' !LBrm 513 495 zqair(:,:) = q_sat( sf(jp_humi)%fnow(:,:,1), sf(jp_slp)%fnow(:,:,1) ) 514 496 CASE( np_humi_rlh ) 515 IF (lwp) WRITE(numout,*) ' *** blk_oce => computing q_air out of RH, t_air and slp !' !LBrm497 !IF (lwp) WRITE(numout,*) ' *** blk_oce => computing q_air out of RH, t_air and slp !' !LBrm 516 498 zqair(:,:) = q_air_rh( 0.01_wp*sf(jp_humi)%fnow(:,:,1), sf(jp_tair)%fnow(:,:,1), sf(jp_slp)%fnow(:,:,1) ) !LB: 0.01 => RH is % percent in file 517 499 END SELECT … … 529 511 530 512 CASE( np_COARE_3p0 ) 531 IF (lwp) WRITE(numout,*) ' *** blk_oce => calling "turb_coare3p0" WITH CSWL options!!!, gdept_1d(1)=', gdept_1d(1) !LBrm 532 CALL turb_coare3p0 ( kt, rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl,& ! COARE v3.0 533 & Cd_atm, Ch_atm, Ce_atm, t_zu, q_zu, zU_zu, cdn_oce, chn_oce, cen_oce, & 513 CALL turb_coare3p0 ( kt, rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl, & 514 & Cd_atm, Ch_atm, Ce_atm, t_zu, q_zu, zU_zu, cdn_oce, chn_oce, cen_oce, & 534 515 & Qsw=qsr(:,:), rad_lw=sf(jp_qlw)%fnow(:,:,1), slp=sf(jp_slp)%fnow(:,:,1) ) 535 516 536 517 CASE( np_COARE_3p6 ) 537 IF (lwp) WRITE(numout,*) ' *** blk_oce => calling "turb_coare3p6" WITH CSWL options!!!, gdept_1d(1)=', gdept_1d(1) !LBrm 538 CALL turb_coare3p6 ( kt, rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl,& ! COARE v3.6 539 & Cd_atm, Ch_atm, Ce_atm, t_zu, q_zu, zU_zu, cdn_oce, chn_oce, cen_oce, & 518 CALL turb_coare3p6 ( kt, rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl, & 519 & Cd_atm, Ch_atm, Ce_atm, t_zu, q_zu, zU_zu, cdn_oce, chn_oce, cen_oce, & 540 520 & Qsw=qsr(:,:), rad_lw=sf(jp_qlw)%fnow(:,:,1), slp=sf(jp_slp)%fnow(:,:,1) ) 541 521 542 522 CASE( np_ECMWF ) 543 IF (lwp) WRITE(numout,*) ' *** blk_oce => calling "turb_ecmwf" WITH CSWL options!!!, gdept_1d(1)=', gdept_1d(1) !LBrm 544 CALL turb_ecmwf ( kt, rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl, & ! ECMWF 545 & Cd_atm, Ch_atm, Ce_atm, t_zu, q_zu, zU_zu, cdn_oce, chn_oce, cen_oce, & 523 CALL turb_ecmwf ( kt, rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl, & 524 & Cd_atm, Ch_atm, Ce_atm, t_zu, q_zu, zU_zu, cdn_oce, chn_oce, cen_oce, & 546 525 & Qsw=qsr(:,:), rad_lw=sf(jp_qlw)%fnow(:,:,1), slp=sf(jp_slp)%fnow(:,:,1) ) 547 526 … … 566 545 ELSE !IF ( ln_skin_cs .OR. ln_skin_wl ) 567 546 568 569 547 SELECT CASE( nblk ) !== transfer coefficients ==! Cd, Ch, Ce at T-point 570 548 ! 571 549 CASE( np_NCAR ) 572 CALL turb_ncar ( rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, & ! NCAR-COREv2550 CALL turb_ncar ( rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, & 573 551 & Cd_atm, Ch_atm, Ce_atm, t_zu, q_zu, zU_zu, cdn_oce, chn_oce, cen_oce ) 574 552 575 553 CASE( np_COARE_3p0 ) 576 IF (lwp) WRITE(numout,*) ' *** blk_oce => calling "turb_coare3p0" WITHOUT CSWL optional arrays!!!' !LBrm 577 CALL turb_coare3p0 ( kt, rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl,& ! COARE v3.0 554 CALL turb_coare3p0 ( kt, rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl, & 578 555 & Cd_atm, Ch_atm, Ce_atm, t_zu, q_zu, zU_zu, cdn_oce, chn_oce, cen_oce ) 579 556 580 557 CASE( np_COARE_3p6 ) 581 IF (lwp) WRITE(numout,*) ' *** blk_oce => calling "turb_coare3p6" WITHOUT CSWL optional arrays!!!' !LBrm 582 CALL turb_coare3p6 ( kt, rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl,& ! COARE v3.6 558 CALL turb_coare3p6 ( kt, rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl, & 583 559 & Cd_atm, Ch_atm, Ce_atm, t_zu, q_zu, zU_zu, cdn_oce, chn_oce, cen_oce ) 584 560 585 561 CASE( np_ECMWF ) 586 IF (lwp) WRITE(numout,*) ' *** blk_oce => calling "turb_ecmwf" WITHOUT CSWL optional arrays!!!' !LBrm 587 CALL turb_ecmwf ( kt, rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl, & ! ECMWF 562 CALL turb_ecmwf ( kt, rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl, & 588 563 & Cd_atm, Ch_atm, Ce_atm, t_zu, q_zu, zU_zu, cdn_oce, chn_oce, cen_oce ) 589 564 … … 596 571 !! CALL iom_put( "Cd_oce", Cd_atm) ! output value of pure ocean-atm. transfer coef. 597 572 !! CALL iom_put( "Ch_oce", Ch_atm) ! output value of pure ocean-atm. transfer coef. 598 573 599 574 IF( ABS(rn_zu - rn_zqt) < 0.1_wp ) THEN 600 575 !! If zu == zt, then ensuring once for all that: … … 602 577 q_zu(:,:) = zqair(:,:) 603 578 END IF 604 579 605 580 606 581 ! Turbulent fluxes over ocean => BULK_FORMULA @ sbcblk_phy.F90 607 582 ! ------------------------------------------------------------- 608 609 CALL BULK_FORMULA( rn_zu, zst(:,:), zsq(:,:), t_zu(:,:), q_zu(:,:), Cd_atm(:,:), Ch_atm(:,:), Ce_atm(:,:), & 610 & wndm(:,:), zU_zu(:,:), sf(jp_slp)%fnow(:,:,1), & 611 & taum(:,:), zqsb(:,:), zqla(:,:), & 612 & pEvap=zevap(:,:), prhoa=rhoa(:,:) ) 613 583 584 CALL BULK_FORMULA( rn_zu, zst(:,:), zsq(:,:), t_zu(:,:), q_zu(:,:), & 585 & Cd_atm(:,:), Ch_atm(:,:), Ce_atm(:,:), & 586 & wndm(:,:), zU_zu(:,:), sf(jp_slp)%fnow(:,:,1), & 587 & taum(:,:), zqsb(:,:), zqla(:,:), & 588 & pEvap=zevap(:,:), prhoa=rhoa(:,:) ) 589 614 590 zqla(:,:) = zqla(:,:) * tmask(:,:,1) 615 591 zqsb(:,:) = zqsb(:,:) * tmask(:,:,1) 616 592 taum(:,:) = taum(:,:) * tmask(:,:,1) 617 593 zevap(:,:) = zevap(:,:) * tmask(:,:,1) 618 594 619 595 ! Tau i and j component on T-grid points, using array "Cd_atm" as a temporary array... 620 596 Cd_atm = 0._wp 621 597 WHERE ( wndm > 0._wp ) Cd_atm = taum / wndm 622 598 zwnd_i = Cd_atm * zwnd_i 623 zwnd_j = Cd_atm * zwnd_j 624 !DO jj = 1, jpj ! tau i and j component on T-grid points 625 ! DO ji = 1, jpi 626 ! zztmp = taum(ji,jj) / MAX( wndm(ji,jj) , 0.01_wp ) 627 ! zwnd_i(ji,jj) = zztmp * zwnd_i(ji,jj) 628 ! zwnd_j(ji,jj) = zztmp * zwnd_j(ji,jj) 629 ! END DO 630 !END DO 599 zwnd_j = Cd_atm * zwnd_j 600 631 601 ! ! add the HF tau contribution to the wind stress module 632 602 IF( lhftau ) taum(:,:) = taum(:,:) + sf(jp_tdif)%fnow(:,:,1) … … 652 622 ! ----------------------------------------------------------------------------- ! 653 623 654 !! LB: now moved after Turbulent fluxes because must use the skin temperature rather that the SST ! (zst is skin temperature if ln_skin_cs==.TRUE. .OR. ln_skin_wl==.TRUE.) 624 !! LB: now moved after Turbulent fluxes because must use the skin temperature rather that the SST 625 !! (zst is skin temperature if ln_skin_cs==.TRUE. .OR. ln_skin_wl==.TRUE.) 655 626 zqlw(:,:) = emiss_w * ( sf(jp_qlw)%fnow(:,:,1) - stefan*zst(:,:)*zst(:,:)*zst(:,:)*zst(:,:) ) * tmask(:,:,1) ! Net radiative longwave flux 656 627 … … 848 819 zqair(:,:) = sf(jp_humi)%fnow(:,:,1) ! what we read in file is already a spec. humidity! 849 820 CASE( np_humi_dpt ) 850 IF (lwp) WRITE(numout,*) ' *** blk_ice_flx => computing q_air out of d_air and slp !' !LBrm821 !IF (lwp) WRITE(numout,*) ' *** blk_ice_flx => computing q_air out of d_air and slp !' !LBrm 851 822 zqair(:,:) = q_sat( sf(jp_humi)%fnow(:,:,1), sf(jp_slp)%fnow(:,:,1) ) 852 823 CASE( np_humi_rlh ) 853 IF (lwp) WRITE(numout,*) ' *** blk_ice_flx => computing q_air out of RH, t_air and slp !' !LBrm824 !IF (lwp) WRITE(numout,*) ' *** blk_ice_flx => computing q_air out of RH, t_air and slp !' !LBrm 854 825 zqair(:,:) = q_air_rh( 0.01_wp*sf(jp_humi)%fnow(:,:,1), sf(jp_tair)%fnow(:,:,1), sf(jp_slp)%fnow(:,:,1) ) !LB: 0.01 => RH is % percent in file 855 826 END SELECT -
NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk_algo_coare3p0.F90
r11851 r11962 43 43 PRIVATE 44 44 45 PUBLIC :: COARE3P0_INIT, TURB_COARE3P045 PUBLIC :: SBCBLK_ALGO_COARE3P0_INIT, TURB_COARE3P0 46 46 47 47 !! COARE own values for given constants: … … 55 55 56 56 57 SUBROUTINE coare3p0_init(l_use_cs, l_use_wl)57 SUBROUTINE sbcblk_algo_coare3p0_init(l_use_cs, l_use_wl) 58 58 !!--------------------------------------------------------------------- 59 !! *** FUNCTION coare3p0_init ***59 !! *** FUNCTION sbcblk_algo_coare3p0_init *** 60 60 !! 61 61 !! INPUT : … … 71 71 ierr = 0 72 72 ALLOCATE ( Tau_ac(jpi,jpj) , Qnt_ac(jpi,jpj), dT_wl(jpi,jpj), Hz_wl(jpi,jpj), STAT=ierr ) 73 IF( ierr > 0 ) CALL ctl_stop( ' COARE3P0_INIT => allocation of Tau_ac, Qnt_ac, dT_wl & Hz_wl failed!' )73 IF( ierr > 0 ) CALL ctl_stop( ' SBCBLK_ALGO_COARE3P0_INIT => allocation of Tau_ac, Qnt_ac, dT_wl & Hz_wl failed!' ) 74 74 Tau_ac(:,:) = 0._wp 75 75 Qnt_ac(:,:) = 0._wp … … 80 80 ierr = 0 81 81 ALLOCATE ( dT_cs(jpi,jpj), STAT=ierr ) 82 IF( ierr > 0 ) CALL ctl_stop( ' COARE3P0_INIT => allocation of dT_cs failed!' )82 IF( ierr > 0 ) CALL ctl_stop( ' SBCBLK_ALGO_COARE3P0_INIT => allocation of dT_cs failed!' ) 83 83 dT_cs(:,:) = -0.25_wp ! First guess of skin correction 84 84 END IF 85 END SUBROUTINE coare3p0_init85 END SUBROUTINE sbcblk_algo_coare3p0_init 86 86 87 87 … … 185 185 REAL(wp), DIMENSION(jpi,jpj) :: ztmp0, ztmp1, ztmp2 186 186 ! 187 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: & 188 & zeta_t, & ! stability parameter at height zt 189 & zsst, & ! to back up the initial bulk SST 190 & pdTc, & ! SST increment "dT" for cool-skin correction [K] 191 & pdTw, & ! SST increment "dT" for warm layer correction [K] 192 & zHwl ! depth of warm-layer [m] 187 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: zeta_t ! stability parameter at height zt 188 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: zsst ! to back up the initial bulk SST 189 ! 193 190 CHARACTER(len=40), PARAMETER :: crtnm = 'turb_coare3p0@sbcblk_algo_coare3p0' 194 191 !!---------------------------------------------------------------------------------- 195 IF ( kt == nit000 ) CALL COARE3P0_INIT(l_use_cs, l_use_wl)192 IF ( kt == nit000 ) CALL SBCBLK_ALGO_COARE3P0_INIT(l_use_cs, l_use_wl) 196 193 197 194 l_zt_equal_zu = .FALSE. -
NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk_algo_coare3p6.F90
r11851 r11962 43 43 PRIVATE 44 44 45 PUBLIC :: COARE3P6_INIT, TURB_COARE3P645 PUBLIC :: SBCBLK_ALGO_COARE3P6_INIT, TURB_COARE3P6 46 46 47 47 !! COARE own values for given constants: … … 55 55 56 56 57 SUBROUTINE coare3p6_init(l_use_cs, l_use_wl)57 SUBROUTINE sbcblk_algo_coare3p6_init(l_use_cs, l_use_wl) 58 58 !!--------------------------------------------------------------------- 59 !! *** FUNCTION coare3p6_init ***59 !! *** FUNCTION sbcblk_algo_coare3p6_init *** 60 60 !! 61 61 !! INPUT : … … 71 71 ierr = 0 72 72 ALLOCATE ( Tau_ac(jpi,jpj) , Qnt_ac(jpi,jpj), dT_wl(jpi,jpj), Hz_wl(jpi,jpj), STAT=ierr ) 73 IF( ierr > 0 ) CALL ctl_stop( ' COARE3P6_INIT => allocation of Tau_ac, Qnt_ac, dT_wl & Hz_wl failed!' )73 IF( ierr > 0 ) CALL ctl_stop( ' SBCBLK_ALGO_COARE3P6_INIT => allocation of Tau_ac, Qnt_ac, dT_wl & Hz_wl failed!' ) 74 74 Tau_ac(:,:) = 0._wp 75 75 Qnt_ac(:,:) = 0._wp … … 80 80 ierr = 0 81 81 ALLOCATE ( dT_cs(jpi,jpj), STAT=ierr ) 82 IF( ierr > 0 ) CALL ctl_stop( ' COARE3P6_INIT => allocation of dT_cs failed!' )82 IF( ierr > 0 ) CALL ctl_stop( ' SBCBLK_ALGO_COARE3P6_INIT => allocation of dT_cs failed!' ) 83 83 dT_cs(:,:) = -0.25_wp ! First guess of skin correction 84 84 END IF 85 END SUBROUTINE coare3p6_init85 END SUBROUTINE sbcblk_algo_coare3p6_init 86 86 87 87 … … 185 185 REAL(wp), DIMENSION(jpi,jpj) :: ztmp0, ztmp1, ztmp2 186 186 ! 187 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: & 188 & zeta_t, & ! stability parameter at height zt 189 & zsst, & ! to back up the initial bulk SST 190 & pdTc, & ! SST increment "dT" for cool-skin correction [K] 191 & pdTw, & ! SST increment "dT" for warm layer correction [K] 192 & zHwl ! depth of warm-layer [m] 187 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: zeta_t ! stability parameter at height zt 188 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: zsst ! to back up the initial bulk SST 189 ! 193 190 CHARACTER(len=40), PARAMETER :: crtnm = 'turb_coare3p6@sbcblk_algo_coare3p6' 194 191 !!---------------------------------------------------------------------------------- 195 IF ( kt == nit000 ) CALL COARE3P6_INIT(l_use_cs, l_use_wl)192 IF ( kt == nit000 ) CALL SBCBLK_ALGO_COARE3P6_INIT(l_use_cs, l_use_wl) 196 193 197 194 l_zt_equal_zu = .FALSE. -
NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk_algo_ecmwf.F90
r11845 r11962 44 44 PRIVATE 45 45 46 PUBLIC :: ECMWF_INIT, TURB_ECMWF46 PUBLIC :: SBCBLK_ALGO_ECMWF_INIT, TURB_ECMWF 47 47 48 48 !! ECMWF own values for given constants, taken form IFS documentation... … … 61 61 62 62 63 SUBROUTINE ecmwf_init(l_use_cs, l_use_wl)63 SUBROUTINE sbcblk_algo_ecmwf_init(l_use_cs, l_use_wl) 64 64 !!--------------------------------------------------------------------- 65 !! *** FUNCTION ecmwf_init ***65 !! *** FUNCTION sbcblk_algo_ecmwf_init *** 66 66 !! 67 67 !! INPUT : … … 77 77 ierr = 0 78 78 ALLOCATE ( dT_wl(jpi,jpj), Hz_wl(jpi,jpj), STAT=ierr ) 79 IF( ierr > 0 ) CALL ctl_stop( ' ECMWF_INIT => allocation of dT_wl & Hz_wl failed!' )79 IF( ierr > 0 ) CALL ctl_stop( ' SBCBLK_ALGO_ECMWF_INIT => allocation of dT_wl & Hz_wl failed!' ) 80 80 dT_wl(:,:) = 0._wp 81 81 Hz_wl(:,:) = rd0 ! (rd0, constant, = 3m is default for Zeng & Beljaars) … … 84 84 ierr = 0 85 85 ALLOCATE ( dT_cs(jpi,jpj), STAT=ierr ) 86 IF( ierr > 0 ) CALL ctl_stop( ' ECMWF_INIT => allocation of dT_cs failed!' )86 IF( ierr > 0 ) CALL ctl_stop( ' SBCBLK_ALGO_ECMWF_INIT => allocation of dT_cs failed!' ) 87 87 dT_cs(:,:) = -0.25_wp ! First guess of skin correction 88 88 END IF 89 END SUBROUTINE ecmwf_init89 END SUBROUTINE sbcblk_algo_ecmwf_init 90 90 91 91 … … 182 182 LOGICAL :: l_zt_equal_zu = .FALSE. ! if q and t are given at same height as U 183 183 ! 184 REAL(wp), DIMENSION(jpi,jpj) :: & 185 & u_star, t_star, q_star, & 186 & dt_zu, dq_zu, & 187 & znu_a, & !: Nu_air, Viscosity of air 188 & Linv, & !: 1/L (inverse of Monin Obukhov length... 189 & z0, z0t, z0q 190 ! 191 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: & 192 & zsst, & ! to back up the initial bulk SST 193 & pdTc, & ! SST increment "dT" for cool-skin correction [K] 194 & pdTw ! SST increment "dT" for warm layer correction [K] 195 ! 196 REAL(wp), DIMENSION(jpi,jpj) :: func_m, func_h 197 REAL(wp), DIMENSION(jpi,jpj) :: ztmp0, ztmp1, ztmp2 184 REAL(wp), DIMENSION(jpi,jpj) :: u_star, t_star, q_star 185 REAL(wp), DIMENSION(jpi,jpj) :: dt_zu, dq_zu 186 REAL(wp), DIMENSION(jpi,jpj) :: znu_a !: Nu_air, Viscosity of air 187 REAL(wp), DIMENSION(jpi,jpj) :: Linv !: 1/L (inverse of Monin Obukhov length... 188 REAL(wp), DIMENSION(jpi,jpj) :: z0, z0t, z0q 189 ! 190 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: zsst ! to back up the initial bulk SST 191 ! 192 REAL(wp), DIMENSION(jpi,jpj) :: func_m, func_h 193 REAL(wp), DIMENSION(jpi,jpj) :: ztmp0, ztmp1, ztmp2 198 194 CHARACTER(len=40), PARAMETER :: crtnm = 'turb_ecmwf@sbcblk_algo_ecmwf.F90' 199 195 !!---------------------------------------------------------------------------------- 200 196 201 IF ( kt == nit000 ) CALL ECMWF_INIT(l_use_cs, l_use_wl)197 IF ( kt == nit000 ) CALL SBCBLK_ALGO_ECMWF_INIT(l_use_cs, l_use_wl) 202 198 203 199 l_zt_equal_zu = .FALSE. … … 367 363 CALL WL_ECMWF( Qsw, ztmp1, u_star, zsst ) 368 364 !! Updating T_s and q_s !!! 369 T_s(:,:) = zsst(:,:) + dT_wl(:,:)*tmask(:,:,1) 365 T_s(:,:) = zsst(:,:) + dT_wl(:,:)*tmask(:,:,1) ! 370 366 IF( l_use_cs ) T_s(:,:) = T_s(:,:) + dT_cs(:,:)*tmask(:,:,1) 371 367 q_s(:,:) = rdct_qsat_salt*q_sat(MAX(T_s(:,:), 200._wp), slp(:,:)) -
NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk_phy.F90
r11845 r11962 531 531 zCe = zz0*pqst(ji,jj)/zdq 532 532 533 CALL BULK_FORMULA( pzu, pTs(ji,jj), pqs(ji,jj), pTa(ji,jj), pqa(ji,jj), zCd, zCh, zCe, & 534 & pwnd(ji,jj), pUb(ji,jj), pslp(ji,jj), & 533 CALL BULK_FORMULA( pzu, pTs(ji,jj), pqs(ji,jj), pTa(ji,jj), pqa(ji,jj), & 534 & zCd, zCh, zCe, & 535 & pwnd(ji,jj), pUb(ji,jj), pslp(ji,jj), & 535 536 & pTau(ji,jj), zQsen, zQlat ) 536 537 537 538 zTs2 = pTs(ji,jj)*pTs(ji,jj) 538 539 zQlw = emiss_w*(prlw(ji,jj) - stefan*zTs2*zTs2) ! Net longwave flux … … 545 546 END SUBROUTINE UPDATE_QNSOL_TAU 546 547 547 548 SUBROUTINE BULK_FORMULA_VCTR( pzu, pTs, pqs, pTa, pqa, pCd, pCh, pCe, pwnd, pUb, pslp, & 549 & pTau, pQsen, pQlat, pEvap, prhoa ) 548 549 SUBROUTINE BULK_FORMULA_VCTR( pzu, pTs, pqs, pTa, pqa, & 550 & pCd, pCh, pCe, & 551 & pwnd, pUb, pslp, & 552 & pTau, pQsen, pQlat, pEvap, prhoa ) 550 553 !!---------------------------------------------------------------------------------- 551 554 REAL(wp), INTENT(in) :: pzu ! height above the sea-level where all this takes place (normally 10m) … … 598 601 END SUBROUTINE BULK_FORMULA_VCTR 599 602 600 601 SUBROUTINE BULK_FORMULA_SCLR( pzu, pTs, pqs, pTa, pqa, pCd, pCh, pCe, pwnd, pUb, pslp, & 602 & pTau, pQsen, pQlat, pEvap, prhoa ) 603 604 SUBROUTINE BULK_FORMULA_SCLR( pzu, pTs, pqs, pTa, pqa, & 605 & pCd, pCh, pCe, & 606 & pwnd, pUb, pslp, & 607 & pTau, pQsen, pQlat, pEvap, prhoa ) 603 608 !!---------------------------------------------------------------------------------- 604 609 REAL(wp), INTENT(in) :: pzu ! height above the sea-level where all this takes place (normally 10m) -
NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk_skin_coare.F90
r11826 r11962 61 61 !!--------------------------------------------------------------------- 62 62 !! 63 !! Cool-skin parameterization, based on Fairall et al., 1996, revisited for COARE 3.6 (Fairall et al., 2019) 63 !! Cool-skin parameterization, based on Fairall et al., 1996, 64 !! revisited for COARE 3.6 (Fairall et al., 2019) 64 65 !! 65 66 !! Fairall, C. W., Bradley, E. F., Godfrey, J. S., Wick, G. A., -
NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/EXPREF/plot_station_asf.py
r11930 r11962 58 58 if narg != 3: 59 59 print 'Usage: '+sys.argv[0]+' <DIR_OUT_SASF> <year>'; sys.exit(0) 60 cdir_ in= sys.argv[1]61 cyyyy = sys.argv[2]60 cdir_data = sys.argv[1] 61 cyyyy = sys.argv[2] 62 62 63 63
Note: See TracChangeset
for help on using the changeset viewer.
