Changeset 3425 for LMDZ6/branches/IPSLCM6.0.15/libf/phylmd/radlwsw_m.F90

Timestamp:

Dec 13, 2018, 1:07:06 PM (6 years ago)

Author:

Laurent Fairhead

Message:

Followup to r3408 modifications for VolMIP
NL

File:

• LMDZ6/branches/IPSLCM6.0.15/libf/phylmd/radlwsw_m.F90 (modified) (1 diff)

LMDZ6/branches/IPSLCM6.0.15/libf/phylmd/radlwsw_m.F90

@@ -1037 +1037 @@
            zcool0(i,k)=(ZLWFT0_i(i,k)-ZLWFT0_i(i,k+1))*RDAY*RG/RCPD/PDP(i,k)
            IF(ok_volcan) THEN
-              zheat_volc(i,k)=(ZSWADAERO(i,k+1)-ZSWADAERO(i,k))*RDAY*RG/RCPD/PDP(i,k) !NL
-              zcool_volc(i,k)=(ZLWADAERO(i,k)-ZLWADAERO(i,k+1))*RDAY*RG/RCPD/PDP(i,k) !NL
+              zheat_volc(i,k)=(ZSWADAERO(i,k+1)-ZSWADAERO(i,k))*RG/RCPD/PDP(i,k) !NL
+              zcool_volc(i,k)=(ZLWADAERO(i,k)-ZLWADAERO(i,k+1))*RG/RCPD/PDP(i,k) !NL
            ENDIF
 !          print *,'heat cool heat0 cool0 ',zheat(i,k),zcool(i,k),zheat0(i,k),zcool0(i,k)