Changeset 3409 for trunk/LMDZ.PLUTO/libf/phypluto/dsolver.F

Timestamp:

Aug 20, 2024, 12:12:40 PM (3 months ago)

Author:

afalco

Message:

Pluto PCM:
Take into account zeros in aerosol optical properties.
AF

File:

• trunk/LMDZ.PLUTO/libf/phypluto/dsolver.F (modified) (4 diffs)

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

@@ -40 +40 @@
 
       L=2*NL
- 
+
 C     ************MIXED COEFFICENTS**********
 C     THIS VERSION AVOIDS SINGULARITIES ASSOC.
@@ -53 +53 @@
 
 C     EVEN TERMS
- 
+
       DO I=2,LM2,2
         N     = N+1
-        AF(I) = (E1(N)+E3(N))*(GAMA(N+1)-1.)       
+        AF(I) = (E1(N)+E3(N))*(GAMA(N+1)-1.)
         BF(I) = (E2(N)+E4(N))*(GAMA(N+1)-1.)
+        IF (BF(I).eq.0) THEN
+            BF(I) = 1e-16
+        END IF
         CF(I) = 2.0*(1.-GAMA(N+1)**2)
         DF(I) = (GAMA(N+1)-1.) * (CPM1(N+1) - CP(N)) +
      *            (1.-GAMA(N+1))* (CM(N)-CMM1(N+1))
       END DO
- 
+
       N   = 0
       LM1 = L-1
@@ -69 +72 @@
         AF(I) = 2.0*(1.-GAMA(N)**2)
         BF(I) = (E1(N)-E3(N))*(1.+GAMA(N+1))
+        IF (BF(I).eq.0) THEN
+            BF(I) = 1e-16
+        END IF
         CF(I) = (E1(N)+E3(N))*(GAMA(N+1)-1.)
         DF(I) = E3(N)*(CPM1(N+1) - CP(N)) + E1(N)*(CM(N) - CMM1(N+1))
       END DO
- 
+
       AF(L) = E1(NL)-RSF*E3(NL)
       BF(L) = E2(NL)-RSF*E4(NL)
+      IF (BF(L).eq.0) THEN
+        BF(L) = 1e-16
+      END IF
       CF(L) = 0.0
       DF(L) = BSURF-CP(NL)+RSF*CM(NL)
- 
+
       CALL DTRIDGL(L,AF,BF,CF,DF,XK)
- 
+
 C     ***UNMIX THE COEFFICIENTS****
 
@@ -98 +107 @@
 
    28 CONTINUE
- 
+
       RETURN
       END