Changeset 3936 for trunk/LMDZ.PLUTO/libf/phypluto/hazecloud.F90

Timestamp:

Oct 25, 2025, 3:21:13 PM (7 weeks ago)

Author:

tbertrand

Message:

PLUTO PCM : correcting a bug in hazecloud (wrong lyman alpha fluxes due to mu0 being negative during nighttime) + cleaning routines
TB

File:

• trunk/LMDZ.PLUTO/libf/phypluto/hazecloud.F90 (modified) (3 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/hazecloud.F90

@@ -164 +164 @@
          flym_ipm(1)=flym_ipm(1)*pfluxuv/5.
       ELSE IF (.not.diurnal) THEN
-         mu_sol(:) = mu0(:)
+         mu_sol(:) = max(mu0(:),0.)
          mu_ipm(:) = max(mu_sol(:), 0.5)
-         flym_ipm(:)= mu0(:)*72.5e10
+         flym_ipm(:)= mu_sol(:)*72.5e10
       ELSE ! case with full fit to Gladstone et al. results
 !     1)  get longitude/latitude (in radian) of anti-subsolar point (max de mu0 - 180)
@@ -182 +182 @@
         DO ig=1,ngrid
         ! calculation of cosinus of incident angle for IPM flux
-          mu_sol(ig) = mu0(ig)
+          mu_sol(ig) = max(mu0(ig),0.)
           mu_ipm(ig) = max(mu_sol(ig), 0.5)
-          IF (mu0(ig).LT.1.e-4) THEN ! Daytime
+          IF (mu_sol(ig).LE.0.) THEN ! Nighttime
            ! Distance to subsolar point
            dist=acos(sin(latitude(ig))*sin(latit)+cos(latitude(ig))*    &
@@ -191 +191 @@
            flym_ipm(ig)=(valmin_dl-valmin)/(90.**puis)*(dist)**puis  &
                                 +valmin
-          ELSE ! Nightime
-           flym_ipm(ig)= mu0(ig)*(valmax-valmin_dl)+valmin_dl
+          ELSE ! Daytime
+           flym_ipm(ig)= mu_sol(ig)*(valmax-valmin_dl)+valmin_dl
           ENDIF
           ! Proportional to lyman alpha solar flux (reference 2015 : 5e11 ph/cm2/s) ?