Changeset 3702 for LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

Timestamp:

Jun 8, 2020, 6:20:04 PM (4 years ago)

Author:

oboucher

Message:

Shorten variable names and make polynom calculation faster for onshort/offshore wind capacity factors

File:

• LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -733 +733 @@
                zx_tmp_fi2d(i)=1.0
              ELSE
-               zx_tmp_fi2d(i)= 1.059e-09*x**10. - 1.351e-07*x**9. + 7.478e-06*x**8. - 0.0002352*x**7. + 0.004627*x**6. & 
-                             - 0.05898*x**5. + 0.4893*x**4. - 2.59*x**3. + 8.339*x**2. - 14.69*x + 10.73
+               zx_tmp_fi2d(i)= 10.73 + x*(-14.69 + x*(8.339 + x*(-2.59 + x*(0.4893 + x*(-0.05898 + x*(0.004627 + & 
+                               x*(-0.0002352 + x*(7.478e-06 + x*(-1.351e-07 + x*(1.059e-09))))))))))
                zx_tmp_fi2d(i)=MIN(MAX(zx_tmp_fi2d(i),0.0),1.0)
              ENDIF
@@ -754 +754 @@
                zx_tmp_fi2d(i)=1.0
              ELSE
-               zx_tmp_fi2d(i)= 3.352e-10*x**10. - 4.959e-08*x**9. + 3.195e-06*x**8. - 0.0001175*x**7. + 0.002716*x**6. &
-                             - 0.04099*x**5. + 0.4065*x**4. - 2.601*x**3. + 10.25*x**2. - 22.39*x + 20.59
+               zx_tmp_fi2d(i)= 20.59 + x*(-22.39 + x*(10.25 + x*(-2.601 + x*(0.4065 + x*(-0.04099 + x*(0.002716 + & 
+                               x*(-0.0001175 + x*(3.195e-06 + x*(-4.959e-08 + x*(3.352e-10))))))))))
                zx_tmp_fi2d(i)=MIN(MAX(zx_tmp_fi2d(i),0.0),1.0)
              ENDIF