diff --git a/src/MNH/default_desfmn.f90 b/src/MNH/default_desfmn.f90
index 0b4de6d6fe7162213d718f2f7598431b4221ab89..fa6a915d120212d58495e7b07a800aeedc389202 100644
--- a/src/MNH/default_desfmn.f90
+++ b/src/MNH/default_desfmn.f90
@@ -209,6 +209,7 @@ END MODULE MODI_DEFAULT_DESFM_n
 !!                   11/2019 C.Lac correction in the drag formula and application to building in addition to tree
 !  P. Wautelet 17/04/2020: move budgets switch values into modd_budget
 !  P. Wautelet 30/06/2020: add NNETURSV, NNEADVSV and NNECONSV variables
+!  PA. Joulin     10/2020: add default variables for wind tubrines (EOL)
 !  F. Auguste, T. Nagel 02/2021: add IBM defaults parameters
 !  T. Nagel       02/2021: add turbulence recycling defaults parameters
 !  P-A Joulin  21/05/2021: add Wind turbines
@@ -223,6 +224,7 @@ END MODULE MODI_DEFAULT_DESFM_n
 !  C. Barthe      03/2022: add CIBU and RDSF options in LIMA
 !  Delbeke/Vie    03/2022: KHKO option in LIMA
 !  P. Wautelet 27/04/2022: add namelist for profilers
+!  PA. Joulin     04/2023: add EOL/ADR
 !-------------------------------------------------------------------------------
 !
 !*       0.    DECLARATIONS
@@ -270,9 +272,9 @@ USE MODD_CONDSAMP
 USE MODD_MEAN_FIELD
 USE MODD_DRAGTREE_n
 USE MODD_DRAGBLDG_n
-USE MODD_EOL_MAIN
-USE MODD_EOL_ADNR
-USE MODD_EOL_ALM
+USE MODD_EOL_MAIN, ONLY: LMAIN_EOL,CMETH_EOL,CSMEAR,NMODEL_EOL
+USE MODD_EOL_ADR,  ONLY: NNB_RADELT, NNB_AZIELT
+USE MODD_EOL_ALM,  ONLY: NNB_BLAELT, LTIMESPLIT
 USE MODD_EOL_SHARED_IO
 USE MODD_ALLPROFILER_n
 USE MODD_ALLSTATION_n
@@ -592,14 +594,20 @@ CBLADE_CSVDATA    = 'data_blade.csv'
 CAIRFOIL_CSVDATA  = 'data_airfoil.csv'
 !
 CINTERP           = 'CLS'
+LTIPLOSSG         = .TRUE.
+LTECOUTPTS        = .FALSE.
+LCSVOUTFRM        = .FALSE.
+!
+!       10d.iii) MODD_EOL_ADR
+!
+NNB_RADELT        =  18
+NNB_AZIELT        =  56
 !
-!       10d.iii) MODD_EOL_ALM
+!       10d.iv) MODD_EOL_ALM
 !
 NNB_BLAELT        =  42
 LTIMESPLIT        = .FALSE.
-LTIPLOSSG         = .TRUE.
-LTECOUTPTS        = .FALSE.
-!
+
 !------------------------------------------------------------------------------
 !*      10.e   SET DEFAULT VALUES FOR MODD_ALLPROFILER_n :
 !             ----------------------------------
diff --git a/src/MNH/ini_mean_field.f90 b/src/MNH/ini_mean_field.f90
index 1abffffd053f5d6ead06c6979420e037acc58249..5334dc0d190a98c4a45a1e9d919479f84a76bff1 100644
--- a/src/MNH/ini_mean_field.f90
+++ b/src/MNH/ini_mean_field.f90
@@ -51,6 +51,7 @@ END MODULE MODI_INI_MEAN_FIELD
 !!                      02/2021 (T.Nagel) add passive scalar (XSVT) and UW wind component
 !!                      05/2021 (PA.Joulin) add wind turbine variables
 !!                      11/2022 (E. Jezequel) add covariances
+!!                      09/2022 (H.Toumi) add EOL/ADR variables
 !-------------------------------------------------------------------------------
 !
 !*       0.    DECLARATIONS
@@ -61,8 +62,9 @@ USE MODD_MEAN_FIELD_n
 USE MODD_MEAN_FIELD
 USE MODD_PARAM_n        
 USE MODD_EOL_MAIN, ONLY: LMAIN_EOL, CMETH_EOL, NMODEL_EOL
-USE MODD_EOL_SHARED_IO, ONLY: XTHRU_SUM, XTORQ_SUM, XPOW_SUM
-USE MODD_EOL_ALM
+USE MODD_EOL_SHARED_IO, ONLY: XTHRU_SUM, XTORQ_SUM, XPOW_SUM, XAOA_SUM
+USE MODD_EOL_ALM, ONLY: XFAERO_RE_SUM
+USE MODD_EOL_ADR, ONLY: XFAERO_RA_SUM
 USE MODE_MODELN_HANDLER
 !
 IMPLICIT NONE
@@ -98,6 +100,12 @@ IF (LMAIN_EOL .AND. IMI==NMODEL_EOL) THEN
  SELECT CASE(CMETH_EOL)
   CASE('ADNR') ! Actuator Disc Non-Rotating
    XTHRU_SUM      = 0.0
+  CASE('ADR') ! Actuator Disc Method
+   XAOA_SUM_ADR   = 0.0
+   XFAERO_RA_SUM  = 0.0
+   XTHRU_SUM      = 0.0
+   XTORQ_SUM      = 0.0
+   XPOW_SUM       = 0.0
   CASE('ALM') ! Actuator Line Method
    XAOA_SUM       = 0.0
    XFAERO_RE_SUM  = 0.0
diff --git a/src/MNH/ini_modeln.f90 b/src/MNH/ini_modeln.f90
index d1c00bae782862fd7ecee6d9dac255887aa27af4..bbcacd53bb99b43e0296d595a70a1f4bcbe6a25d 100644
--- a/src/MNH/ini_modeln.f90
+++ b/src/MNH/ini_modeln.f90
@@ -290,10 +290,12 @@ END MODULE MODI_INI_MODEL_n
 !  P. Wautelet 13/09/2019: budget: simplify and modernize date/time management
 !  C. Lac         11/2019: correction in the drag formula and application to building in addition to tree
 !  S. Riette      04/2020: XHL* fields
+!  PA. Joulin     10/2020: add wind tubrines (EOL)
 !  F. Auguste     02/2021: add IBM
 !  T.Nigel        02/2021: add turbulence recycling
 ! J.L.Redelsperger 06/2011: OCEAN case
 ! A. Costes       12/2021: Blaze fire model
+!  H. Toumi       09/2022: add EOL/ADR
 !---------------------------------------------------------------------------------
 !
 !*       0.    DECLARATIONS
@@ -439,6 +441,7 @@ USE MODI_INI_DRAG
 USE MODI_INI_DYNAMICS
 USE MODI_INI_ELEC_n
 USE MODI_INI_EOL_ADNR
+USE MODI_INI_EOL_ADR
 USE MODI_INI_EOL_ALM
 USE MODI_INI_LES_N
 USE MODI_INI_LG
@@ -2900,6 +2903,8 @@ IF (LMAIN_EOL .AND. KMI == NMODEL_EOL) THEN
  SELECT CASE(CMETH_EOL)
   CASE('ADNR')
    CALL INI_EOL_ADNR
+  CASE('ADR')
+   CALL INI_EOL_ADR(XDXX,XDYY)
   CASE('ALM')
    CALL INI_EOL_ALM(XDXX,XDYY)
  END SELECT
diff --git a/src/MNH/mean_field.f90 b/src/MNH/mean_field.f90
index a94d92633b762071be438dd28d4dd9c6f4bda493..41eb9349fa56e113ef2b2b8e56ff13d0518ac4d6 100644
--- a/src/MNH/mean_field.f90
+++ b/src/MNH/mean_field.f90
@@ -45,10 +45,11 @@ END MODULE MODI_MEAN_FIELD
 !!
 !!    MODIFICATIONS
 !!    -------------
-!!      Original    07/2009
-!!      (C.Lac)     09/2016 Max values
-!!      (PA.Joulin) 12/2020 Wind turbine variables
+!!      Original      07/2009
+!!      (C.Lac)       09/2016 Max values
+!!      (PA.Joulin)   12/2020 Wind turbine variables
 !!      (E. Jezequel) 11/2022 Welford algorithm and covariances
+!!      (H. Toumi)    009/2022: add ADR
 !!---------------------------------------------------------------
 !
 !
@@ -63,10 +64,10 @@ USE MODD_CST
 USE MODD_PASPOL
 !
 USE MODD_EOL_MAIN, ONLY: LMAIN_EOL, CMETH_EOL, NMODEL_EOL
-USE MODD_EOL_SHARED_IO, ONLY: XTHRUT, XTORQT, XPOWT
-USE MODD_EOL_SHARED_IO, ONLY: XTHRU_SUM, XTORQ_SUM, XPOW_SUM
-USE MODD_EOL_ALM
-USE MODD_EOL_ADNR
+USE MODD_EOL_SHARED_IO, ONLY: XTHRUT, XTORQT, XPOWT, XAOA_GLB
+USE MODD_EOL_SHARED_IO, ONLY: XTHRU_SUM, XTORQ_SUM, XPOW_SUM, XAOA_SUM
+USE MODD_EOL_ADR, ONLY: XFAERO_RA_SUM
+USE MODD_EOL_ALM, ONLY: XFAERO_RE_SUM
 USE MODE_MODELN_HANDLER
 USE MODI_UPDATE_WELFORD
 !
@@ -117,6 +118,12 @@ IKE=IKU-JPVEXT
     SELECT CASE(CMETH_EOL)
      CASE('ADNR') ! Actuator Disc Non-Rotating
       XTHRU_SUM       = XTHRUT        + XTHRU_SUM
+     CASE('ADR') ! Actuator Line Method
+      XAOA_SUM        = XAOA_GLB      + XAOA_SUM
+      XFAERO_RA_SUM   = XFAERO_RA_GLB + XFAERO_RA_SUM
+      XTHRU_SUM       = XTHRUT        + XTHRU_SUM
+      XTORQ_SUM       = XTORQT        + XTORQ_SUM
+      XPOW_SUM        = XPOWT         + XPOW_SUM
      CASE('ALM') ! Actuator Line Method
       XAOA_SUM        = XAOA_GLB      + XAOA_SUM
       XFAERO_RE_SUM   = XFAERO_RE_GLB + XFAERO_RE_SUM
diff --git a/src/MNH/read_exsegn.f90 b/src/MNH/read_exsegn.f90
index 70d8d7e98e669b1f0dd856a2789a64ff9491e676..9774fc4b9449248108d77ddcd125b765078f816b 100644
--- a/src/MNH/read_exsegn.f90
+++ b/src/MNH/read_exsegn.f90
@@ -295,6 +295,7 @@ END MODULE MODI_READ_EXSEG_n
 !  P. Wautelet 10/04/2019: replace ABORT and STOP calls by Print_msg
 !  C. Lac         11/2019: correction in the drag formula and application to building in addition to tree
 !  Q. Rodier      03/2020: add abort if use of any LHORELAX and cyclic conditions
+!  PA. Joulin     10/2020: add namelists for wind tubrines (EOL)
 !  F.Auguste      02/2021: add IBM
 !  T.Nagel        02/2021: add turbulence recycling
 !  E.Jezequel     02/2021: add stations read from CSV file
@@ -308,6 +309,8 @@ END MODULE MODI_READ_EXSEG_n
 !  P. Wautelet 24/06/2022: remove check on CSTORAGE_TYPE for restart of ForeFire variables
 !  P. Wautelet 13/07/2022: add namelist for flyers and balloons
 !  P. Wautelet 19/08/2022: add namelist for aircrafts
+!  H. Toumi       09/2022: add EOL/ADR
+!  PA. Joulin     04/2023: update for main
 !------------------------------------------------------------------------------
 !
 !*       0.    DECLARATIONS
@@ -369,6 +372,7 @@ USE MODN_DYN_n      ! to avoid the duplication of this routine for each model.
 USE MODN_ELEC
 USE MODN_EOL
 USE MODN_EOL_ADNR
+USE MODN_EOL_ADR
 USE MODN_EOL_ALM
 USE MODN_FIRE_n
 USE MODN_FLYERS
@@ -565,6 +569,8 @@ CALL POSNAM(ILUSEG,'NAM_EOL',GFOUND,ILUOUT)
 IF (GFOUND) READ(UNIT=ILUSEG,NML=NAM_EOL)
 CALL POSNAM(ILUSEG,'NAM_EOL_ADNR',GFOUND,ILUOUT)
 IF (GFOUND) READ(UNIT=ILUSEG,NML=NAM_EOL_ADNR)
+CALL POSNAM(ILUSEG,'NAM_EOL_ADR',GFOUND,ILUOUT)
+IF (GFOUND) READ(UNIT=ILUSEG,NML=NAM_EOL_ADR)
 CALL POSNAM(ILUSEG,'NAM_EOL_ALM',GFOUND,ILUOUT)
 IF (GFOUND) READ(UNIT=ILUSEG,NML=NAM_EOL_ALM)
 CALL POSNAM(ILUSEG,'NAM_PROFILERN',GFOUND,ILUOUT)
diff --git a/src/MNH/write_lfin.f90 b/src/MNH/write_lfin.f90
index 356692550d70a74414763532bb08f5314cfa10c1..052674bad5762f38dfebc36c73dc3a52bbebb906 100644
--- a/src/MNH/write_lfin.f90
+++ b/src/MNH/write_lfin.f90
@@ -181,6 +181,8 @@ END MODULE MODI_WRITE_LFIFM_n
 !  A. Costes      12/2021: add Blaze fire model
 !  P. Wautelet 04/02/2022: use TSVLIST to manage metadata of scalar variables
 !  E. Jezequel    11/2022: add covariances from MEAN fields
+!  H. Toumi       09/2022: add ADR
+!  PA. Joulin     04/2023: update EOL metadata management
 !-------------------------------------------------------------------------------
 !
 !*       0.    DECLARATIONS
@@ -213,6 +215,7 @@ USE MODD_ELEC_DESCR,      ONLY: LLNOX_EXPLICIT
 USE MODD_ELEC_FLASH
 USE MODD_ELEC_n
 USE MODD_EOL_ADNR
+USE MODD_EOL_ADR
 USE MODD_EOL_ALM
 USE MODD_EOL_MAIN
 USE MODD_EOL_SHARED_IO
@@ -2057,18 +2060,20 @@ SELECT CASE(CMETH_EOL)
       CALL IO_Field_write(TPFILE,TZFIELD,XTHRU_SUM/MEAN_COUNT)
 !
     END IF
-!             iii) Actuator Line Model
 !
-  CASE('ALM') ! Actuator Line Method
+!             iii) Actuator Disc with Rotation Model
+!
+  CASE('ADR') ! Actuator Disc with Rotation
 !
+! * 1D Variables (rotor id)
     TZFIELD = TFIELDMETADATA(                   &
-      CMNHNAME   = 'generic for ALM variables', & !Temporary name to ease identification
+      CMNHNAME   = '1D ADR variables: (rot)',   & !Temporary name to ease identification
       CSTDNAME   = '',                          &
       CDIR       = '--',                        &
       NGRID      = 1,                           &
       NTYPE      = TYPEREAL,                    &
       NDIMS      = 1,                           &
-      LTIMEDEP   = .TRUE.                       )
+      LTIMEDEP   = .TRUE.   
 !
     TZFIELD%CMNHNAME   = 'THRUT'
     TZFIELD%CLONGNAME  = 'THRUSTT_EOL'
@@ -2088,14 +2093,184 @@ SELECT CASE(CMETH_EOL)
     TZFIELD%CCOMMENT   = 'RID instantaneous power (W) of wind turbines'
     CALL IO_Field_write(TPFILE,TZFIELD,XPOWT)
 !
+!
+! * 3D Variables (rotor id, azimuthal id, radial id)
+    TZFIELD = TFIELDMETADATA(                           &
+      CMNHNAME   = '3D ADR variables : (rot,azi,rad)',  & 
+      CSTDNAME   = '',                                  &
+      CDIR       = '--',                                &
+      NGRID      = 1,                                   &
+      NTYPE      = TYPEREAL,                            &
+      NDIMS      = 3,                                   &
+      LTIMEDEP   = .TRUE.                               )
+!
+    TZFIELD%CMNHNAME   = 'ELT_RAD'
+    TZFIELD%CLONGNAME  = 'ELT_RAD'
+    TZFIELD%CUNITS     = 'm'
+    TZFIELD%CCOMMENT   = 'RID_AZI_RAD radius (m) of wind turbine blade elements'
+    CALL IO_Field_write(TPFILE,TZFIELD,XELT_RAD)
+!
+    TZFIELD%CMNHNAME   = 'ELT_AZI'
+    TZFIELD%CLONGNAME  = 'ELT_AZI'
+    TZFIELD%CUNITS     = 'rad'
+    TZFIELD%CCOMMENT   = 'RID_AZI_RAD Azimutal angle (rad) of wind turbine blade elements'
+    CALL IO_Field_write(TPFILE,TZFIELD,XELT_AZI)
+!
+    TZFIELD%CMNHNAME   = 'AOA'
+    TZFIELD%CLONGNAME  = 'ANGLE OF ATTACK'
+    TZFIELD%CUNITS     = 'rad'
+    TZFIELD%CCOMMENT   = 'RID_AZI_RAD instantaneous angle of attack (rad)'
+    CALL IO_Field_write(TPFILE,TZFIELD,XAOA_GLB)
+!
+    TZFIELD%CMNHNAME   = 'FLIFT'
+    TZFIELD%CLONGNAME  = 'LIFT FORCE'
+    TZFIELD%CUNITS     = 'N'
+    TZFIELD%CCOMMENT   = 'RID_AZI_RAD instantaneous lift (N) in relative frame'
+    CALL IO_Field_write(TPFILE,TZFIELD,XFLIFT_GLB)
+!
+    TZFIELD%CMNHNAME   = 'FDRAG'
+    TZFIELD%CLONGNAME  = 'DRAG FORCE'
+    TZFIELD%CUNITS     = 'N'
+    TZFIELD%CCOMMENT   = 'RID_AZI_RAD instantaneous drag (N) in relative frame'
+    CALL IO_Field_write(TPFILE,TZFIELD,XFDRAG_GLB)
+!
+! * 4D Variables (rotor id, azimuthal id, radial id, xyz)
+    TZFIELD = TFIELDMETADATA(                               &
+      CMNHNAME   = '4D ADR variables : (rot,azi,rad,xyz)',  & 
+      CSTDNAME   = '',                                      &
+      CDIR       = '--',                                    &
+      NGRID      = 1,                                       &
+      NTYPE      = TYPEREAL,                                &
+      NDIMS      = 4,                                       &
+      LTIMEDEP   = .TRUE.                                   )
+!
+    TZFIELD%CMNHNAME   = 'FAERO_RA'
+    TZFIELD%CLONGNAME  = 'AERODYNAMIC FORCE RA'
+    TZFIELD%CUNITS     = 'N'
+    TZFIELD%CCOMMENT   = 'RID_AZI_RAD_XYZ instantaneous forces (N) in RA'
+    CALL IO_Field_write(TPFILE,TZFIELD,XFAERO_RA_GLB)
+!
+    TZFIELD%CMNHNAME   = 'FAERO_RG'
+    TZFIELD%CLONGNAME  = 'AERODYNAMIC FORCE RG'
+    TZFIELD%CUNITS     = 'N'
+    TZFIELD%CCOMMENT   = 'RID_AZI_RAD_XYZ instantaneous forces (N) in RG'
+    CALL IO_Field_write(TPFILE,TZFIELD,XFAERO_RG_GLB)
+!
+! * Blade Equivalent Variables (rotor id, radial id, xyz)
+    TZFIELD = TFIELDMETADATA(                           &
+      CMNHNAME   = '3D ADR variables : (rot,rad,xyz)',  & 
+      CSTDNAME   = '',                                  &
+      CDIR       = '--',                                &
+      NGRID      = 1,                                   &
+      NTYPE      = TYPEREAL,                            &
+      NDIMS      = 3,                                   &
+      LTIMEDEP   = .TRUE.                               )
+!
+    TZFIELD%CMNHNAME   = 'FAERO_BLA_RA'
+    TZFIELD%CLONGNAME  = 'AERODYNAMIC FORCE ON ONE BLADE IN RA'
+    TZFIELD%CUNITS     = 'N'
+    TZFIELD%CCOMMENT   = 'RID_RAD_XYZ forces (N) in RA'
+    CALL IO_Field_write(TPFILE,TZFIELD,XFAERO_RA_BLA)
+!
+    IF (MEAN_COUNT /= 0) THEN
+!
+! * 1D Variables (rotor id)
+      TZFIELD = TFIELDMETADATA(                        &
+        CMNHNAME   = '1D ADR mean variables: (rot)',   & !Temporary name to ease identification
+        CSTDNAME   = '',                               &
+        CDIR       = '--',                             &
+        NGRID      = 1,                                &
+        NTYPE      = TYPEREAL,                         &
+        NDIMS      = 1,                                &
+        LTIMEDEP   = .TRUE.   
+!
+      TZFIELD%CMNHNAME   = 'THRUMME'
+      TZFIELD%CLONGNAME  = 'MEAN_THRUST_EOL'
+      TZFIELD%CUNITS     = 'N'
+      TZFIELD%CCOMMENT   = 'RID mean thrust of the wind turbines (N)'
+      CALL IO_Field_write(TPFILE,TZFIELD,XTHRU_SUM/MEAN_COUNT)
+!
+      TZFIELD%CMNHNAME   = 'TORQMME'
+      TZFIELD%CLONGNAME  = 'MEAN_TORQUE_EOL'
+      TZFIELD%CUNITS     = 'Nm'
+      TZFIELD%CCOMMENT   = 'RID mean torque of the wind turbines (Nm)'
+      CALL IO_Field_write(TPFILE,TZFIELD,XTORQ_SUM/MEAN_COUNT)
+!
+      TZFIELD%CMNHNAME   = 'POWMME'
+      TZFIELD%CLONGNAME  = 'MEAN_POWER_EOL'
+      TZFIELD%CUNITS     = 'W'
+      TZFIELD%CCOMMENT   = 'RID mean power of the wind turbines (W)'
+      CALL IO_Field_write(TPFILE,TZFIELD,XPOW_SUM/MEAN_COUNT)
+!
+      TZFIELD = TFIELDMETADATA(                                &
+        CMNHNAME   = 'AOAMME',                                 &
+        CSTDNAME   = '',                                       &
+        CLONGNAME  = 'MEAN_ANGLE_OF_ATTACK',                   &
+        CUNITS     = 'rad',                                    &
+        CDIR       = '--',                                     &
+        CCOMMENT   = 'RID_AZI_RAD mean angle of attack (rad)', &
+        NGRID      = 1,                                        &
+        NTYPE      = TYPEREAL,                                 &
+        NDIMS      = 3,                                        &
+        LTIMEDEP   = .TRUE.                                    )
+      CALL IO_Field_write(TPFILE,TZFIELD,XAOA_SUM/MEAN_COUNT)
+!
+      TZFIELD = TFIELDMETADATA(                               &
+        CMNHNAME   = 'FAEROMME_RE',                           &
+        CSTDNAME   = '',                                      &
+        CLONGNAME  = 'MEAN_AERODYNAMIC_FORCE_RA',             &
+        CUNITS     = 'N',                                     &
+        CDIR       = '--',                                    &
+        CCOMMENT   = 'RID_AZI_RAD_XYZ mean forces (N) in RA', &
+        NGRID      = 1,                                       &
+        NTYPE      = TYPEREAL,                                &
+        NDIMS      = 4,                                       &
+        LTIMEDEP   = .TRUE.                                   )
+      CALL IO_Field_write(TPFILE,TZFIELD,XFAERO_RA_SUM/MEAN_COUNT)
+!
+    END IF
+!
+!             iv) Actuator Line Model
+!
+  CASE('ALM') ! Actuator Line Method
+!
+! * 1D Variables (rotor id)
     TZFIELD = TFIELDMETADATA(                   &
-      CMNHNAME   = 'generic for ALM variables', & !Temporary name to ease identification
+      CMNHNAME   = '1D ALM variables: (rot)',   & !Temporary name to ease identification
       CSTDNAME   = '',                          &
       CDIR       = '--',                        &
       NGRID      = 1,                           &
       NTYPE      = TYPEREAL,                    &
-      NDIMS      = 3,                           &
+      NDIMS      = 1,                           &
       LTIMEDEP   = .TRUE.                       )
+!
+    TZFIELD%CMNHNAME   = 'THRUT'
+    TZFIELD%CLONGNAME  = 'THRUSTT_EOL'
+    TZFIELD%CUNITS     = 'N'
+    TZFIELD%CCOMMENT   = 'RID instantaneous thrust (N) of wind turbines'
+    CALL IO_Field_write(TPFILE,TZFIELD,XTHRUT)
+!
+    TZFIELD%CMNHNAME   = 'TORQT'
+    TZFIELD%CLONGNAME  = 'TORQUET_EOL'
+    TZFIELD%CUNITS     = 'Nm'
+    TZFIELD%CCOMMENT   = 'RID instantaneous torque (Nm) of wind turbines'
+    CALL IO_Field_write(TPFILE,TZFIELD,XTORQT)
+!
+    TZFIELD%CMNHNAME   = 'POWT'
+    TZFIELD%CLONGNAME  = 'POWERT_EOL'
+    TZFIELD%CUNITS     = 'W'
+    TZFIELD%CCOMMENT   = 'RID instantaneous power (W) of wind turbines'
+    CALL IO_Field_write(TPFILE,TZFIELD,XPOWT)
+!
+! * 3D Variables (rotor id, blade id, radial id)
+    TZFIELD = TFIELDMETADATA(                             &
+      CMNHNAME   = '3D ALM variables : (rot,blade,rad)',  & 
+      CSTDNAME   = '',                                    &
+      CDIR       = '--',                                  &
+      NGRID      = 1,                                     &
+      NTYPE      = TYPEREAL,                              &
+      NDIMS      = 3,                                     &
+      LTIMEDEP   = .TRUE.                                 )
 !
     TZFIELD%CMNHNAME   = 'ELT_RAD'
     TZFIELD%CLONGNAME  = 'ELT_RAD'
@@ -2121,14 +2296,15 @@ SELECT CASE(CMETH_EOL)
     TZFIELD%CCOMMENT   = 'RID_BID_EID instantaneous drag (N) in relative frame'
     CALL IO_Field_write(TPFILE,TZFIELD,XFDRAG_GLB)
 !
-    TZFIELD = TFIELDMETADATA(                   &
-      CMNHNAME   = 'generic for ALM variables', & !Temporary name to ease identification
-      CSTDNAME   = '',                          &
-      CDIR       = '--',                        &
-      NGRID      = 1,                           &
-      NTYPE      = TYPEREAL,                    &
-      NDIMS      = 4,                           &
-      LTIMEDEP   = .TRUE.                       )
+! * 4D Variables (rotor id, azimuthal id, radial id, xyz)
+    TZFIELD = TFIELDMETADATA(                                 &
+      CMNHNAME   = '4D ALM variables : (rot,blade,rad,xyz)',  & 
+      CSTDNAME   = '',                                        &
+      CDIR       = '--',                                      &
+      NGRID      = 1,                                         &
+      NTYPE      = TYPEREAL,                                  &
+      NDIMS      = 4,                                         &
+      LTIMEDEP   = .TRUE.                                     )
 !
     TZFIELD%CMNHNAME   = 'FAERO_RE'
     TZFIELD%CLONGNAME  = 'AERODYNAMIC FORCE RE'
@@ -2144,8 +2320,9 @@ SELECT CASE(CMETH_EOL)
 !
     IF (MEAN_COUNT /= 0) THEN
 !
+! * 1D Variables (rotor id)
       TZFIELD = TFIELDMETADATA(                        &
-        CMNHNAME   = 'generic for ALM mean variables', & !Temporary name to ease identification
+        CMNHNAME   = '1D ALM mean variables: (rot)',   & !Temporary name to ease identification
         CSTDNAME   = '',                               &
         CDIR       = '--',                             &
         NGRID      = 1,                                &