Changeset 393 for LMDZ.3.3/branches/rel-LF/libf/phylmd/radlwsw.F

Timestamp:

Jul 19, 2002, 11:51:34 AM (22 years ago)

Author:

lmdzadmin

Message:

Modifications de JLD sur la conservation de l'energie
On supprime les modifs de Pascale sur le cdrag, elles refroidissaient trop
l'atmosphere
LF

File:

• LMDZ.3.3/branches/rel-LF/libf/phylmd/radlwsw.F (modified) (3 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/radlwsw.F

@@ -39 +39 @@
 #include "dimphy.h"
 #include "raddim.h"
+#include "YOETHF.h"
 c
       real rmu0(klon), fract(klon), dist
@@ -88 +89 @@
       REAL*8 ztopsw0(kdlon), ztoplw0(kdlon)
       REAL*8 zsolsw0(kdlon), zsollw0(kdlon)
-
+      REAL*8 zznormcp
 c
 c-------------------------------------------
@@ -215 +216 @@
       ENDDO
       DO k = 1, kflev
+c      DO i = 1, kdlon
+c         heat(iof+i,k) = zheat(i,k)
+c         cool(iof+i,k) = zcool(i,k)
+c         heat0(iof+i,k) = zheat0(i,k)
+c         cool0(iof+i,k) = zcool0(i,k)
+c      ENDDO
       DO i = 1, kdlon
-         heat(iof+i,k) = zheat(i,k)
-         cool(iof+i,k) = zcool(i,k)
-         heat0(iof+i,k) = zheat0(i,k)
-         cool0(iof+i,k) = zcool0(i,k)
+C        scale factor to take into account the difference between
+C        dry air and watter vapour scpecific heat capacity
+         zznormcp=1.0+RVTMP2*PWV(i,k)
+         heat(iof+i,k) = zheat(i,k)/zznormcp
+         cool(iof+i,k) = zcool(i,k)/zznormcp
+         heat0(iof+i,k) = zheat0(i,k)/zznormcp
+         cool0(iof+i,k) = zcool0(i,k)/zznormcp
       ENDDO
       ENDDO