Changeset 2879

Timestamp:

Jan 25, 2023, 1:51:23 PM (22 months ago)

Author:

llange

Message:

MARS PCM
The condition to call param_slope has been completed to prevent a
division by zero (bug that appeared during a twisted case with the 1D
GCM)
LL

File:

• trunk/LMDZ.MARS/libf/phymars/physiq_mod.F (modified) (4 diffs)

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -888 +888 @@
 c       Local Solar zenith angle
 c       ~~~~~~~~~~~~~~~~~~~~~~~~
+        write(*,*) 'mu0=',mu0(:) 
         CALL orbite(zls,dist_sol,declin)
-
+        write(*,*) 'mu0=',mu0(:) 
         IF (diurnal) THEN
             ztim1=SIN(declin)
@@ -901 +902 @@
             CALL mucorr(ngrid,declin,latitude,mu0,fract,10000.,rad)
         ENDIF ! of IF (diurnal)
-
+        write(*,*) 'mu0=',mu0(:),zls,declin
          IF( MOD(icount-1,iradia).EQ.0) THEN
 
@@ -1066 +1067 @@
 c          ---------------------------------------------------------
            IF(callslope) THEN
-
+          write(*,*) 'Flux=',fluxsurf_dn_sw(:,1),fluxsurf_dn_sw(:,2)
+          write(*,*) 'mu0=',mu0(:),zls,declin
             print *, 'Slope scheme is on and computing...'
             DO ig=1,ngrid  
@@ -1081 +1083 @@
                  sl_fl0 = fluxsurf_dn_sw(ig,l)
                  sl_di0 = 0.
-                 if (mu0(ig) .gt. 0.) then
+                 if ((mu0(ig) .gt. 0.).and.(ztim1.gt.0.)) then
                   sl_di0 = mu0(ig)*(exp(-sl_tau/mu0(ig)))
                   sl_di0 = sl_di0*flux_1AU/dist_sol/dist_sol