Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real, | intent(in) | :: | PZ_O_LMO |
FUNCTION PAULSON_PSIH_0D(PZ_O_LMO)
IMPLICIT NONE
REAL, INTENT(IN) :: PZ_O_LMO
REAL :: PAULSON_PSIH_0D
!
REAL :: ZY
REAL(KIND=JPRB) :: ZHOOK_HANDLE
IF (LHOOK) CALL DR_HOOK('MODE_SBLS:PAULSON_PSIH_0D',0,ZHOOK_HANDLE)
ZY=1.
IF ( PZ_O_LMO < 0. ) THEN
ZY=(1.-9.*PZ_O_LMO)**(0.5)
PAULSON_PSIH_0D = LOG( (1.+ZY)/2. )
ELSE
PAULSON_PSIH_0D = - 4.7 * PZ_O_LMO / 0.74
END IF
IF (LHOOK) CALL DR_HOOK('MODE_SBLS:PAULSON_PSIH_0D',1,ZHOOK_HANDLE)
END FUNCTION PAULSON_PSIH_0D