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

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

@@ -37 +37 @@
       USE comgeomphy
       USE phys_state_var_mod, only: falbe,heat,cool,radsol,
-     .      topsw,toplw,solsw,sollw,sollwdown,lwnet,swnet
+     .      topsw,toplw,solsw,sollw,sollwdown,lwnet,swnet,
+     .      lwup,lwdn,swup,swdn
       USE write_field_phy
        IMPLICIT none
@@ -60 +61 @@
       real zheatc(klon,klev), zcoolc(klon,klev)
       real zheatp(klon,klev), zcoolp(klon,klev)
-      REAL zswnetc(klon,klev+1),zlwnetp(klon,klev+1)
-      REAL zswnetp(klon,klev+1),zlwnetc(klon,klev+1)
+      REAL zswupc(klon,klev+1),zlwupc(klon,klev+1)
+      REAL zswupp(klon,klev+1),zlwupp(klon,klev+1)
+      REAL zswdnc(klon,klev+1),zlwdnc(klon,klev+1)
+      REAL zswdnp(klon,klev+1),zlwdnp(klon,klev+1)
       REAL zsollwdownc(klon),zsollwdownp(klon)
       INTEGER icld
@@ -124 +127 @@
        IF (clouds.eq.1) THEN
          ICLD = 1   ! colonne avec nuages
-         CALL heating(dist,rmu0,fract,falbe,zheatc,zswnetc,icld)
+         CALL heating(dist,rmu0,fract,falbe,zheatc,zswupc,zswdnc,icld)
        ELSE
          zheatc  = 0.
-         zswnetc = 0.
+         zswupc = 0.
+         zswdnc = 0.
        ENDIF 
        ICLD = 0   ! colonne sans nuages
-       CALL heating(dist,rmu0,fract,falbe,zheatp,zswnetp,icld)
+       CALL heating(dist,rmu0,fract,falbe,zheatp,zswupp,zswdnp,icld)
 
 c inversion de l'axe vertical
@@ -141 +145 @@
        do l=1,klev+1
          do i=1,klon
-           swnet(i,l)=zswnetc(i,klev+2-l)*xnuf +
-     &                zswnetp(i,klev+2-l)*(1.-xnuf)
+           swup(i,l) =zswupc(i,klev+2-l)*xnuf +
+     &                zswupp(i,klev+2-l)*(1.-xnuf)
+           swdn(i,l) =zswdnc(i,klev+2-l)*xnuf +
+     &                zswdnp(i,klev+2-l)*(1.-xnuf)
+           swnet(i,l)=swdn(i,l)-swup(i,l)
          enddo
        enddo
@@ -157 +164 @@
        IF (clouds.eq.1) THEN
          ICLD = 1
-         CALL cooling(klon,klev+1,zp,zt,zz,zcoolc,zlwnetc,zsollwdownc,
-     &   icld)
+         CALL cooling(klon,klev+1,zp,zt,zz,zcoolc,zlwupc,zlwdnc,
+     &   zsollwdownc,icld)
        ELSE
          zcoolc      = 0.
-         zlwnetc     = 0.
+         zlwupc      = 0.
+         zlwdnc      = 0.
          zsollwdownc = 0.
        ENDIF
        ICLD = 0
-       CALL cooling(klon,klev+1,zp,zt,zz,zcoolp,zlwnetp,zsollwdownp,
-     & icld)
+       CALL cooling(klon,klev+1,zp,zt,zz,zcoolp,zlwupp,zlwdnp,
+     & zsollwdownp,icld)
 
 c inversion de l'axe vertical
@@ -177 +185 @@
        do l=1,klev+1
          do i=1,klon
-           lwnet(i,l)=zlwnetc(i,klev+2-l)*xnuf +
-     &                zlwnetp(i,klev+2-l)*(1.-xnuf)
+           lwup(i,l) =zlwupc(i,klev+2-l)*xnuf +
+     &                zlwupp(i,klev+2-l)*(1.-xnuf)
+           lwdn(i,l) =zlwdnc(i,klev+2-l)*xnuf +
+     &                zlwdnp(i,klev+2-l)*(1.-xnuf)
+           lwnet(i,l)=lwup(i,l)-lwdn(i,l)
          enddo
        enddo