- Timestamp:
- 2019-12-02T10:02:07+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_ASINTER-01-05_merged/src/OCE/SBC/sbcblk_phy.F90
r12015 r12021 20 20 21 21 USE dom_oce ! ocean space and time domain 22 USE sbc_oce, ONLY : ln_abl !!GS23 22 USE phycst ! physical constants 24 23 … … 575 574 INTEGER :: ji, jj, jq ! dummy loop indices 576 575 !!---------------------------------------------------------------------------------- 577 578 !== ABL formulation ==! multiplication by rho_air and turbulent fluxes computation done in ablstp579 IF (ln_abl) THEN580 DO jj = 1, jpj581 DO ji = 1, jpi582 583 zUrho = pUb(ji,jj)584 585 pTau(ji,jj) = zUrho * pCd(ji,jj) !!GS: store cd_du in wind stress586 587 zevap = zUrho * pCe(ji,jj)588 pQsen(ji,jj) = zUrho * pCh(ji,jj)589 pQlat(ji,jj) = zevap !!GS: store evap in latent590 591 END DO592 END DO593 ELSE594 576 DO jj = 1, jpj 595 577 DO ji = 1, jpi … … 617 599 END DO 618 600 END DO 619 END IF620 621 601 END SUBROUTINE BULK_FORMULA_VCTR 622 602 … … 650 630 !!---------------------------------------------------------------------------------- 651 631 652 !== ABL formulation ==! multiplication by rho_air and turbulent fluxes computation done in ablstp653 IF (ln_abl) THEN654 655 zUrho = pUb656 657 pTau = zUrho * pCd658 659 zevap = zUrho * pCe660 pQsen = zUrho * pCh661 pQlat = zevap662 663 ELSE664 632 !! Need ztaa, absolute temperature at pzu (formula to estimate rho_air needs absolute temperature, not the potential temperature "pTa") 665 633 ztaa = pTa ! first guess... … … 681 649 IF( PRESENT(pEvap) ) pEvap = - zevap 682 650 IF( PRESENT(prhoa) ) prhoa = zrho 683 END IF684 651 685 652 END SUBROUTINE BULK_FORMULA_SCLR
Note: See TracChangeset
for help on using the changeset viewer.