Changeset 2151 for trunk/LMDZ.MARS/libf/phymars/nlte_aux.F

Timestamp:

Sep 5, 2019, 3:33:16 PM (5 years ago)

Author:

aslmd

Message:

changed old functions dexp dlog in generic exp and log

File:

• trunk/LMDZ.MARS/libf/phymars/nlte_aux.F (modified) (5 diffs)

trunk/LMDZ.MARS/libf/phymars/nlte_aux.F

@@ -354 +354 @@
       factor =  (h-zy(k)) /  (zy(k+1)-zy(k))
       ap = dble( exp( log(py(k)) + log(py(k+1)/py(k)) * factor ) )
-      aco2 = dlog(con(k)) + dlog( con(k+1)/con(k) ) * dble(factor)
+      aco2 = log(con(k)) + log( con(k+1)/con(k) ) * dble(factor)
       aco2 = exp( aco2 )
       at = dble( ty(k) + (ty(k+1)-ty(k)) * factor )
@@ -392 +392 @@
       ap = dble( exp( log(py_cts(k)) +
      @     log(py_cts(k+1)/py_cts(k)) * factor ) )
-      aco2 = dlog(con(k)) + dlog( con(k+1)/con(k) ) * dble(factor)
+      aco2 = log(con(k)) + log( con(k+1)/con(k) ) * dble(factor)
       aco2 = exp( aco2 )
       at = dble( ty_cts(k) + (ty_cts(k+1)-ty_cts(k)) * factor )
@@ -1168 +1168 @@
             y3(i) = x3(j) + (x3(j+1)-x3(j)) * factor
          else
-            y1(i) = dexp( dlog(x1(j)) + dlog(x1(j+1)/x1(j)) * factor )
-            y2(i) = dexp( dlog(x2(j)) + dlog(x2(j+1)/x2(j)) * factor )
-            y3(i) = dexp( dlog(x3(j)) + dlog(x3(j+1)/x3(j)) * factor )
+            y1(i) = exp( log(x1(j)) + log(x1(j+1)/x1(j)) * factor )
+            y2(i) = exp( log(x2(j)) + log(x2(j+1)/x2(j)) * factor )
+            y3(i) = exp( log(x3(j)) + log(x3(j+1)/x3(j)) * factor )
          end if
 
@@ -1235 +1235 @@
             y4(i) = x4(j) + (x4(j+1)-x4(j)) * factor
          else
-            y1(i) = dexp( dlog(x1(j)) + dlog(x1(j+1)/x1(j)) * factor )
-            y2(i) = dexp( dlog(x2(j)) + dlog(x2(j+1)/x2(j)) * factor )
-            y3(i) = dexp( dlog(x3(j)) + dlog(x3(j+1)/x3(j)) * factor )
-            y4(i) = dexp( dlog(x4(j)) + dlog(x4(j+1)/x4(j)) * factor )
+            y1(i) = exp( log(x1(j)) + log(x1(j+1)/x1(j)) * factor )
+            y2(i) = exp( log(x2(j)) + log(x2(j+1)/x2(j)) * factor )
+            y3(i) = exp( log(x3(j)) + log(x3(j+1)/x3(j)) * factor )
+            y4(i) = exp( log(x4(j)) + log(x4(j+1)/x4(j)) * factor )
          end if
 
@@ -1300 +1300 @@
             y1(i) = x1(j) + (x1(j+1)-x1(j)) * factor 
          else
-            y1(i) = dexp( dlog(x1(j)) + dlog(x1(j+1)/x1(j)) * factor )
+            y1(i) = exp( log(x1(j)) + log(x1(j+1)/x1(j)) * factor )
          end if