Changeset 3936 for trunk/LMDZ.PLUTO/libf/phypluto/stelang.F

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/stelang.F (modified) (3 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/stelang.F

@@ -85 +85 @@
 c----- AF24: replaced flat param by 0
 
-      rap = 1./((1.-0)**2)
+c      rap = 1./((1.-0)**2)
 
  100  CONTINUE
@@ -98 +98 @@
 C
       DO jl=1,kgrid
-        ztim1=psilat(jl)*ptim1*rap
+        ztim1=psilat(jl)*ptim1 !*rap
         ztim2=pcolat(jl)*ptim2
         ztim3=pcolat(jl)*ptim3
         pmu0(jl)=ztim1+ztim2*pcolon(jl)+ztim3*psilon(jl)
-        pmu0(jl)=pmu0(jl)/SQRT(pcolat(jl)**2+(rap**2)*(psilat(jl)**2))
+c       pmu0(jl)=pmu0(jl)/SQRT(pcolat(jl)**2+(rap**2)*(psilat(jl)**2))
 
       ENDDO
@@ -112 +112 @@
         IF (pmu0(jl).gt.0.) THEN
           pfract(jl)=1.
-c       pmu0(jl)=sqrt(1224.*pmu0(jl)*pmu0(jl)+1.)/35.
-      ELSE
-c       pmu0(jl)=0.
-        pfract(jl)=0.
+c          pmu0(jl)=sqrt(1224.*pmu0(jl)*pmu0(jl)+1.)/35.
+        ELSE
+c          pmu0(jl)=0.
+          pfract(jl)=0.
         ENDIF
       ENDDO