Ignore:
Timestamp:
Jul 21, 2020, 11:08:53 AM (4 years ago)
Author:
adurocher
Message:

Offload post-processing to XIOS

use expressions in .xml files to offload computations previously done in phys_output_write to XIOS

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/Optimisation_LMDZ/DefLists/field_def_lmdz.xml

    r3656 r3768  
    66<!--        <field id="Bp"   long_name="Hybrid coordinate: sigma part"    unit="-" />  -->
    77<!--        <field id="Alt"  long_name="Pseudo-altitude (for an 8km scale height)"    unit="km" />  -->
    8 <!--    </field_group>  -->
    9 
     8<!--    </field_group>  -->   
    109    <field_group id="remap_expr" level="_reset_" prec="_reset_" operation="_reset_" freq_op="_reset_" enabled="_reset_" default_value="_reset_">
    1110      <field_group id="remap_1ts" />
     
    3837        <field id="phis"    long_name="Surface geop.height"     unit="m2/s2" />
    3938        <field id="aire"    long_name="Grid area"               unit="m2"    />
    40         <field id="contfracATM"     long_name="% sfce ter+lic"  unit="-"   />
    41         <field id="contfracOR"      long_name="% sfce terre OR" unit="-" />
     39        <field id="contfracATM"     long_name="% sfce ter+lic"  unit="-"   > fract_ter + fract_lic </field>
     40        <field id="contfracOR"      long_name="% sfce terre OR" unit="-"   > fract_ter </field>
    4241        <field id="sza"     long_name="Solar zenithal angle "   unit="deg" />
    4342        <field id="flat"    long_name="Latent heat flux"        unit="W/m2" />
     
    4948        <field id="tasmax"  field_ref="t2m"  operation="maximum" />
    5049        <field id="tasmin"  field_ref="t2m"  operation="minimum" />
    51         <field id="t2m_min"    long_name="Temp 2m min"    unit="K" operation="minimum" />
    52         <field id="t2m_max"    long_name="Temp 2m max"    unit="K" operation="maximum" />
     50        <field id="t2m_min"  field_ref="t2m" long_name="Temp 2m min" operation="minimum" />
     51        <field id="t2m_max"  field_ref="t2m" long_name="Temp 2m max" operation="maximum" />
    5352        <field id="t2m_ter"    long_name="Temp 2m ter"    unit="K" />
    5453        <field id="t2m_max_mon"   long_name="Monthly average max 2m temperature" unit="K" />
     
    6160        <field id="t2m_probsup30"   field_ref="t2m" long_name="Prob. t2m exceeds 30 degC" unit="-"> t2m &gt; 303.15 </field>
    6261        <field id="gusts"    long_name="surface gustiness"    unit="m2/s2" />
    63         <field id="wind10m"    long_name="10-m wind speed"      unit="m/s" />
    64         <field id="wind10max"  long_name="10m wind speed max"   unit="m/s" operation="maximum"/>
     62        <field id="wind10m"    long_name="10-m wind speed"      unit="m/s" >  sqrt(v10m*v10m + u10m*u10m) </field>
     63        <field id="wind10max"  long_name="10m wind speed max"   unit="m/s" operation="maximum"> wind10m </field>
    6564        <field id="wind100m"   long_name="100-m wind speed"      unit="m/s" />
    6665        <field id="sfcWindmax" field_ref="wind10m" operation="maximum" />
    67         <field id="sicf"    long_name="Sea-ice fraction"        unit="-" />
     66        <field id="sicf"    long_name="Sea-ice fraction"        unit="-" > fract_sic </field>
    6867        <field id="q2m"     long_name="Specific humidity 2m"    unit="kg/kg" />
    6968        <field id="ustar"   long_name="Friction velocity"       unit="m/s" />
     
    9190        <field id="qsol"        long_name="Soil watter content"         unit="mm" />
    9291        <field id="JrNt"        long_name="1 if Day 0 if Night"         unit="-" />
    93         <field id="cldhjn"      long_name="High-level cloudiness Day"         unit="-" />
    94         <field id="cldmjn"      long_name="Mid-level cloudiness Day"         unit="-" />
    95         <field id="cldljn"      long_name="Low-level cloudiness Day"         unit="-" />
    96         <field id="cldtjn"      long_name="Total cloudiness Day"         unit="-" />
     92        <field id="cldhjn"      long_name="High-level cloudiness Day"         unit="-" > cldh/JrNt </field>
     93        <field id="cldmjn"      long_name="Mid-level cloudiness Day"         unit="-" > cldm/JrNt </field>
     94        <field id="cldljn"      long_name="Low-level cloudiness Day"         unit="-" > cldl/JrNt </field>
     95        <field id="cldtjn"      long_name="Total cloudiness Day"         unit="-"  > cldt/JrNt </field>
    9796        <field id="rnebjn"      long_name="Cloud fraction in day"         unit="-" />
    9897        <field id="ndayrain"    long_name="Number of dayrain(liq+sol)"  unit="-" operation="instant"/>
    99         <field id="precip"      long_name="Precip Totale liq+sol"       unit="kg/(s*m2)" />
     98        <field id="precip"      long_name="Precip Totale liq+sol"       unit="kg/(s*m2)" > rain_fall+snow </field>
    10099        <field id="prhmax"      field_ref="precip" operation="maximum"           />
    101100        <field id="rain_fall"   long_name="Precip Totale liq"   unit="kg/(s*m2)" />
    102101        <field id="rain_con"    long_name="Precip convective liq"   unit="kg/(s*m2)" />
    103         <field id="plul"    long_name="Large-scale Precip."     unit="kg/(s*m2)" />
     102        <field id="plul"    long_name="Large-scale Precip."     unit="kg/(s*m2)" > this+snowl </field>
    104103        <field id="plun"    long_name="Numerical Precip."       unit="kg/(s*m2)" />
    105         <field id="pluc"    long_name="Convective Precip."      unit="kg/(s*m2)" />
     104        <field id="pluc"    long_name="Convective Precip."      unit="kg/(s*m2)" > rain_con + this </field>
    106105        <field id="snow"    long_name="Snow fall"               unit="kg/(s*m2)" />
    107106        <field id="evap"    long_name="Evaporat"                unit="kg/(s*m2)" />
     
    112111        <field id="msnow"    long_name="Surface snow amount"    unit="kg/m2" />
    113112        <field id="fsnow"    long_name="Surface snow area fraction"    unit="-" />
    114         <field id="tops"    long_name="Solar rad. at TOA"    unit="W/m2" />
    115         <field id="tops0"    long_name="CS Solar rad. at TOA"    unit="W/m2" />
     113        <field id="swradcorr"    long_name="Temporary" />
     114        <field id="tops"    long_name="Solar rad. at TOA"    unit="W/m2" > this*swradcorr </field>
     115        <field id="tops0"    long_name="CS Solar rad. at TOA"    unit="W/m2" > this*swradcorr </field>
    116116        <field id="topl"    long_name="IR rad. at TOA"    unit="W/m2" />
    117117        <field id="topl0"    long_name="IR rad. at TOA"    unit="W/m2" />
    118         <field id="SWupTOA"    long_name="SWup at TOA"    unit="W/m2" />
    119         <field id="SWupTOAclr"    long_name="SWup clear sky at TOA"    unit="W/m2" />
    120         <field id="SWupTOAcleanclr"    long_name="SWup clean (no aerosol) clear sky at TOA"    unit="W/m2" />
    121         <field id="SWdnTOA"    long_name="SWdn at TOA"    unit="W/m2" />
    122         <field id="SWdnTOAclr"    long_name="SWdn clear sky at TOA"    unit="W/m2" />
    123         <field id="nettop"    long_name="Net dn radiatif flux at TOA"    unit="W/m2" />
    124         <field id="SWup200"    long_name="SWup at 200mb"    unit="W/m2" />
    125         <field id="SWup200clr"    long_name="SWup clear sky at 200mb"    unit="W/m2" />
    126         <field id="SWdn200"    long_name="SWdn at 200mb"    unit="W/m2" />
    127         <field id="SWdn200clr"    long_name="SWdn clear sky at 200mb"    unit="W/m2" />
     118        <field id="SWupTOA"    long_name="SWup at TOA"    unit="W/m2" > this*swradcorr </field>
     119        <field id="SWupTOAclr"    long_name="SWup clear sky at TOA"    unit="W/m2"  > this*swradcorr </field>
     120        <field id="SWupTOAcleanclr"    long_name="SWup clean (no aerosol) clear sky at TOA"    unit="W/m2" > this*swradcorr </field>
     121        <field id="SWdnTOA"    long_name="SWdn at TOA"    unit="W/m2" > this*swradcorr </field>
     122        <field id="SWdnTOAclr"    long_name="SWdn clear sky at TOA"    unit="W/m2"> this*swradcorr </field>
     123        <field id="nettop"    long_name="Net dn radiatif flux at TOA"    unit="W/m2" > tops - topl </field>
     124        <field id="SWup200"    long_name="SWup at 200mb"    unit="W/m2" > this*swradcorr </field>
     125        <field id="SWup200clr"    long_name="SWup clear sky at 200mb"    unit="W/m2" > this*swradcorr </field>
     126        <field id="SWdn200"    long_name="SWdn at 200mb"    unit="W/m2" > this*swradcorr </field>
     127        <field id="SWdn200clr"    long_name="SWdn clear sky at 200mb"    unit="W/m2" > this*swradcorr </field>
    128128        <field id="LWup200"    long_name="LWup at 200mb"    unit="W/m2" />
    129129        <field id="LWup200clr"    long_name="LWup clear sky at 200mb"    unit="W/m2" />
    130130        <field id="LWdn200"    long_name="LWdn at 200mb"    unit="W/m2" />
    131131        <field id="LWdn200clr"    long_name="LWdn clear sky at 200mb"    unit="W/m2" />
    132         <field id="sols"    long_name="Solar rad. at surf."    unit="W/m2" />
    133         <field id="sols0"    long_name="Solar rad. at surf."    unit="W/m2" />
     132        <field id="sols"    long_name="Solar rad. at surf."    unit="W/m2" > this*swradcorr </field>
     133        <field id="sols0"    long_name="Solar rad. at surf."    unit="W/m2" > this*swradcorr </field>
    134134        <field id="soll"    long_name="IR rad. at surface"    unit="W/m2" />
    135135        <field id="soll0"    long_name="IR rad. at surface"    unit="W/m2" />
    136136        <field id="radsol"    long_name="Rayonnement au sol"    unit="W/m2" />
    137         <field id="SWupSFC"    long_name="SWup at surface"    unit="W/m2" />
    138         <field id="SWupSFCclr"    long_name="SWup clear sky at surface"    unit="W/m2" />
    139         <field id="SWupSFCcleanclr"    long_name="SWup clean (no aerosol) clear sky at surface"    unit="W/m2" />
    140         <field id="SWdnSFC"    long_name="SWdn at surface"    unit="W/m2" />
    141         <field id="SWdnSFCclr"    long_name="SWdn clear sky at surface"    unit="W/m2" />
    142         <field id="SWdnSFCcleanclr"    long_name="SWdn clean (no aerosol) clear sky at surface"    unit="W/m2" />
    143         <field id="LWupSFC"    long_name="Upwd. IR rad. at surface"    unit="W/m2" />
    144         <field id="LWupSFCclr"    long_name="CS Upwd. IR rad. at surface"    unit="W/m2" />
     137 
     138        <field id="SWupSFC"    long_name="SWup at surface"    unit="W/m2" > this*swradcorr </field>
     139        <field id="SWupSFCclr"    long_name="SWup clear sky at surface"    unit="W/m2" > this*swradcorr </field>
     140        <field id="SWupSFCcleanclr"    long_name="SWup clean (no aerosol) clear sky at surface"    unit="W/m2" > this*swradcorr </field>
     141        <field id="SWdnSFC"    long_name="SWdn at surface"    unit="W/m2" > this*swradcorr </field>
     142        <field id="SWdnSFCclr"    long_name="SWdn clear sky at surface"    unit="W/m2" > this*swradcorr </field>
     143        <field id="SWdnSFCcleanclr"    long_name="SWdn clean (no aerosol) clear sky at surface"    unit="W/m2" > this*swradcorr </field>
     144        <field id="LWupSFC"    long_name="Upwd. IR rad. at surface"    unit="W/m2"> LWdnSFC - soll </field>
     145       
     146        <field id="LWupSFCclr"    long_name="CS Upwd. IR rad. at surface"    unit="W/m2" > LWdnSFCclr - soll0  </field>
    145147        <field id="LWupTOAcleanclr"    long_name="CS clean (no aerosol) Upwd. IR rad. at TOA"    unit="W/m2" />
    146         <field id="LWdnSFC"    long_name="Down. IR rad. at surface"    unit="W/m2" />
    147         <field id="LWdnSFCclr"    long_name="Down. CS IR rad. at surface"    unit="W/m2" />
    148         <field id="LWdnSFCcleanclr"    long_name="Down. CS clean (no aerosol) IR rad. at surface"    unit="W/m2" />
     148        <field id="LWdnSFC"    long_name="Down. IR rad. at surface"    unit="W/m2"  />
     149        <field id="LWdnSFCclr"    long_name="Down. CS IR rad. at surface"    unit="W/m2" > -this </field>
     150        <field id="LWdnSFCcleanclr"    long_name="Down. CS clean (no aerosol) IR rad. at surface"    unit="W/m2" > -this </field>
    149151        <field id="bils"    long_name="Surf. total heat flux"    unit="W/m2" />
    150152        <field id="bils_tke"    long_name="Surf. total heat flux"    unit="W/m2" />
     
    161163        <field id="ustr_gwd_rando"    long_name="zonal wind stress random gravity waves"    unit="Pa" />
    162164        <field id="vstr_gwd_rando"    long_name="meridinal wind stress random gravity waves"    unit="Pa" />
    163         <field id="sens"    long_name="Sensible heat flux"    unit="W/m2" />
     165        <field id="sens"    long_name="Sensible heat flux"    unit="W/m2" > -this </field>
    164166        <field id="sens_x"    long_name="sens off-wake"    unit="W/m2" />
    165167        <field id="sens_w"    long_name="sens within-wake"    unit="W/m2" />
    166168        <field id="flat_x"    long_name="flat off-wake"    unit="W/m2" />
    167169        <field id="flat_w"    long_name="flat within-wake"    unit="W/m2" />
    168         <field id="delta_tsurf"    long_name="Temperature difference (w-x)"    unit="K" enabled="false" /> <!-- ym --field disabled-- delta_tsurf is not conform to the associated grid (the is an arry containing all sub_surface) -->
     170        <field id="delta_tsurf"    long_name="Temperature difference (w-x)"    unit="K" enabled="false" /> <!-- ym field disabled delta_tsurf is not conform to the associated grid (the is an arry containing all sub_surface) -->
    169171        <field id="cdragh_x"    long_name="cdragh off-wake"    unit="-" />
    170172        <field id="cdragh_w"    long_name="cdragh within-wake"    unit="-" />
     
    179181        <field id="fqfonte"    long_name="Land ice melt"    unit="kg/m2/s" />
    180182        <field id="runofflic"  long_name="Land ice melt to ocean" unit="kg/m2/s" />
    181         <field id="mrroli"  long_name="Runoff flux over land ice" unit="kg/m2/s" />
    182         <field id="taux"    long_name="Zonal wind stress"    unit="Pa" />
    183         <field id="tauy"    long_name="Meridional wind stress"    unit="Pa" />
     183        <field id="mrroli"  long_name="Runoff flux over land ice" unit="kg/m2/s"> (fqfonte + rain_fall)*fract_lic </field>
     184        <field id="taux"    long_name="Zonal wind stress"    unit="Pa" >      fract_ter*taux_ter + fract_lic*taux_lic + fract_oce*taux_oce +fract_sic*taux_sic </field>
     185        <field id="tauy"    long_name="Meridional wind stress"    unit="Pa" > fract_ter*tauy_ter + fract_lic*tauy_lic + fract_oce*tauy_oce +fract_sic*tauy_sic </field>
    184186        <field id="taux_ter"    long_name="Zonal wind stress ter"    unit="Pa" />
    185187        <field id="taux_lic"    long_name="Zonal wind stress lic"    unit="Pa" />
     
    190192        <field id="tauy_oce"    long_name="Meridional wind stress oce"    unit="Pa" />
    191193        <field id="tauy_sic"    long_name="Meridional wind stress sic"    unit="Pa" />
    192         <field id="pourc_ter"    long_name="% ter"    unit="%" />
    193         <field id="pourc_lic"    long_name="% lic"    unit="%" />
    194         <field id="pourc_oce"    long_name="% oce"    unit="%" />
    195         <field id="pourc_sic"    long_name="% sic"    unit="%" />
     194        <field id="pourc_ter"    long_name="% ter"    unit="%" > 100*fract_ter </field>
     195        <field id="pourc_lic"    long_name="% lic"    unit="%" > 100*fract_lic </field>
     196        <field id="pourc_oce"    long_name="% oce"    unit="%" > 100*fract_oce </field>
     197        <field id="pourc_sic"    long_name="% sic"    unit="%" > 100*fract_sic </field>
    196198        <field id="fract_ter"    long_name="Fraction ter"    unit="1" />
    197199        <field id="fract_lic"    long_name="Fraction lic"    unit="1" />
     
    335337        <field id="z50"    long_name="Geopotential height 50hPa"    unit="m" detect_missing_value=".true." />
    336338        <field id="z10"    long_name="Geopotential height 10hPa"    unit="m" detect_missing_value=".true." />
    337         <field id="t_oce_sic"    long_name="Temp mixte oce-sic"    unit="K" />
     339        <field id="t_oce_sic"    long_name="Temp mixte oce-sic"    unit="K" > ((fract_oce > 1.0E-05)+(fract_sic > 1.0E-05))?(tsol_oce*fract_oce+tsol_sic*fract_sic)/(fract_oce+fract_sic):273.15 </field>
     340        <!--<field id="t_oce_sic"    long_name="Temp mixte oce-sic"    unit="K" > (fract_oce > 1.0E-05 || fract_sic > 1.0E-05)?(tsol_oce*fract_oce+tsol_sic*fract_sic)/(fract_oce+fract_sic):273.15 </field > -->
    338341        <field id="weakinv"    long_name="Weak inversion"    unit="-" />
    339342        <field id="dthmin"    long_name="dTheta mini"    unit="K/m" />
     
    346349        <field id="v10_oce"    long_name="-"    unit="-" />
    347350        <field id="v10_sic"    long_name="-"    unit="-" />
    348         <field id="rh2m"    long_name="Relative humidity at 2m"    unit="%" />
     351        <field id="rh2m"    long_name="Relative humidity at 2m"    unit="%" > (this > 1.) ? 100. : this*100. </field>
    349352        <field id="rh2m_max"    field_ref="rh2m" operation="maximum" enabled="false"/> <!-- ym Problem : field is disabled check rh2m_min has a reference (field_ref) to rh2 but is also sent by model which is not compatible-->
    350353        <field id="rh2m_min"    field_ref="rh2m" operation="minimum" enabled="false"/> <!-- ym Problem : field is disabled check rh2m_max has a reference (field_ref) to rh2 but is also sent by model which is not compatible-->
     
    354357        <field id="tpot"    long_name="Surface air potential temperature"    unit="K" />
    355358        <field id="tpote"    long_name="Surface air equivalent potential temperature"    unit="K" />
    356         <field id="SWnetOR"    long_name="Sfce net SW radiation OR"    unit="W/m2" />
     359        <field id="SWnetOR"   field_ref="fsw_ter" long_name="Sfce net SW radiation OR"    unit="W/m2" />
    357360        <field id="SWdownOR"    long_name="Sfce incident SW radiation OR"    unit="W/m2" />
    358         <field id="LWdownOR"    long_name="Sfce incident LW radiation OR"    unit="W/m2" />
     361        <field id="LWdownOR"   field_ref="LWdnSFC" long_name="Sfce incident LW radiation OR"    unit="W/m2" />
    359362        <field id="snowl"    long_name="Solid Large-scale Precip."    unit="kg/(m2*s)" />
    360         <field id="cape_max"    long_name="CAPE max."    unit="J/kg" operation="maximum"/>
    361         <field id="solldown"    long_name="Down. IR rad. at surface"    unit="W/m2" />
     363        <field id="cape_max"  field_ref="cape"   long_name="CAPE max." operation="maximum"/>
     364        <field id="solldown"   field_ref="LWdnSFC" long_name="Down. IR rad. at surface"    unit="W/m2" />
    362365        <field id="dtsvdfo"    long_name="Boundary-layer dTs(o)"    unit="K/s" />
    363366        <field id="dtsvdft"    long_name="Boundary-layer dTs(t)"    unit="K/s" />
     
    369372        <field id="topswad"    long_name="ADE at TOA"    unit="W/m2" />
    370373        <field id="topswad0"    long_name="ADE clear-sky at TOA"    unit="W/m2" />
    371         <field id="topswai"    long_name="AIE at TOA"    unit="W/m2" />
    372         <field id="solswad"    long_name="ADE at SRF"    unit="W/m2" />
     374        <field id="topswai"    long_name="AIE at TOA"    unit="W/m2" > this/swradcorr </field>
     375        <field id="solswad"    long_name="ADE at SRF"    unit="W/m2" > this/swradcorr </field>
    373376        <field id="solswad0"    long_name="ADE clear-sky at SRF"    unit="W/m2" />
    374         <field id="solswai"    long_name="AIE at SFR"    unit="W/m2" />
    375         <field id="toplwad"    long_name="LW-ADE at TOA" unit="W/m2" />
    376         <field id="toplwad0"    long_name="LW-ADE clear-sky at TOA"    unit="W/m2" />
    377         <field id="toplwai"    long_name="LW-AIE clear-sky at TOA"    unit="W/m2" />
     377        <field id="solswai"    long_name="AIE at SFR"    unit="W/m2" > this/swradcorr </field>
     378        <field id="toplwad"    long_name="LW-ADE at TOA" unit="W/m2"  > this/swradcorr </field>
     379        <field id="toplwad0"    long_name="LW-ADE clear-sky at TOA"    unit="W/m2"  > this/swradcorr </field>
     380        <field id="toplwai"    long_name="LW-AIE clear-sky at TOA"    unit="W/m2" > this/swradcorr </field>
    378381        <field id="sollwad"    long_name="LW-ADE at SRF" unit="W/m2" />
    379382        <field id="sollwad0"    long_name="LW-ADE clear-sky at SRF"    unit="W/m2" />
    380         <field id="sollwai"    long_name="LW-AIE at SRF" unit="W/m2" />
     383        <field id="sollwai"    long_name="LW-AIE at SRF" unit="W/m2" > this/swradcorr </field>
    381384        <field id="od550_ASBCM"    long_name="Aerosol Optical depth at 550 nm "    unit="1" />
    382385        <field id="od550_ASPOMM"   long_name="Aerosol Optical depth at 550 nm "    unit="1" />
     
    425428        <field id="loaddust"  long_name="Column Load of Dust "    unit="kg/m2" />
    426429        <field id="loadno3"   long_name="Column Load of Nitrate " unit="kg/m2" />
    427         <field id="swtoaas_nat"    long_name="Natural aerosol radiative forcing all-sky at TOA"    unit="W/m2" />
    428         <field id="swsrfas_nat"    long_name="Natural aerosol radiative forcing all-sky at SRF"    unit="W/m2" />
    429         <field id="swtoacs_nat"    long_name="Natural aerosol radiative forcing clear-sky at TOA"    unit="W/m2" />
    430         <field id="swsrfcs_nat"    long_name="Natural aerosol radiative forcing clear-sky at SRF"    unit="W/m2" />
    431         <field id="swtoaas_ant"    long_name="Anthropogenic aerosol radiative forcing all-sky at TOA"    unit="W/m2" />
    432         <field id="swsrfas_ant"    long_name="Anthropogenic aerosol radiative forcing all-sky at SRF"    unit="W/m2" />
    433         <field id="swtoacs_ant"    long_name="Anthropogenic aerosol radiative forcing clear-sky at TOA"    unit="W/m2" />
    434         <field id="swsrfcs_ant"    long_name="Anthropogenic aerosol radiative forcing clear-sky at SRF"    unit="W/m2" />
    435         <field id="swtoacf_nat"    long_name="Natural aerosol impact on cloud radiative forcing at TOA"    unit="W/m2" />
    436         <field id="swsrfcf_nat"    long_name="Natural aerosol impact on cloud radiative forcing  at SRF"    unit="W/m2" />
    437         <field id="swtoacf_ant"    long_name="Anthropogenic aerosol impact on cloud radiative forcing at TOA"    unit="W/m2" />
    438         <field id="swsrfcf_ant"    long_name="Anthropogenic aerosol impact on cloud radiative forcing at SRF"    unit="W/m2" />
    439         <field id="swtoacf_zero"    long_name="Cloud radiative forcing (allsky-clearsky fluxes) at TOA"    unit="W/m2" />
    440         <field id="swsrfcf_zero"    long_name="Cloud radiative forcing (allsky-clearsky fluxes) at SRF"    unit="W/m2" />
     430        <field id="swtoaas_nat"    long_name="Natural aerosol radiative forcing all-sky at TOA"    unit="W/m2" > this/swradcorr </field>
     431        <field id="swsrfas_nat"    long_name="Natural aerosol radiative forcing all-sky at SRF"    unit="W/m2" > this/swradcorr </field>
     432        <field id="swtoacs_nat"    long_name="Natural aerosol radiative forcing clear-sky at TOA"    unit="W/m2" > this/swradcorr </field>
     433        <field id="swsrfcs_nat"    long_name="Natural aerosol radiative forcing clear-sky at SRF"    unit="W/m2" > this/swradcorr </field>
     434        <field id="swtoaas_ant"    long_name="Anthropogenic aerosol radiative forcing all-sky at TOA"    unit="W/m2" > this/swradcorr </field>
     435        <field id="swsrfas_ant"    long_name="Anthropogenic aerosol radiative forcing all-sky at SRF"    unit="W/m2"  > this/swradcorr </field>
     436        <field id="swtoacs_ant"    long_name="Anthropogenic aerosol radiative forcing clear-sky at TOA"    unit="W/m2"  > this/swradcorr </field>
     437        <field id="swsrfcs_ant"    long_name="Anthropogenic aerosol radiative forcing clear-sky at SRF"    unit="W/m2" > this/swradcorr </field>
     438        <field id="swtoacf_nat"    long_name="Natural aerosol impact on cloud radiative forcing at TOA"    unit="W/m2" > this/swradcorr </field>
     439        <field id="swsrfcf_nat"    long_name="Natural aerosol impact on cloud radiative forcing  at SRF"    unit="W/m2" > this/swradcorr </field>
     440        <field id="swtoacf_ant"    long_name="Anthropogenic aerosol impact on cloud radiative forcing at TOA"    unit="W/m2" > this/swradcorr </field>
     441        <field id="swsrfcf_ant"    long_name="Anthropogenic aerosol impact on cloud radiative forcing at SRF"    unit="W/m2" > this/swradcorr </field>
     442        <field id="swtoacf_zero"    long_name="Cloud radiative forcing (allsky-clearsky fluxes) at TOA"    unit="W/m2" > this/swradcorr </field>
     443        <field id="swsrfcf_zero"    long_name="Cloud radiative forcing (allsky-clearsky fluxes) at SRF"    unit="W/m2" > this/swradcorr </field>
    441444        <field id="cldncl"    long_name="CDNC at top of liquid water cloud"    unit="m-3" />
    442445        <field id="reffclwtop"    long_name="Droplet effective radius at top of liquid water cloud"    unit="m" />
     
    478481        <field id="alb1"    long_name="Surface VIS albedo"    unit="-" />
    479482        <field id="alb2"    long_name="Surface Near IR albedo"    unit="-" />
    480         <field id="ftime_con"    long_name="Fraction of time convection Occurs"    unit="-" operation="instant" />
     483        <field id="ftime_con"    long_name="Fraction of time convection Occurs"    unit="-" operation="instant" > (ftime_th>ftime_deepcv)?ftime_th:ftime_deepcv </field>
    481484        <field id="ftime_deepcv"    long_name="Fraction of time deep convection Occurs"    unit="-" />
    482485        <field id="wake_h"    long_name="wake height"    unit="m" />
     
    490493        <field id="f0_th"    long_name="Thermal closure mass flux"    unit="K/s" />
    491494        <field id="zmax_th"    long_name="Thermal plume height"    unit="K/s" />
    492         <field id="rsut4co2"    long_name="TOA Out SW in 4xCO2 atmosphere"    unit="W/m2" />
     495        <field id="rsut4co2"    long_name="TOA Out SW in 4xCO2 atmosphere"    unit="W/m2" > this*swradcorr </field>
    493496        <field id="rlut4co2"    long_name="TOA Out LW in 4xCO2 atmosphere"    unit="W/m2" />
    494         <field id="rsutcs4co2"    long_name="TOA Out CS SW in 4xCO2 atmosphere"    unit="W/m2" />
     497        <field id="rsutcs4co2"    long_name="TOA Out CS SW in 4xCO2 atmosphere"    unit="W/m2" > this*swradcorr </field>
    495498        <field id="rlutcs4co2"    long_name="TOA Out CS LW in 4xCO2 atmosphere"    unit="W/m2" />
    496499        <field id="dqphy2d"    long_name="Physics dQ"    unit="(kg/m2)/s" />
     
    527530        <field id="tke_oce"    long_name="Max Turb. Kinetic Energy oce"    unit="m2/s2" />
    528531        <field id="tke_sic"    long_name="Max Turb. Kinetic Energy sic"    unit="m2/s2" />
    529         <field id="tke_max_ter"    long_name="Max Turb. Kinetic Energy ter"    unit="m2/s2" operation="maximum"/>
    530         <field id="tke_max_lic"    long_name="Max Turb. Kinetic Energy lic"    unit="m2/s2" operation="maximum"/>
    531         <field id="tke_max_oce"    long_name="Max Turb. Kinetic Energy oce"    unit="m2/s2" operation="maximum"/>
    532         <field id="tke_max_sic"    long_name="Max Turb. Kinetic Energy sic"    unit="m2/s2" operation="maximum"/>
     532        <field id="tke_max_ter"  field_ref="tke_ter" operation="maximum"/>
     533        <field id="tke_max_lic"  field_ref="tke_lic" operation="maximum"/>
     534        <field id="tke_max_oce"  field_ref="tke_oce" operation="maximum"/>
     535        <field id="tke_max_sic"  field_ref="tke_sic" operation="maximum"/>
    533536        <field id="l_mix_ter"    long_name="PBL mixing length ter"    unit="m" />
    534537        <field id="l_mix_lic"    long_name="PBL mixing length lic"    unit="m" />
     
    549552        <field id="cldtau"    long_name="Cloud optical thickness"    unit="1" />
    550553        <field id="cldemi"    long_name="Cloud optical emissivity"    unit="1" />
    551         <field id="tke_max"    long_name="TKE max"    unit="m2/s2" operation="maximum"/>
     554        <field id="tke_max"  long_name="TKE max" field_ref="tke" operation="maximum"/>
    552555        <field id="kz"    long_name="Kz melange"    unit="m2/s" enabled="FALSE" />
    553         <field id="kz_max"    long_name="Kz melange max"    unit="m2/s" enabled="FALSE" operation="maximum" />
     556        <field id="kz_max"    field_ref="kz" long_name="Kz melange max" enabled="FALSE" operation="maximum" />
    554557        <field id="concso4"    long_name="Concentration of Sulfate "    unit="kg/m3" />
    555558        <field id="concno3"    long_name="Concentration of Nitrate "    unit="kg/m3" />
     
    566569        <field id="ovapinit"    long_name="Specific humidity (begin of timestep)"    unit="kg/kg" />
    567570        <field id="oliq"    long_name="Liquid water"    unit="kg/kg" />
    568         <field id="ocond"    long_name="Condensed water"    unit="kg/kg" />
     571        <field id="ocond"    long_name="Condensed water"    unit="kg/kg" >oliq+this</field>
    569572        <field id="geop"    long_name="Geopotential height"    unit="m2/s2" />
    570573        <field id="geoph" field_ref="geop" long_name="Geopotential height" unit="m" > geop/9.80665 </field>
     
    575578        <field id="paprs"    long_name="Air pressure Inter-Couches"    unit="Pa" />
    576579        <field id="mass"    long_name="Masse Couches"    unit="kg/m2" />
    577         <field id="zfull"    long_name="Altitude of full pressure levels"    unit="m" />
     580        <field id="zfull"    long_name="Altitude of full pressure levels"    unit="m"  > geop / $RG </field>
    578581        <field id="zhalf"    long_name="Altitude of half pressure levels"    unit="m" />
    579582        <field id="rneb"    long_name="Cloud fraction"    unit="-" />
     
    582585        <field id="rneblsvol" long_name="LS Cloud fraction by volume"    unit="-" />
    583586        <field id="rhum"    long_name="Relative humidity"    unit="-" />
    584         <field id="ozone"    long_name="Ozone mole fraction"    unit="-" />
    585         <field id="ozone_daylight"    long_name="Daylight ozone mole fraction"    unit="-" />
     587        <field id="wo_1"    long_name="Temporary"/>
     588        <field id="ozone"    long_name="Ozone mole fraction"    unit="-" > wo_1 * $dobson_u * 1E3 / mass / $rmo3 * $rmd </field>
     589        <field id="ozone_daylight"    long_name="Daylight ozone mole fraction"    unit="-" > this * $dobson_u * 1E3 / mass / $rmo3 * $rmd </field>
    586590        <field id="upwd"    long_name="saturated updraft"    unit="kg/m2/s" />
    587591        <field id="ep"    long_name="ep"    unit="su" />
     
    603607        <field id="lcc3dcon"    long_name="Convective cloud liquid fraction"    unit="1" />
    604608        <field id="lcc3dstra"    long_name="Stratiform cloud liquid fraction"    unit="1" />
    605         <field id="cldwatmxrat"    long_name="Cloud Water Mixing Ratio"    unit="1" />
     609        <field id="cldwatmxrat"    long_name="Cloud Water Mixing Ratio"    unit="1" > 1-this </field>
    606610        <field id="icc3dcon"    long_name="Mass Fraction of Convective Cloud Ice"    unit="kg/kg" />
    607611        <field id="icc3dstra"    long_name="Mass Fraction of Stratiform Cloud Ice "    unit="kg/kg" />
     
    612616        <field id="dnwd"    long_name="saturated downdraft"    unit="kg/m2/s" />
    613617        <field id="dnwd0"    long_name="unsat. downdraft"    unit="kg/m2/s" />
    614         <field id="mc"    long_name="Convective mass flux"    unit="kg/m2/s" />
     618        <field id="mc"    long_name="Convective mass flux"    unit="kg/m2/s" > -dnwd + dnwd0 + upwd + f_th </field>
    615619        <field id="dtdyn"    long_name="Dynamics dT"    unit="K/s" />
    616620        <field id="dqdyn"    long_name="Dynamics dQ"    unit="(kg/kg)/s" />
     
    619623        <field id="dudyn"    long_name="Dynamics dU"    unit="m/s2" />
    620624        <field id="dvdyn"    long_name="Dynamics dV"    unit="m/s2" />
    621         <field id="dtcon"    long_name="Convection dT"    unit="K/s" />
    622         <field id="ducon"    long_name="Convection du"    unit="m/s2" />
    623         <field id="dvcon"    long_name="Convection dv"    unit="m/s2" />
    624         <field id="dqcon"    long_name="Convection dQ"    unit="(kg/kg)/s" />
    625         <field id="dtwak"    long_name="Wake dT"    unit="K/s" />
    626         <field id="dqwak"    long_name="Wake dQ"    unit="(kg/kg)/s" />
     625        <field id="dtcon"    long_name="Convection dT"    unit="K/s" > this/$pdtphys </field>
     626        <field id="ducon"    long_name="Convection du"    unit="m/s2" > this/$pdtphys </field>
     627        <field id="dvcon"    long_name="Convection dv"    unit="m/s2" > this/$pdtphys </field>
     628        <field id="dqcon"    long_name="Convection dQ"    unit="(kg/kg)/s" > this/$pdtphys </field>
     629        <field id="dtwak"    long_name="Wake dT"    unit="K/s" > this/$pdtphys </field>
     630        <field id="dqwak"    long_name="Wake dQ"    unit="(kg/kg)/s" > this/$pdtphys </field>
    627631        <field id="wake_deltat"    long_name="wake_deltat"    unit="-" />
    628632        <field id="wake_deltaq"    long_name="wake_deltaq"    unit="-" />
     
    640644        <field id="dqlsc"    long_name="Condensation dQ"    unit="(kg/kg)/s" />
    641645        <field id="beta_prec"    long_name="LS Conversion rate to prec"    unit="(kg/kg)/s" />
    642         <field id="dtvdf"    long_name="Boundary-layer dT"    unit="K/s" />
    643         <field id="dtvdf_x"  long_name="dtvdf off_wake"    unit="K/s" />
    644         <field id="dtvdf_w"  long_name="dtvdf within_wake"    unit="K/s" />
    645         <field id="dtdis"    long_name="TKE dissipation dT"    unit="K/s" />
    646         <field id="dqvdf"    long_name="Boundary-layer dQ"    unit="(kg/kg)/s" />
    647         <field id="dqvdf_x" long_name="dqvdf off_wake" unit="(kg/kg)/s" />
    648         <field id="dqvdf_w" long_name="dqvdf within_wake" unit="(kg/kg)/s" />
    649         <field id="dteva"    long_name="Reevaporation dT"    unit="K/s" />
    650         <field id="dqeva"    long_name="Reevaporation dQ"    unit="(kg/kg)/s" />
    651         <field id="dqlscth"    long_name="dQ therm."    unit="(kg/kg)/s" />
    652         <field id="dqlscst"    long_name="dQ strat."    unit="(kg/kg)/s" />
    653         <field id="dtlscth"    long_name="dQ therm."    unit="K/s" />
    654         <field id="dtlscst"    long_name="dQ strat."    unit="K/s" />
    655         <field id="ptconvth"    long_name="POINTS CONVECTIFS therm."    unit="-" />
     646        <field id="dtvdf"    long_name="Boundary-layer dT"    unit="K/s" > this/$pdtphys </field>
     647        <field id="dtvdf_x"  long_name="dtvdf off_wake"    unit="K/s"> this/$pdtphys </field>
     648        <field id="dtvdf_w"  long_name="dtvdf within_wake"    unit="K/s"> this/$pdtphys </field>
     649        <field id="dtdis"    long_name="TKE dissipation dT"    unit="K/s" > this/$pdtphys </field>
     650        <field id="dqvdf"    long_name="Boundary-layer dQ"    unit="(kg/kg)/s" > this/$pdtphys </field>
     651        <field id="dqvdf_x" long_name="dqvdf off_wake" unit="(kg/kg)/s" > this/$pdtphys </field>
     652        <field id="dqvdf_w" long_name="dqvdf within_wake" unit="(kg/kg)/s" > this/$pdtphys </field>
     653        <field id="dteva"    long_name="Reevaporation dT"    unit="K/s" > this/$pdtphys </field>
     654        <field id="dqeva"    long_name="Reevaporation dQ"    unit="(kg/kg)/s" > this/$pdtphys </field>
     655        <field id="dqlscth"    long_name="dQ therm."    unit="(kg/kg)/s" > this/$pdtphys </field>
     656        <field id="dqlscst"    long_name="dQ strat."    unit="(kg/kg)/s" > this/$pdtphys </field>
     657        <field id="dtlscth"    long_name="dQ therm."    unit="K/s" > this/$pdtphys </field>
     658        <field id="dtlscst"    long_name="dQ strat."    unit="K/s" > this/$pdtphys </field>
     659        <field id="ptconvth"    long_name="POINTS CONVECTIFS therm."    unit="-" > this/=0 </field>
    656660        <field id="ptconv"    long_name="POINTS CONVECTIFS"    unit="-" />
    657661        <field id="ratqs"    long_name="RATQS"    unit="-" />
    658         <field id="dtthe"    long_name="Thermal dT"    unit="K/s" />
    659         <field id="duthe"    long_name="Thermal du"    unit="m/s2" />
    660         <field id="dvthe"    long_name="Thermal dv"    unit="m/s2" />
     662        <field id="dtthe"    long_name="Thermal dT"    unit="K/s" > this/$pdtphys - dtajs </field>
     663        <field id="duthe"    long_name="Thermal du"    unit="m/s2" > this/$pdtphys </field>
     664        <field id="dvthe"    long_name="Thermal dv"    unit="m/s2" > this/$pdtphys </field>
    661665        <field id="f_th"    long_name="Thermal plume mass flux"    unit="kg/(m2*s)" />
    662666        <field id="e_th"    long_name="Thermal plume entrainment"    unit="K/s" />
     
    671675        <field id="d_th"    long_name="Thermal plume detrainment"    unit="K/s" />
    672676        <field id="dqthe"    long_name="Thermal dQ"    unit="(kg/kg)/s" />
    673         <field id="dtajs"    long_name="Dry adjust. dT"    unit="K/s" />
    674         <field id="dqajs"    long_name="Dry adjust. dQ"    unit="(kg/kg)/s" />
    675         <field id="dtswr"    long_name="SW radiation dT"    unit="K/s" />
    676         <field id="dtsw0"    long_name="CS SW radiation dT"    unit="K/s" />
    677         <field id="dtlwr"    long_name="LW radiation dT"    unit="K/s" />
    678         <field id="dtlw0"    long_name="CS LW radiation dT"    unit="K/s" />
    679         <field id="dtec"    long_name="Cinetic dissip dT"    unit="K/s" />
    680         <field id="duvdf"    long_name="Boundary-layer dU"    unit="m/s2" />
    681         <field id="dvvdf"    long_name="Boundary-layer dV"    unit="m/s2" />
    682         <field id="duoro"    long_name="Orography dU"    unit="m/s2" />
    683         <field id="dvoro"    long_name="Orography dV"    unit="m/s2" />
    684         <field id="dulif"    long_name="Orography dU"    unit="m/s2" />
    685         <field id="dvlif"    long_name="Orography dV"    unit="m/s2" />
    686         <field id="du_gwd_rando"    long_name="peu importe"    unit="W/m2" />
    687         <field id="dv_gwd_rando"    long_name="peu importe"    unit="W/m2" />
    688         <field id="du_gwd_hines" long_name="Hines GWD dU" unit="m/s2" />
    689         <field id="dv_gwd_hines" long_name="Hines GWD dV" unit="m/s2" />
    690         <field id="du_gwd_front" long_name="Fronts GWD dU" unit="m/s2" />
    691         <field id="dv_gwd_front" long_name="Fronts GWD dV" unit="m/s2" />
     677        <field id="dtajs"    long_name="Dry adjust. dT"    unit="K/s" > this/$pdtphys </field>
     678        <field id="dqajs"    long_name="Dry adjust. dQ"    unit="(kg/kg)/s" > this/$pdtphys </field>
     679        <field id="dtswr"    long_name="SW radiation dT"    unit="K/s" > this/$pdtphys </field>
     680        <field id="dtsw0"    long_name="CS SW radiation dT"    unit="K/s" > this/$pdtphys </field>
     681        <field id="dtlwr"    long_name="LW radiation dT"    unit="K/s" > this/$pdtphys </field>
     682        <field id="dtlw0"    long_name="CS LW radiation dT"    unit="K/s" > this/$pdtphys </field>
     683        <field id="dtec"    long_name="Cinetic dissip dT"    unit="K/s" > this/$pdtphys </field>
     684        <field id="duvdf"    long_name="Boundary-layer dU"    unit="m/s2" > this/$pdtphys </field>
     685        <field id="dvvdf"    long_name="Boundary-layer dV"    unit="m/s2" > this/$pdtphys </field>
     686        <field id="duoro"    long_name="Orography dU"    unit="m/s2" > this/$pdtphys </field>
     687        <field id="dvoro"    long_name="Orography dV"    unit="m/s2" > this/$pdtphys </field>
     688        <field id="dulif"    long_name="Orography dU"    unit="m/s2" > this/$pdtphys </field>
     689        <field id="dvlif"    long_name="Orography dV"    unit="m/s2" > this/$pdtphys </field>
     690        <field id="du_gwd_rando"    long_name="peu importe"    unit="W/m2" > this/$pdtphys </field>
     691        <field id="dv_gwd_rando"    long_name="peu importe"    unit="W/m2" > this/$pdtphys </field>
     692        <field id="du_gwd_hines" long_name="Hines GWD dU" unit="m/s2" > this/$pdtphys </field>
     693        <field id="dv_gwd_hines" long_name="Hines GWD dV" unit="m/s2" > this/$pdtphys </field>
     694        <field id="du_gwd_front" long_name="Fronts GWD dU" unit="m/s2" > this/$pdtphys </field>
     695        <field id="dv_gwd_front" long_name="Fronts GWD dV" unit="m/s2" > this/$pdtphys </field>
    692696        <field id="east_gwstress" long_name="Eastward GW Stress" unit="Pa" />
    693697        <field id="west_gwstress" long_name="Westward GW Stress" unit="Pa" />
    694         <field id="dtoro"    long_name="Orography dT"    unit="K/s" />
    695         <field id="dtlif"    long_name="Orography dT"    unit="K/s" />
    696         <field id="dthin"    long_name="Hines GWD dT"    unit="K/s" />
     698        <field id="dtoro"    long_name="Orography dT"    unit="K/s" > this/$pdtphys </field>
     699        <field id="dtlif"    long_name="Orography dT"    unit="K/s" > this/$pdtphys </field>
     700        <field id="dthin"    long_name="Hines GWD dT"    unit="K/s" > this/$pdtphys </field>
    697701        <field id="rsu"    long_name="SW upward radiation"    unit="W m-2" />
    698702        <field id="rsd"    long_name="SW downward radiation"    unit="W m-2" />
     
    705709        <field id="rlucs"    long_name="LW CS upward radiation"    unit="W m-2" />
    706710        <field id="rldcs"    long_name="LW CS downward radiation"    unit="W m-2" />
    707         <field id="tnt"    long_name="Tendency of air temperature"    unit="K s-1" />
     711        <field id="tnt"    long_name="Tendency of air temperature"    unit="K s-1" > dtphy + dtdyn </field>
    708712        <field id="tntc"    long_name="Tendency of air temperature due to Moist Convection"    unit="K s-1" />
    709         <field id="tntr"    long_name="Air temperature tendency due to Radiative heating"    unit="K s-1" />
    710         <field id="tntscpbl"    long_name="Air temperature tendency due to St cloud and precipitation and BL mixing"    unit="K s-1" />
    711         <field id="tnhus"    long_name="Tendency of specific humidity"    unit="s-1" />
     713        <field id="tntr"    long_name="Air temperature tendency due to Radiative heating"    unit="K s-1" > dtswr + dtlwr </field>
     714        <field id="tntscpbl"    long_name="Air temperature tendency due to St cloud and precipitation and BL mixing"    unit="K s-1" > dtlsc +  dteva + dtvdf </field>
     715        <field id="tnhus"    long_name="Tendency of specific humidity"    unit="s-1" > dqphy + dqdyn </field>
    712716        <field id="tnhusc"    long_name="Tendency of specific humidity due to convection"    unit="s-1" />
    713         <field id="tnhusscpbl"    long_name="Tendency of Specific humidity due to ST cl, precip and BL mixing"    unit="s-1" />
     717        <field id="tnhusscpbl"    long_name="Tendency of Specific humidity due to ST cl, precip and BL mixing"    unit="s-1" > dqlsc + dqeva </field>
    714718        <field id="evu"    long_name="Eddy viscosity coefficient for Momentum Variables"    unit="m2 s-1" />
    715         <field id="h2o"    long_name="Mass Fraction of Water"    unit="1" />
     719        <field id="h2o"    long_name="Mass Fraction of Water"    unit="1" > ovap + oliq </field>
    716720        <field id="mcd"    long_name="Downdraft COnvective Mass Flux"    unit="kg/(m2*s)" />
    717721        <field id="dmc"    long_name="Deep COnvective Mass Flux"    unit="kg/(m2*s)" />
     
    783787      <field id="wap" long_name="Lagrangian tendency of air pressure" unit="Pa s-1" detect_missing_value=".true." />
    784788<!--  <field id="psbg" long_name="Pressure sfce below ground" unit="%" detect_missing_value=".true." /> -->
    785       <field id="tro3" long_name="Ozone mole fraction" unit="1e-9" detect_missing_value=".true." />
    786       <field id="tro3_daylight" long_name="Daylight ozone mole fraction" unit="1e-9" detect_missing_value=".true." />
     789      <field id="tro3" long_name="Ozone mole fraction" unit="1e-9" detect_missing_value=".true." > (this/=$missing_val)?this*1.e+9:$missing_val </field>
     790      <field id="tro3_daylight" long_name="Daylight ozone mole fraction" unit="1e-9" detect_missing_value=".true." > (this/=$missing_val)?this*1.e+9:$missing_val </field>
    787791      <field id="uv" long_name="uv" unit="m2/s2" detect_missing_value=".true." />
    788792      <field id="vxq" long_name="vxq" unit="m/s * (kg/kg)" detect_missing_value=".true." />
Note: See TracChangeset for help on using the changeset viewer.