how to modify initial conditions manually in a single point case?

[Tanja]

I am running clm in single point mode.

I would like to modify the initial conditions manually, especially H2OSOI_LIQ, in the finidat-file.

H2OSOI_LIQ has following dimensions:

H2OSOI_LIQ(column,levtot) where columns=46813 and levtot=37.

What does columns refer to in the case?

How do I know which column I have to modify to change the initial H2OSOI_LIQ in my single point CLM model?

 

[oleson]

Columns are a unit in the sub-grid hierarchy of CLM. Columns types can be vegetated, crop, landice, lake, wetland, urban, etc. Column type numbering can be found in the global attributes of a history file:

:ctype_vegetated_or_bare_soil = 1 ;
:ctype_crop = 2 ;
:ctype_crop_noncompete = "2*100+m, m=cft_lb,cft_ub" ;
:ctype_landice = 3 ;
:ctype_landice_multiple_elevation_classes = "4*100+m, m=1,glcnec" ;
:ctype_deep_lake = 5 ;
:ctype_wetland = 6 ;
:ctype_urban_roof = 71 ;
:ctype_urban_sunwall = 72 ;
:ctype_urban_shadewall = 73 ;
:ctype_urban_impervious_road = 74 ;
:ctype_urban_pervious_road = 75 ;

There is column information in the restart file that can help you find which column(s) you are actually using in the single-point simulation:

cols1d_lon, cols1d_lat, cols1d_ixy, cols1d_jxy, cols1d_gridcell_index, cols1d_landunit_index, cols1d_ityplun, etc.

You could use cols1d_lon and cols1d_lat to find the specific column(s) in the single-point simulation, and change H2OSOI_LIQ in the finidat (restart) file.

Alternatively, you could run your single-point simulation for a couple of time steps, modify the resulting restart file, and then use that in your final simulation.