Changeset 1992 for LMDZ5/trunk/libf/phylmd/transp.F90
- Timestamp:
- Mar 5, 2014, 2:19:12 PM (11 years ago)
- File:
-
- 1 moved
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/transp.F90
r1988 r1992 1 ! 1 2 2 ! $Header$ 3 ! 4 SUBROUTINE transp (paprs,tsol, 5 e t, q, u, v, geom, 6 s vtran_e, vtran_q, utran_e, utran_q) 7 c 8 USE dimphy 9 IMPLICIT none 10 c====================================================================== 11 c Auteur(s): Z.X.Li (LMD/CNRS) 12 c Date: le 25 avril 1994 13 c Objet: Calculer le transport de l'energie et de la vapeur d'eau 14 c====================================================================== 15 c 16 cym#include "dimensions.h" 17 cym#include "dimphy.h" 18 #include "YOMCST.h" 19 c 20 REAL paprs(klon,klev+1), tsol(klon) 21 REAL t(klon,klev), q(klon,klev), u(klon,klev), v(klon,klev) 22 REAL utran_e(klon), utran_q(klon), vtran_e(klon), vtran_q(klon) 23 c 24 INTEGER i, l 25 c ------------------------------------------------------------------ 26 REAL geom(klon,klev), e 27 c ------------------------------------------------------------------ 28 DO i = 1, klon 29 utran_e(i) = 0.0 30 utran_q(i) = 0.0 31 vtran_e(i) = 0.0 32 vtran_q(i) = 0.0 33 ENDDO 34 c 35 DO l = 1, klev 36 DO i = 1, klon 37 e = RCPD*t(i,l) + RLVTT*q(i,l) + geom(i,l) 38 utran_e(i)=utran_e(i)+ u(i,l)*e*(paprs(i,l)-paprs(i,l+1))/RG 39 utran_q(i)=utran_q(i)+ u(i,l)*q(i,l) 40 . *(paprs(i,l)-paprs(i,l+1))/RG 41 vtran_e(i)=vtran_e(i)+ v(i,l)*e*(paprs(i,l)-paprs(i,l+1))/RG 42 vtran_q(i)=vtran_q(i)+ v(i,l)*q(i,l) 43 . *(paprs(i,l)-paprs(i,l+1))/RG 44 ENDDO 45 ENDDO 46 c 47 RETURN 48 END 3 4 SUBROUTINE transp(paprs, tsol, t, q, u, v, geom, vtran_e, vtran_q, utran_e, & 5 utran_q) 6 7 USE dimphy 8 IMPLICIT NONE 9 ! ====================================================================== 10 ! Auteur(s): Z.X.Li (LMD/CNRS) 11 ! Date: le 25 avril 1994 12 ! Objet: Calculer le transport de l'energie et de la vapeur d'eau 13 ! ====================================================================== 14 15 ! ym#include "dimensions.h" 16 ! ym#include "dimphy.h" 17 include "YOMCST.h" 18 19 REAL paprs(klon, klev+1), tsol(klon) 20 REAL t(klon, klev), q(klon, klev), u(klon, klev), v(klon, klev) 21 REAL utran_e(klon), utran_q(klon), vtran_e(klon), vtran_q(klon) 22 23 INTEGER i, l 24 ! ------------------------------------------------------------------ 25 REAL geom(klon, klev), e 26 ! ------------------------------------------------------------------ 27 DO i = 1, klon 28 utran_e(i) = 0.0 29 utran_q(i) = 0.0 30 vtran_e(i) = 0.0 31 vtran_q(i) = 0.0 32 END DO 33 34 DO l = 1, klev 35 DO i = 1, klon 36 e = rcpd*t(i, l) + rlvtt*q(i, l) + geom(i, l) 37 utran_e(i) = utran_e(i) + u(i, l)*e*(paprs(i,l)-paprs(i,l+1))/rg 38 utran_q(i) = utran_q(i) + u(i, l)*q(i, l)*(paprs(i,l)-paprs(i,l+1))/rg 39 vtran_e(i) = vtran_e(i) + v(i, l)*e*(paprs(i,l)-paprs(i,l+1))/rg 40 vtran_q(i) = vtran_q(i) + v(i, l)*q(i, l)*(paprs(i,l)-paprs(i,l+1))/rg 41 END DO 42 END DO 43 44 RETURN 45 END SUBROUTINE transp
Note: See TracChangeset
for help on using the changeset viewer.