Changeset 2151 for trunk/LMDZ.MARS/libf/phymars/nlte_calc.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_calc.F (modified) (11 diffs)

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

@@ -80 +80 @@
          tyd(i) = dble(ty_cts(i))
          con(i) =  dble( co2y_cts(i) * imr(isot) )
-         correc = 2.d0 * dexp( -ee*dble(elow(isot,2))/tyd(i) )
+         correc = 2.d0 * exp( -ee*dble(elow(isot,2))/tyd(i) )
          con(i) = con(i) * ( 1.d0 - correc )
          mr_cts(i) = dble(co2y_cts(i)/nty_cts(i))
@@ -201 +201 @@
 
          argumento = eqw / deltanudbl
-         tauinf(i) = dexp( - argumento )
+         tauinf(i) = exp( - argumento )
          if (i.eq.nl_cts_real) then
             taustar11_cts(i) = 0.0d0
@@ -282 +282 @@
       do i=1,nzy
          con(i) =  dble( co2y(i) * imr(isot) )
-         correc = 2.d0 * dexp( -ee*dble(elow(isot,2))/tvtbs(i) )
+         correc = 2.d0 * exp( -ee*dble(elow(isot,2))/tvtbs(i) )
          con(i) = con(i) * ( 1.d0 - correc )
          mr(i) = dble(co2y(i)/nty(i))
@@ -424 +424 @@
 
          argumento = eqw / deltanudbl
-         tauinf(i) = dexp( - argumento )
+         tauinf(i) = exp( - argumento )
 
       end do
@@ -978 +978 @@
 
          argumento = eqw / deltanudbl
-         tauinf(i) = dexp( - argumento )
+         tauinf(i) = exp( - argumento )
 
 
@@ -1056 +1056 @@
 
             argumento = eqw / deltanudbl
-            tau(in,ir) = dexp( - argumento )
+            tau(in,ir) = exp( - argumento )
 
  2       continue
@@ -1127 +1127 @@
 
          argumento = eqw / deltanudbl
-         tau(in,ir) = dexp( - argumento )
+         tau(in,ir) = exp( - argumento )
 
  4    continue
@@ -1375 +1375 @@
       do i=1,nzy
          con(i) =  dble( co2y(i) * imr(isot) )
-         correc = 2.d0 * dexp( -ee*dble(elow(isot,2))/tvtbs(i) )
+         correc = 2.d0 * exp( -ee*dble(elow(isot,2))/tvtbs(i) )
          con(i) = con(i) * ( 1.d0 - correc )
          mr(i) = dble(co2y(i)/nty(i))
@@ -1492 +1492 @@
             eqw = eqw + no(kr)*w
          end do
-         tauinf(i) = dexp( - eqw / deltanudbl )
+         tauinf(i) = exp( - eqw / deltanudbl )
          if (tauinf(i).lt.0.d0) tauinf(i) = 0.0d0
 
@@ -1645 +1645 @@
       do i=1,nzy
          con(i) =  dble( co2y(i) * imr(isot) )
-         correc = 2.d0 * dexp( -ee*dble(elow(isot,2))/tvtbs(i) )
+         correc = 2.d0 * exp( -ee*dble(elow(isot,2))/tvtbs(i) )
          con(i) = con(i) * ( 1.d0 - correc )
          mr(i) = dble(co2y(i)/nty(i))
@@ -1759 +1759 @@
 
          argumento = eqw / deltanudbl
-         tau(in,ir) = dexp( - argumento )
+         tau(in,ir) = exp( - argumento )
 
  2    continue