"ERRWAT is positive" - "Water budget problem in NOAHMP LSM"

[bhanu-magotra]

Bug Description

I am facing an error " The model is gaining water (ERRWAT is positive)" while running LISF 7.4 with NoahMP3.6 and HyMAP2 in restart configuration. I am using ESMWF's ERA5 with MERRA2 in overlay at 1 hour model timestep with outputs at daily timestep. LIS works fine in "coldstart" mode but in "restart" mode, it runs for a few days and then stops abruptly. I have attached the relevant logs, config and restart files below. Please help me to resolve the issue.

The mpirun error log records the following :

application called MPI_Abort(comm=0x84000002, 1) - process 10
 The model is gaining water (ERRWAT is positive)
  ERRWAT =  0.2526955     kg m{-2} timestep{-1}
     I      J     END_WB     BEG_WB       PRCP       ECAN       EDIR      ETRAN 
      RUNSRF     RUNSUB
     31    200        5778.792       5778.929    0.08638   -0.05000    0.00164  
   0.00000    0.52024    0.00440    2.76236
real 1097.61
user 11431.53
sys 624.29

And the mpirun output log records :

csky093.hpc.iitd.ac.in,csky095.hpc.iitd.ac.in,vsky010.hpc.iitd.ac.in
csky093: Attaching 18226 to 3244556.pbshpc
csky095: Attaching 12888 to 3244556.pbshpc
vsky010: Attaching 506 to 3244556.pbshpc
Attaching 576 to 3244556.pbshpc
Attaching 12953 to 3244556.pbshpc
 Water budget problem in NOAHMP LSM

The lislog.xxxx appears like this :

  [INFO] *** NASA Land Information System (LIS) ***
 [INFO] starting time:            4 /          30 /        1983
 [INFO] ending time:            1 /           1 /        1985
   
 [INFO] Output start time:   1983.32602739726             1983           4
          30           0           0           0
  [INFO] ************************************************
  [INFO] Start date (ymd tod):             1983           4          30
  [INFO] Stop date (ymd tod):              1985           1           1
  [INFO] Current date (ymd tod):           1983           4          30
  [INFO] ************************************************
 [INFO] DOMAIN details:
 [INFO] local domain:(          50          54 )
 [INFO] local domain without halo:(          50          54 )
 [INFO] running domain:(         300         320 )
 [INFO] local domain    5.550000       10.85000       68.05000    
   72.95000     d:            1
 [INFO] --------------------Domain            1 ----------------------
 [INFO] (           1 ,  0.0000000E+00 )
 [INFO] (           2 ,   50.00000     )
 [INFO] (           3 ,   54.00000     )
 [INFO] (           4 ,   5.550000     )
 [INFO] (           5 ,   68.05000     )
 [INFO] (           6 ,   128.0000     )
 [INFO] (           7 ,   10.85000     )
 [INFO] (           8 ,   72.95000     )
 [INFO] (           9 ,  0.1000000     )
 [INFO] (          10 ,  0.1000000     )
 [INFO] (          11 ,   64.00000     )
 [INFO] (          12 ,  0.0000000E+00 )
 [INFO] (          13 ,  0.0000000E+00 )
 [INFO] --------------------------------------------------------------
 [INFO] Reading mask from ./lis_input_files/lis_input_010.nc
 [INFO] Saving global land mask...
 [INFO] Reading landcover map from ./lis_input_files/lis_input_010.nc
 [INFO] Reading surfacetype map from ./lis_input_files/lis_input_010.nc
 [INFO] Finished reading landmask and landcover data
 [INFO] Reading elevation map from ./lis_input_files/lis_input_010.nc
 [INFO] Reading soiltexture map from ./lis_input_files/lis_input_010.nc
 [INFO] Landmask LAND points =            0
 [INFO] Total number of surface tiles           0
 [INFO] Total number of tiles for LSM           0
 [INFO] Total number of grid points           0
 [INFO] Opening Model Output Attributes File
 ./NoahMP_InputFiles/MODEL_OUTPUT_LIST_FINAL.TBL                                
                                                                                
                                                                                
                                                                                
                                                                                
                                                                                
                                                                                
                                                                                
                                                                                
                                                                                
                                                                                
                                                                                
                                                                             
 [INFO] Reading greenness map for month            4 from 
 ./lis_input_files/lis_input_010.nc
 [INFO] Reading greenness map for month            5 from 
 ./lis_input_files/lis_input_010.nc
 [INFO] Reading albedo map for month            4
 [INFO] Reading albedo map for month            5
 [INFO] Reading max snow albedo map from ./lis_input_files/lis_input_010.nc
 Start reading LIS configuration file for Noah-MP.3.6 model
 Finish reading LIS configuration file for Noah-MP.3.6 model
 [INFO] Opening Forcing Variables list file 
 ./NoahMP_InputFiles/forcing_variables.txt                                      
                                                                                
                                                                                
                                                                                
                                                                                
                                                                                
                                                                                
                                                                                
                                                                                
                                                                                
                                                                                
                                                                                
                                                                             
 [INFO] FORCING: Near Surface Air Temperature           1 K         
 [INFO] FORCING: Near Surface Specific Humidity           1 kg/kg     
 [INFO] FORCING: Incident Shortwave Radiation           1 W/m2      
 [INFO] FORCING: Incident Longwave Radiation           1 W/m2      
 [INFO] FORCING: Eastward Wind           1 m/s       
 [INFO] FORCING: Northward Wind           1 m/s       
 [INFO] FORCING: Surface Pressure           1 Pa        
 [INFO] FORCING: Rainfall Rate           1 kg/m2s    
 [INFO] Using MERRA2 forcing
 [INFO] MERRA2 forcing directory: ../../datasets/MERRA2
 [INFO] Using ERA5 forcing
 [INFO] ERA5 forcing directory: ../../datasets/strERA5
 [INFO] HYMAP2 2-way coupling: deactivated
 [INFO] Initializing HYMAP....
 [INFO] columns and rows          50          54
 [INFO] Get number of HYMAP2 grid cells
 [get_vector_size] number of cells           0
 [INFO] Get HYMAP2 cell vector sequence
 [INFO] Remove zeros from groundwater time delay
 [INFO] Processing data before running HYMAP
 [INFO] Calculate maximum river storage
 [INFO] HYMAP2 restart file used: 
 ./restart_files/ERA5/LIS_RST_HYMAP2_router_198304302300.d01.nc
 NoahMP36: retrieve parameter VEGETYPE from LIS
 NoahMP36: retrieve parameter SOILTYPE from LIS
 NoahMP36: reading parameter SLOPETYPE from ./lis_input_files/lis_input_010.nc
 NoahMP36: reading parameter TBOT from ./lis_input_files/lis_input_010.nc
 NoahMP36: reading parameter PBLH from ./lis_input_files/lis_input_010.nc
 NoahMP36: reading parameter SHDFAC_MONTHLY from 
 ./lis_input_files/lis_input_010.nc
 Reading GREENNESS map for level            1
 Reading GREENNESS map for level            2
 Reading GREENNESS map for level            3
 Reading GREENNESS map for level            4
 Reading GREENNESS map for level            5
 Reading GREENNESS map for level            6
 Reading GREENNESS map for level            7
 Reading GREENNESS map for level            8
 Reading GREENNESS map for level            9
 Reading GREENNESS map for level           10
 Reading GREENNESS map for level           11
 Reading GREENNESS map for level           12
  Skipping over LUTYPE = USGS
  LANDUSE TYPE = MODIFIED_IGBP_MODIS_NOAH FOUND          20  CATEGORIES
  INPUT SOIL TEXTURE CLASSIFICATION = STAS
  SOIL TEXTURE CLASSIFICATION = STAS FOUND          19  CATEGORIES
 NoahMP36 restart file used: 
 ./restart_files/ERA5/LIS_RST_NOAHMP36_198304302300.d01.nc                      
                                                                                
                                                                                
                                                                                
                                                                                
                                                                                
                                                                                
                                                                                
                                                                                
                                                                                
                                                                                
                                                                                
                                                                             
 [INFO]: LIS timestep    3600.000    
[INFO] LIS cycle time: 04/30/1983 01:00:00
 [INFO] Reading greenness map for month            4 from 
 ./lis_input_files/lis_input_010.nc
 [INFO] Reading greenness map for month            5 from 
 ./lis_input_files/lis_input_010.nc
 [INFO] Reading albedo map for month            4
 [INFO] Reading albedo map for month            5
 [INFO] Reading MERRA-2 file (bookend,           1  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M04/MERRA2_100.tavg1_2d_slv_Nx.19830430.
 nc4
 [INFO] Reading MERRA-2 file (bookend,           1  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M04/MERRA2_100.tavg1_2d_flx_Nx.19830430.
 nc4 (for lowest model level fields)
 [INFO] Reading MERRA-2 file (bookend,           1  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M04/MERRA2_100.tavg1_2d_flx_Nx.19830430.
 nc4 (for corrected precipitation fields)
 [INFO] Reading MERRA-2 file (bookend,           1  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M04/MERRA2_100.tavg1_2d_rad_Nx.19830430.
 nc4
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_slv_Nx.19830501.
 nc4
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830501.
 nc4 (for lowest model level fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830501.
 nc4 (for corrected precipitation fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_rad_Nx.19830501.
 nc4
 [INFO] Reading ERA5 file (bookend,           1  -
 ../../datasets/strERA5/FORCING_198304.nc)
 [INFO] Reading ERA5 file (bookend,           2  -
 ../../datasets/strERA5/FORCING_198304.nc)
[INFO] LIS cycle time: 04/30/1983 02:00:00
[INFO] LIS cycle time: 04/30/1983 03:00:00
[INFO] LIS cycle time: 04/30/1983 04:00:00
[INFO] LIS cycle time: 04/30/1983 05:00:00
[INFO] LIS cycle time: 04/30/1983 06:00:00
[INFO] LIS cycle time: 04/30/1983 07:00:00
[INFO] LIS cycle time: 04/30/1983 08:00:00
[INFO] LIS cycle time: 04/30/1983 09:00:00
[INFO] LIS cycle time: 04/30/1983 10:00:00
[INFO] LIS cycle time: 04/30/1983 11:00:00
[INFO] LIS cycle time: 04/30/1983 12:00:00
[INFO] LIS cycle time: 04/30/1983 13:00:00
[INFO] LIS cycle time: 04/30/1983 14:00:00
[INFO] LIS cycle time: 04/30/1983 15:00:00
[INFO] LIS cycle time: 04/30/1983 16:00:00
[INFO] LIS cycle time: 04/30/1983 17:00:00
[INFO] LIS cycle time: 04/30/1983 18:00:00
[INFO] LIS cycle time: 04/30/1983 19:00:00
[INFO] LIS cycle time: 04/30/1983 20:00:00
[INFO] LIS cycle time: 04/30/1983 21:00:00
[INFO] LIS cycle time: 04/30/1983 22:00:00
[INFO] LIS cycle time: 04/30/1983 23:00:00
 NoahMP36 archive restart written: 
 ./output/LIS_ERA5/SURFACEMODEL/198304/LIS_RST_NOAHMP36_198304302300.d01.nc
 [INFO] Writing routing restart 
 ./output/LIS_ERA5/ROUTING/198304/LIS_RST_HYMAP2_router_198304302300.d01.nc
[INFO] LIS cycle time: 05/01/1983 00:00:00
 [INFO] Reading ERA5 file (bookend,           2  -
 ../../datasets/strERA5/FORCING_198305.nc)
 [INFO] Writing surface model output to:  
 ./output/LIS_ERA5/SURFACEMODEL/198305/LIS_HIST_198305010000.d01.nc
 [INFO] Writing routing model output to:  
 ./output/LIS_ERA5/ROUTING/198305/LIS_HIST_198305010000.d01.nc
[INFO] LIS cycle time: 05/01/1983 01:00:00
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_slv_Nx.19830502.
 nc4
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830502.
 nc4 (for lowest model level fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830502.
 nc4 (for corrected precipitation fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_rad_Nx.19830502.
 nc4
[INFO] LIS cycle time: 05/01/1983 02:00:00
[INFO] LIS cycle time: 05/01/1983 03:00:00
[INFO] LIS cycle time: 05/01/1983 04:00:00
[INFO] LIS cycle time: 05/01/1983 05:00:00
[INFO] LIS cycle time: 05/01/1983 06:00:00
[INFO] LIS cycle time: 05/01/1983 07:00:00
[INFO] LIS cycle time: 05/01/1983 08:00:00
[INFO] LIS cycle time: 05/01/1983 09:00:00
[INFO] LIS cycle time: 05/01/1983 10:00:00
[INFO] LIS cycle time: 05/01/1983 11:00:00
[INFO] LIS cycle time: 05/01/1983 12:00:00
[INFO] LIS cycle time: 05/01/1983 13:00:00
[INFO] LIS cycle time: 05/01/1983 14:00:00
[INFO] LIS cycle time: 05/01/1983 15:00:00
[INFO] LIS cycle time: 05/01/1983 16:00:00
[INFO] LIS cycle time: 05/01/1983 17:00:00
[INFO] LIS cycle time: 05/01/1983 18:00:00
[INFO] LIS cycle time: 05/01/1983 19:00:00
[INFO] LIS cycle time: 05/01/1983 20:00:00
[INFO] LIS cycle time: 05/01/1983 21:00:00
[INFO] LIS cycle time: 05/01/1983 22:00:00
[INFO] LIS cycle time: 05/01/1983 23:00:00
[INFO] LIS cycle time: 05/02/1983 00:00:00
 [INFO] Writing surface model output to:  
 ./output/LIS_ERA5/SURFACEMODEL/198305/LIS_HIST_198305020000.d01.nc
 [INFO] Writing routing model output to:  
 ./output/LIS_ERA5/ROUTING/198305/LIS_HIST_198305020000.d01.nc
[INFO] LIS cycle time: 05/02/1983 01:00:00
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_slv_Nx.19830503.
 nc4
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830503.
 nc4 (for lowest model level fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830503.
 nc4 (for corrected precipitation fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_rad_Nx.19830503.
 nc4
[INFO] LIS cycle time: 05/02/1983 02:00:00
[INFO] LIS cycle time: 05/02/1983 03:00:00
[INFO] LIS cycle time: 05/02/1983 04:00:00
[INFO] LIS cycle time: 05/02/1983 05:00:00
[INFO] LIS cycle time: 05/02/1983 06:00:00
[INFO] LIS cycle time: 05/02/1983 07:00:00
[INFO] LIS cycle time: 05/02/1983 08:00:00
[INFO] LIS cycle time: 05/02/1983 09:00:00
[INFO] LIS cycle time: 05/02/1983 10:00:00
[INFO] LIS cycle time: 05/02/1983 11:00:00
[INFO] LIS cycle time: 05/02/1983 12:00:00
[INFO] LIS cycle time: 05/02/1983 13:00:00
[INFO] LIS cycle time: 05/02/1983 14:00:00
[INFO] LIS cycle time: 05/02/1983 15:00:00
[INFO] LIS cycle time: 05/02/1983 16:00:00
[INFO] LIS cycle time: 05/02/1983 17:00:00
[INFO] LIS cycle time: 05/02/1983 18:00:00
[INFO] LIS cycle time: 05/02/1983 19:00:00
[INFO] LIS cycle time: 05/02/1983 20:00:00
[INFO] LIS cycle time: 05/02/1983 21:00:00
[INFO] LIS cycle time: 05/02/1983 22:00:00
[INFO] LIS cycle time: 05/02/1983 23:00:00
[INFO] LIS cycle time: 05/03/1983 00:00:00
 [INFO] Writing surface model output to:  
 ./output/LIS_ERA5/SURFACEMODEL/198305/LIS_HIST_198305030000.d01.nc
 [INFO] Writing routing model output to:  
 ./output/LIS_ERA5/ROUTING/198305/LIS_HIST_198305030000.d01.nc
[INFO] LIS cycle time: 05/03/1983 01:00:00
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_slv_Nx.19830504.
 nc4
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830504.
 nc4 (for lowest model level fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830504.
 nc4 (for corrected precipitation fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_rad_Nx.19830504.
 nc4
[INFO] LIS cycle time: 05/03/1983 02:00:00
[INFO] LIS cycle time: 05/03/1983 03:00:00
[INFO] LIS cycle time: 05/03/1983 04:00:00
[INFO] LIS cycle time: 05/03/1983 05:00:00
[INFO] LIS cycle time: 05/03/1983 06:00:00
[INFO] LIS cycle time: 05/03/1983 07:00:00
[INFO] LIS cycle time: 05/03/1983 08:00:00
[INFO] LIS cycle time: 05/03/1983 09:00:00
[INFO] LIS cycle time: 05/03/1983 10:00:00
[INFO] LIS cycle time: 05/03/1983 11:00:00
[INFO] LIS cycle time: 05/03/1983 12:00:00
[INFO] LIS cycle time: 05/03/1983 13:00:00
[INFO] LIS cycle time: 05/03/1983 14:00:00
[INFO] LIS cycle time: 05/03/1983 15:00:00
[INFO] LIS cycle time: 05/03/1983 16:00:00
[INFO] LIS cycle time: 05/03/1983 17:00:00
[INFO] LIS cycle time: 05/03/1983 18:00:00
[INFO] LIS cycle time: 05/03/1983 19:00:00
[INFO] LIS cycle time: 05/03/1983 20:00:00
[INFO] LIS cycle time: 05/03/1983 21:00:00
[INFO] LIS cycle time: 05/03/1983 22:00:00
[INFO] LIS cycle time: 05/03/1983 23:00:00
[INFO] LIS cycle time: 05/04/1983 00:00:00
 [INFO] Writing surface model output to:  
 ./output/LIS_ERA5/SURFACEMODEL/198305/LIS_HIST_198305040000.d01.nc
 [INFO] Writing routing model output to:  
 ./output/LIS_ERA5/ROUTING/198305/LIS_HIST_198305040000.d01.nc
[INFO] LIS cycle time: 05/04/1983 01:00:00
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_slv_Nx.19830505.
 nc4
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830505.
 nc4 (for lowest model level fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830505.
 nc4 (for corrected precipitation fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_rad_Nx.19830505.
 nc4
[INFO] LIS cycle time: 05/04/1983 02:00:00
[INFO] LIS cycle time: 05/04/1983 03:00:00
[INFO] LIS cycle time: 05/04/1983 04:00:00
[INFO] LIS cycle time: 05/04/1983 05:00:00
[INFO] LIS cycle time: 05/04/1983 06:00:00
[INFO] LIS cycle time: 05/04/1983 07:00:00
[INFO] LIS cycle time: 05/04/1983 08:00:00
[INFO] LIS cycle time: 05/04/1983 09:00:00
[INFO] LIS cycle time: 05/04/1983 10:00:00
[INFO] LIS cycle time: 05/04/1983 11:00:00
[INFO] LIS cycle time: 05/04/1983 12:00:00
[INFO] LIS cycle time: 05/04/1983 13:00:00
[INFO] LIS cycle time: 05/04/1983 14:00:00
[INFO] LIS cycle time: 05/04/1983 15:00:00
[INFO] LIS cycle time: 05/04/1983 16:00:00
[INFO] LIS cycle time: 05/04/1983 17:00:00
[INFO] LIS cycle time: 05/04/1983 18:00:00
[INFO] LIS cycle time: 05/04/1983 19:00:00
[INFO] LIS cycle time: 05/04/1983 20:00:00
[INFO] LIS cycle time: 05/04/1983 21:00:00
[INFO] LIS cycle time: 05/04/1983 22:00:00
[INFO] LIS cycle time: 05/04/1983 23:00:00
[INFO] LIS cycle time: 05/05/1983 00:00:00
 [INFO] Writing surface model output to:  
 ./output/LIS_ERA5/SURFACEMODEL/198305/LIS_HIST_198305050000.d01.nc
 [INFO] Writing routing model output to:  
 ./output/LIS_ERA5/ROUTING/198305/LIS_HIST_198305050000.d01.nc
[INFO] LIS cycle time: 05/05/1983 01:00:00
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_slv_Nx.19830506.
 nc4
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830506.
 nc4 (for lowest model level fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830506.
 nc4 (for corrected precipitation fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_rad_Nx.19830506.
 nc4
[INFO] LIS cycle time: 05/05/1983 02:00:00
[INFO] LIS cycle time: 05/05/1983 03:00:00
[INFO] LIS cycle time: 05/05/1983 04:00:00
[INFO] LIS cycle time: 05/05/1983 05:00:00
[INFO] LIS cycle time: 05/05/1983 06:00:00
[INFO] LIS cycle time: 05/05/1983 07:00:00
[INFO] LIS cycle time: 05/05/1983 08:00:00
[INFO] LIS cycle time: 05/05/1983 09:00:00
[INFO] LIS cycle time: 05/05/1983 10:00:00
[INFO] LIS cycle time: 05/05/1983 11:00:00
[INFO] LIS cycle time: 05/05/1983 12:00:00
[INFO] LIS cycle time: 05/05/1983 13:00:00
[INFO] LIS cycle time: 05/05/1983 14:00:00
[INFO] LIS cycle time: 05/05/1983 15:00:00
[INFO] LIS cycle time: 05/05/1983 16:00:00
[INFO] LIS cycle time: 05/05/1983 17:00:00
[INFO] LIS cycle time: 05/05/1983 18:00:00
[INFO] LIS cycle time: 05/05/1983 19:00:00
[INFO] LIS cycle time: 05/05/1983 20:00:00
[INFO] LIS cycle time: 05/05/1983 21:00:00
[INFO] LIS cycle time: 05/05/1983 22:00:00
[INFO] LIS cycle time: 05/05/1983 23:00:00
[INFO] LIS cycle time: 05/06/1983 00:00:00
 [INFO] Writing surface model output to:  
 ./output/LIS_ERA5/SURFACEMODEL/198305/LIS_HIST_198305060000.d01.nc
 [INFO] Writing routing model output to:  
 ./output/LIS_ERA5/ROUTING/198305/LIS_HIST_198305060000.d01.nc
[INFO] LIS cycle time: 05/06/1983 01:00:00
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_slv_Nx.19830507.
 nc4
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830507.
 nc4 (for lowest model level fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830507.
 nc4 (for corrected precipitation fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_rad_Nx.19830507.
 nc4
[INFO] LIS cycle time: 05/06/1983 02:00:00
[INFO] LIS cycle time: 05/06/1983 03:00:00
[INFO] LIS cycle time: 05/06/1983 04:00:00
[INFO] LIS cycle time: 05/06/1983 05:00:00
[INFO] LIS cycle time: 05/06/1983 06:00:00
[INFO] LIS cycle time: 05/06/1983 07:00:00
[INFO] LIS cycle time: 05/06/1983 08:00:00
[INFO] LIS cycle time: 05/06/1983 09:00:00
[INFO] LIS cycle time: 05/06/1983 10:00:00
[INFO] LIS cycle time: 05/06/1983 11:00:00
[INFO] LIS cycle time: 05/06/1983 12:00:00
[INFO] LIS cycle time: 05/06/1983 13:00:00
[INFO] LIS cycle time: 05/06/1983 14:00:00
[INFO] LIS cycle time: 05/06/1983 15:00:00
[INFO] LIS cycle time: 05/06/1983 16:00:00
[INFO] LIS cycle time: 05/06/1983 17:00:00
[INFO] LIS cycle time: 05/06/1983 18:00:00
[INFO] LIS cycle time: 05/06/1983 19:00:00
[INFO] LIS cycle time: 05/06/1983 20:00:00
[INFO] LIS cycle time: 05/06/1983 21:00:00
[INFO] LIS cycle time: 05/06/1983 22:00:00
[INFO] LIS cycle time: 05/06/1983 23:00:00
[INFO] LIS cycle time: 05/07/1983 00:00:00
 [INFO] Writing surface model output to:  
 ./output/LIS_ERA5/SURFACEMODEL/198305/LIS_HIST_198305070000.d01.nc
 [INFO] Writing routing model output to:  
 ./output/LIS_ERA5/ROUTING/198305/LIS_HIST_198305070000.d01.nc
[INFO] LIS cycle time: 05/07/1983 01:00:00
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_slv_Nx.19830508.
 nc4
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830508.
 nc4 (for lowest model level fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830508.
 nc4 (for corrected precipitation fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_rad_Nx.19830508.
 nc4
[INFO] LIS cycle time: 05/07/1983 02:00:00
[INFO] LIS cycle time: 05/07/1983 03:00:00
[INFO] LIS cycle time: 05/07/1983 04:00:00
[INFO] LIS cycle time: 05/07/1983 05:00:00
[INFO] LIS cycle time: 05/07/1983 06:00:00
[INFO] LIS cycle time: 05/07/1983 07:00:00
[INFO] LIS cycle time: 05/07/1983 08:00:00
[INFO] LIS cycle time: 05/07/1983 09:00:00
[INFO] LIS cycle time: 05/07/1983 10:00:00
[INFO] LIS cycle time: 05/07/1983 11:00:00
[INFO] LIS cycle time: 05/07/1983 12:00:00
[INFO] LIS cycle time: 05/07/1983 13:00:00
[INFO] LIS cycle time: 05/07/1983 14:00:00
[INFO] LIS cycle time: 05/07/1983 15:00:00
[INFO] LIS cycle time: 05/07/1983 16:00:00
[INFO] LIS cycle time: 05/07/1983 17:00:00
[INFO] LIS cycle time: 05/07/1983 18:00:00
[INFO] LIS cycle time: 05/07/1983 19:00:00
[INFO] LIS cycle time: 05/07/1983 20:00:00
[INFO] LIS cycle time: 05/07/1983 21:00:00
[INFO] LIS cycle time: 05/07/1983 22:00:00
[INFO] LIS cycle time: 05/07/1983 23:00:00
[INFO] LIS cycle time: 05/08/1983 00:00:00
 [INFO] Writing surface model output to:  
 ./output/LIS_ERA5/SURFACEMODEL/198305/LIS_HIST_198305080000.d01.nc
 [INFO] Writing routing model output to:  
 ./output/LIS_ERA5/ROUTING/198305/LIS_HIST_198305080000.d01.nc
[INFO] LIS cycle time: 05/08/1983 01:00:00
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_slv_Nx.19830509.
 nc4
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830509.
 nc4 (for lowest model level fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830509.
 nc4 (for corrected precipitation fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_rad_Nx.19830509.
 nc4
[INFO] LIS cycle time: 05/08/1983 02:00:00
[INFO] LIS cycle time: 05/08/1983 03:00:00
[INFO] LIS cycle time: 05/08/1983 04:00:00
[INFO] LIS cycle time: 05/08/1983 05:00:00
[INFO] LIS cycle time: 05/08/1983 06:00:00
[INFO] LIS cycle time: 05/08/1983 07:00:00
[INFO] LIS cycle time: 05/08/1983 08:00:00
[INFO] LIS cycle time: 05/08/1983 09:00:00
[INFO] LIS cycle time: 05/08/1983 10:00:00
[INFO] LIS cycle time: 05/08/1983 11:00:00
[INFO] LIS cycle time: 05/08/1983 12:00:00
[INFO] LIS cycle time: 05/08/1983 13:00:00
[INFO] LIS cycle time: 05/08/1983 14:00:00
[INFO] LIS cycle time: 05/08/1983 15:00:00
[INFO] LIS cycle time: 05/08/1983 16:00:00
[INFO] LIS cycle time: 05/08/1983 17:00:00
[INFO] LIS cycle time: 05/08/1983 18:00:00
[INFO] LIS cycle time: 05/08/1983 19:00:00
[INFO] LIS cycle time: 05/08/1983 20:00:00
[INFO] LIS cycle time: 05/08/1983 21:00:00
[INFO] LIS cycle time: 05/08/1983 22:00:00
[INFO] LIS cycle time: 05/08/1983 23:00:00
[INFO] LIS cycle time: 05/09/1983 00:00:00
 [INFO] Writing surface model output to:  
 ./output/LIS_ERA5/SURFACEMODEL/198305/LIS_HIST_198305090000.d01.nc
 [INFO] Writing routing model output to:  
 ./output/LIS_ERA5/ROUTING/198305/LIS_HIST_198305090000.d01.nc
[INFO] LIS cycle time: 05/09/1983 01:00:00
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_slv_Nx.19830510.
 nc4
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830510.
 nc4 (for lowest model level fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830510.
 nc4 (for corrected precipitation fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_rad_Nx.19830510.
 nc4
[INFO] LIS cycle time: 05/09/1983 02:00:00
[INFO] LIS cycle time: 05/09/1983 03:00:00
[INFO] LIS cycle time: 05/09/1983 04:00:00
[INFO] LIS cycle time: 05/09/1983 05:00:00
[INFO] LIS cycle time: 05/09/1983 06:00:00
[INFO] LIS cycle time: 05/09/1983 07:00:00
[INFO] LIS cycle time: 05/09/1983 08:00:00
[INFO] LIS cycle time: 05/09/1983 09:00:00
[INFO] LIS cycle time: 05/09/1983 10:00:00
[INFO] LIS cycle time: 05/09/1983 11:00:00
[INFO] LIS cycle time: 05/09/1983 12:00:00
[INFO] LIS cycle time: 05/09/1983 13:00:00
[INFO] LIS cycle time: 05/09/1983 14:00:00
[INFO] LIS cycle time: 05/09/1983 15:00:00
[INFO] LIS cycle time: 05/09/1983 16:00:00
[INFO] LIS cycle time: 05/09/1983 17:00:00
[INFO] LIS cycle time: 05/09/1983 18:00:00
[INFO] LIS cycle time: 05/09/1983 19:00:00
[INFO] LIS cycle time: 05/09/1983 20:00:00
[INFO] LIS cycle time: 05/09/1983 21:00:00
[INFO] LIS cycle time: 05/09/1983 22:00:00
[INFO] LIS cycle time: 05/09/1983 23:00:00
[INFO] LIS cycle time: 05/10/1983 00:00:00
 [INFO] Writing surface model output to:  
 ./output/LIS_ERA5/SURFACEMODEL/198305/LIS_HIST_198305100000.d01.nc
 [INFO] Writing routing model output to:  
 ./output/LIS_ERA5/ROUTING/198305/LIS_HIST_198305100000.d01.nc
[INFO] LIS cycle time: 05/10/1983 01:00:00
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_slv_Nx.19830511.
 nc4
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830511.
 nc4 (for lowest model level fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830511.
 nc4 (for corrected precipitation fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_rad_Nx.19830511.
 nc4
[INFO] LIS cycle time: 05/10/1983 02:00:00
[INFO] LIS cycle time: 05/10/1983 03:00:00
[INFO] LIS cycle time: 05/10/1983 04:00:00
[INFO] LIS cycle time: 05/10/1983 05:00:00
[INFO] LIS cycle time: 05/10/1983 06:00:00
[INFO] LIS cycle time: 05/10/1983 07:00:00
[INFO] LIS cycle time: 05/10/1983 08:00:00
[INFO] LIS cycle time: 05/10/1983 09:00:00
[INFO] LIS cycle time: 05/10/1983 10:00:00
[INFO] LIS cycle time: 05/10/1983 11:00:00
[INFO] LIS cycle time: 05/10/1983 12:00:00
[INFO] LIS cycle time: 05/10/1983 13:00:00
[INFO] LIS cycle time: 05/10/1983 14:00:00
[INFO] LIS cycle time: 05/10/1983 15:00:00
[INFO] LIS cycle time: 05/10/1983 16:00:00
[INFO] LIS cycle time: 05/10/1983 17:00:00
[INFO] LIS cycle time: 05/10/1983 18:00:00
[INFO] LIS cycle time: 05/10/1983 19:00:00
[INFO] LIS cycle time: 05/10/1983 20:00:00
[INFO] LIS cycle time: 05/10/1983 21:00:00
[INFO] LIS cycle time: 05/10/1983 22:00:00
[INFO] LIS cycle time: 05/10/1983 23:00:00
[INFO] LIS cycle time: 05/11/1983 00:00:00
 [INFO] Writing surface model output to:  
 ./output/LIS_ERA5/SURFACEMODEL/198305/LIS_HIST_198305110000.d01.nc
 [INFO] Writing routing model output to:  
 ./output/LIS_ERA5/ROUTING/198305/LIS_HIST_198305110000.d01.nc
[INFO] LIS cycle time: 05/11/1983 01:00:00
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_slv_Nx.19830512.
 nc4
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830512.
 nc4 (for lowest model level fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830512.
 nc4 (for corrected precipitation fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_rad_Nx.19830512.
 nc4
[INFO] LIS cycle time: 05/11/1983 02:00:00
[INFO] LIS cycle time: 05/11/1983 03:00:00
[INFO] LIS cycle time: 05/11/1983 04:00:00
[INFO] LIS cycle time: 05/11/1983 05:00:00
[INFO] LIS cycle time: 05/11/1983 06:00:00
[INFO] LIS cycle time: 05/11/1983 07:00:00
[INFO] LIS cycle time: 05/11/1983 08:00:00
[INFO] LIS cycle time: 05/11/1983 09:00:00
[INFO] LIS cycle time: 05/11/1983 10:00:00
[INFO] LIS cycle time: 05/11/1983 11:00:00
[INFO] LIS cycle time: 05/11/1983 12:00:00
[INFO] LIS cycle time: 05/11/1983 13:00:00
[INFO] LIS cycle time: 05/11/1983 14:00:00
[INFO] LIS cycle time: 05/11/1983 15:00:00
[INFO] LIS cycle time: 05/11/1983 16:00:00
[INFO] LIS cycle time: 05/11/1983 17:00:00
[INFO] LIS cycle time: 05/11/1983 18:00:00
[INFO] LIS cycle time: 05/11/1983 19:00:00
[INFO] LIS cycle time: 05/11/1983 20:00:00
[INFO] LIS cycle time: 05/11/1983 21:00:00
[INFO] LIS cycle time: 05/11/1983 22:00:00
[INFO] LIS cycle time: 05/11/1983 23:00:00
[INFO] LIS cycle time: 05/12/1983 00:00:00
 [INFO] Writing surface model output to:  
 ./output/LIS_ERA5/SURFACEMODEL/198305/LIS_HIST_198305120000.d01.nc
 [INFO] Writing routing model output to:  
 ./output/LIS_ERA5/ROUTING/198305/LIS_HIST_198305120000.d01.nc
[INFO] LIS cycle time: 05/12/1983 01:00:00
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_slv_Nx.19830513.
 nc4
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830513.
 nc4 (for lowest model level fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830513.
 nc4 (for corrected precipitation fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_rad_Nx.19830513.
 nc4
[INFO] LIS cycle time: 05/12/1983 02:00:00
[INFO] LIS cycle time: 05/12/1983 03:00:00
[INFO] LIS cycle time: 05/12/1983 04:00:00
[INFO] LIS cycle time: 05/12/1983 05:00:00
[INFO] LIS cycle time: 05/12/1983 06:00:00
[INFO] LIS cycle time: 05/12/1983 07:00:00
[INFO] LIS cycle time: 05/12/1983 08:00:00
[INFO] LIS cycle time: 05/12/1983 09:00:00
[INFO] LIS cycle time: 05/12/1983 10:00:00
[INFO] LIS cycle time: 05/12/1983 11:00:00
[INFO] LIS cycle time: 05/12/1983 12:00:00
[INFO] LIS cycle time: 05/12/1983 13:00:00
[INFO] LIS cycle time: 05/12/1983 14:00:00
[INFO] LIS cycle time: 05/12/1983 15:00:00
[INFO] LIS cycle time: 05/12/1983 16:00:00
[INFO] LIS cycle time: 05/12/1983 17:00:00
[INFO] LIS cycle time: 05/12/1983 18:00:00
[INFO] LIS cycle time: 05/12/1983 19:00:00
[INFO] LIS cycle time: 05/12/1983 20:00:00
[INFO] LIS cycle time: 05/12/1983 21:00:00
[INFO] LIS cycle time: 05/12/1983 22:00:00
[INFO] LIS cycle time: 05/12/1983 23:00:00
[INFO] LIS cycle time: 05/13/1983 00:00:00
 [INFO] Writing surface model output to:  
 ./output/LIS_ERA5/SURFACEMODEL/198305/LIS_HIST_198305130000.d01.nc
 [INFO] Writing routing model output to:  
 ./output/LIS_ERA5/ROUTING/198305/LIS_HIST_198305130000.d01.nc
[INFO] LIS cycle time: 05/13/1983 01:00:00
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_slv_Nx.19830514.
 nc4
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830514.
 nc4 (for lowest model level fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830514.
 nc4 (for corrected precipitation fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_rad_Nx.19830514.
 nc4
[INFO] LIS cycle time: 05/13/1983 02:00:00
[INFO] LIS cycle time: 05/13/1983 03:00:00
[INFO] LIS cycle time: 05/13/1983 04:00:00
[INFO] LIS cycle time: 05/13/1983 05:00:00
[INFO] LIS cycle time: 05/13/1983 06:00:00
[INFO] LIS cycle time: 05/13/1983 07:00:00
[INFO] LIS cycle time: 05/13/1983 08:00:00
[INFO] LIS cycle time: 05/13/1983 09:00:00
[INFO] LIS cycle time: 05/13/1983 10:00:00
[INFO] LIS cycle time: 05/13/1983 11:00:00
[INFO] LIS cycle time: 05/13/1983 12:00:00
[INFO] LIS cycle time: 05/13/1983 13:00:00
[INFO] LIS cycle time: 05/13/1983 14:00:00
[INFO] LIS cycle time: 05/13/1983 15:00:00
[INFO] LIS cycle time: 05/13/1983 16:00:00
[INFO] LIS cycle time: 05/13/1983 17:00:00
[INFO] LIS cycle time: 05/13/1983 18:00:00
[INFO] LIS cycle time: 05/13/1983 19:00:00
[INFO] LIS cycle time: 05/13/1983 20:00:00
[INFO] LIS cycle time: 05/13/1983 21:00:00
[INFO] LIS cycle time: 05/13/1983 22:00:00
[INFO] LIS cycle time: 05/13/1983 23:00:00
[INFO] LIS cycle time: 05/14/1983 00:00:00
 [INFO] Writing surface model output to:  
 ./output/LIS_ERA5/SURFACEMODEL/198305/LIS_HIST_198305140000.d01.nc
 [INFO] Writing routing model output to:  
 ./output/LIS_ERA5/ROUTING/198305/LIS_HIST_198305140000.d01.nc
[INFO] LIS cycle time: 05/14/1983 01:00:00
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_slv_Nx.19830515.
 nc4
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830515.
 nc4 (for lowest model level fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830515.
 nc4 (for corrected precipitation fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_rad_Nx.19830515.
 nc4
[INFO] LIS cycle time: 05/14/1983 02:00:00
[INFO] LIS cycle time: 05/14/1983 03:00:00
[INFO] LIS cycle time: 05/14/1983 04:00:00
[INFO] LIS cycle time: 05/14/1983 05:00:00
[INFO] LIS cycle time: 05/14/1983 06:00:00
[INFO] LIS cycle time: 05/14/1983 07:00:00
[INFO] LIS cycle time: 05/14/1983 08:00:00
[INFO] LIS cycle time: 05/14/1983 09:00:00
[INFO] LIS cycle time: 05/14/1983 10:00:00
[INFO] LIS cycle time: 05/14/1983 11:00:00
[INFO] LIS cycle time: 05/14/1983 12:00:00
[INFO] LIS cycle time: 05/14/1983 13:00:00
[INFO] LIS cycle time: 05/14/1983 14:00:00
[INFO] LIS cycle time: 05/14/1983 15:00:00
[INFO] LIS cycle time: 05/14/1983 16:00:00
[INFO] LIS cycle time: 05/14/1983 17:00:00
[INFO] LIS cycle time: 05/14/1983 18:00:00
[INFO] LIS cycle time: 05/14/1983 19:00:00
[INFO] LIS cycle time: 05/14/1983 20:00:00
[INFO] LIS cycle time: 05/14/1983 21:00:00
[INFO] LIS cycle time: 05/14/1983 22:00:00
[INFO] LIS cycle time: 05/14/1983 23:00:00
[INFO] LIS cycle time: 05/15/1983 00:00:00
 [INFO] Writing surface model output to:  
 ./output/LIS_ERA5/SURFACEMODEL/198305/LIS_HIST_198305150000.d01.nc
 [INFO] Writing routing model output to:  
 ./output/LIS_ERA5/ROUTING/198305/LIS_HIST_198305150000.d01.nc
[INFO] LIS cycle time: 05/15/1983 01:00:00
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_slv_Nx.19830516.
 nc4
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830516.
 nc4 (for lowest model level fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830516.
 nc4 (for corrected precipitation fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_rad_Nx.19830516.
 nc4
[INFO] LIS cycle time: 05/15/1983 02:00:00
[INFO] LIS cycle time: 05/15/1983 03:00:00
[INFO] LIS cycle time: 05/15/1983 04:00:00
[INFO] LIS cycle time: 05/15/1983 05:00:00
[INFO] LIS cycle time: 05/15/1983 06:00:00
[INFO] LIS cycle time: 05/15/1983 07:00:00
[INFO] LIS cycle time: 05/15/1983 08:00:00
[INFO] LIS cycle time: 05/15/1983 09:00:00
[INFO] LIS cycle time: 05/15/1983 10:00:00
[INFO] LIS cycle time: 05/15/1983 11:00:00
[INFO] LIS cycle time: 05/15/1983 12:00:00
[INFO] LIS cycle time: 05/15/1983 13:00:00
[INFO] LIS cycle time: 05/15/1983 14:00:00
[INFO] LIS cycle time: 05/15/1983 15:00:00
[INFO] LIS cycle time: 05/15/1983 16:00:00
[INFO] LIS cycle time: 05/15/1983 17:00:00
[INFO] LIS cycle time: 05/15/1983 18:00:00
[INFO] LIS cycle time: 05/15/1983 19:00:00
[INFO] LIS cycle time: 05/15/1983 20:00:00
[INFO] LIS cycle time: 05/15/1983 21:00:00
[INFO] LIS cycle time: 05/15/1983 22:00:00
[INFO] LIS cycle time: 05/15/1983 23:00:00
[INFO] LIS cycle time: 05/16/1983 00:00:00
 [INFO] Reading greenness map for month            5 from 
 ./lis_input_files/lis_input_010.nc
 [INFO] Reading greenness map for month            6 from 
 ./lis_input_files/lis_input_010.nc
 [INFO] Reading albedo map for month            5
 [INFO] Reading albedo map for month            6
 [INFO] Writing surface model output to:  
 ./output/LIS_ERA5/SURFACEMODEL/198305/LIS_HIST_198305160000.d01.nc
 [INFO] Writing routing model output to:  
 ./output/LIS_ERA5/ROUTING/198305/LIS_HIST_198305160000.d01.nc
[INFO] LIS cycle time: 05/16/1983 01:00:00
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_slv_Nx.19830517.
 nc4
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830517.
 nc4 (for lowest model level fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830517.
 nc4 (for corrected precipitation fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_rad_Nx.19830517.
 nc4
[INFO] LIS cycle time: 05/16/1983 02:00:00
[INFO] LIS cycle time: 05/16/1983 03:00:00
[INFO] LIS cycle time: 05/16/1983 04:00:00
[INFO] LIS cycle time: 05/16/1983 05:00:00
[INFO] LIS cycle time: 05/16/1983 06:00:00
[INFO] LIS cycle time: 05/16/1983 07:00:00
[INFO] LIS cycle time: 05/16/1983 08:00:00
[INFO] LIS cycle time: 05/16/1983 09:00:00
[INFO] LIS cycle time: 05/16/1983 10:00:00
[INFO] LIS cycle time: 05/16/1983 11:00:00
[INFO] LIS cycle time: 05/16/1983 12:00:00
[INFO] LIS cycle time: 05/16/1983 13:00:00
[INFO] LIS cycle time: 05/16/1983 14:00:00
[INFO] LIS cycle time: 05/16/1983 15:00:00
[INFO] LIS cycle time: 05/16/1983 16:00:00
[INFO] LIS cycle time: 05/16/1983 17:00:00
[INFO] LIS cycle time: 05/16/1983 18:00:00
[INFO] LIS cycle time: 05/16/1983 19:00:00
[INFO] LIS cycle time: 05/16/1983 20:00:00
[INFO] LIS cycle time: 05/16/1983 21:00:00
[INFO] LIS cycle time: 05/16/1983 22:00:00
[INFO] LIS cycle time: 05/16/1983 23:00:00
[INFO] LIS cycle time: 05/17/1983 00:00:00
 [INFO] Writing surface model output to:  
 ./output/LIS_ERA5/SURFACEMODEL/198305/LIS_HIST_198305170000.d01.nc
 [INFO] Writing routing model output to:  
 ./output/LIS_ERA5/ROUTING/198305/LIS_HIST_198305170000.d01.nc
[INFO] LIS cycle time: 05/17/1983 01:00:00
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_slv_Nx.19830518.
 nc4
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830518.
 nc4 (for lowest model level fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830518.
 nc4 (for corrected precipitation fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_rad_Nx.19830518.
 nc4
[INFO] LIS cycle time: 05/17/1983 02:00:00
[INFO] LIS cycle time: 05/17/1983 03:00:00
[INFO] LIS cycle time: 05/17/1983 04:00:00
[INFO] LIS cycle time: 05/17/1983 05:00:00
[INFO] LIS cycle time: 05/17/1983 06:00:00
[INFO] LIS cycle time: 05/17/1983 07:00:00
[INFO] LIS cycle time: 05/17/1983 08:00:00
[INFO] LIS cycle time: 05/17/1983 09:00:00
[INFO] LIS cycle time: 05/17/1983 10:00:00
[INFO] LIS cycle time: 05/17/1983 11:00:00
[INFO] LIS cycle time: 05/17/1983 12:00:00
[INFO] LIS cycle time: 05/17/1983 13:00:00
[INFO] LIS cycle time: 05/17/1983 14:00:00
[INFO] LIS cycle time: 05/17/1983 15:00:00
[INFO] LIS cycle time: 05/17/1983 16:00:00
[INFO] LIS cycle time: 05/17/1983 17:00:00
[INFO] LIS cycle time: 05/17/1983 18:00:00
[INFO] LIS cycle time: 05/17/1983 19:00:00
[INFO] LIS cycle time: 05/17/1983 20:00:00
[INFO] LIS cycle time: 05/17/1983 21:00:00
[INFO] LIS cycle time: 05/17/1983 22:00:00
[INFO] LIS cycle time: 05/17/1983 23:00:00
[INFO] LIS cycle time: 05/18/1983 00:00:00
 [INFO] Writing surface model output to:  
 ./output/LIS_ERA5/SURFACEMODEL/198305/LIS_HIST_198305180000.d01.nc
 [INFO] Writing routing model output to:  
 ./output/LIS_ERA5/ROUTING/198305/LIS_HIST_198305180000.d01.nc
[INFO] LIS cycle time: 05/18/1983 01:00:00
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_slv_Nx.19830519.
 nc4
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830519.
 nc4 (for lowest model level fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830519.
 nc4 (for corrected precipitation fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_rad_Nx.19830519.
 nc4
[INFO] LIS cycle time: 05/18/1983 02:00:00
[INFO] LIS cycle time: 05/18/1983 03:00:00
[INFO] LIS cycle time: 05/18/1983 04:00:00
[INFO] LIS cycle time: 05/18/1983 05:00:00
[INFO] LIS cycle time: 05/18/1983 06:00:00
[INFO] LIS cycle time: 05/18/1983 07:00:00
[INFO] LIS cycle time: 05/18/1983 08:00:00
[INFO] LIS cycle time: 05/18/1983 09:00:00
[INFO] LIS cycle time: 05/18/1983 10:00:00
[INFO] LIS cycle time: 05/18/1983 11:00:00
[INFO] LIS cycle time: 05/18/1983 12:00:00
[INFO] LIS cycle time: 05/18/1983 13:00:00
[INFO] LIS cycle time: 05/18/1983 14:00:00
[INFO] LIS cycle time: 05/18/1983 15:00:00
[INFO] LIS cycle time: 05/18/1983 16:00:00
[INFO] LIS cycle time: 05/18/1983 17:00:00
[INFO] LIS cycle time: 05/18/1983 18:00:00
[INFO] LIS cycle time: 05/18/1983 19:00:00
[INFO] LIS cycle time: 05/18/1983 20:00:00
[INFO] LIS cycle time: 05/18/1983 21:00:00
[INFO] LIS cycle time: 05/18/1983 22:00:00
[INFO] LIS cycle time: 05/18/1983 23:00:00
[INFO] LIS cycle time: 05/19/1983 00:00:00
 [INFO] Writing surface model output to:  
 ./output/LIS_ERA5/SURFACEMODEL/198305/LIS_HIST_198305190000.d01.nc
 [INFO] Writing routing model output to:  
 ./output/LIS_ERA5/ROUTING/198305/LIS_HIST_198305190000.d01.nc
[INFO] LIS cycle time: 05/19/1983 01:00:00
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_slv_Nx.19830520.
 nc4
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830520.
 nc4 (for lowest model level fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_flx_Nx.19830520.
 nc4 (for corrected precipitation fields)
 [INFO] Reading MERRA-2 file (bookend,           2  ... 
 ../../datasets/MERRA2/MERRA2_100/Y1983/M05/MERRA2_100.tavg1_2d_rad_Nx.19830520.
 nc4
[INFO] LIS cycle time: 05/19/1983 02:00:00
[INFO] LIS cycle time: 05/19/1983 03:00:00
[INFO] LIS cycle time: 05/19/1983 04:00:00
[INFO] LIS cycle time: 05/19/1983 05:00:00
[INFO] LIS cycle time: 05/19/1983 06:00:00

Following zip file contains restart and LIS config files.
lis_errwat.zip

[zhengzhenk]

Hello, I also encountered the problem "The model is gaining water (ERRWAT is positive)" when running LISF 7.4 with NoahMP3.6 as you had before. I am glad to see that you have solved this problem by increasing the NoahMP3.6's water balance check threshold to 10. May I ask where you modify this threshold parameter？Thank you very much.

[dmocko]

Hi @zhengzhenk

It's generally not advised to increase this water balance check threshold, but instead to identify the other source of your problem.

But if you insist, for diagnostic purposes only, you can increase the value from "0.1" on this line:

LISF/lis/surfacemodels/land/noahmp.3.6/module_sf_noahmplsm_36.F90

Line 1328 in 23a60bc

IF(ABS(ERRWAT) > 0.1) THEN

[zhengzhenk]

Hi @dmocko
Thank you for your valuable advice and I will try other ways to solve this problem. I found that shortening the NoahMP3.6's time step in the lis.config file would solve this problem to some extent, but if the time series of running data is too long, the problem still exists. Thanks again for your advice.

[sudarsan3010]

I am using LISF 7.4.3 with Noah-MP4.0.1. Even after changing the threshold to 10, the error remains the same. I even reduced the timestep to 1 mn, still it dint solve the error.

I am using GDAS forcing (6 hr) and GPM IMERG (at 30 mins).

The error values when running at 1mn:
The model is gaining water (ERRWAT is positive)
ERRWAT = 1.45071098E+26 kg m{-2} timestep{-1}
Water budget problem in NOAHMP LSM
Note: The following floating-point exceptions are signalling: IEEE_UNDERFLOW_FLAG IEEE_DENORMAL

I overcome this error by downloading the GPM IMERG data without checking "Refine Region" checkbox, which I have done before.

[dmocko]

Hi @sudarsan3010

The LIS reader for the GPM IMERG data assumes that the files have NOT been spatially sub-setted. If you try to read in a file that was spatially sub-setted, then the code will not perform correctly, as you have seen.

You will need to download the entire region of the data, but NOT clicking on "refine region". I realize that this will take more disk space and longer to download.