Changeset 5099 for LMDZ6/branches/Amaury_dev/libf/phylmd/qsat_seawater.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/qsat_seawater.F90
r5093 r5099 1 1 !------------------------------------------------------------------------------- 2 ! 2 3 3 ! ###################################### 4 4 REAL FUNCTION QSAT_SEAWATER(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 67 67 USE indice_sol_mod 68 ! 68 69 69 IMPLICIT NONE 70 ! 70 71 71 !* 0.1 Declarations of arguments and results 72 ! 73 ! 72 74 73 75 74 REAL, DIMENSION(klon), INTENT(IN) :: PT ! Temperature … … 82 81 ! water (kg/kg) 83 82 84 !85 83 !* 0.2 Declarations of local variables 86 ! 84 87 85 REAL, DIMENSION(SIZE(PT)) :: ZFOES ! saturation vapor 88 86 ! pressure 89 87 ! (Pascal) 90 ! 88 91 89 REAL, DIMENSION(SIZE(PT)) :: ZWORK1 92 90 REAL :: ZWORK2 93 91 !REAL(KIND=JPRB) :: ZHOOK_HANDLE 94 92 !------------------------------------------------------------------------------- 95 ! 93 96 94 !IF (LHOOK) CALL DR_HOOK('MODE_THERMOS:QSATSEAW_1D',0,ZHOOK_HANDLE) 97 ! 95 98 96 !ZFOES = 1 !PSAT(PT(:)) 99 97 !ZFOES = 0.98*ZFOES … … 103 101 ! impact on the computation of surface latent heat flux under strong wind 104 102 ! conditions (Zeng et al, 1998). 105 ! 103 106 104 ZWORK1 = ZFOES/PP 107 105 ZWORK2 = XRD/XRV 108 ! 106 109 107 !* 2. COMPUTE SATURATION HUMIDITY 110 108 ! --------------------------- 111 ! 109 112 110 PQSAT = ZWORK2*ZWORK1 / (1.+(ZWORK2-1.)*ZWORK1) 113 ! 111 114 112 !IF (LHOOK) CALL DR_HOOK('MODE_THERMOS:QSATSEAW_1D',1,ZHOOK_HANDLE) 115 113 !------------------------------------------------------------------------------- 116 ! 114 117 115 END FUNCTION QSAT_SEAWATER 118 ! 116 119 117 !-------------------------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.