Changeset 3794 for trunk/LMDZ.VENUS/libf/phyvenus/sfluxv.F

Timestamp:

Jun 5, 2025, 3:39:30 PM (3 days ago)

Author:

slebonnois

Message:

SL: VenusPCM - update of CIA in solar module (adapted from Generic PCM) + potential output of band-resolved net flux

File:

• trunk/LMDZ.VENUS/libf/phyvenus/sfluxv.F (modified) (5 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/sfluxv.F

@@ -1 +1 @@
       SUBROUTINE SFLUXV(DTAUV,TAUV,TAUCUMV,RSFV,DWNV,WBARV,COSBV,
      *                  UBAR0,STEL,NFLUXTOPV,FLUXTOPVDN,
-     *                  NFLUXOUTV_nu,NFLUXGNDV_nu,
+     *                  NFLUXV_nu,NFLUXGNDV_nu,
      *                  FMNETV,FLUXUPV,FLUXDNV,FZEROV,taugsurf)
 
@@ -19 +19 @@
       real*8 FLUXUPV(L_NLAYRAD), FLUXDNV(L_NLAYRAD)
       real*8 NFLUXTOPV, FLUXUP, FLUXDN,FLUXTOPVDN
-      real*8 NFLUXOUTV_nu(L_NSPECTV)
+      real*8 NFLUXV_nu(L_NLAYRAD,L_NSPECTV)
       real*8 NFLUXGNDV_nu(L_NSPECTV)
 
@@ -40 +40 @@
 
       DO NW=1,L_NSPECTV
-         NFLUXOUTV_nu(NW)=0.0
+         NFLUXV_nu(:,NW)=0.0
          NFLUXGNDV_nu(NW)=0.0
       END DO
@@ -115 +115 @@
           END DO
 
-c     band-resolved flux leaving TOA (RDW)
-          NFLUXOUTV_nu(NW) = NFLUXOUTV_nu(NW)
-     *      +FLUXUP*GWEIGHT(NG)*(1.0-FZEROV(NW))
+c     band-resolved net flux 
+          NFLUXV_nu(:,NW) = NFLUXV_nu(:,NW)
+     *      +(FMUPV-FMDV)*GWEIGHT(NG)*(1.0-FZEROV(NW))
 
 c     band-resolved flux at ground (RDW)
@@ -173 +173 @@
         END DO
 
-c     band-resolved flux leaving TOA (RDW)
-          NFLUXOUTV_nu(NW) = NFLUXOUTV_nu(NW)
-     *      +FLUXUP*FZERO
+c     band-resolved net flux 
+          NFLUXV_nu(:,NW) = NFLUXV_nu(:,NW)
+     *      +(FMUPV-FMDV)*FZERO
 
 c     band-resolved flux at ground (RDW)