diff --git a/lis/arch/Config.pl b/lis/arch/Config.pl
index c26605256..1d954cc1e 100755
--- a/lis/arch/Config.pl
+++ b/lis/arch/Config.pl
@@ -165,9 +165,11 @@
        $sys_c_opt .= " -Wunused-function -Wunused-parameter";
        $sys_c_opt .= " -Wunused-variable -Wwrite-strings";
        # Run-time flags
-       $sys_c_opt .= " -check=conversions,stack,uninit";
+       #EMK 20231109...Disabled several flags that are rejected by the new ICX
+       #compiler on Narwhal.
+       #$sys_c_opt .= " -check=conversions,stack,uninit";
        $sys_c_opt .= " -fp-stack-check -fp-trap=common -fp-trap-all=common";
-       $sys_c_opt .= " -ftrapuv";
+       #$sys_c_opt .= " -ftrapuv";
    }
    elsif($sys_arch eq "linux_pgi") {
       print "Optimization level $opt_lev is not defined for $sys_arch.\n";
diff --git a/lis/surfacemodels/land/rdhm.3.5.6/frz/gen_frzst.c b/lis/surfacemodels/land/rdhm.3.5.6/frz/gen_frzst.c
index ca6f22a08..295d8ead3 100755
--- a/lis/surfacemodels/land/rdhm.3.5.6/frz/gen_frzst.c
+++ b/lis/surfacemodels/land/rdhm.3.5.6/frz/gen_frzst.c
@@ -6,6 +6,10 @@
 #include "linux.h"
 /* #include "models.h" */
 
+extern void fstfg1_(int *, int *, int *, int *, float *, float *, float *, float *, float *, float *, float *, float *, float *, float *, float *, float *, float *);
+
+extern void frzind1_(float *, float *, float *, float *, int *, float *, float *, float *);
+
 int gen_frzst(int npix, int *nsoil, int *nupl, float *sac_parx, 
                float *frz_parx, float *frz_stx, float *sac_stx,
                float *sacst_prv, float *smc, float *sh2o, float dthr,
diff --git a/lis/surfacemodels/land/rdhm.3.5.6/frz/scale_frzst.c b/lis/surfacemodels/land/rdhm.3.5.6/frz/scale_frzst.c
index a3b7c3062..a15c23eb8 100755
--- a/lis/surfacemodels/land/rdhm.3.5.6/frz/scale_frzst.c
+++ b/lis/surfacemodels/land/rdhm.3.5.6/frz/scale_frzst.c
@@ -4,6 +4,14 @@
 #include "com_header.h"
 #include "linux.h"
 
+extern void fstfg1_(int *, int *, int *, int *, float *, float *, float *, float *, float *, float *, float *, float *, float *, float *, float *, float *, float *);
+
+extern void soil_int1_(float *,int *,float *,float *,int *,float *);
+
+extern void scale_sh2o_(float *, float *, float *, int *, int *, float *, float *, float *, float *, float *, float *, float *);
+
+extern void frzind1_(float *, float *, float *, float *, int *, float *, float *, float *);
+
 int scale_frzst(int npix, int *nsoil, int *nupl, float *sac_parx, 
                float *frz_parx, float *frz_stx, float *sac_stx,
                float *sacst_prv, float *smc, float *sh2o, float dthr,
diff --git a/lis/surfacemodels/land/rdhm.3.5.6/rdhm_get_ped.c b/lis/surfacemodels/land/rdhm.3.5.6/rdhm_get_ped.c
index b0a3358b1..ad97bf88e 100755
--- a/lis/surfacemodels/land/rdhm.3.5.6/rdhm_get_ped.c
+++ b/lis/surfacemodels/land/rdhm.3.5.6/rdhm_get_ped.c
@@ -9,6 +9,9 @@
 //-------------------------END NOTICE -- DO NOT EDIT-----------------------
  /* wrapper for calling get_ped 
     Shugong Wang  10/25/2013  */
+
+extern void get_ped(int, int, int, int, int, float *, float *, float *);
+
  void rdhm_get_ped_(int *year,            /* */
                     int *month,           /* */
                     int *day,             /* */
diff --git a/lis/surfacemodels/land/rdhm.3.5.6/sac/days_of_month.c b/lis/surfacemodels/land/rdhm.3.5.6/sac/days_of_month.c
index cef3f38c3..b109acf20 100755
--- a/lis/surfacemodels/land/rdhm.3.5.6/sac/days_of_month.c
+++ b/lis/surfacemodels/land/rdhm.3.5.6/sac/days_of_month.c
@@ -13,6 +13,8 @@
 #include "com_header.h"
 #include "models.h"
 
+extern int is_leap_year(int);
+
 int days_of_month(int m, int y)
 {
   int num_days[12]={31,28,31,30,31,30,31,31,30,31,30,31};
diff --git a/lis/surfacemodels/land/rdhm.3.5.6/sac/get_daily_pe.c b/lis/surfacemodels/land/rdhm.3.5.6/sac/get_daily_pe.c
index e6f2aa2c1..d30f3a0fa 100755
--- a/lis/surfacemodels/land/rdhm.3.5.6/sac/get_daily_pe.c
+++ b/lis/surfacemodels/land/rdhm.3.5.6/sac/get_daily_pe.c
@@ -15,6 +15,8 @@
 
 #include "models.h"
 
+extern int days_of_month(int, int);
+
 void get_daily_pe(int yr, int mon, int day, int npix,
 		  float *pe, float *pe_adj, float *pe_day)
 {
diff --git a/lis/surfacemodels/land/rdhm.3.5.6/sac/get_ped.c b/lis/surfacemodels/land/rdhm.3.5.6/sac/get_ped.c
index 65dca8f2b..80af81203 100755
--- a/lis/surfacemodels/land/rdhm.3.5.6/sac/get_ped.c
+++ b/lis/surfacemodels/land/rdhm.3.5.6/sac/get_ped.c
@@ -16,6 +16,8 @@
 #include "com_header.h"
 #include "models.h"
 
+extern void get_daily_pe(int, int, int, int, float *, float *, float *);
+
 void get_ped(int y, int m, int d, int npix, int dtm, 
 	     float *pe, float *pe_adj, float *ped)
 {
diff --git a/lis/surfacemodels/land/vic.4.1.1/vic411_read_vegparam.c b/lis/surfacemodels/land/vic.4.1.1/physics/vic411_read_vegparam.c
similarity index 76%
rename from lis/surfacemodels/land/vic.4.1.1/vic411_read_vegparam.c
rename to lis/surfacemodels/land/vic.4.1.1/physics/vic411_read_vegparam.c
index c89ba26df..a87371872 100644
--- a/lis/surfacemodels/land/vic.4.1.1/vic411_read_vegparam.c
+++ b/lis/surfacemodels/land/vic.4.1.1/physics/vic411_read_vegparam.c
@@ -1,25 +1,13 @@
-//-----------------------BEGIN NOTICE -- DO NOT EDIT-----------------------
-// NASA Goddard Space Flight Center
-// Land Information System Framework (LISF)
-// Version 7.4
-//
-// Copyright (c) 2022 United States Government as represented by the
-// Administrator of the National Aeronautics and Space Administration.
-// All Rights Reserved.
-//-------------------------END NOTICE -- DO NOT EDIT-----------------------
-
 #include <stdio.h>
 #include <stdlib.h>
-#include "vic411_vicNl.h" 
+#include <vic411_vicNl.h>
 #include <string.h>
 
 static char vcid[] = "$Id: vic411_read_vegparam.c,v 4.5.2.11 2009/10/01 21:12:16 vicadmin Exp $";
-// shugong: if cannot find entry for the gridid+vegetation class, return 0, other return 1
-int vic411_read_vegparam(FILE *vegparam,
-                         int   gridcel,
-                         int   Nveg_type,
-                         int   lis_veg_class,
-                         vic411_veg_con_struct *lis_temp)
+
+vic411_veg_con_struct *vic411_read_vegparam(FILE *vegparam,
+                              int   gridcel,
+                              int   Nveg_type)
 /**********************************************************************
   vic411_read_vegparam.c    Keith Cherkauer and Dag Lohmann       1997
 
@@ -58,7 +46,7 @@ int vic411_read_vegparam(FILE *vegparam,
 **********************************************************************/
 {
 
-  void vic411_ttrim1( char *string );
+  void vic411_ttrim( char *string );
   extern vic411_veg_lib_struct *vic411_veg_lib;
   extern vic411_option_struct   vic411_options;
 #if LINK_DEBUG
@@ -80,18 +68,6 @@ int vic411_read_vegparam(FILE *vegparam,
   char            *token;
   char            *vegarr[500];
   size_t	  length;
-  int vic411_flag;
-  vic411_flag = 0; // defualt: not found  
-  // added by Shugong Wang.
-  /* vic_veg_clas is defined to save the vegetation class in vegetation parameter file (called real vegetation class after). 
-     temp is a vic411_veg_con_struct. Firstly, the real vegetation class will be read in and save in the veg_class of temp. Then, 
-     vic411_read_vegparam will find the index number of the real vegetation class in vegetation parameter library (vic411_veg_lib). 
-     In the end, the veg_class of temp will be replced with the index number. In other word, the vegetation class in vic411_veg_con_struct
-     will not be the real vegetatoin class. In stead, it will be the corresponding index number in vegetatoin parameter library.
-     
-     In this function (vic411_read_vegparam), lis_veg_class is real vegetation class. In order to update the LAI and related values in
-     vic411_veg_lib correctly, lis_veg_class should be compared with the real vegetation class (vic_veg_class) */
-  int vic_veg_class;
 
   if(vic411_options.GLOBAL_LAI) skip=2;
   else skip=1;
@@ -121,7 +97,7 @@ int vic411_read_vegparam(FILE *vegparam,
     exit(99);
   }
   if(vegetat_type_num >= MAX_VEG) {
-    sprintf(ErrStr,"Vegetation parameter file wants more vegetation tiles in grid cell %i (%i) than are allowed by MAX_VEG (%i) [NOTE: bare soil class is assumed].  Edit vicNl_def.h and recompile.",gridcel,vegetat_type_num+1,MAX_VEG);
+    sprintf(ErrStr,"Vegetation parameter file wants more vegetation tiles in grid cell %i (%i) than are allowed by MAX_VEG (%i) [NOTE: bare soil class is assumed].  Edit vic411_vicNl_def.h and recompile.",gridcel,vegetat_type_num+1,MAX_VEG);
     vic411_nrerror(ErrStr);
   }
 
@@ -146,7 +122,7 @@ int vic411_read_vegparam(FILE *vegparam,
       vic411_nrerror(ErrStr);
     }
     strcpy(tmpline, line);
-    vic411_ttrim1( tmpline );
+    vic411_ttrim( tmpline );
     token = strtok (tmpline, delimiters);    /*  token => veg_class, move 'line' pointer to next field */  
     Nfields = 0;
     vegarr[Nfields] = calloc( 500, sizeof(char));
@@ -174,20 +150,15 @@ int vic411_read_vegparam(FILE *vegparam,
 
     temp[i].LAKE = 0;
     temp[i].veg_class = atoi( vegarr[0] );
-    /* added by Shugong Wang */
-    vic_veg_class = temp[i].veg_class; // this is the real vegetation class in vegetation parameter file. 
-
     temp[i].Cv = atof( vegarr[1] );
     depth_sum = 0;
     sum = 0.;
-    /* the following VIC code read root zoon information from vegetation parameter file*/
     for(j=0;j<vic411_options.ROOT_ZONES;j++) {
       temp[i].zone_depth[j] = atof( vegarr[2 + j*2] );
       temp[i].zone_fract[j] = atof( vegarr[3 + j*2] );
       depth_sum += temp[i].zone_depth[j];
       sum += temp[i].zone_fract[j];
     }
-
     if(depth_sum <= 0) {
       sprintf(str,"Root zone depths must sum to a value greater than 0.");
       vic411_nrerror(str);
@@ -199,21 +170,6 @@ int vic411_read_vegparam(FILE *vegparam,
       }
     }
 
-    /* added by Shugong Wang to read root zone information for LIS-tiling. 
-       if current vic_veg_class is the same as the lis_veg_class, then set the root zone
-       information to lis_temp, which points to vic411_lis_veg_con.                            */
-    if(vic_veg_class == lis_veg_class)
-    {
-         // added by Shugong Wang. lis_temp is a pointer of vic411_lis_veg_con. In the following lines, lis_temp will retrieve
-         // the root zone depth and root zone fraction from VIC vegetation parameter file
-         for(j=0; j<vic411_options.ROOT_ZONES; j++)
-         {
-           lis_temp->zone_depth[j] = temp[i].zone_depth[j];
-           lis_temp->zone_fract[j] = temp[i].zone_fract[j];     
-         }
-    }
-    
-
     if(vic411_options.BLOWING) {
       j = 2 * vic411_options.ROOT_ZONES;
       temp[i].sigma_slope = atof( vegarr[2 + j] );
@@ -238,16 +194,14 @@ int vic411_read_vegparam(FILE *vegparam,
       vic411_nrerror(ErrStr);
     }
     else
-     {
       temp[i].veg_class = veg_class;
-      }
 
     temp[0].Cv_sum += temp[i].Cv;
 
     for(k=0; k<Nfields; k++)
       free(vegarr[k]);
 
-    if ( vic411_options.GLOBAL_LAI) {
+    if ( vic411_options.GLOBAL_LAI ) {
       // Read the LAI line
       if ( fgets( line, MAXSTRING, vegparam ) == NULL ){
         sprintf(ErrStr,"ERROR unexpected EOF for cell %i while reading LAI for vegetat_type_num %d\n",vegcel,vegetat_type_num);
@@ -256,7 +210,7 @@ int vic411_read_vegparam(FILE *vegparam,
       Nfields = 0;
       vegarr[Nfields] = calloc( 500, sizeof(char));      
       strcpy(tmpline, line);
-      vic411_ttrim1( tmpline );
+      vic411_ttrim( tmpline );
       token = strtok (tmpline, delimiters); 
       strcpy(vegarr[Nfields],token);
       Nfields++;
@@ -272,20 +226,14 @@ int vic411_read_vegparam(FILE *vegparam,
         vic411_nrerror(ErrStr);
       }
 
-      if(vic_veg_class == lis_veg_class) // added by shugong for LIS tiling. Only if the current tile has the same vegetation class of lis tile
-      { 
-        // set vic411_flag = 1 indicating find the parameters for the vegetion class
-        vic411_flag = 1; 
-        for ( j = 0; j < 12; j++ ) 
-        {
-           vic411_veg_lib[temp[i].veg_class].LAI[j] = atof( vegarr[j] );
-           if (vic411_veg_lib[temp[i].veg_class].overstory && vic411_veg_lib[temp[i].veg_class].LAI[j] == 0) 
-           {
-               sprintf(ErrStr,"ERROR: cell %d, veg tile %d: the specified veg class (%d) is listed as an overstory class in the veg LIBRARY, but the LAI given in the veg PARAM FILE for this tile for month %d is 0.\n",gridcel, i+1, temp[i].veg_class+1, j+1);
-               vic411_nrerror(ErrStr);
-           }
-           vic411_veg_lib[temp[i].veg_class].Wdmax[j] =  LAI_WATER_FACTOR * vic411_veg_lib[temp[i].veg_class].LAI[j];
-         }
+      for ( j = 0; j < 12; j++ ) {
+        vic411_veg_lib[temp[i].veg_class].LAI[j] = atof( vegarr[j] );
+        if (vic411_veg_lib[temp[i].veg_class].overstory && vic411_veg_lib[temp[i].veg_class].LAI[j] == 0) {
+          sprintf(ErrStr,"ERROR: cell %d, veg tile %d: the specified veg class (%d) is listed as an overstory class in the veg LIBRARY, but the LAI given in the veg PARAM FILE for this tile for month %d is 0.\n",gridcel, i+1, temp[i].veg_class+1, j+1);
+          vic411_nrerror(ErrStr);
+        }
+        vic411_veg_lib[temp[i].veg_class].Wdmax[j] = 
+	  LAI_WATER_FACTOR * vic411_veg_lib[temp[i].veg_class].LAI[j];
       }
       for(k=0; k<Nfields; k++)
         free(vegarr[k]);
@@ -404,28 +352,19 @@ int vic411_read_vegparam(FILE *vegparam,
     }
   }
 
-  vic411_free_vegcon(&temp);
-  return(vic411_flag);
+  return temp;
 } 
 
 
 
 
 
-
+#if 0
 /* trim trailing newlines */
 
 #define END '\0'
 #define NEW '\n'
 
-void vic411_ttrim1( char *c ) 
-{
-  while( (*c++ != END) );
-    --c;
-  for( ; *--c == NEW; *c = END );
-
-}
-
 void vic411_ttrim( char *c ) 
 {
   while( (*c++ != END) );
@@ -433,4 +372,4 @@ void vic411_ttrim( char *c )
   for( ; *--c == NEW; *c = END );
 
 }
-
+#endif
diff --git a/lis/surfacemodels/land/vic.4.1.1/physics/vic411_vicNl.h b/lis/surfacemodels/land/vic.4.1.1/physics/vic411_vicNl.h
index 72507a592..e04da3738 100644
--- a/lis/surfacemodels/land/vic.4.1.1/physics/vic411_vicNl.h
+++ b/lis/surfacemodels/land/vic.4.1.1/physics/vic411_vicNl.h
@@ -395,7 +395,7 @@ void   read_snowmodel(vic411_atmos_data_struct *, FILE *, int, int, int, int);
 vic411_soil_con_struct vic411_read_soilparam(FILE *, int);
 //vic411_soil_con_struct vic411_read_soilparam_arc(FILE *, char *, int *, int *, int);
 vic411_veg_lib_struct *vic411_read_veglib(FILE *, int *);
-// vic411_veg_con_struct *vic411_read_vegparam(FILE *, int, int);
+vic411_veg_con_struct *vic411_read_vegparam(FILE *, int, int);
 int    vic411_redistribute_during_storm(vic411_cell_data_struct ***, vic411_veg_var_struct ***,
 				 int, int, int, double, double, double, 
 				 double *);
diff --git a/lis/surfacemodels/land/vic.4.1.1/setup_vic411.c b/lis/surfacemodels/land/vic.4.1.1/setup_vic411.c
index eaec5d847..b7fabf560 100644
--- a/lis/surfacemodels/land/vic.4.1.1/setup_vic411.c
+++ b/lis/surfacemodels/land/vic.4.1.1/setup_vic411.c
@@ -17,6 +17,10 @@
 #define VIC_BASED ( 0 )
 #define LIS_BASED ( 1 )
 
+extern void vic411_lis_scan_soilparam_flat_ascii(FILE *, float *, float *);
+
+extern void vic411_compute_treeline_lis(double, double *, char *);
+
 vic411_soil_con_struct vic411_lis_read_soilparam_arc(FILE *, char *, float *, float *);
 vic411_veg_con_struct *vic411_set_vegparam(int tile_idx,       // index of current tile 
                              int vegclass,       // vegetation class determined by LIS
diff --git a/lis/surfacemodels/land/vic.4.1.1/physics/vic411_read_vegparam.c.origin b/lis/surfacemodels/land/vic.4.1.1/vic411_read_vegparam_lis.c
similarity index 75%
rename from lis/surfacemodels/land/vic.4.1.1/physics/vic411_read_vegparam.c.origin
rename to lis/surfacemodels/land/vic.4.1.1/vic411_read_vegparam_lis.c
index b50986ec4..e02ca8bfb 100644
--- a/lis/surfacemodels/land/vic.4.1.1/physics/vic411_read_vegparam.c.origin
+++ b/lis/surfacemodels/land/vic.4.1.1/vic411_read_vegparam_lis.c
@@ -1,13 +1,25 @@
+//-----------------------BEGIN NOTICE -- DO NOT EDIT-----------------------
+// NASA Goddard Space Flight Center
+// Land Information System Framework (LISF)
+// Version 7.4
+//
+// Copyright (c) 2022 United States Government as represented by the
+// Administrator of the National Aeronautics and Space Administration.
+// All Rights Reserved.
+//-------------------------END NOTICE -- DO NOT EDIT-----------------------
+
 #include <stdio.h>
 #include <stdlib.h>
-#include <vic411_vicNl.h>
+#include "vic411_vicNl.h" 
 #include <string.h>
 
-static char vcid[] = "$Id: vic411_read_vegparam.c,v 4.5.2.11 2009/10/01 21:12:16 vicadmin Exp $";
-
-vic411_veg_con_struct *vic411_read_vegparam(FILE *vegparam,
-                              int   gridcel,
-                              int   Nveg_type)
+static char vcid[] = "$Id: vic411_read_vegparam_lis.c,v 4.5.2.11 2009/10/01 21:12:16 vicadmin Exp $";
+// shugong: if cannot find entry for the gridid+vegetation class, return 0, other return 1
+int vic411_read_vegparam_lis(FILE *vegparam,
+                         int   gridcel,
+                         int   Nveg_type,
+                         int   lis_veg_class,
+                         vic411_veg_con_struct *lis_temp)
 /**********************************************************************
   vic411_read_vegparam.c    Keith Cherkauer and Dag Lohmann       1997
 
@@ -46,7 +58,7 @@ vic411_veg_con_struct *vic411_read_vegparam(FILE *vegparam,
 **********************************************************************/
 {
 
-  void vic411_ttrim( char *string );
+  void vic411_ttrim1( char *string );
   extern vic411_veg_lib_struct *vic411_veg_lib;
   extern vic411_option_struct   vic411_options;
 #if LINK_DEBUG
@@ -68,6 +80,18 @@ vic411_veg_con_struct *vic411_read_vegparam(FILE *vegparam,
   char            *token;
   char            *vegarr[500];
   size_t	  length;
+  int vic411_flag;
+  vic411_flag = 0; // defualt: not found  
+  // added by Shugong Wang.
+  /* vic_veg_clas is defined to save the vegetation class in vegetation parameter file (called real vegetation class after). 
+     temp is a vic411_veg_con_struct. Firstly, the real vegetation class will be read in and save in the veg_class of temp. Then, 
+     vic411_read_vegparam will find the index number of the real vegetation class in vegetation parameter library (vic411_veg_lib). 
+     In the end, the veg_class of temp will be replced with the index number. In other word, the vegetation class in vic411_veg_con_struct
+     will not be the real vegetatoin class. In stead, it will be the corresponding index number in vegetatoin parameter library.
+     
+     In this function (vic411_read_vegparam), lis_veg_class is real vegetation class. In order to update the LAI and related values in
+     vic411_veg_lib correctly, lis_veg_class should be compared with the real vegetation class (vic_veg_class) */
+  int vic_veg_class;
 
   if(vic411_options.GLOBAL_LAI) skip=2;
   else skip=1;
@@ -97,7 +121,7 @@ vic411_veg_con_struct *vic411_read_vegparam(FILE *vegparam,
     exit(99);
   }
   if(vegetat_type_num >= MAX_VEG) {
-    sprintf(ErrStr,"Vegetation parameter file wants more vegetation tiles in grid cell %i (%i) than are allowed by MAX_VEG (%i) [NOTE: bare soil class is assumed].  Edit vic411_vicNl_def.h and recompile.",gridcel,vegetat_type_num+1,MAX_VEG);
+    sprintf(ErrStr,"Vegetation parameter file wants more vegetation tiles in grid cell %i (%i) than are allowed by MAX_VEG (%i) [NOTE: bare soil class is assumed].  Edit vicNl_def.h and recompile.",gridcel,vegetat_type_num+1,MAX_VEG);
     vic411_nrerror(ErrStr);
   }
 
@@ -122,7 +146,7 @@ vic411_veg_con_struct *vic411_read_vegparam(FILE *vegparam,
       vic411_nrerror(ErrStr);
     }
     strcpy(tmpline, line);
-    vic411_ttrim( tmpline );
+    vic411_ttrim1( tmpline );
     token = strtok (tmpline, delimiters);    /*  token => veg_class, move 'line' pointer to next field */  
     Nfields = 0;
     vegarr[Nfields] = calloc( 500, sizeof(char));
@@ -150,15 +174,20 @@ vic411_veg_con_struct *vic411_read_vegparam(FILE *vegparam,
 
     temp[i].LAKE = 0;
     temp[i].veg_class = atoi( vegarr[0] );
+    /* added by Shugong Wang */
+    vic_veg_class = temp[i].veg_class; // this is the real vegetation class in vegetation parameter file. 
+
     temp[i].Cv = atof( vegarr[1] );
     depth_sum = 0;
     sum = 0.;
+    /* the following VIC code read root zoon information from vegetation parameter file*/
     for(j=0;j<vic411_options.ROOT_ZONES;j++) {
       temp[i].zone_depth[j] = atof( vegarr[2 + j*2] );
       temp[i].zone_fract[j] = atof( vegarr[3 + j*2] );
       depth_sum += temp[i].zone_depth[j];
       sum += temp[i].zone_fract[j];
     }
+
     if(depth_sum <= 0) {
       sprintf(str,"Root zone depths must sum to a value greater than 0.");
       vic411_nrerror(str);
@@ -170,6 +199,21 @@ vic411_veg_con_struct *vic411_read_vegparam(FILE *vegparam,
       }
     }
 
+    /* added by Shugong Wang to read root zone information for LIS-tiling. 
+       if current vic_veg_class is the same as the lis_veg_class, then set the root zone
+       information to lis_temp, which points to vic411_lis_veg_con.                            */
+    if(vic_veg_class == lis_veg_class)
+    {
+         // added by Shugong Wang. lis_temp is a pointer of vic411_lis_veg_con. In the following lines, lis_temp will retrieve
+         // the root zone depth and root zone fraction from VIC vegetation parameter file
+         for(j=0; j<vic411_options.ROOT_ZONES; j++)
+         {
+           lis_temp->zone_depth[j] = temp[i].zone_depth[j];
+           lis_temp->zone_fract[j] = temp[i].zone_fract[j];     
+         }
+    }
+    
+
     if(vic411_options.BLOWING) {
       j = 2 * vic411_options.ROOT_ZONES;
       temp[i].sigma_slope = atof( vegarr[2 + j] );
@@ -194,14 +238,16 @@ vic411_veg_con_struct *vic411_read_vegparam(FILE *vegparam,
       vic411_nrerror(ErrStr);
     }
     else
+     {
       temp[i].veg_class = veg_class;
+      }
 
     temp[0].Cv_sum += temp[i].Cv;
 
     for(k=0; k<Nfields; k++)
       free(vegarr[k]);
 
-    if ( vic411_options.GLOBAL_LAI ) {
+    if ( vic411_options.GLOBAL_LAI) {
       // Read the LAI line
       if ( fgets( line, MAXSTRING, vegparam ) == NULL ){
         sprintf(ErrStr,"ERROR unexpected EOF for cell %i while reading LAI for vegetat_type_num %d\n",vegcel,vegetat_type_num);
@@ -210,7 +256,7 @@ vic411_veg_con_struct *vic411_read_vegparam(FILE *vegparam,
       Nfields = 0;
       vegarr[Nfields] = calloc( 500, sizeof(char));      
       strcpy(tmpline, line);
-      vic411_ttrim( tmpline );
+      vic411_ttrim1( tmpline );
       token = strtok (tmpline, delimiters); 
       strcpy(vegarr[Nfields],token);
       Nfields++;
@@ -226,14 +272,20 @@ vic411_veg_con_struct *vic411_read_vegparam(FILE *vegparam,
         vic411_nrerror(ErrStr);
       }
 
-      for ( j = 0; j < 12; j++ ) {
-        vic411_veg_lib[temp[i].veg_class].LAI[j] = atof( vegarr[j] );
-        if (vic411_veg_lib[temp[i].veg_class].overstory && vic411_veg_lib[temp[i].veg_class].LAI[j] == 0) {
-          sprintf(ErrStr,"ERROR: cell %d, veg tile %d: the specified veg class (%d) is listed as an overstory class in the veg LIBRARY, but the LAI given in the veg PARAM FILE for this tile for month %d is 0.\n",gridcel, i+1, temp[i].veg_class+1, j+1);
-          vic411_nrerror(ErrStr);
-        }
-        vic411_veg_lib[temp[i].veg_class].Wdmax[j] = 
-	  LAI_WATER_FACTOR * vic411_veg_lib[temp[i].veg_class].LAI[j];
+      if(vic_veg_class == lis_veg_class) // added by shugong for LIS tiling. Only if the current tile has the same vegetation class of lis tile
+      { 
+        // set vic411_flag = 1 indicating find the parameters for the vegetion class
+        vic411_flag = 1; 
+        for ( j = 0; j < 12; j++ ) 
+        {
+           vic411_veg_lib[temp[i].veg_class].LAI[j] = atof( vegarr[j] );
+           if (vic411_veg_lib[temp[i].veg_class].overstory && vic411_veg_lib[temp[i].veg_class].LAI[j] == 0) 
+           {
+               sprintf(ErrStr,"ERROR: cell %d, veg tile %d: the specified veg class (%d) is listed as an overstory class in the veg LIBRARY, but the LAI given in the veg PARAM FILE for this tile for month %d is 0.\n",gridcel, i+1, temp[i].veg_class+1, j+1);
+               vic411_nrerror(ErrStr);
+           }
+           vic411_veg_lib[temp[i].veg_class].Wdmax[j] =  LAI_WATER_FACTOR * vic411_veg_lib[temp[i].veg_class].LAI[j];
+         }
       }
       for(k=0; k<Nfields; k++)
         free(vegarr[k]);
@@ -352,7 +404,8 @@ vic411_veg_con_struct *vic411_read_vegparam(FILE *vegparam,
     }
   }
 
-  return temp;
+  vic411_free_vegcon(&temp);
+  return(vic411_flag);
 } 
 
 
@@ -365,6 +418,14 @@ vic411_veg_con_struct *vic411_read_vegparam(FILE *vegparam,
 #define END '\0'
 #define NEW '\n'
 
+void vic411_ttrim1( char *c ) 
+{
+  while( (*c++ != END) );
+    --c;
+  for( ; *--c == NEW; *c = END );
+
+}
+
 void vic411_ttrim( char *c ) 
 {
   while( (*c++ != END) );
diff --git a/lis/surfacemodels/land/vic.4.1.1/vic411_run.c b/lis/surfacemodels/land/vic.4.1.1/vic411_run.c
index f30fa281b..c41114f99 100644
--- a/lis/surfacemodels/land/vic.4.1.1/vic411_run.c
+++ b/lis/surfacemodels/land/vic.4.1.1/vic411_run.c
@@ -18,8 +18,7 @@
 #define VIC_BASED 0  /* VIC based tiling */
 #define LIS_BASED 1  /* LIS based tiling */
 
-//extern void FTN(vic411_diagnoseoutputvar)(int, int, int, int, float, 
-//                                          char *, char *);
+extern void FTN(vic411_diagnoseoutputvar)(int *, int *, int *, int *, double *, char *, int *, char *, int *);
 
 //BOP
 //
diff --git a/lis/surfacemodels/land/vic.4.1.1/vic411_set_vegparam.c b/lis/surfacemodels/land/vic.4.1.1/vic411_set_vegparam.c
index 0c0845779..cfe7e8c04 100644
--- a/lis/surfacemodels/land/vic.4.1.1/vic411_set_vegparam.c
+++ b/lis/surfacemodels/land/vic.4.1.1/vic411_set_vegparam.c
@@ -14,6 +14,8 @@
 
 #define BARETYPE ( 12 ) // UMD -- HARDCODED!
 
+extern int vic411_read_vegparam_lis(FILE *, int, int, int, vic411_veg_con_struct *);
+
 int vic411_real_2_vic(int real_type, int Nveg_type)
 {
    extern vic411_veg_lib_struct *vic411_veg_lib;
@@ -359,7 +361,7 @@ vic411_veg_con_struct *vic411_set_vegparam(int tile_idx,
          depth_sum = 0;
 
          // added by Shugong Wang to read in root zone information and update vic411_lis_veg_lib
-         vic411_flag = vic411_read_vegparam(fp_vegparam, gridcel , Nveg_type, vegclass, &(temp[i]));
+         vic411_flag = vic411_read_vegparam_lis(fp_vegparam, gridcel , Nveg_type, vegclass, &(temp[i]));
          if(vic411_flag==0) // cannot find veg parameters
          {
             // Since root zones are not defined they are copied from the last
diff --git a/lis/surfacemodels/land/vic.4.1.2.l/physics/initialize_model_state.c b/lis/surfacemodels/land/vic.4.1.2.l/physics/initialize_model_state.c
index a79475e00..d9a4e01e0 100644
--- a/lis/surfacemodels/land/vic.4.1.2.l/physics/initialize_model_state.c
+++ b/lis/surfacemodels/land/vic.4.1.2.l/physics/initialize_model_state.c
@@ -5,6 +5,8 @@
 
 static char vcid[] = "$Id$";
 
+extern void unpack_model_state(dist_prcp_struct *, global_param_struct *, int, int, int, soil_con_struct *, int, char *, int *, lake_con_struct, float *);
+
 int initialize_model_state(dist_prcp_struct    *prcp,
                            dmy_struct           dmy,
                            global_param_struct *global_param,
diff --git a/lis/surfacemodels/land/vic.4.1.2.l/set_vegparam.c b/lis/surfacemodels/land/vic.4.1.2.l/set_vegparam.c
index f1d3923ce..043efd7ea 100755
--- a/lis/surfacemodels/land/vic.4.1.2.l/set_vegparam.c
+++ b/lis/surfacemodels/land/vic.4.1.2.l/set_vegparam.c
@@ -14,6 +14,8 @@
 
 #define BARETYPE ( 12 ) // UMD -- HARDCODED!
 
+extern int vic412_read_vegparam(FILE *, int, int, int, veg_con_struct *);
+
 int real_2_vic(int real_type, int Nveg_type)
 {
    extern veg_lib_struct *veg_lib;
diff --git a/lis/surfacemodels/land/vic.4.1.2.l/setup_vic412.c b/lis/surfacemodels/land/vic.4.1.2.l/setup_vic412.c
index 305cf4586..0d11b4f81 100755
--- a/lis/surfacemodels/land/vic.4.1.2.l/setup_vic412.c
+++ b/lis/surfacemodels/land/vic.4.1.2.l/setup_vic412.c
@@ -17,6 +17,10 @@
 #define VIC_BASED ( 0 )
 #define LIS_BASED ( 1 )
 
+extern void lis_scan_soilparam_flat_ascii(FILE *, float *, float *);
+
+extern void vic412_compute_treeline(double, double *, char *);
+
 soil_con_struct         *lis_soil_con;
 dmy_struct              *dmy;
 atmos_data_struct      **lis_atmos;
diff --git a/lis/surfacemodels/land/vic.4.1.2.l/vic412_run.c b/lis/surfacemodels/land/vic.4.1.2.l/vic412_run.c
index eaa138800..f6e322776 100644
--- a/lis/surfacemodels/land/vic.4.1.2.l/vic412_run.c
+++ b/lis/surfacemodels/land/vic.4.1.2.l/vic412_run.c
@@ -18,8 +18,7 @@
 #define VIC_BASED 0  /* VIC based tiling */
 #define LIS_BASED 1  /* LIS based tiling */
 
-//extern void FTN(vic412_diagnoseoutputvar)(int, int, int, int, float, 
-//                                          char *, char *);
+extern void FTN(vic412_diagnoseoutputvar)(int *, int *, int *, int *, double *, char *, int *, char *, int *);
 
 //BOP
 //
diff --git a/lis/surfacemodels/land/vic.4.1.2.l/vic412_write_atmos.c b/lis/surfacemodels/land/vic.4.1.2.l/vic412_write_atmos.c
index 81874bb5d..9830394a3 100755
--- a/lis/surfacemodels/land/vic.4.1.2.l/vic412_write_atmos.c
+++ b/lis/surfacemodels/land/vic.4.1.2.l/vic412_write_atmos.c
@@ -11,6 +11,8 @@
 #include <vicNl.h>
 #include "ftn.h"
 
+extern void write_atmosdata(atmos_data_struct *, int);
+
 //BOP
 //
 // !ROUTINE: vic412_write_atmos
diff --git a/lis/surfacemodels/land/vic.4.1.2.l/vic412_write_state.c b/lis/surfacemodels/land/vic.4.1.2.l/vic412_write_state.c
index 3a57c3e89..4fbd36c50 100755
--- a/lis/surfacemodels/land/vic.4.1.2.l/vic412_write_state.c
+++ b/lis/surfacemodels/land/vic.4.1.2.l/vic412_write_state.c
@@ -13,6 +13,9 @@
 #include "ftn.h"
 #define VIC_BASED 0  /** VIC based vegetation tiling  Added by Shugong Wang on 05/01/2012 **/
 #define LIS_BASED 1  /** LIS based vegetation tiling  Added by Shugong Wang on 05/01/2012 **/
+
+extern void pack_model_state(dist_prcp_struct *, global_param_struct *, int, int, filep_struct *, soil_con_struct *, char *, int *, float *, int *);
+
 //BOP
 //
 // !ROUTINE: vic412_write_state