Changeset 1356 for trunk/LMDZ.TITAN/libf/phytitan/heating.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/heating.F (modified) (5 diffs)

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

@@ -1 @@
-       SUBROUTINE heating(dist,rmu0,fract,falbe,sol_htg,swnet,icld)
+       SUBROUTINE heating(dist,rmu0,fract,falbe,sol_htg,zswup,zswdn,
+     .                    icld)
 
 
@@ -23 +24 @@
 c      -------
 c sol_htg-----output-R- echauffement atmospherique (visible) (K/s)
-c swnet-------output-R- flux solaire net (+ vers le bas)     (W/m2)
+c zswup-------output-R- flux solaire upward  (+ vers le haut)     (W/m2)
+c zswdn-------output-R- flux solaire downward (+ vers le bas)     (W/m2)
 c
 c=======================================================================
@@ -51 +53 @@
 
       real sol_htg(klon,klev)
-      real swnet(klon,klev+1)
+      real zswup(klon,klev+1)
+      real zswdn(klon,klev+1)
       
 c   Local:
@@ -87 +90 @@
          fnetv  = 0.0
          sol_htg= 0.0
-         swnet  = 0.0
+         zswup  = 0.0
+         zswdn  = 0.0
 c pour sorties dans gfluxv...
          iprint = 0
@@ -101 +105 @@
                CALL sfluxv(iprint,ig,dist,falbe,icld)      ! #3
 
-               fnetv(ig,:) = fnetv(ig,:) *fract(ig)   ! >0 vers le haut
-c >0 vers le bas + conversion en W/m2:
-               swnet(ig,:) = -1.e-3*fnetv(ig,:)          
+               do K=1,NSPECV
+                 zswup(ig,:) = zswup(ig,:)+FUPV(ig,:,K)*fract(ig) ! >0 up
+                 zswdn(ig,:) = zswdn(ig,:)+FDV(ig,:,K) *fract(ig) ! >0 down
+               enddo
+               fnetv(ig,:) = fnetv(ig,:) *fract(ig)   ! >0 up
+
+c conversion en W/m2:
+               zswup(ig,:) = 1.e-3*zswup(ig,:)
+               zswdn(ig,:) = 1.e-3*zswdn(ig,:)
                
                DO j=1,nlayer