Changeset 3702


Ignore:
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:
3 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/DefLists/field_def_lmdz.xml

    r3691 r3702  
    1717      <field_group id="remap_1mo" />
    1818    </field_group>
    19    
    2019   
    2120    <field_group id="coordinates" grid_ref="grid_glo">
     
    6463        <field id="wind10max"  long_name="10m wind speed max"   unit="m/s" operation="maximum"/>
    6564        <field id="wind100m"   long_name="100-m wind speed"      unit="m/s" />
    66         <field id="loadfactor_wind_onshore"   long_name="Load factor for onshore windmill"      unit="-" />
    67         <field id="loadfactor_wind_offshore"  long_name="Load factor for offshore windmill"     unit="-" />
     65        <field id="woncfr"   long_name="Onshore Wind Capacity factor"      unit="kW/kW_installed" />
     66        <field id="wofcfr"   long_name="Offshore Wind Capacity factor"     unit="kW/kW_installed" />
    6867        <field id="sfcWindmax" field_ref="wind10m" operation="maximum" />
    6968        <field id="sicf"    long_name="Sea-ice fraction"        unit="-" />
  • LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

    r3691 r3702  
    280280    'wind100m', '100-m wind speed', 'm/s', (/ ('', i=1, 10) /))
    281281  TYPE(ctrl_out), SAVE :: o_loadfactor_wind_onshore = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
    282     'load_factor_wind_onshore', 'Load factor for onshore windmill', '-', (/ ('', i=1, 10) /))
     282    'woncfr', 'Onshore Wind Capacity factor', 'kW/kW_installed', (/ ('', i=1, 10) /))
    283283  TYPE(ctrl_out), SAVE :: o_loadfactor_wind_offshore = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
    284     'load_factor_wind_offshore', 'Load factor for offshore windmill', '-', (/ ('', i=1, 10) /))
     284    'wofcfr', 'Offshore Wind Capacity factor', 'kW/kW_installed', (/ ('', i=1, 10) /))
    285285  TYPE(ctrl_out), SAVE :: o_wind10max = ctrl_out((/ 10, 1, 10, 10, 10, 10, 11, 11, 11, 11/), &
    286286    'wind10max', '10m wind speed max', 'm/s', &
  • LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

    r3692 r3702  
    733733               zx_tmp_fi2d(i)=1.0
    734734             ELSE
    735                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. &
    736                              - 0.05898*x**5. + 0.4893*x**4. - 2.59*x**3. + 8.339*x**2. - 14.69*x + 10.73
     735               zx_tmp_fi2d(i)= 10.73 + x*(-14.69 + x*(8.339 + x*(-2.59 + x*(0.4893 + x*(-0.05898 + x*(0.004627 + &
     736                               x*(-0.0002352 + x*(7.478e-06 + x*(-1.351e-07 + x*(1.059e-09))))))))))
    737737               zx_tmp_fi2d(i)=MIN(MAX(zx_tmp_fi2d(i),0.0),1.0)
    738738             ENDIF
     
    754754               zx_tmp_fi2d(i)=1.0
    755755             ELSE
    756                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. &
    757                              - 0.04099*x**5. + 0.4065*x**4. - 2.601*x**3. + 10.25*x**2. - 22.39*x + 20.59
     756               zx_tmp_fi2d(i)= 20.59 + x*(-22.39 + x*(10.25 + x*(-2.601 + x*(0.4065 + x*(-0.04099 + x*(0.002716 + &
     757                               x*(-0.0001175 + x*(3.195e-06 + x*(-4.959e-08 + x*(3.352e-10))))))))))
    758758               zx_tmp_fi2d(i)=MIN(MAX(zx_tmp_fi2d(i),0.0),1.0)
    759759             ENDIF
Note: See TracChangeset for help on using the changeset viewer.