Changeset 2151 for trunk/LMDZ.MARS/libf/aeronomars

Timestamp:

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

Author:

aslmd

Message:

changed old functions dexp dlog in generic exp and log

File:

• trunk/LMDZ.MARS/libf/aeronomars/paramfoto_compact.F (modified) (16 diffs)

trunk/LMDZ.MARS/libf/aeronomars/paramfoto_compact.F

@@ -4707 +4707 @@
             
             if (o1d_eq(i).eq.'Y') then
-               log1 = dlog10(o1dxpares(1))
-               log2 = dlog10(o1dxpares(2))
-               log3 = dlog10(o1dxpares(3))
+               log1 = log10(o1dxpares(1))
+               log2 = log10(o1dxpares(2))
+               log3 = log10(o1dxpares(3))
                avg_pares = avg(log1,log2,log3)
                dif_pares = dif(log1,log2,log3,avg_pares)
-               log1 = dlog10(o1dximpar(1))
-               log2 = dlog10(o1dximpar(2))
-               log3 = dlog10(o1dximpar(3))
+               log1 = log10(o1dximpar(1))
+               log2 = log10(o1dximpar(2))
+               log3 = log10(o1dximpar(3))
                avg_impar = avg(log1,log2,log3)
                dif_impar = dif(log1,log2,log3,avg_impar)
@@ -4726 +4726 @@
             
             if (oh_eq(i).eq.'Y') then
-               log1 = dlog10(ohxpares(1))
-               log2 = dlog10(ohxpares(2))
-               log3 = dlog10(ohxpares(3))
+               log1 = log10(ohxpares(1))
+               log2 = log10(ohxpares(2))
+               log3 = log10(ohxpares(3))
                avg_pares = avg(log1,log2,log3)
                dif_pares = dif(log1,log2,log3,avg_pares)
-               log1 = dlog10(ohximpar(1))
-               log2 = dlog10(ohximpar(2))
-               log3 = dlog10(ohximpar(3))
+               log1 = log10(ohximpar(1))
+               log2 = log10(ohximpar(2))
+               log3 = log10(ohximpar(3))
                avg_impar = avg(log1,log2,log3)
                dif_impar = dif(log1,log2,log3,avg_impar)
@@ -4746 +4746 @@
             
             if (ho2_eq(i).eq.'Y') then
-               log1 = dlog10(ho2xpares(1))
-               log2 = dlog10(ho2xpares(2))
-               log3 = dlog10(ho2xpares(3))
+               log1 = log10(ho2xpares(1))
+               log2 = log10(ho2xpares(2))
+               log3 = log10(ho2xpares(3))
                avg_pares = avg(log1,log2,log3)
                dif_pares = dif(log1,log2,log3,avg_pares)
-               log1 = dlog10(ho2ximpar(1))
-               log2 = dlog10(ho2ximpar(2))
-               log3 = dlog10(ho2ximpar(3))
+               log1 = log10(ho2ximpar(1))
+               log2 = log10(ho2ximpar(2))
+               log3 = log10(ho2ximpar(3))
                avg_impar = avg(log1,log2,log3)
                dif_impar = dif(log1,log2,log3,avg_impar)
@@ -4765 +4765 @@
             
             if (h_eq(i).eq.'Y') then
-               log1 = dlog10(hxpares(1))
-               log2 = dlog10(hxpares(2))
-               log3 = dlog10(hxpares(3))
+               log1 = log10(hxpares(1))
+               log2 = log10(hxpares(2))
+               log3 = log10(hxpares(3))
                avg_pares = avg(log1,log2,log3)
                dif_pares = dif(log1,log2,log3,avg_pares)
-               log1 = dlog10(hximpar(1))
-               log2 = dlog10(hximpar(2))
-               log3 = dlog10(hximpar(3))
+               log1 = log10(hximpar(1))
+               log2 = log10(hximpar(2))
+               log3 = log10(hximpar(3))
                avg_impar = avg(log1,log2,log3)
                dif_impar = dif(log1,log2,log3,avg_impar)
@@ -4786 +4786 @@
             if(chemthermod.ge.2) then
                if (n2d_eq(i).eq.'Y') then
-                  log1 = dlog10(n2dxpares(1))
-                  log2 = dlog10(n2dxpares(2))
-                  log3 = dlog10(n2dxpares(3))
+                  log1 = log10(n2dxpares(1))
+                  log2 = log10(n2dxpares(2))
+                  log3 = log10(n2dxpares(3))
                   avg_pares = avg(log1,log2,log3)
                   dif_pares = dif(log1,log2,log3,avg_pares)
-                  log1 = dlog10(n2dximpar(1))
-                  log2 = dlog10(n2dximpar(2))
-                  log3 = dlog10(n2dximpar(3))
+                  log1 = log10(n2dximpar(1))
+                  log2 = log10(n2dximpar(2))
+                  log3 = log10(n2dximpar(3))
                   avg_impar = avg(log1,log2,log3)
                   dif_impar = dif(log1,log2,log3,avg_impar)
@@ -4805 +4805 @@
             
                if (no2_eq(i).eq.'Y') then
-                  log1 = dlog10(no2xpares(1))
-                  log2 = dlog10(no2xpares(2))
-                  log3 = dlog10(no2xpares(3))
+                  log1 = log10(no2xpares(1))
+                  log2 = log10(no2xpares(2))
+                  log3 = log10(no2xpares(3))
                   avg_pares = avg(log1,log2,log3)
                   dif_pares = dif(log1,log2,log3,avg_pares)
-                  log1 = dlog10(no2ximpar(1))
-                  log2 = dlog10(no2ximpar(2))
-                  log3 = dlog10(no2ximpar(3))
+                  log1 = log10(no2ximpar(1))
+                  log2 = log10(no2ximpar(2))
+                  log3 = log10(no2ximpar(3))
                   avg_impar = avg(log1,log2,log3)
                   dif_impar = dif(log1,log2,log3,avg_impar)
@@ -4829 +4829 @@
             if(chemthermod.eq.3) then
                if (cplus_eq(i).eq.'Y') then
-                  log1 = dlog10(cplusxpares(1))
-                  log2 = dlog10(cplusxpares(2))
-                  log3 = dlog10(cplusxpares(3))
+                  log1 = log10(cplusxpares(1))
+                  log2 = log10(cplusxpares(2))
+                  log3 = log10(cplusxpares(3))
                   avg_pares = avg(log1,log2,log3)
                   dif_pares = dif(log1,log2,log3,avg_pares)
-                  log1 = dlog10(cplusximpar(1))
-                  log2 = dlog10(cplusximpar(2))
-                  log3 = dlog10(cplusximpar(3))
+                  log1 = log10(cplusximpar(1))
+                  log2 = log10(cplusximpar(2))
+                  log3 = log10(cplusximpar(3))
                   avg_impar = avg(log1,log2,log3)
                   dif_impar = dif(log1,log2,log3,avg_impar)
@@ -4848 +4848 @@
                
                if (coplus_eq(i).eq.'Y') then
-                  log1 = dlog10(coplusxpares(1))
-                  log2 = dlog10(coplusxpares(2))
-                  log3 = dlog10(coplusxpares(3))
+                  log1 = log10(coplusxpares(1))
+                  log2 = log10(coplusxpares(2))
+                  log3 = log10(coplusxpares(3))
                   avg_pares = avg(log1,log2,log3)
                   dif_pares = dif(log1,log2,log3,avg_pares)
-                  log1 = dlog10(coplusximpar(1))
-                  log2 = dlog10(coplusximpar(2))
-                  log3 = dlog10(coplusximpar(3))
+                  log1 = log10(coplusximpar(1))
+                  log2 = log10(coplusximpar(2))
+                  log3 = log10(coplusximpar(3))
                   avg_impar = avg(log1,log2,log3)
                   dif_impar = dif(log1,log2,log3,avg_impar)
@@ -4867 +4867 @@
                
                if (oplus_eq(i).eq.'Y') then
-                  log1 = dlog10(oplusxpares(1))
-                  log2 = dlog10(oplusxpares(2))
-                  log3 = dlog10(oplusxpares(3))
+                  log1 = log10(oplusxpares(1))
+                  log2 = log10(oplusxpares(2))
+                  log3 = log10(oplusxpares(3))
                   avg_pares = avg(log1,log2,log3)
                   dif_pares = dif(log1,log2,log3,avg_pares)
-                  log1 = dlog10(oplusximpar(1))
-                  log2 = dlog10(oplusximpar(2))
-                  log3 = dlog10(oplusximpar(3))
+                  log1 = log10(oplusximpar(1))
+                  log2 = log10(oplusximpar(2))
+                  log3 = log10(oplusximpar(3))
                   avg_impar = avg(log1,log2,log3)
                   dif_impar = dif(log1,log2,log3,avg_impar)
@@ -4886 +4886 @@
                
                if (n2plus_eq(i).eq.'Y') then
-                  log1 = dlog10(n2plusxpares(1))
-                  log2 = dlog10(n2plusxpares(2))
-                  log3 = dlog10(n2plusxpares(3))
+                  log1 = log10(n2plusxpares(1))
+                  log2 = log10(n2plusxpares(2))
+                  log3 = log10(n2plusxpares(3))
                   avg_pares = avg(log1,log2,log3)
                   dif_pares = dif(log1,log2,log3,avg_pares)
-                  log1 = dlog10(n2plusximpar(1))
-                  log2 = dlog10(n2plusximpar(2))
-                  log3 = dlog10(n2plusximpar(3))
+                  log1 = log10(n2plusximpar(1))
+                  log2 = log10(n2plusximpar(2))
+                  log3 = log10(n2plusximpar(3))
                   avg_impar = avg(log1,log2,log3)
                   dif_impar = dif(log1,log2,log3,avg_impar)
@@ -4905 +4905 @@
                
                if (hplus_eq(i).eq.'Y') then
-                  log1 = dlog10(hplusxpares(1))
-                  log2 = dlog10(hplusxpares(2))
-                  log3 = dlog10(hplusxpares(3))
+                  log1 = log10(hplusxpares(1))
+                  log2 = log10(hplusxpares(2))
+                  log3 = log10(hplusxpares(3))
                   avg_pares = avg(log1,log2,log3)
                   dif_pares = dif(log1,log2,log3,avg_pares)
-                  log1 = dlog10(hplusximpar(1))
-                  log2 = dlog10(hplusximpar(2))
-                  log3 = dlog10(hplusximpar(3))
+                  log1 = log10(hplusximpar(1))
+                  log2 = log10(hplusximpar(2))
+                  log3 = log10(hplusximpar(3))
                   avg_impar = avg(log1,log2,log3)
                   dif_impar = dif(log1,log2,log3,avg_impar)
@@ -4924 +4924 @@
                
                if (co2plus_eq(i).eq.'Y') then
-                  log1 = dlog10(co2plusxpares(1))
-                  log2 = dlog10(co2plusxpares(2))
-                  log3 = dlog10(co2plusxpares(3))
+                  log1 = log10(co2plusxpares(1))
+                  log2 = log10(co2plusxpares(2))
+                  log3 = log10(co2plusxpares(3))
                   avg_pares = avg(log1,log2,log3)
                   dif_pares = dif(log1,log2,log3,avg_pares)
-                  log1 = dlog10(co2plusximpar(1))
-                  log2 = dlog10(co2plusximpar(2))
-                  log3 = dlog10(co2plusximpar(3))
+                  log1 = log10(co2plusximpar(1))
+                  log2 = log10(co2plusximpar(2))
+                  log3 = log10(co2plusximpar(3))
                   avg_impar = avg(log1,log2,log3)
                   dif_impar = dif(log1,log2,log3,avg_impar)
@@ -4943 +4943 @@
                
                if (o2plus_eq(i).eq.'Y') then
-                  log1 = dlog10(o2plusxpares(1))
-                  log2 = dlog10(o2plusxpares(2))
-                  log3 = dlog10(o2plusxpares(3))
+                  log1 = log10(o2plusxpares(1))
+                  log2 = log10(o2plusxpares(2))
+                  log3 = log10(o2plusxpares(3))
                   avg_pares = avg(log1,log2,log3)
                   dif_pares = dif(log1,log2,log3,avg_pares)
-                  log1 = dlog10(o2plusximpar(1))
-                  log2 = dlog10(o2plusximpar(2))
-                  log3 = dlog10(o2plusximpar(3))
+                  log1 = log10(o2plusximpar(1))
+                  log2 = log10(o2plusximpar(2))
+                  log3 = log10(o2plusximpar(3))
                   avg_impar = avg(log1,log2,log3)
                   dif_impar = dif(log1,log2,log3,avg_impar)
@@ -4962 +4962 @@
                
                if (noplus_eq(i).eq.'Y') then
-                  log1 = dlog10(noplusxpares(1))
-                  log2 = dlog10(noplusxpares(2))
-                  log3 = dlog10(noplusxpares(3))
+                  log1 = log10(noplusxpares(1))
+                  log2 = log10(noplusxpares(2))
+                  log3 = log10(noplusxpares(3))
                   avg_pares = avg(log1,log2,log3)
                   dif_pares = dif(log1,log2,log3,avg_pares)
-                  log1 = dlog10(noplusximpar(1))
-                  log2 = dlog10(noplusximpar(2))
-                  log3 = dlog10(noplusximpar(3))
+                  log1 = log10(noplusximpar(1))
+                  log2 = log10(noplusximpar(2))
+                  log3 = log10(noplusximpar(3))
                   avg_impar = avg(log1,log2,log3)
                   dif_impar = dif(log1,log2,log3,avg_impar)
@@ -4981 +4981 @@
                
                if (nplus_eq(i).eq.'Y') then
-                  log1 = dlog10(nplusxpares(1))
-                  log2 = dlog10(nplusxpares(2))
-                  log3 = dlog10(nplusxpares(3))
+                  log1 = log10(nplusxpares(1))
+                  log2 = log10(nplusxpares(2))
+                  log3 = log10(nplusxpares(3))
                   avg_pares = avg(log1,log2,log3)
                   dif_pares = dif(log1,log2,log3,avg_pares)
-                  log1 = dlog10(nplusximpar(1))
-                  log2 = dlog10(nplusximpar(2))
-                  log3 = dlog10(nplusximpar(3))
+                  log1 = log10(nplusximpar(1))
+                  log2 = log10(nplusximpar(2))
+                  log3 = log10(nplusximpar(3))
                   avg_impar = avg(log1,log2,log3)
                   dif_impar = dif(log1,log2,log3,avg_impar)
@@ -5000 +5000 @@
 
                if (hco2plus_eq(i).eq.'Y') then
-                 log1 = dlog10(hco2plusxpares(1))
-                 log2 = dlog10(hco2plusxpares(2))
-                 log3 = dlog10(hco2plusxpares(3))
+                 log1 = log10(hco2plusxpares(1))
+                 log2 = log10(hco2plusxpares(2))
+                 log3 = log10(hco2plusxpares(3))
                  avg_pares = avg(log1,log2,log3)
                  dif_pares = dif(log1,log2,log3,avg_pares)
-                 log1 = dlog10(hco2plusximpar(1))
-                 log2 = dlog10(hco2plusximpar(2))
-                 log3 = dlog10(hco2plusximpar(3))
+                 log1 = log10(hco2plusximpar(1))
+                 log2 = log10(hco2plusximpar(2))
+                 log3 = log10(hco2plusximpar(3))
                  avg_impar = avg(log1,log2,log3)
                  dif_impar = dif(log1,log2,log3,avg_impar)