Changeset 1356 for trunk/LMDZ.TITAN/libf/phytitan/cooling.F

Timestamp:

Oct 8, 2014, 9:26:28 AM (10 years ago)

Author:

slebonnois

Message:

SL: update to newstart/start2archive tools in Venus+Titan / additional diagnostics in radiative fluxes for Titan

File:

• trunk/LMDZ.TITAN/libf/phytitan/cooling.F (modified) (6 diffs)

trunk/LMDZ.TITAN/libf/phytitan/cooling.F

@@ -1 @@
-      SUBROUTINE COOLING(NG,NL,PRESS,TEMP,Z,Q0,lwnet,pfluxi,icld)
+      SUBROUTINE COOLING(NG,NL,PRESS,TEMP,Z,Q0,zlwup,zlwdn,pfluxi,icld)
 
 c=======================================================================
@@ -36 +36 @@
 c
 c      q0(nl-1)         radiative cooling in K/sec
-c      lwnet(nl)        net fluxes, (+) upward
+c      zlwup(nl)         up fluxes,   (+) upward
+c      zlwdn(nl)         down fluxes, (+) downward
 c      pfluxi          IR descendant a la surface (+ vers le bas)
 c
@@ -73 +74 @@
       REAL PRESS(NG,NL),TEMP(NG,NL)
       REAL Z(NG,NL),Q0(NG,NL-1)
-      REAL lwnet(NG,NL),UBARI2
+      REAL zlwup(NG,NL),zlwdn(NG,NL),UBARI2
       real pfluxi(NG)
 
@@ -147 +148 @@
       UBARI2=UBARI
 
-C ZERO THE NET FLUXES
+C ZERO THE FLUXES
          Q0    = 0.0
-         lwnet = 0.0 
+         zlwup = 0.0 
+         zlwdn = 0.0 
 
 c-----------------------------------------------------------------------
@@ -275 +277 @@
 3520      CONTINUE
 
-c compute the net IR flux, (+) upward:
+c compute the up (+ upward) and down (+ downward) IR fluxes:
 c 
           DO J=1,NL
           DO ig=1,NG
-             lwnet(ig,J)= lwnet(ig,J)+ DWNI(K)*(FUPI(ig,J)-FDI(ig,J))
+             zlwup(ig,J)= zlwup(ig,J)+ DWNI(K)*FUPI(ig,J)
+             zlwdn(ig,J)= zlwdn(ig,J)+ DWNI(K)*FDI(ig,J)
           ENDDO
           ENDDO
@@ -312 +315 @@
       DO 3550 ig=1,NG
          pfluxi(ig)  = 1.e-3*pfluxi(ig)
-         lwnet(ig,:) = 1.e-3*lwnet(ig,:)
+         zlwup(ig,:) = 1.e-3*zlwup(ig,:)
+         zlwdn(ig,:) = 1.e-3*zlwdn(ig,:)
 3550  CONTINUE