Changeset 3702 for LMDZ6

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

Location:

LMDZ6/trunk

Files:

• DefLists/field_def_lmdz.xml (modified) (2 diffs)
• libf/phylmd/phys_output_ctrlout_mod.F90 (modified) (1 diff)
• libf/phylmd/phys_output_write_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/DefLists/field_def_lmdz.xml

@@ -17 +17 @@
       <field_group id="remap_1mo" />
     </field_group>
-    
     
     <field_group id="coordinates" grid_ref="grid_glo">
@@ -64 +63 @@
         <field id="wind10max"  long_name="10m wind speed max"   unit="m/s" operation="maximum"/>
         <field id="wind100m"   long_name="100-m wind speed"      unit="m/s" />
-        <field id="loadfactor_wind_onshore"   long_name="Load factor for onshore windmill"      unit="-" />
-        <field id="loadfactor_wind_offshore"  long_name="Load factor for offshore windmill"     unit="-" />
+        <field id="woncfr"   long_name="Onshore Wind Capacity factor"      unit="kW/kW_installed" />
+        <field id="wofcfr"   long_name="Offshore Wind Capacity factor"     unit="kW/kW_installed" />
         <field id="sfcWindmax" field_ref="wind10m" operation="maximum" />
         <field id="sicf"    long_name="Sea-ice fraction"        unit="-" />

LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -280 +280 @@
     'wind100m', '100-m wind speed', 'm/s', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_loadfactor_wind_onshore = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
-    'load_factor_wind_onshore', 'Load factor for onshore windmill', '-', (/ ('', i=1, 10) /))
+    'woncfr', 'Onshore Wind Capacity factor', 'kW/kW_installed', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_loadfactor_wind_offshore = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
-    'load_factor_wind_offshore', 'Load factor for offshore windmill', '-', (/ ('', i=1, 10) /))
+    'wofcfr', 'Offshore Wind Capacity factor', 'kW/kW_installed', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_wind10max = ctrl_out((/ 10, 1, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'wind10max', '10m wind speed max', 'm/s', &

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