Changeset 3259

Timestamp:

Mar 11, 2024, 6:43:35 PM (9 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.

Location:

trunk/LMDZ.GENERIC

Files:

• changelog.txt (modified) (1 diff)
• libf/phystd/callcorrk.F90 (modified) (2 diffs)

trunk/LMDZ.GENERIC/changelog.txt

@@ -1874 +1874 @@
 Removed the outdated 'How_to_add_an_aerosol_to_the_GCM.pdf' file, the information is now in the wiki page : https://lmdz-forge.lmd.jussieu.fr/mediawiki/Planets/index.php/Advanced_Use_of_the_GCM#How_to_Change_the_Aerosols_Optical_Properties
 Add 'optpropgen' routines in the SVN (code written by JB Madeleine and Franck Montmessin, to generate aerosol optical properties)
+
+== 11/03/2024 == GM
+
+Fix a bug with global1d in callcorrk. 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. 

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.