search

EC-Earth / ece2cmor3

Post-processing and cmorization of ec-earth output
Apache License 2.0
13 stars 6 forks source link

Adding the IFS variables which are identified in step 2 #143

Closed treerink closed 5 years ago

treerink commented 6 years ago

These IFS variables have to be added to the ifspar.json with their expressions / unit conversions:

LImon       tasIs       T2M as in IFS (grib code 128.167)
LImon       tsIs        We use soil layer 1 temeperature as the proxy (grib 128.139)
LImon       tsnIs       This is Grib 128.238. But we never test it
LImon       prsnIs      grib 128.144
LImon       prraIs      grib 128.(142+143-144)
LImon       sblIs       grib 128.182
LImon       snmIs       grib 128.45 
LImon       mrroIs      grib 128.205
LImon       hflsIs      grib 128.147
LImon       hfssIs      grib 128.146
LImon       rsdsIs      grib 128.169
LImon       rsusIs      grib 128.(169-176)
LImon       rldsIs      grib 128.175
LImon       rlusIs      grib 128.(175-177)
CFmon       tntr        Available in IFS: T-tendency from radiation: grib 128.95
CFmon       tntc        Available in IFS: T-tendency from convection : grib 128.105
CFmon       tnhus       Adding all the q-tendencies, thus: grib 128.94 + 128.99 + 128.106 + 128.110.  Alternatively, in IFS: just estimating the delta q per month. So far no direct grib code for the totoal q-tendency found
CFmon       tnhusc      Available in IFS: q-tendency from convection: grib 128.106
CFmon       tnhusmp     Adding all the q-tendencies without advection, thus: grib 128.99 + 128.106 + 128.110.
EmonZ       tntc        Available in IFS: T-tendency from convection : grib 128.105
6hrPlev     wsgmax10m   Available in IFS: 10 metre wind gust in the 6 hours: grib 128.123
AERmon      cltc        Available in IFS: convective cloud cover ccc: grib 128.185
AERmon      lwp     Available in IFS: Total column liquid water: grib 128.78
AERmon      wa      Available in IFS:Vertical velocity:w: grib 128.135 requires unit conversion from Pa.s-1 to m/s
Eday        zmla        Available in IFS: Boundary layer height: grib 128.159
Emon        rls     Available in IFS: Surface net thermal radiation: grib 128.177
Emon        rss     Available in IFS: Surface net solar radiation: grib 128.176

Efx     ksat        H-TESSEL parameters inside the model that is prescribed for each dominant soil texture  (not in the output; no grib code associated).
Efx     rootdsl     H-TESSEL parameters inside the model that is prescribed for each dominant soil texture (not in the output, no grib code associated).
Efx     slthick     Available in IFS: independent of lon, lat. Top soil layer: 0-7 cm, Soil layer 2: 7-28 cm, Soil layer 3: 28-100 cm, Soil layer 4: 100-289 cm, added as issue #126
LImon       sncIs       Not available in IFS. Although it could be calculated from tile fractions and written out as extra output
LImon       sftgif      This is the land ice mask and will be an extra variable in IFS (thomas: via PEXTRA?)
LImon       sftgrf      For Greenland this is the same as above sftgif. We do not have Antarctic ice sheet.
LImon       orogIs      This is grib 128.129 at surface. But be aware not missing up with the grib 128.129 for free atmospehre.
CFmon       tnt     Alternatively, just estimating the delta T per month. No direct grib code for the totoal T-tendency found. In IFS from Cycle 39R1: add all the T-Tendencies: grib 128.93 + 128.95 + 128.98 + 128.102 + 128.105 + 128.109. But with IFS cycle 36 the T-tendency of gravity wave drag grib 128.102 is bugged until Cycle 39R1. This has been checked by Gijs with ECMWF: https://software.ecmwf.int/wiki/pages/viewpage.action?pageId=97384581
CFsubhr     longitude   This site lon,lat is specifed by the experiment we guess, and should be available in the netcdf file. So far it doesn't seem to be specified in the data request: in the CMIP6_coordinate.json table file the requested": "" is empty for site.
CFsubhr     latitude    This site lon,lat is specifed by the experiment we guess, and should be available in the netcdf file. So far it doesn't seem to be specified in the data request: in the CMIP6_coordinate.json table file the requested": "" is empty for site.
treerink commented 6 years ago

Two points here:

  1. In ece2cmor3 several of the PEXTRA grib codes have been used as codes for derived variables, this needs to be changed such that derived variables have their own table, see #99.
  2. The tendency related variables here are based on grib codes as provided on the ecmwf grib code page for OpenIFS and this matches the latest cylcle, while we use cycle 36. The page mentioned that these grib codes might differ depending on the cycly, so this needs verification
treerink commented 6 years ago

I think those ones could be make directly available as well: CFmon tnhus CFmon tnhusmp Omon prra

treerink commented 6 years ago 
Implemented:               LImon         tasIs        T2M as in IFS (grib code 128.167)
Implemented:               LImon         tsIs         We use soil layer 1 temeperature as the proxy (grib 128.139)
Implemented:               LImon         tsnIs        This is Grib 128.238. But we never test it
Implemented:               LImon         prsnIs       grib 128.144
Implemented:               LImon         prraIs       grib 128.(142+143-144) ["source": "122.129"]  ["var122=var142+var143-var144"]
Implemented:               LImon         sblIs        grib 128.182
Implemented:               LImon         snmIs        grib 128.45
Implemented:               LImon         mrroIs       grib 128.205
Implemented:               LImon         hflsIs       grib 128.147
Implemented:               LImon         hfssIs       grib 128.146
Implemented:               LImon         rsdsIs       grib 128.169
Implemented:               LImon         rldsIs       grib 128.175
Implemented:               CFmon         tntr         Available in IFS: T-tendency from radiation: grib 128.95
Implemented:               CFmon         tntc         Available in IFS: T-tendency from convection : grib 128.105
Implemented:               CFmon         tnhusc       Available in IFS: q-tendency from convection: grib 128.106
Implemented:               EmonZ         tntc         Available in IFS: T-tendency from convection : grib 128.105
Implemented:               AERmon        cltc         Available in IFS: convective cloud cover ccc: grib 128.185
Implemented:               AERmon        lwp          Available in IFS: Total column liquid water: grib 128.78
Implemented:               AERmon        wa           Available in IFS:Vertical velocity:w: grib 128.135 requires unit conversion from Pa.s-1 to m/s  ["source": "123.129"]  ["var123=-287*var130*(1+0.61*var133)*var135/(9.81*var54)"]
Implemented:               Eday          zmla         Available in IFS: Boundary layer height: grib 128.159
Implemented:               Emon          rls          Available in IFS: Surface net thermal radiation: grib 128.177
Implemented:               Emon          rss          Available in IFS: Surface net solar radiation: grib 128.176

Implemented:               6hrPlev       wsgmax10m    Available in IFS: 10 metre wind gust in the 6 hours: grib 128.123  ["source": "49.128"] Note : Changed grib code identification
Implemented:               LImon         rsusIs       grib 128.(169-176)  ["source": "95.129"] ["var95=var176-var169"]  CHECK MINUS SIGN
Implemented:               LImon         rlusIs       grib 128.(175-177)  ["source": "96.129"] ["var96=var177-var175"]  CHECK MINUS SIGN

Implemented:               CFmon         tnhus        Adding all the q-tendencies, thus: grib 128.94 + 128.99 + 128.106 + 128.110.  Alternatively, in IFS: just estimating the delta q per month. So far no direct grib code for the totoal q-tendency found
Implemented:               CFmon         tnhusmp      Adding all the q-tendencies without advection, thus: grib 128.99 + 128.106 + 128.110.
Implemented:               ImonGre,Omon  prra         Large-scale precipitation (rain+snow) + Convective precipitation (rain+snow) - Snowfall: grib 128.142 + 128.143 - 128.144
                           ImonGre       mrroLi       IFS Surface runoff grib 128.8 but for EC_Earth-GrIs additional melt etc is included
                           E3hr          prrc         Availlable in IFS: Precip. flux from convection liquid grib 128.107, this is a 3D field so the surface field has to be extracted from this. So only level 91 needs to be outputted.

Implemented:               LImon         orogIs       This is grib 128.129 at surface. But be aware not missing up with the grib 128.129 for free atmospehre.

                           LImon         sncIs        Not available in IFS. Although it could be calculated from tile fractions and written out as extra output
                           LImon         sftgif       This is the land ice mask and will be an extra variable in IFS (thomas: via PEXTRA?)
                           LImon         sftgrf       For Greenland this is the same as above sftgif. We do not have Antarctic ice sheet.
Unclear yet                Efx           wilt         Wilting point is defined in HTESSEL as a fraction (different for each soil texture in HTESSEL; variable slt grib code 43 - Table 128) of field capacity. (Alessandri A.: an Effective Wilting Point can be defined ... same as below discussion for the field capacity)
Unclear yet                Efx           fldcapacity  Maximum amount of water the soil can hold (different for each soil texture in HTESSEL; variable slt grib code 43 - Table 128). It is constant map in time in H-TESSEL (Alessandri A.: an Effective Field Capacity can be defined for coupling analysis pourpouses that will change with the effective vegetation cover; as soon as we have daily effective vegetation cover in the outputs we can compute as well the Effective Field capacity).
Unclear yet                Efx           clayfrac     soil texture (could be derived by soil HTESSEL soil types [1-6]: variable slt grib code 43 - Table 128): very fine (& fine?) in H-TESSEL boundary condition constant in time. It is available in a kind of mask: for each type a number 1-6 if this type at a lon,lat location the number corresponding with this type is in the mask. (No grib code we think.)
Unclear yet                Efx           siltfrac     soil texture (could be derived by soil HTESSEL soil types [1-6]: variable slt grib code 43 - Table 128): medium (& medium fine?) in H-TESSEL boundary condition constant in time. It is available in a kind of mask: for each type a number 1-6 if this type at a lon,lat location the number corresponding with this type is in the mask. (No grib code we think.)
Unclear yet                Efx           sandfrac     soil texture (could be derived by soil HTESSEL soil types [1-6]: variable slt grib code 43 - Table 128): coarse in H-TESSEL boundary condition constant in time. It is available in a kind of mask: for each type a number 1-6 if this type at a lon,lat location the number corresponding with this type is in the mask. (No grib code we think.)
H-TESSEL parameter         Efx           ksat         H-TESSEL parameters inside the model that is prescribed for each dominant soil texture  (not in the output; no grib code associated).
H-TESSEL parameter         Efx           rootdsl      H-TESSEL parameters inside the model that is prescribed for each dominant soil texture (not in the output, no grib code associated).
                           Efx           slthick      Available in IFS: independent of lon, lat. Top soil layer: 0-7 cm, Soil layer 2: 7-28 cm, Soil layer 3: 28-100 cm, Soil layer 4: 100-289 cm, added as issue #126

delta field                CFmon         tnt          Alternatively, just estimating the delta T per month. No direct grib code for the totoal T-tendency found. In IFS from Cycle 39R1: add all the T-Tendencies: grib 128.93 + 128.95 + 128.98 + 128.102 + 128.105 + 128.109. But with IFS cycle 36 the T-tendency of gravity wave drag grib 128.102 is bugged until Cycle 39R1. This has been checked by Gijs with ECMWF: https://software.ecmwf.int/wiki/pages/viewpage.action?pageId=97384581
                           CFsubhr       longitude    This site lon,lat is specifed by the experiment we guess, and should be available in the netcdf file. So far it doesn't seem to be specified in the data request: in the CMIP6_coordinate.json table file the requested": "" is empty for site.
                           CFsubhr       latitude     This site lon,lat is specifed by the experiment we guess, and should be available in the netcdf file. So far it doesn't seem to be specified in the data request: in the CMIP6_coordinate.json table file the requested": "" is empty for site.

PEXTRA => grib code issue  CFmon         cltisccp     COSP: ISCCP total cloud area, CVEXTRA(5)='94 ISCCP_TOTALCLDAREA'
PEXTRA => grib code issue  CFmon         albisccp     COSP: ISCCP cloud albedo, CVEXTRA(7)='96 ISCCP_MEANALBEDOCLD'
PEXTRA => grib code issue  CFmon         pctisccp     COSP: ISCCP cloud top pressure, CVEXTRA(6)='95 ISCCP_MEANPTOP'
PEXTRA => grib code issue  CFmon         cltcalipso   COSP: CALIPSO total cloud cover, CVEXTRA(4)='93 CALIPSO_CLDLAYER TOTAL'
PEXTRA => grib code issue  CFmon         cllcalipso   COSP: CALIPSO low cloud cover, CVEXTRA(1)='90 CALIPSO_CLDLAYER LOW'
PEXTRA => grib code issue  CFmon         clmcalipso   COSP: CALIPSO mid cloud cover, CVEXTRA(2)='91 CALIPSO_CLDLAYER MID'
PEXTRA => grib code issue  CFmon         clhcalipso   COSP: CALIPSO high cloud cover, CVEXTRA(3)='92 CALIPSO_CLDLAYER HIGH'

                           AERmon        cdnc         In runtime/classic/ctrl/namelist.ifs.cloudact+diag.sh CVEXTRA(1)='CDNC' which is a PEXTRA variable.

Need identification:       AERmon        rlutaf       Available from double radiation call in IFS. See also PEXTRA issue #403
Need identification:       AERmon        rlutcsaf     Available from double radiation call in IFS. See also PEXTRA issue #403
Need identification:       AERmon        rsutaf       Available from double radiation call in IFS. See also PEXTRA issue #403   aerosol free
Need identification:       AERmon,E3hrPt rsutcsaf     Available from double radiation call in IFS. See also PEXTRA issue #403   aerosol free

Need identification:       E3hrPt        rsdscsaf     Available from double radiation call in IFS. See also PEXTRA issue #403
Need identification:       E3hrPt        rsuscsaf     Available from double radiation call in IFS. See also PEXTRA issue #403
Need identification:       E3hrPt        sza          In the IFS code we found: A CALL TO SUBROUTINE *SOLANG* GIVES FIELDS OF SOLAR ZENITH. Best would be to copy the IFS routine/formula to ece2cmor3 and produce this field off-line.

Needs check                Eday          tauupbl      Available in IFS: Eastward turbulent surface stress (grib 128.180) plus the Eastward gravity wave surface stress (grib 128.195).  Check with data request or DynVar people if this is indeed what is asked here.
Needs check                Eday          tauvpbl      Available in IFS: Northward turbulent surface stress (grib 128.181) plus the Northward gravity wave surface stress (grib 128.196).  Check with data request or DynVar people if this is indeed what is asked here.

Implemented LS3MIP step 3: Eday          hfdsl        Available in IFS: grib code 146 + 147 : convert to a time step average
Implemented LS3MIP step 3: Eday          hfmlt        Available in IFS: grib code 45 is amount of melt in mweq, so convert to melt energy
Implemented LS3MIP step 3: Eday          hfsbl        Available in IFS: grib code 44 is amount of snow evaporation in mweq, so convert to sublimation energy
Implemented LS3MIP step 3: Eday          tau          Available in IFS: grib code 180 and 181: squareroot of the squares of these two components
Implemented LS3MIP step 3: Eday          tr           Available in IFS: Skin temperature could be used code 235
Implemented LS3MIP step 3: Eday          albsn        Available in IFS: grib code 32
Implemented LS3MIP step 3: Eday          esn          Available in IFS: code 44
Implemented LS3MIP step 3: Eday          cw           Available in IFS: Same as interception reservoir code: 198
treerink commented 6 years ago

The tnhus, tnhusmp, prra variables have been implemented now. This status is adjusted in the list in the comment above.

treerink commented 6 years ago

Added in the overview above a block of LS3MIP variables coming from step 3.

treerink commented 6 years ago

For tauupbl and tauvpbl see also #191

goord commented 6 years ago

LS3MIP variables were added in 0c630a1f039386abdd4f45c04f17f85562c5504c

goord commented 6 years ago

Link to the document with extra EC-Earth PEXTRA variables: https://docs.google.com/spreadsheets/d/17ipCYh48U0N0or8b0qG-_ALGmeD79wwK-s4W38r_YtM/edit#gid=1285118949

treerink commented 5 years ago

Starting to implement the PEXTRA variables on the ec-earth table 126 in the grib_table_126 branch.

treerink commented 5 years ago

To test the DynVar MIP pextra variables in this branch:

cp list-of-ignored-cmpi6-requested-variables-for-DynVar-pextra-testing.xlsx list-of-ignored-cmpi6-requested-variables.xlsx
goord commented 5 years ago

Tendency variables were added in rev. 830473adbb1410abcf6a779cda38b81b051ceacf

treerink commented 5 years ago

You mean all pextra variables added? Ready for genecec to test?

goord commented 5 years ago

Yes you can test

treerink commented 5 years ago

Looks like genecec is fine with the changes. No crash encountered, and see some additional 126 table grib codes in the ppt files and increased IFS volume estimates, as expected.

The full test has to be done now, run some experiment with EC-Earth with this additions (that however requires also the specific PEXTRA IFS namelist configuration, and these are different for the different PEXTRA sets of variables, and unfortunately just from the grib codes (table 126) it is not directly visible which sets needs to be addressed, maybe we can add comments about that in the ppt file with drq2ppt to make it easier and to prevent incorrect usage). Thereafter the cmorization needs to be tested on those EC-Earth results.

treerink commented 5 years ago 
Implemented:                LImon         tasIs        T2M as in IFS (grib code 128.167)
Implemented:                LImon         tsIs         We use soil layer 1 temeperature as the proxy (grib 128.139)
Implemented:                LImon         tsnIs        This is Grib 128.238. But we never test it
Implemented:                LImon         prsnIs       grib 128.144
Implemented:                LImon         prraIs       grib 128.(142+143-144) ["source": "122.129"]  ["var122=var142+var143-var144"]
Implemented:                LImon         sblIs        grib 128.182
Implemented:                LImon         snmIs        grib 128.45
Implemented:                LImon         mrroIs       grib 128.205
Implemented:                LImon         hflsIs       grib 128.147
Implemented:                LImon         hfssIs       grib 128.146
Implemented:                LImon         rsdsIs       grib 128.169
Implemented:                LImon         rldsIs       grib 128.175
Implemented:                CFmon         tntr         Available in IFS: T-tendency from radiation: grib 128.95
Implemented:                CFmon         tntc         Available in IFS: T-tendency from convection : grib 128.105
Implemented:                CFmon         tnhusc       Available in IFS: q-tendency from convection: grib 128.106
Implemented:                EmonZ         tntc         Available in IFS: T-tendency from convection : grib 128.105
Implemented:                AERmon        cltc         Available in IFS: convective cloud cover ccc: grib 128.185
Implemented:                AERmon        lwp          Available in IFS: Total column liquid water: grib 128.78
Implemented:                AERmon        wa           Available in IFS:Vertical velocity:w: grib 128.135 requires unit conversion from Pa.s-1 to m/s  ["source": "123.129"]  ["var123=-287*var130*(1+0.61*var133)*var135/(9.81*var54)"]
Implemented:                Eday          zmla         Available in IFS: Boundary layer height: grib 128.159
Implemented:                Emon          rls          Available in IFS: Surface net thermal radiation: grib 128.177
Implemented:                Emon          rss          Available in IFS: Surface net solar radiation: grib 128.176

Implemented:                6hrPlev       wsgmax10m    Available in IFS: 10 metre wind gust in the 6 hours: grib 128.123  ["source": "49.128"] Note : Changed grib code identification
Implemented:                LImon         rsusIs       grib 128.(169-176)  ["source": "95.129"] ["var95=var176-var169"]  CHECK MINUS SIGN
Implemented:                LImon         rlusIs       grib 128.(175-177)  ["source": "96.129"] ["var96=var177-var175"]  CHECK MINUS SIGN

Implemented:                CFmon         tnhus        Adding all the q-tendencies, thus: grib 128.94 + 128.99 + 128.106 + 128.110.  Alternatively, in IFS: just estimating the delta q per month. So far no direct grib code for the totoal q-tendency found
Implemented:                CFmon         tnhusmp      Adding all the q-tendencies without advection, thus: grib 128.99 + 128.106 + 128.110.
Implemented:                ImonGre,Omon  prra         Large-scale precipitation (rain+snow) + Convective precipitation (rain+snow) - Snowfall: grib 128.142 + 128.143 - 128.144
To do (PISM related):       ImonGre       mrroLi       IFS Surface runoff grib 128.8 but for EC_Earth-GrIs additional melt etc is included
To do:                      E3hr          prrc         Availlable in IFS: Precip. flux from convection liquid grib 128.107, this is a 3D field so the surface field has to be extracted from this. So only level 91 needs to be outputted.

Implemented:                LImon         orogIs       This is grib 128.129 at surface. But be aware not missing up with the grib 128.129 for free atmospehre.

??                          LImon         sncIs        Not available in IFS. Although it could be calculated from tile fractions and written out as extra output
To do (PISM related):       LImon         sftgif       This is the land ice mask and will be an extra variable in IFS (thomas: via PEXTRA?)
To do (PISM related):       LImon         sftgrf       For Greenland this is the same as above sftgif. We do not have Antarctic ice sheet.
Unclear yet                 Efx           wilt         Wilting point is defined in HTESSEL as a fraction (different for each soil texture in HTESSEL; variable slt grib code 43 - Table 128) of field capacity. (Alessandri A.: an Effective Wilting Point can be defined ... same as below discussion for the field capacity)
Unclear yet                 Efx           fldcapacity  Maximum amount of water the soil can hold (different for each soil texture in HTESSEL; variable slt grib code 43 - Table 128). It is constant map in time in H-TESSEL (Alessandri A.: an Effective Field Capacity can be defined for coupling analysis pourpouses that will change with the effective vegetation cover; as soon as we have daily effective vegetation cover in the outputs we can compute as well the Effective Field capacity).
Unclear yet                 Efx           clayfrac     soil texture (could be derived by soil HTESSEL soil types [1-6]: variable slt grib code 43 - Table 128): very fine (& fine?) in H-TESSEL boundary condition constant in time. It is available in a kind of mask: for each type a number 1-6 if this type at a lon,lat location the number corresponding with this type is in the mask. (No grib code we think.)
Unclear yet                 Efx           siltfrac     soil texture (could be derived by soil HTESSEL soil types [1-6]: variable slt grib code 43 - Table 128): medium (& medium fine?) in H-TESSEL boundary condition constant in time. It is available in a kind of mask: for each type a number 1-6 if this type at a lon,lat location the number corresponding with this type is in the mask. (No grib code we think.)
Unclear yet                 Efx           sandfrac     soil texture (could be derived by soil HTESSEL soil types [1-6]: variable slt grib code 43 - Table 128): coarse in H-TESSEL boundary condition constant in time. It is available in a kind of mask: for each type a number 1-6 if this type at a lon,lat location the number corresponding with this type is in the mask. (No grib code we think.)
H-TESSEL parameter          Efx           ksat         H-TESSEL parameters inside the model that is prescribed for each dominant soil texture  (not in the output; no grib code associated).
H-TESSEL parameter          Efx           rootdsl      H-TESSEL parameters inside the model that is prescribed for each dominant soil texture (not in the output, no grib code associated).
To do #126                  Efx           slthick      Available in IFS: independent of lon, lat. Top soil layer: 0-7 cm, Soil layer 2: 7-28 cm, Soil layer 3: 28-100 cm, Soil layer 4: 100-289 cm, added as issue #126

delta field                 CFmon         tnt          Alternatively, just estimating the delta T per month. No direct grib code for the totoal T-tendency found. In IFS from Cycle 39R1: add all the T-Tendencies: grib 128.93 + 128.95 + 128.98 + 128.102 + 128.105 + 128.109. But with IFS cycle 36 the T-tendency of gravity wave drag grib 128.102 is bugged until Cycle 39R1. This has been checked by Gijs with ECMWF: https://software.ecmwf.int/wiki/pages/viewpage.action?pageId=97384581
Ignore subhr tables:        CFsubhr       longitude    This site lon,lat is specifed by the experiment we guess, and should be available in the netcdf file. So far it doesn't seem to be specified in the data request: in the CMIP6_coordinate.json table file the requested": "" is empty for site.
Ignore subhr tables:        CFsubhr       latitude     This site lon,lat is specifed by the experiment we guess, and should be available in the netcdf file. So far it doesn't seem to be specified in the data request: in the CMIP6_coordinate.json table file the requested": "" is empty for site.

Implemented:                CFmon         cltisccp     COSP: ISCCP total cloud area, CVEXTRA(5)='94 ISCCP_TOTALCLDAREA'
Implemented:                CFmon         albisccp     COSP: ISCCP cloud albedo, CVEXTRA(7)='96 ISCCP_MEANALBEDOCLD'
Implemented:                CFmon         pctisccp     COSP: ISCCP cloud top pressure, CVEXTRA(6)='95 ISCCP_MEANPTOP'
Implemented:                CFmon         cltcalipso   COSP: CALIPSO total cloud cover, CVEXTRA(4)='93 CALIPSO_CLDLAYER TOTAL'
Implemented:                CFmon         cllcalipso   COSP: CALIPSO low cloud cover, CVEXTRA(1)='90 CALIPSO_CLDLAYER LOW'
Implemented:                CFmon         clmcalipso   COSP: CALIPSO mid cloud cover, CVEXTRA(2)='91 CALIPSO_CLDLAYER MID'
Implemented:                CFmon         clhcalipso   COSP: CALIPSO high cloud cover, CVEXTRA(3)='92 CALIPSO_CLDLAYER HIGH'

Implemented:                AERmon        cdnc         In runtime/classic/ctrl/namelist.ifs.cloudact+diag.sh CVEXTRA(1)='CDNC' which is a PEXTRA variable.

Implemented #368:           AERmon        rlutaf       Available from double radiation call in IFS. See also PEXTRA issue #403
Implemented #368:           AERmon        rlutcsaf     Available from double radiation call in IFS. See also PEXTRA issue #403
Implemented #368:           AERmon        rsutaf       Available from double radiation call in IFS. See also PEXTRA issue #403   aerosol free
Moved to ignored file #368: AERmon,E3hrPt rsutcsaf     Available from double radiation call in IFS. See also PEXTRA issue #403   aerosol free
Moved to ignored file #368: E3hrPt        rsdscsaf     Available from double radiation call in IFS. See also PEXTRA issue #403
Moved to ignored file #368: E3hrPt        rsuscsaf     Available from double radiation call in IFS. See also PEXTRA issue #403

Moved to ignored file #191: Eday          tauupbl      Available in IFS: Eastward  turbulent surface stress (grib 128.180) plus the Eastward  gravity wave surface stress (grib 128.195).  Check with data request or DynVar people if this is indeed what is asked here.
Moved to ignored file #191: Eday          tauvpbl      Available in IFS: Northward turbulent surface stress (grib 128.181) plus the Northward gravity wave surface stress (grib 128.196).  Check with data request or DynVar people if this is indeed what is asked here.

Implemented LS3MIP step 3:  Eday          hfdsl        Available in IFS: grib code 146 + 147 : convert to a time step average
Implemented LS3MIP step 3:  Eday          hfmlt        Available in IFS: grib code 45 is amount of melt in mweq, so convert to melt energy
Implemented LS3MIP step 3:  Eday          hfsbl        Available in IFS: grib code 44 is amount of snow evaporation in mweq, so convert to sublimation energy
Implemented LS3MIP step 3:  Eday          tau          Available in IFS: grib code 180 and 181: squareroot of the squares of these two components
Implemented LS3MIP step 3:  Eday          tr           Available in IFS: Skin temperature could be used code 235
Implemented LS3MIP step 3:  Eday          albsn        Available in IFS: grib code 32
Implemented LS3MIP step 3:  Eday          esn          Available in IFS: code 44
Implemented LS3MIP step 3:  Eday          cw           Available in IFS: Same as interception reservoir code: 198
treerink commented 5 years ago

Updated the last overview.

treerink commented 5 years ago

Overview of variables which still need further action:

To do:                 E3hr     prrc         Availlable in IFS: Precip. flux from convection liquid grib 128.107, this is a 3D field so the surface field has to be extracted from this. So only level 91 needs to be outputted.

Unclear yet            Efx      wilt         Wilting point is defined in HTESSEL as a fraction (different for each soil texture in HTESSEL; variable slt grib code 43 - Table 128) of field capacity. (Alessandri A.: an Effective Wilting Point can be defined ... same as below discussion for the field capacity)
Unclear yet            Efx      fldcapacity  Maximum amount of water the soil can hold (different for each soil texture in HTESSEL; variable slt grib code 43 - Table 128). It is constant map in time in H-TESSEL (Alessandri A.: an Effective Field Capacity can be defined for coupling analysis pourpouses that will change with the effective vegetation cover; as soon as we have daily effective vegetation cover in the outputs we can compute as well the Effective Field capacity).
Unclear yet            Efx      clayfrac     soil texture (could be derived by soil HTESSEL soil types [1-6]: variable slt grib code 43 - Table 128): very fine (& fine?) in H-TESSEL boundary condition constant in time. It is available in a kind of mask: for each type a number 1-6 if this type at a lon,lat location the number corresponding with this type is in the mask. (No grib code we think.)
Unclear yet            Efx      siltfrac     soil texture (could be derived by soil HTESSEL soil types [1-6]: variable slt grib code 43 - Table 128): medium (& medium fine?) in H-TESSEL boundary condition constant in time. It is available in a kind of mask: for each type a number 1-6 if this type at a lon,lat location the number corresponding with this type is in the mask. (No grib code we think.)
Unclear yet            Efx      sandfrac     soil texture (could be derived by soil HTESSEL soil types [1-6]: variable slt grib code 43 - Table 128): coarse in H-TESSEL boundary condition constant in time. It is available in a kind of mask: for each type a number 1-6 if this type at a lon,lat location the number corresponding with this type is in the mask. (No grib code we think.)
H-TESSEL parameter     Efx      ksat         H-TESSEL parameters inside the model that is prescribed for each dominant soil texture  (not in the output; no grib code associated).
H-TESSEL parameter     Efx      rootdsl      H-TESSEL parameters inside the model that is prescribed for each dominant soil texture (not in the output, no grib code associated).
To do #126             Efx      slthick      Available in IFS: independent of lon, lat. Top soil layer: 0-7 cm, Soil layer 2: 7-28 cm, Soil layer 3: 28-100 cm, Soil layer 4: 100-289 cm, added as issue #126

moved to ignored       CFmon    tnt          Alternatively, just estimating the delta T per month. No direct grib code for the totoal T-tendency found. In IFS from Cycle 39R1: add all the T-Tendencies: grib 128.93 + 128.95 + 128.98 + 128.102 + 128.105 + 128.109. But with IFS cycle 36 the T-tendency of gravity wave drag grib 128.102 is bugged until Cycle 39R1. This has been checked by Gijs with ECMWF: https://software.ecmwf.int/wiki/pages/viewpage.action?pageId=97384581

To do (PISM related):  ImonGre  mrroLi       IFS Surface runoff grib 128.8 but for EC_Earth-GrIs additional melt etc is included
To do (PISM related):  LImon    sftgif       This is the land ice mask and will be an extra variable in IFS (thomas: via PEXTRA?)
To do (PISM related):  LImon    sftgrf       For Greenland this is the same as above sftgif. We do not have Antarctic ice sheet.
??                     LImon    sncIs        Not available in IFS. Although it could be calculated from tile fractions and written out as extra output
treerink commented 5 years ago

I asked Andrea for the HTESSEL parameter values for ksat & rootdsl. Also asked whether specific expressions can be provided for the Efx table variables: wilt, fldcapacity, clayfrac, siltfrac, and sandfrac. The fall back workaround which was discussed in a previous mail conversation is that we take for each of these variables the following same expression: 66.128 + 67.128 + 27.128 + 28.128 + 43.128as a work around. Then the LS3MIP people can later decide which exact combination they want to use. However, we would appreciate very much that the most likely expressions will be provided, an eventual later fix can be always applied in a later stage.

treerink commented 5 years ago

CFmon tnt has been moved to the (pre) ignored lists (and taken off from the (pre) identified missing lists in dfd7caa.

treerink commented 5 years ago

I have made separate issues #444, #443, #439 and #442 for the remaining ones in this issues:

#444   To do:                 E3hr     prrc         Availlable in IFS: Precip. flux from convection liquid grib 128.107, this is a 3D field so the surface field has to be extracted from this. So only level 91 needs to be outputted.

#443   Unclear yet            Efx      wilt         Wilting point is defined in HTESSEL as a fraction (different for each soil texture in HTESSEL; variable slt grib code 43 - Table 128) of field capacity. (Alessandri A.: an Effective Wilting Point can be defined ... same as below discussion for the field capacity)
#443   Unclear yet            Efx      fldcapacity  Maximum amount of water the soil can hold (different for each soil texture in HTESSEL; variable slt grib code 43 - Table 128). It is constant map in time in H-TESSEL (Alessandri A.: an Effective Field Capacity can be defined for coupling analysis pourpouses that will change with the effective vegetation cover; as soon as we have daily effective vegetation cover in the outputs we can compute as well the Effective Field capacity).
#443   Unclear yet            Efx      clayfrac     soil texture (could be derived by soil HTESSEL soil types [1-6]: variable slt grib code 43 - Table 128): very fine (& fine?) in H-TESSEL boundary condition constant in time. It is available in a kind of mask: for each type a number 1-6 if this type at a lon,lat location the number corresponding with this type is in the mask. (No grib code we think.)
#443   Unclear yet            Efx      siltfrac     soil texture (could be derived by soil HTESSEL soil types [1-6]: variable slt grib code 43 - Table 128): medium (& medium fine?) in H-TESSEL boundary condition constant in time. It is available in a kind of mask: for each type a number 1-6 if this type at a lon,lat location the number corresponding with this type is in the mask. (No grib code we think.)
#443   Unclear yet            Efx      sandfrac     soil texture (could be derived by soil HTESSEL soil types [1-6]: variable slt grib code 43 - Table 128): coarse in H-TESSEL boundary condition constant in time. It is available in a kind of mask: for each type a number 1-6 if this type at a lon,lat location the number corresponding with this type is in the mask. (No grib code we think.)
#443   H-TESSEL parameter     Efx      ksat         H-TESSEL parameters inside the model that is prescribed for each dominant soil texture  (not in the output; no grib code associated).
#443   H-TESSEL parameter     Efx      rootdsl      H-TESSEL parameters inside the model that is prescribed for each dominant soil texture (not in the output, no grib code associated).

#439   To do #126             Efx      slthick      Available in IFS: independent of lon, lat. Top soil layer: 0-7 cm, Soil layer 2: 7-28 cm, Soil layer 3: 28-100 cm, Soil layer 4: 100-289 cm, added as issue #126

#442   To do (PISM related):  ImonGre  mrroLi       IFS Surface runoff grib 128.8 but for EC_Earth-GrIs additional melt etc is included
#442   To do (PISM related):  LImon    sftgif       This is the land ice mask and will be an extra variable in IFS (thomas: via PEXTRA?)
#442   To do (PISM related):  LImon    sftgrf       For Greenland this is the same as above sftgif. We do not have Antarctic ice sheet.
#442   ??                     LImon    sncIs        Not available in IFS. Although it could be calculated from tile fractions and written out as extra output

So closing this issue.