CRHM .int file, Initial_State and Final_State

[loganxingfang]

I used the following .prj to demonstrate the Initial_State and Final_State from Borland CRHM.

1. I ran the "MarmotCreek_05-06test_modified_25Jul18_basinaggregationscheme3autoa_generateInt.prj" to generate the final state file "MarmotInitialstate06.int" at the end of run.

This project can be found here:
MarmotCreek_05-06test_modified_25Jul18_basinaggregationscheme3autoa_generateInt.zip

The obs file for the prj can be found here:
Marmot_Hourly_ArrayMetData_withT_g_1Oct05-30Sept20_update_22Jan2021.zip

In this prj, I saved the the final state file name "MarmotInitialstate06.int" under Final_State,

Final_State

MarmotInitialstate06.int

This .int file stores all the state variables in this prj. A state variable is declared using declstatvar(...); the state variables include variables such as Albedo, SWE, and so on.

The purpose of .int file is to show information of these state variables at the end of run, but also it can provide initial values of these state variables used by the subsequent model run.

2. For the subsequent model run that uses the "MarmotInitialstate06.int", I ran the "MarmotCreek_06-07test_modified_25Jul18_basinaggregationscheme3autoa_useInt.prj".

This project can be found here:
MarmotCreek_06-07test_modified_25Jul18_basinaggregationscheme3autoa_useInt.zip

This prj reads the "MarmotInitialstate06.int" file under Initial_State:

Initial_State

MarmotInitialstate06.int

It should be noted that value for Time in the "MarmotInitialstate06.int" and Start value for Date in the "MarmotCreek_06-07test_modified_25Jul18_basinaggregationscheme3autoa_useInt.prj" must be the same:

Time:
2006 10 1

Dates:

2006 10 1
2007 10 1

If these do not match, CRHM will have complaint message.

Note that in the "MarmotCreek_06-07test_modified_25Jul18_basinaggregationscheme3autoa_useInt.prj", I also saved the final state file "MarmotInitialstate07.int" so that it could be used by the subsequent model run starting 2007 10 1.

Note that in above prj files, there is no file directory specified for the .int files under Final_State and Initial_State, this saves .int file in same file directory as where .prj files are and assumes both prj files are in the same file directory.
On the other hand, .int file can be placed at specified file directory that is different from where .prj files are, such as C:\Research\MarmotInitialstate06.int. As long as whereabout of the int file under Final_State and Initial_State are specified correctly, the prj file generating and reading the .int file can run without problem.

[jhs507]

Thank you Logan,

I do not think this feature ever worked in the crhmcode version of crhm.

How high of a priority would you say this feature is? Just trying to get a feel for the urgency.

[loganxingfang]

This can be useful for CRHM users when they have information about state variables from field measurements such as SWE, soil_moist, soil_rechr, z_s, z_s_0, z_s_l, and so on, then the users can use field measured values from a specified date, e.g. 2006 4 15, and save that date in Time value for .int file and start to run a .prj using .int from 2006 4 15. This is quite useful especially when users test some new process algorithms developed from field study.

If this does not work in gcc crhm, would this feature work in vcc crhm?

[jhs507]

Thanks,

I will look into this and let you know what my estimate for how long it will take to get working is.

[loganxingfang]

Thank you Justin.

[jhs507]

Hello Logan,

I have found the relevant code for saving a final state and added it into CRHMcode gcc.

I have it on branch initial_and_final_state_files

Could you pull that branch, make clean, make, and then check if final state files are being saved as you would expect?

I have not yet implemented reading an initial state from such a file.

[loganxingfang]

Thanks Justin.
I compiled crhm using branch initial_and_final_state_files and ran a project. I confirm that gcc crhm now can generate the final state file from a project.

[jhs507]

I took a closer look and reading in an initial state looks like it should work. Could you try a project with a initial state and let me know if it looks like everything is correct?

[loganxingfang]

I took a closer look and reading in an initial state looks like it should work. Could you try a project with a initial state and let me know if it looks like everything is correct?

I tried a project "MarmotCreek_06-07test_modified_25Jul18_basinaggregationscheme3autoa_useInt.prj"; it reads "MarmotInitialstate06.int" file from Initial_State. This is generated by "MarmotCreek_05-06test_modified_25Jul18_basinaggregationscheme3autoa_generateInt.prj".

However, it appears that "MarmotCreek_06-07test_modified_25Jul18_basinaggregationscheme3autoa_useInt.prj" can read the .int file but cannot finish the run. I have the screen shot below:

C:\Research\GitHub\crhm\10May2021>crhm MarmotCreek_06-07test_modified_25Jul18_basinaggregationscheme3autoa_useInt.prj
Project file: MarmotCreek_06-07test_modified_25Jul18_basinaggregationscheme3autoa_useInt.prj

Observation file: E:\Research\Model\CRHM\CRHM Tests\CRHM Marmot Creek\Marmot Creek Hydraulic Conductivity Macro Test\CRHM_OBS\Marmot_Hourly_ArrayMetData_withT_g_1Oct05-30Sept20_update_22Jan2021.obs

State file: MarmotInitialstate06.int

Time of model run: 44326.5 2021 5 10. Program Version: 4.7_16

C:\Research\GitHub\crhm\10May2021>

[jhs507]

Where is the MarmotInitialstate06.int file located on your system? To me it looks like it wouldn't be found unless it was in the current working directory.

[loganxingfang]

The MarmotInitialstate06.int file is in the same file directory as the project "MarmotCreek_06-07test_modified_25Jul18_basinaggregationscheme3autoa_useInt.prj".

So, under Initial_State, I did not include the file directory.

Initial_State

MarmotInitialstate06.int

Anyway, to be more sure, I specify the file directory now in a project "MarmotCreek_06-07test_modified_25Jul18_basinaggregationscheme3autoa_useInt1.prj"

Initial_State

C:\Research\GitHub\crhm\10May2021\MarmotInitialstate06.int

gcc crhm had same outcome; see the screen shot below,

C:\Research\GitHub\crhm\10May2021>crhm MarmotCreek_06-07test_modified_25Jul18_basinaggregationscheme3autoa_useInt1.prj
Project file: MarmotCreek_06-07test_modified_25Jul18_basinaggregationscheme3autoa_useInt1.prj

Observation file: E:\Research\Model\CRHM\CRHM Tests\CRHM Marmot Creek\Marmot Creek Hydraulic Conductivity Macro Test\CRHM_OBS\Marmot_Hourly_ArrayMetData_withT_g_1Oct05-30Sept20_update_22Jan2021.obs

State file: C:\Research\GitHub\crhm\10May2021\MarmotInitialstate06.int

Time of model run: 44326.6 2021 5 10. Program Version: 4.7_16

C:\Research\GitHub\crhm\10May2021>

Borland CRHM has no problem to locate the .int file and run both projects above.

[jhs507]

So to be clear for you it crashes when you run a project that tries to set an initial state?

I have been able to run MarmotCreek_06-07test_modified_25Jul18_basinaggregationscheme3autoa_useInt.prj that you provided above and I get some error messages but the run completes.

[loganxingfang]

So to be clear for you it crashes when you run a project that tries to set an initial state?

MarmotCreek_06-07test_modified_25Jul18_basinaggregationscheme3autoa_useInt.prj is the project that reads the initial state file "MarmotInitialstate06.int". In above, I just tried two options for the project to locate and read this initial state file. In both cases, the screen shot information shows that gcc crhm could read the initial state file but did not finish run.

[loganxingfang]

Hi Justin,
I am also providing some information here about how the initial state file works with the project that reads it.

First, for the initial state file "MarmotInitialstate06.int", there is Time on line 3, value is 2006 10 1 on line4; this must match the Start value under Dates: 2006 10 1 on line 131 in "MarmotCreek_06-07test_modified_25Jul18_basinaggregationscheme3autoa_useInt.prj".

Then, values of those state variables in "MarmotInitialstate06.int" must be assigned to those variables in "MarmotCreek_06-07test_modified_25Jul18_basinaggregationscheme3autoa_useInt.prj" on the first time step of run.