New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 12021 for NEMO/branches/2019/dev_ASINTER-01-05_merged/src/OCE/SBC/sbcblk_phy.F90 – NEMO

Ignore:
Timestamp:
2019-12-02T10:02:07+01:00 (4 years ago)
Author:
gsamson
Message:

dev_ASINTER-01-05_merged: update sbcblk, bugfix with latent heat sign and update orca2_ice_abl namelist_cfg (tickets #2159 and #2131)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_ASINTER-01-05_merged/src/OCE/SBC/sbcblk_phy.F90

    r12015 r12021  
    2020 
    2121   USE dom_oce        ! ocean space and time domain 
    22    USE sbc_oce, ONLY  : ln_abl !!GS 
    2322   USE phycst         ! physical constants 
    2423 
     
    575574      INTEGER  :: ji, jj, jq     ! dummy loop indices 
    576575      !!---------------------------------------------------------------------------------- 
    577  
    578       !==  ABL formulation  ==!   multiplication by rho_air and turbulent fluxes computation done in ablstp 
    579       IF (ln_abl) THEN 
    580          DO jj = 1, jpj 
    581             DO ji = 1, jpi 
    582     
    583                zUrho = pUb(ji,jj) 
    584     
    585                pTau(ji,jj)  = zUrho * pCd(ji,jj)   !!GS: store cd_du in wind stress 
    586     
    587                zevap        = zUrho * pCe(ji,jj) 
    588                pQsen(ji,jj) = zUrho * pCh(ji,jj) 
    589                pQlat(ji,jj) = zevap                !!GS: store evap in latent 
    590     
    591             END DO 
    592          END DO 
    593       ELSE 
    594576         DO jj = 1, jpj 
    595577            DO ji = 1, jpi 
     
    617599            END DO 
    618600         END DO 
    619       END IF 
    620  
    621601   END SUBROUTINE BULK_FORMULA_VCTR 
    622602 
     
    650630      !!---------------------------------------------------------------------------------- 
    651631 
    652       !==  ABL formulation  ==!   multiplication by rho_air and turbulent fluxes computation done in ablstp 
    653       IF (ln_abl) THEN 
    654     
    655          zUrho = pUb 
    656     
    657          pTau  = zUrho * pCd 
    658     
    659          zevap = zUrho * pCe 
    660          pQsen = zUrho * pCh 
    661          pQlat = zevap 
    662  
    663       ELSE 
    664632         !! Need ztaa, absolute temperature at pzu (formula to estimate rho_air needs absolute temperature, not the potential temperature "pTa") 
    665633         ztaa = pTa ! first guess... 
     
    681649         IF( PRESENT(pEvap) ) pEvap = - zevap 
    682650         IF( PRESENT(prhoa) ) prhoa = zrho 
    683      END IF 
    684651 
    685652   END SUBROUTINE BULK_FORMULA_SCLR 
Note: See TracChangeset for help on using the changeset viewer.