Changeset 5099 for LMDZ6/branches/Amaury_dev/libf/phylmd/qsatseaw_1D.F90
- Timestamp:
- Jul 22, 2024, 9:29:09 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/qsatseaw_1D.F90
r5093 r5099 1 1 !------------------------------------------------------------------------------- 2 ! 2 3 3 ! ###################################### 4 4 REAL FUNCTION QSATSEAW_1D(PT,PP) 5 5 ! ###################################### 6 ! 6 7 7 !!**** *QSATW * - function to compute saturation vapor humidity from 8 8 !! temperature … … 12 12 ! The purpose of this function is to compute the saturation vapor 13 13 ! pressure from temperature over saline seawater 14 ! 15 ! 14 15 16 16 !!** METHOD 17 17 !! ------ … … 59 59 !! Original 6/04/2005 60 60 !------------------------------------------------------------------------------- 61 ! 61 62 62 !* 0. DECLARATIONS 63 63 ! ------------ 64 ! 64 65 65 USE MODD_CSTS 66 66 USE dimphy … … 70 70 71 71 !* 0.1 Declarations of arguments and results 72 ! 73 ! 72 73 74 74 REAL, DIMENSION(klon), INTENT(IN) :: PT ! Temperature 75 75 ! (Kelvin) … … 80 80 ! with respect to 81 81 ! water (kg/kg) 82 ! 82 83 83 !* 0.2 Declarations of local variables 84 ! 84 85 85 REAL, DIMENSION(SIZE(PT)) :: ZFOES ! saturation vapor 86 86 ! pressure 87 87 ! (Pascal) 88 ! 88 89 89 INTEGER :: JJ ! loop index 90 90 !REAL(KIND=JPRB) :: ZHOOK_HANDLE 91 91 !------------------------------------------------------------------------------- 92 ! 92 93 93 !IF (LHOOK) CALL DR_HOOK('MODE_THERMOS:QSATSEAW_1D',0,ZHOOK_HANDLE) 94 94 !DO JJ = 1, SIZE(PT) 95 95 !* 1. COMPUTE SATURATION VAPOR PRESSURE 96 96 ! --------------------------------- 97 ! 97 98 98 ZFOES = 0.98*EXP( XALPW - XBETAW/PT - XGAMW*LOG(PT) ) 99 99 ! vapor pressure reduction of 2% over saline seawater could have a significant 100 100 ! impact on the computation of surface latent heat flux under strong wind 101 101 ! conditions (Zeng et al, 1998). 102 ! 102 103 103 !* 2. COMPUTE SATURATION HUMIDITY 104 104 ! --------------------------- 105 ! 105 106 106 PQSAT = XRD/XRV*ZFOES/PP /(1.+(XRD/XRV-1.)*ZFOES/PP) 107 ! 107 108 108 !ENDDO 109 109 !IF (LHOOK) CALL DR_HOOK('MODE_THERMOS:QSATSEAW_1D',1,ZHOOK_HANDLE) 110 110 !------------------------------------------------------------------------------- 111 ! 111 112 112 END FUNCTION QSATSEAW_1D 113 113 !------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.