Changeset 2503 for trunk/LMDZ.VENUS

Timestamp:

Apr 23, 2021, 4:41:19 PM (4 years ago)

Author:

emillour

Message:

Venus GCM:
Small correction in the computation of SW fluxes for the
very specific case when running without a diurnal cycle.
VB+EM

File:

• trunk/LMDZ.VENUS/libf/phyvenus/sw_venus_rh.F (modified) (4 diffs)

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

@@ -4 +4 @@
      S              PTOPSW,PSOLSW,ZFSNET)
       
-      use dimphy
+      use dimphy, only: klev
       use cpdet_phy_mod, only: cpdet
       IMPLICIT none
 
-#include "YOMCST.h"
+      include "YOMCST.h"
+      include "clesphys.h"
 C
 C     ------------------------------------------------------------------
@@ -32 +33 @@
 c inputs
 
-      REAL   PRMU0  ! COSINE OF ZENITHAL ANGLE
-      REAL   PFRAC  ! fraction de la journee
-      REAL   latdeg ! |latitude| (in degrees)
-      REAL   PPB(klev+1)  ! inter-couches PRESSURE (bar)
-      REAL   PPA(klev)
-      REAL   pt(klev)     ! mid-layer temperature
+      REAL,INTENT(IN) ::   PRMU0  ! COSINE OF ZENITHAL ANGLE
+      REAL,INTENT(IN) ::   PFRAC  ! fraction de la journee
+      REAL,INTENT(IN) ::   latdeg ! |latitude| (in degrees)
+      REAL,INTENT(IN) ::   PPB(klev+1)  ! inter-couches PRESSURE (bar)
+      REAL,INTENT(IN) ::   PPA(klev)
+      REAL,INTENT(IN) ::   pt(klev)     ! mid-layer temperature
 C
 c output
 
-      REAL   PHEAT(klev)  ! SHORTWAVE HEATING (K/s) within each layer
+      REAL,INTENT(OUT) ::   PHEAT(klev)  ! SHORTWAVE HEATING (K/s) within each layer
       REAL   PHEATPPA(klev)
-      REAL   PTOPSW       ! SHORTWAVE FLUX AT T.O.A. (net)
-      REAL   PSOLSW       ! SHORTWAVE FLUX AT SURFACE (net)
-      REAL   ZFSNET(klev+1) ! net solar flux at ppb levels
+      REAL,INTENT(OUT) ::   PTOPSW       ! SHORTWAVE FLUX AT T.O.A. (net)
+      REAL,INTENT(OUT) ::   PSOLSW       ! SHORTWAVE FLUX AT SURFACE (net)
+      REAL,INTENT(OUT) ::   ZFSNET(klev+1) ! net solar flux at ppb levels
 
 C
@@ -136 +137 @@
       nsza0(:)=2
 
-      do nsza=1,nszarh
-         if (solza(nsza,nlat0-1).le.sza0) then
+      if (.not.cycle_diurne) then
+        ! without a diurnal cycle, no need for any elaborate weights of sza
+        factsza(1)=1
+        factsza(2)=0
+      else
+        ! standard case with diurnal cycle
+        do nsza=1,nszarh
+          if (solza(nsza,nlat0-1).le.sza0) then
               nsza0(1) = nsza+1
-         endif
-      enddo
-      if (nsza0(1).ne.nszarh+1) then
-          factsza(1) = (sza0-solza(nsza0(1)-1,nlat0-1))/
-     .        (solza(nsza0(1),nlat0-1)-solza(nsza0(1)-1,nlat0-1))
-      else
-          factsza(1) = min((sza0-solza(nszarh,nlat0-1))/
-     .         (90.-solza(nszarh,nlat0-1)), 1.)
-      endif
-      if (nlat0.ne.nlatrh+1) then
-       do nsza=1,nszarh
-         if (solza(nsza,nlat0).le.sza0) then
+          endif
+        enddo
+        if (nsza0(1).ne.nszarh+1) then
+            factsza(1) = (sza0-solza(nsza0(1)-1,nlat0-1))/
+     .          (solza(nsza0(1),nlat0-1)-solza(nsza0(1)-1,nlat0-1))
+        else
+            factsza(1) = min((sza0-solza(nszarh,nlat0-1))/
+     .           (90.-solza(nszarh,nlat0-1)), 1.)
+        endif
+        if (nlat0.ne.nlatrh+1) then
+          do nsza=1,nszarh
+            if (solza(nsza,nlat0).le.sza0) then
               nsza0(2) = nsza+1
-         endif
-       enddo
-       if (nsza0(2).eq.nszarh+1) then
-          factsza(2) = min((sza0-solza(nszarh,nlat0))/
-     .         (90.-solza(nszarh,nlat0)), 1.)
-       elseif ((nsza0(2).eq.2).and.(solza(1,nlat0).gt.sza0)) then
-          factsza(2) = 0.
-       else
-          factsza(2) = (sza0-solza(nsza0(2)-1,nlat0))/
-     .        (solza(nsza0(2),nlat0)-solza(nsza0(2)-1,nlat0))
-       endif
-      else
-        nsza0(2)   = nszarh+1
-        factsza(2) = 1.
-      endif
+            endif
+          enddo
+          if (nsza0(2).eq.nszarh+1) then
+             factsza(2) = min((sza0-solza(nszarh,nlat0))/
+     .           (90.-solza(nszarh,nlat0)), 1.)
+          elseif ((nsza0(2).eq.2).and.(solza(1,nlat0).gt.sza0)) then
+            factsza(2) = 0.
+          else
+            factsza(2) = (sza0-solza(nsza0(2)-1,nlat0))/
+     .          (solza(nsza0(2),nlat0)-solza(nsza0(2)-1,nlat0))
+          endif
+        else
+          nsza0(2)   = nszarh+1
+          factsza(2) = 1.
+        endif ! of if (nlat0.ne.nlatrh+1)
+      endif ! of if (.not.cycle_diurne)
+
 c Pressure levels
 c ---------------
@@ -319 +328 @@
      
 
-      return
       end