Changeset 3259 for trunk/LMDZ.GENERIC/libf/phystd

Timestamp:

Mar 11, 2024, 6:43:35 PM (11 months ago)

Author:

gmilcareck

Message:

Fix a bug with global1d. When using global1d, the latitude used when global1d=false was still taken into account at global1d=true in the conditions for the presence of solar flux.

File:

• trunk/LMDZ.GENERIC/libf/phystd/callcorrk.F90 (modified) (2 diffs)

trunk/LMDZ.GENERIC/libf/phystd/callcorrk.F90

@@ -1137 +1137 @@
 !-----------------------------------------------------------------------
 
-         if(fract(ig) .ge. 1.0e-4) then ! Only during daylight.
+         if((fract(ig) .ge. 1.0e-4).or.(global1d)) then ! Only during daylight.
             if((ngrid.eq.1).and.(global1d))then
                do nw=1,L_NSPECTV
@@ -1171 +1171 @@
 
          ! Equivalent Albedo Calculation (for OUTPUT). MT2015
-         if(fract(ig) .ge. 1.0e-4) then ! equivalent albedo makes sense only during daylight.       
+         if((fract(ig) .ge. 1.0e-4).or.(global1d)) then ! equivalent albedo makes sense only during daylight.       
             surface_stellar_flux=sum(nfluxgndv_nu(1:L_NSPECTV))      
             if(surface_stellar_flux .gt. 1.0e-3) then ! equivalent albedo makes sense only if the stellar flux received by the surface is positive.