&ED_NL
!------------------------------------------------------------------------------------------!
! Initial configuration                                                                    !
!------------------------------------------------------------------------------------------!

   NL%EXPNME = 'ED version 2.1 test'   ! Simulation name

   NL%RUNTYPE  = 'INITIAL'   ! Runtype: either INITIAL or HISTORY
                             ! 'INITIAL' - Starts a new run, that can be based on a previous
                             !             run (restart/history), but then it will use.
                             !             the information as a simple initial condition
                             ! 'HISTORY' - Resumes a simulation from the last history.
                             !             This is different from initial in the sense that
                             !             exactly the same information written in the 
                             !             history will be used here.

!----- Start of simulation (RUNTYPE=initial only) -----------------------------------------!
   NL%IMONTHA  = 01        ! Month
   NL%IDATEA   = 01        ! Day
   NL%IYEARA   = 1993      ! Year
   NL%ITIMEA   = 0000      ! UTC

!----- End of simulation ------------------------------------------------------------------!

   NL%IMONTHZ  = 01       ! Month
   NL%IDATEZ   = 01       ! Day
   NL%IYEARZ   = 2004     ! Year
   NL%ITIMEZ   = 0000     ! UTC
!------------------------------------------------------------------------------------------!



!------------------------------------------------------------------------------------------!
! Main time step parameters                                                                !
!------------------------------------------------------------------------------------------!

   NL%DTLSM   = 600.     ! Basic time step [s].  In the Euler scheme, this
                         ! is the length of the Euler time step.  Recommended
                         ! values are between 60-180.  In the Runge Kutta
                         ! scheme, this is the time between photosynthesis
                         ! updates.  Recommended value is 900.

   NL%RADFRQ  = 600.     ! Frequency at which to update radiation. [s]
!------------------------------------------------------------------------------------------!


 
!------------------------------------------------------------------------------------------!
! Specify the grid cells for which you want to run ED.                                     !
!------------------------------------------------------------------------------------------!

!---- ED Region dimensions ----------------------------------------------------------------!
   NL%N_ED_REGION   =  0  ! number of regions for which you want to run ED. 
                          ! This can be zero provided that N_POI is not...
                          
   NL%GRID_TYPE     =  0  !  0 corresponds to a lat/lon grid.  
                          !  1 corresponds to the polar-stereographic

   !---------------------------------------------------------------------------------------!
   !   The following variables are used only if it is a regional run and using a lat/lon   !
   ! grid.                                                                                 !
   !---------------------------------------------------------------------------------------!
   NL%GRID_RES      =  1.0          ! This is the grid resolution scale in degrees. [\*/]
                                    !
   NL%ED_REG_LATMIN  =  39.0 ! List of minimum latitudes;
   NL%ED_REG_LATMAX  =  47.0 ! List of maximum latitudes;
   NL%ED_REG_LONMIN  = -76.0 ! List of minimum longitudes;
   NL%ED_REG_LONMAX  = -68.0! List of maximum latitudes;
   !---------------------------------------------------------------------------------------!


   !---------------------------------------------------------------------------------------!
   !     The following variables are used only if it is a regional run and using the       !
   ! polar-stereographic grid.                                                             !
   !---------------------------------------------------------------------------------------!
   NL%NNXP = 110                  ! Number of points in the X direction
   NL%NNYP = 70                   ! Number of points in the Y direction
                                  !
   NL%DELTAX   = 60000.           ! X and Y grid spacing scale
   NL%DELTAY   = 60000.           ! in metres (dimension 1)
                                  !
   NL%POLELAT  =  42.50           ! Latitude of pole point   (dimension 1)
   NL%POLELON  = -72.50           ! Longitude of pole point  (dimension 1)
                                  !
   NL%CENTLAT  =  42.50           ! Latitude of central point  (n_ed_regions)
   NL%CENTLON  = -72.50           ! Longitude of central point (n_ed_regions)
   !---------------------------------------------------------------------------------------!


   !---------------------------------------------------------------------------------------!
   !     The following variables are used by both grid types of regional run grid.         !
   !---------------------------------------------------------------------------------------!
   NL%NSTRATX = 1,4 ! Fraction of the sizes given above that I should use for these
   NL%NSTRATY = 1,4 ! regions. The variables affected are denoted by
   !---------------------------------------------------------------------------------------!


   !---------------------------------------------------------------------------------------!
   !     Instead of a regular grid, you can also define "polygons of interest"  (POIs),    !
   ! which are individual locations.                                                       !
   !---------------------------------------------------------------------------------------!
   NL%N_POI    =    1   ! Number of POIs.  This can be zero provided that N_ED_REGION
                        ! is not.
                
   NL%POI_LAT  =  42.50 ! list of the latitudes of the POIs (degrees north)
   NL%POI_LON  = -72.50 ! list of the longitudes of the POIs (degrees east)
   !---------------------------------------------------------------------------------------!


!------------------------------------------------------------------------------------------!



!------------------------------------------------------------------------------------------!
!     Load balancing method.  This is only used by regional runs, when RUNTYPE is set to   !
! HISTORY, or when IED_INIT_NODE is set to 4 or 5.                                         !
!------------------------------------------------------------------------------------------!
   NL%LOADMETH  = 0    ! Typical either 0 or 2, unecessary to set 1 or 3
                       ! 0 - Let the loadmaster determine
                       !        what is best (good for npoly > nmach )
                       ! 1 - Load (defined by the polygon vector "workload"
                       !     will be evenly distributed accross 
                       !        the nodes.
                       ! 2 - Try to load an equal number of polygons
                       !        per node (good for npoly = nmach)
                       ! 3 - Similar to 1, but here we round the index 
                       !        instead of finding the next integer
!------------------------------------------------------------------------------------------!



!------------------------------------------------------------------------------------------!
! Analysis/history file output                                                             !
!------------------------------------------------------------------------------------------!
 
   NL%IFOUTPUT  =  0            ! Instantaneous analysis (site average)
   NL%IDOUTPUT  =  0            ! Daily means (site average)
   NL%IMOUTPUT  =  3            ! Monthly means (site average)
   NL%IYOUTPUT  =  0            ! Annual output
   NL%ITOUTPUT  =  0            ! Instantaneous fluxes (site average)
                                ! single annual file, no patch/cohort, for optimization or 
                                !     flux tower comparision 
   NL%ISOUTPUT  =  3            ! History files
                                ! 0 - no file;
                                ! 1 - save ascii (deprecated);
                                ! 2 - save RAMS-vfm file (deprecated);
                                ! 3 - save HDF5;
   
   
   NL%ATTACH_METADATA = 0       ! Flag for attaching metadata to HDF datasets
                                ! Attaching metadata will aid new users
                                ! in quickly identifying dataset descriptions
                                ! but will compromise I/O performance significantly
                                ! 0 = no metadata, 1 = yes metadata

   NL%UNITFAST   = 0            ! Units for FRQFAST/OUTFAST
   NL%UNITSTATE  = 3            ! Units for FRQSTATE/OUTSTATE
                                ! 0 - Seconds;
                                ! 1 - Days;
                                ! 2 - Calendar months (variable)
                                ! 3 - Calendar years  (variable)
                                ! PS: If OUTFAST/OUTSTATE are set to special flags
                                !     (-1 or -2) unitfast/unitstate will be ignored for
                                !     them

   NL%OUTFAST    = 0            ! Time interval between analysis files
   NL%OUTSTATE   = 0            ! Time interval between history files
                                ! Values are in seconds and must be <= 1 month
                                ! Values < FRQFAST/FRQSTATE are set to those values
                                ! Special values: -1 = daily; -2 = monthly
                                ! If OUTFAST/STATE is greater than FRQFAST/STATE then
                                !    multiple timepoints can be recorded to a single file
                                !    reducing the number of files and i/o time in post-
                                !    processing.  Multiple timepoints should be avoided
                                !    in the history files if you intend to use these for
                                !    HISTORY runs as it hasn't been tested.

   NL%ICLOBBER  = 1             ! 0 = stop if files exist, 1 = overwite files
   NL%FRQFAST   = 21600.        ! Time interval between analysis/history files
   NL%FRQSTATE  = 1.            ! Time interval between history files
!------------------------------------------------------------------------------------------!



!------------------------------------------------------------------------------------------!
! Analysis/history/restart file: path, prefix, settings                                    !
!------------------------------------------------------------------------------------------!
   NL%FFILOUT = 'analy/harvard'       ! Analysis output prefix
   NL%SFILOUT = 'histo/harvard'       ! History output prefix

   !---------------------------------------------------------------------------------------!
   ! IED_INIT_MODE determines how the ecosystem should be initialized.                     !
   ! 0 - Start from near-bare ground (only a few seedling from every PFT included in this  !
   !     run).                                                                             !
   ! 1 - This will use history files from ED-1.0 to initialise the ecosystem state (like   !
   !     biomass, LAI, plant density, etc.), but starting the thermodynamic state as a     !
   !     new simulation.                                                                   !
   ! 2 - Same as 1, but using history files from ED-2.0 without multiple sites.            !
   ! 3 - Same as 1, but using history files from ED-2.0 with multiple sites.               !
   ! 4 - Same as 1, but using ED2.1 H5 history/state files that take the form:             !
   !     'dir/prefix-gxx.h5'                                                               !
   !     Initialization files MUST end with -gxx.h5 where xx is a two digit integer grid   !
   !     number.  Each grid has its own initialization file.  As an example, if a user has !
   !     two files to initialize their grids with:                                         !
   !     example_file_init-g01.h5 and example_file_init-g02.h5                             !
   !     NL%SFILIN = 'example_file_init'                                                   !
   !                                                                                       !
   ! 5 - This is similar to option 4, except that you may provide several files (including !
   !     a mix of regional and POI runs, each file ending at a different date).  This will !
   !     not check date nor grid structure, it will simply read all polygons and match the !
   !     nearest neighbour to each polygon of your run.  SFILIN must have the directory    !
   !     common to all history files that are sought to be used, up to the last character  !
   !     the files have in common.  For example if your files are                          !
   !     /mypath/P0001-S-2000-01-01-000000-g01.h5,                                         !
   !     /mypath/P0002-S-1966-01-01-000000-g02.h5,                                         !
   !     ...                                                                               !
   !     /mypath/P1000-S-1687-01-01-000000-g01.h5:                                         !
   !     NL%SFILIN = '/mypath/P'                                                           !
   !                                                                                       !
   ! 6 - Initialize with ED-2 style files without multiple sites, exactly like option 2,   !
   !     except that the PFT types are preserved.                                          !
   !                                                                                       !
   ! -1 - use a mix of ED-1 and ED-2 files. This should                                    !
   !      be avoided by all means, unless you are                                          !
   !      absolute sure of what you're doing;                                              !
   !---------------------------------------------------------------------------------------!
   NL%IED_INIT_MODE   = 2

   NL%EDRES  = 1.0                     ! ED2/ED1 restart grid resolution
                                       !    (Ignored unless ied_init_mode is 3.

   NL%SFILIN   = '../input_data/restart_files/tow.'

!----- Time of the history start ----------------------------------------------------------!
   NL%ITIMEH   = 0000     ! UTC
   NL%IDATEH   = 01       ! Day
   NL%IMONTHH  = 01       ! Month
   NL%IYEARH   = 1999     ! Year
!------------------------------------------------------------------------------------------!



!------------------------------------------------------------------------------------------!
!  Soil/Surface water variables                                                            !
!------------------------------------------------------------------------------------------!

   NL%NZG = 14          ! Number of soil layers
   NL%NZS = 4           ! Maximum number of snowcover layers

   NL%ISOILFLG = 2,2    ! Options are: 1. Read from LAT/LON
                        !              2. Constant default soil
                        ! This variable must be defined for all grids. The order is
                        ! all ED regions first followed by ED POIs. 
   
   NL%NSLCON   = 3      ! Default soil type if ISOILFLG = 2 (choices below:)
   !-----------------------------------------------------------------------!
   ! ED Soil classes:                                                      !
   !-----------------------------------------------------------------------!
   !   1 -- sand            |  2 -- loamy sand    |  3 -- sandy loam       !
   !   4 -- silt loam       |  5 -- loam          |  6 -- sandy clay loam  !
   !   7 -- silty clay loam |  8 -- clay loam     |  9 -- sandy clay       !
   !  10 -- silty clay      | 11 -- clay          | 12 -- peat             !
   !-----------------------------------------------------------------------!

   !---------------------------------------------------------------------------------------!
   !    Prescribed fraction of clay and sand.  This is used only when ISOIFLG is set to 2, !
   ! and slxclay and slxsand are defined between 0 and 1, with the sum of these variables  !
   ! not exceeding 1.  Otherwise, the default ED parameters will be used.                  !
   !---------------------------------------------------------------------------------------!
   NL%SLXCLAY = -1.0
   NL%SLXSAND = -1.0

   !----- Depth (m) of bottom of soil model grid levels -----------------------------------!
   NL%SLZ     = -4.50, -3.50, -2.75, -2.00, -1.50, -1.10, -0.80,
                -0.60, -0.45, -0.30, -0.20, -0.13, -0.08, -0.04
   !----- Initial soil moisture (fraction of saturation) ----------------------------------!
   NL%SLMSTR  =  0.60,  0.60,  0.60,  0.60,  0.60,  0.60,  0.60,
                 0.60,  0.60,  0.60,  0.60,  0.60,  0.60,  0.60
   !----- Initial soil temperature offset (with respect to the atmospheric temperature). --!
   NL%STGOFF  =   0.0,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0,
                  0.0,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0
!------------------------------------------------------------------------------------------!



!------------------------------------------------------------------------------------------!
!  Input databases                                                                         !
!------------------------------------------------------------------------------------------!

   NL%VEG_DATABASE    = '../input_data/veg_OGE2/OGE2_'
   NL%SOIL_DATABASE   = '../input_data/soil/type_FAO/FAO_'
   NL%ED_INPUTS_DIR   = '../input_data/chd+dgd+glu/'
                        ! Directory where chd, dgd and glu data are kept.

   !----- File containing information for meteorological drivers --------------------------!
   NL%ED_MET_DRIVER_DB = './HARVARD_MET'

   !---------------------------------------------------------------------------------------!
   ! SOILSTATE_DB - Soil state (temperature and moisture) database, used only when         !
   !                ISOILSTATEINIT is set to 1.                                            !
   ! SOILDEPTH_DB -  Soil depth database, used only when ISOILDEPTHFLG is set to 1.        !
   !---------------------------------------------------------------------------------------!
   NL%SOILSTATE_DB     = '../input_data/soil/depth/soil_depths.dat' 
   NL%SOILDEPTH_DB     = '../input_data/soil/tempmoist/STW1996OCT.dat'

   NL%ISOILSTATEINIT = 0  ! 0 -- use slmstr
                          ! 1 -- read from file (SOILSTATE_DB)

   NL%ISOILDEPTHFLG  = 0  ! 0 -- use NZG, SLZ for each grid cell.
                          ! 1 -- read file (SOILDEPTH_DB) to get soil depth for each grid 
                          !      cell.

   NL%ISOILBC        = 0  ! Soil boundary condition at lowest soil level
                          ! 0 - Bedrock;
                          ! 1 - Free drainage.
                          ! 2 - Half drainage.
!------------------------------------------------------------------------------------------!



!------------------------------------------------------------------------------------------!
! ED model options                                                                         !
!------------------------------------------------------------------------------------------!

   NL%INTEGRATION_SCHEME = 1  ! This determines how the fast time step
                              ! is integrated:
                              !  0  --  Euler step (not supported)
                              !  1  --  Fourth-order Runge Kutta

   NL%RK4_TOLERANCE  = 0.01   ! This is the relative tolerance for the Runge-Kutta 
                              !   integration.  Larger numbers will make runs go faster, 
                              !   at the expense of being less accurate.  Currently the 
                              !   valid range is between 1.e-7 and 1.e-1, but recommended
                              !   values are between 1.e-4 and 1.e-2

   NL%IBRANCH_THERMO = 0      ! This determines whether branches should be included in
                              !    the vegetation thermodynamics or not.  This is under
                              !    development, not fully tested.
                              !  0  --  No branches in energy/radiation;
                              !  1  --  Use Conijn (1995) numbers;
                              !  2  --  Use parameterization from Jarvela (2004);

   NL%ISTOMA_SCHEME   = 0     ! This is the scheme used to compute stomatal 
                              ! conductance and photosynthesis.  Set to 0 to do 
                              ! the full, exact calculation each time.  Set to 1 
                              ! to do the 'small perturbation' scheme.


   NL%IPHEN_SCHEME    = 1     ! This specifies the phenology scheme:
                              !  0 -- the cold-deciuous (Botta et al.)
                              !       and old drought-deciduous (same with ED2.0);
                              !  1 -- presecribed phenology based on observations
                              !  2 -- the cold-deciduous and new drought-deciduous;
                              !  3 -- the cold-deciduous, new drought-deciduous 
                              !       and light phenology (not fully tested);
  
   NL%REPRO_SCHEME    = 1     ! This specifies the reproduction scheme:
                              !  0 -- reproduction off (useful for short runs)
                              !  1 -- original ED1 default

   NL%LAPSE_SCHEME    = 0     ! This specifies the met lapse rate scheme:
                              !  0 -- lapse rates off
                              !  1 -- phenomenological, global

   NL%CROWN_MOD       = 0     ! Specifies how tree crowns are represent in the 
                              ! canopy radiation model
                              !  0  - ED1 default (continuous horizontal)
                              !  1  - finite radius mixing model (Dietze et al 2008)

   NL%DECOMP_SCHEME   = 0     ! Specifies soil decomposition temperature dependence
                              ! 0 -- original exponential
                              ! 1 -- Lloyd and Taylor (1994) model
                              !      [[requires parameters to be set in xml]]

   NL%N_PLANT_LIM     = 0     ! Determines whether (1) or not (0) plants can
                              !    be limited by nitrogen

   NL%N_DECOMP_LIM    = 0     ! Determines whether (1) or not (0) decomposition
                              !    can be limited by nitrogen

   NL%INCLUDE_FIRE    = 0     ! Determines how to apply fire in the model:
                              ! 0. No fires;
                              ! 1. (deprecated) Fire will be triggered with enough biomass 
                              !    and integrated ground water depth less than a threshold.
                              !    Based on ED-1, the threshold assumes that the soil is 
                              !    1 m, so deeper soils will need to be much drier
                              !    to allow fires to happen and often will never allow 
                              !    fires. 
                              ! 2. Fire will be triggered with enough biomass and the
                              !    total soil water at the top 75 cm falls below a 
                              !    threshold. 
                             

   NL%IANTH_DISTURB   = 0     ! Flag specifying whether (1) or not (0) to include
                              ! anthropogenic disturbances such as agriculture, 
                              ! abandonment, and harvesting.


   NL%ICANTURB        = 2     ! Canopy turbulence scheme
                              ! -1. (very deprecated) This is the original ED-2.0 
                              !     scheme, added here for backward compability.  
                              !     THIS SHOULD BE AVOIDED, since it gives strange 
                              !     results. 
                              !  0. is the default scheme
                              !  1. will rescale reference wind-speed if it
                              !     the reference height is (inapropriately)
                              !     below the canopy
                              !  2. (recommended) uses the method of Massman 
                              !     1997 and the bulk Richardson number of 
                              !     instability.  This method will not work 
                              !     when zref<h


   NL%ISFCLYRM        = 3     ! Surface layer model (used to compute ustar, tstar,...)
                              ! This variable controls the SLM for both LEAF and ED.
                              !  1 - BRAMS default, based on Louis (1979)
                              !      It uses empirical relations to estimate the flux
                              !      based on the bulk Richardson number
                              !  2 - Oncley and Dudhia (1995) model, based on MM5.
                              !      It uses a empirical relation between bulk Richardson
                              !      number and z/L, then computes the flux functions
                              !  3 - (recommended) Beljaars and Holtslag (1991) model.  
                              !      Similar to 2, but it finds z/L using a root-finding 
                              !      method, and it has a different set of equations for 
                              !      stable layers that works better during night time.
                              !  4 - The parameterization of Beljaars and Holtslag, but 
                              !      using z/L estimation as in Oncley and Dudhia. 


   !--------------------------------------------------------------------!
   ! PLANT FUNCTIONAL TYPES (PFTs):                                     !
   !--------------------------------------------------------------------!
   ! 1 - C4 grass                    |  7 - southern pines              !
   ! 2 - early tropical              |  8 - late conifers               !
   ! 3 - mid tropical                |  9 - early temperate deciduous   !
   ! 4 - late tropical               | 10 - mid temperate deciduous     !
   ! 5 - C3 grass                    | 11 - late temperate deciduous    !
   ! 6 - northern pines              |                                  !
   !--------------------------------------------------------------------!
   NL%INCLUDE_THESE_PFT = 5,6,8,9,10,11 ! List of PFTs to be included
   NL%AGRI_STOCK        = 5             ! Agriculture PFT (used only if ianth_disturb=1)
   NL%PLANTATION_STOCK  = 6             ! Plantation PFT  (used only if ianth_disturb=1)


   NL%PFT_1ST_CHECK = 0 ! What should I do when I find a cohort with an unexpected PFT 
                        !    in the history start?
                        ! 0. Stop the run;
                        ! 1. Add the PFT in the INCLUDE_THESE_PFT list;
                        ! 2. Ignore the cohort.
                        !
                        ! This option is ignored in near bare ground runs.

!----- Patch and cohort dimension settings ------------------------------------------------!
   NL%MAXPATCH  =  20 ! If number of patches in a given site exceeds NL%MAXPATCH, 
                      !    force patch fusion.  If NL%MAXPATCH is 0, then
                      !    fusion will never happen.  If MAXPATCH is negative, then the 
                      !    absolute value is used only during the initialization, and  
                      !    fusion will never happen again. 
   NL%MAXCOHORT = 120 ! If number of cohorts in a given patch exceeds NL%MAXCOHORT, 
                      !    force cohort fusion.  If NL%MAXCOHORT is 0, then
                      !    fusion will never happen.  If MAXCOHORT is negative, then the 
                      !    absolute value is used only during the initialization, and  
                      !    fusion will never happen again. 
                      ! 
!------------------------------------------------------------------------------------------!
                                         
                                         
                                         
!------------------------------------------------------------------------------------------!
! Commonly changed parameters                                                              !
!------------------------------------------------------------------------------------------!

   NL%ZROUGH = 0.1                       ! Constant roughness if for all domain

   NL%TREEFALL_DISTURBANCE_RATE  = 0.0   ! - Sign-dependent treefall disturbance rate:
                                         ! > 0. usual disturbance rate, in 1/years;
                                         ! = 0. No treefall disturbance;
                                         ! < 0. Treefall will be added as a mortality
                                         !      rate (it will kill plants, but it won't
                                         !      create a new patch).
                                         !  

   NL%RUNOFF_TIME    = 86400.0           ! the forest could get saturated and 
                                         ! develop a sfcwater pool of meters or 
                                         ! more.  In many cases, this is not realistic.
                                         ! This is the time scale in seconds for 
                                         ! this water to disappear.  If you want no runoff
                                         ! set runoff_time = 0 (not recommended for long 
                                         ! runs).
!------------------------------------------------------------------------------------------!




   !---------------------------------------------------------------------------------------!
   ! Control parameters for printing to standard output                                    !
   !         Any variable can be printed to standard output.                               !
   !         Choose any variables that are defined in the variable table fill              !
   !         routines in ed_state_vars.f90.  Choose the start and end dimension            !
   !         in the vector that you wish to view.  It can be polygon,site,..whatever       !
   !         This 'should' work in parallel ...                                            !
   !         The output will be printed to the screen at frequency frqanl,                 !
   !         the print call is synchronized with instantaneous file writes,                !
   !         so those variables will have been accumulated and normalized                  !
   !         from their rk4state and be meaningful to view.                                !
   !---------------------------------------------------------------------------------------!
   NL%IPRINTPOLYS  =  0        ! 0 = do not print information to screen
                               ! 1 = print polygon arrays to screen, use variables
                               !     described below

   !----- Print these variables ------------------------------------------!
   NL%PRINTVARS    = 'AVG_PCPG','AVG_CAN_TEMP','AVG_VAPOR_AC','AVG_CAN_SHV'

   !---------------------------------------------------------------------------------------!
   ! Use these format strings:                                                             !
   ! The format strings follow standard fortran format notation. i.e. 'f' is for floating  !
   ! point.  First number is total digits, second number is the number of digits to the    !
   ! left of right of the decimal point.                                                   !
   !---------------------------------------------------------------------------------------!
   NL%PFMTSTR      = 'f10.8','f5.1','f7.2','f9.5'

   NL%IPMIN        = 1       ! The first index to print
   NL%IPMAX        = 60      ! The last index to print
  !----------------------------------------------------------------------------------------!


!------------------------------------------------------------------------------------------!
! Met driver control                                                                       !
!------------------------------------------------------------------------------------------!
   NL%IMETTYPE  = 1          ! 0 = ASCII
                             ! 1 = HDF5  (recommended)

   NL%ISHUFFLE  = 0          ! How to choose an year outside the met data range:
                             ! 0. Sequentially cycle over years
                             ! 1. Randomly pick the years, using the same sequence
                             !    (not functional in every OS and compilers.  When it 
                             !     doesn't function, it will behave as option 2).
                             ! 2. Randomly pick the years, choosing a different 
                             !    sequence each time the model is run.

   NL%METCYC1         = 1993 ! First year of met data
   NL%METCYCF         = 2003 ! Last year of met data
   
   NL%INITIAL_CO2     = 370. ! Initial value for CO2 in case no CO2 is provided at the 
                             ! meteorological driver dataset


   !----- For phenology prescribed from observations, we also need to recycle over years. -!
   NL%IPHENYS1 = 1992        ! initial year for spring
   NL%IPHENYSF = 2003        ! final year for spring
   NL%IPHENYF1 = 1992        ! initial year for fall
   NL%IPHENYFF = 2003        ! final year for fall
   !----- Path and prefix for phenology for where to find presecribed phenology driver. ---!
   NL%PHENPATH = '../input_data/phenology/phenology'
   !---------------------------------------------------------------------------------------!



   !---------------------------------------------------------------------------------------!
   !     XML configuration file.  If these files are not found, the default settings will  !
   ! be used.                                                                              !
   !---------------------------------------------------------------------------------------!
   NL%IEDCNFGF   = 'config.xml'
   NL%EVENT_FILE = 'myevents.xml'
   !---------------------------------------------------------------------------------------!



   !---------------------------------------------------------------------------------------!
   ! Optimization configuration. (Currently not used)                                      !
   !---------------------------------------------------------------------------------------!
      NL%IOPTINPT = ''   ! Directory for optimizer inputs
   !---------------------------------------------------------------------------------------!
/ ! end of ED_NL
