Shreve's Hydrology model couples with friction

aisundong

I want to couple Shreve's Hydrological model to an ice sheet simulation I successfully ran beforeI first open the hydrology and set the coupling to 4 in the parameterization file but received the following error message.What I can confirm is that I have given the value of md.hydrology.spcwatercolumn. It seems like there is no variable md.hydrologyshreve.spcwatercolumn in the model, and it is impossible to set the value of it.How can I solve this problem?Can you provide examples of successfully coupling Shreve.

The script to run is as follows:
`[path, filename, extention] = fileparts(pwd);
casename=filename;

time_step = 1;
start_time = 0;
run_length = 1000;
output_num = 2;
output_frequency = run_length/(output_num*time_step);
final_time = run_length;

% create domain
domain = ['./Square.exp'];
hinit = 100000;
hmax = 100000;
hmin = 5000;
gradation = 2;
err = 8;

% Gnerate an inital mesh
md = bamg(model,'domain',domain,'hmax',hinit);

antarData = '/home/mankai/toSun/Map05_PALEOMAP_1deg_middle_Late_Miocene_10Ma.nc';
lat = ncread(antarData,'lat');
lon = ncread(antarData,'lon');
topg = ncread(antarData,'z');
% topg = topg';

[md.mesh.lat md.mesh.long] = xy2ll(md.mesh.x, md.mesh.y, -1);
md.geometry.bed = InterpFromGridToMesh(lat,lon,topg,md.mesh.lat,md.mesh.long,0);

% Parameterization
md = parameterize(md,'./geometry.par');

md = parameterize(md,'./massbalance.par');

md = parameterize(md,'./Antarctica.par');

% Extrude
% md = extrude(md,15,1);
md = extrude(md,30,1);

md.hydrology.spcwatercolumn=ones(md.mesh.numberofvertices,1);

% Parameterization
md = parameterize(md,'./Antarctica2.par');

% Set flow equation
md = setflowequation(md,'HO','all');

md.cluster=generic('name',oshostname,'np',16);

% Transient Run
md.inversion.iscontrol=0;
md.transient.ismasstransport=1;
md.transient.isstressbalance=1;
md.transient.isgroundingline=1;
md.transient.ismovingfront=0;
md.transient.isthermal=1;

md.transient.ishydrology=1;

md.timestepping.time_step=time_step;
md.timestepping.start_time=start_time;
md.settings.output_frequency=output_frequency;
md.timestepping.final_time=final_time;
md.transient.requested_outputs={'default','IceVolume'};

md.settings.solver_residue_threshold = 1e-03;
disp('initial complete');
md=solve(md,'Transient');

%md=solve(md,'Hydrology');`

The error is as follows:
`initial complete
checking model consistency
marshalling file geometry.bin
uploading input file and queueing script
launching solution sequence on remote cluster

Ice Sheet System Model (ISSM) version 4.18
(website: http://issm.jpl.nasa.gov contact: issm@jpl.nasa.gov)

[0] ??? Error using ==> ./classes/IoModel.cpp:2866
[0] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[1] ??? Error using ==> ./classes/IoModel.cpp:2866
[1] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[2] ??? Error using ==> ./classes/IoModel.cpp:2866
[2] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[3] ??? Error using ==> ./classes/IoModel.cpp:2866
[3] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[4] ??? Error using ==> ./classes/IoModel.cpp:2866
[4] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[5] ??? Error using ==> ./classes/IoModel.cpp:2866
[5] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[0] ??? Error using ==> ./classes/IoModel.cpp:2866
[0] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[1] ??? Error using ==> ./classes/IoModel.cpp:2866
[1] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[3] ??? Error using ==> ./classes/IoModel.cpp:2866
[3] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[4] ??? Error using ==> ./classes/IoModel.cpp:2866
[4] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[5] ??? Error using ==> ./classes/IoModel.cpp:2866
[5] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[2] ??? Error using ==> ./classes/IoModel.cpp:2866
[2] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[6] ??? Error using ==> ./classes/IoModel.cpp:2866
[6] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[7] ??? Error using ==> ./classes/IoModel.cpp:2866
[7] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[6] ??? Error using ==> ./classes/IoModel.cpp:2866
[6] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[7] ??? Error using ==> ./classes/IoModel.cpp:2866
[7] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[8] ??? Error using ==> ./classes/IoModel.cpp:2866
[8] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[9] ??? Error using ==> ./classes/IoModel.cpp:2866
[9] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[10] ??? Error using ==> ./classes/IoModel.cpp:2866
[10] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[11] ??? Error using ==> ./classes/IoModel.cpp:2866
[11] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[12] ??? Error using ==> ./classes/IoModel.cpp:2866
[12] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[13] ??? Error using ==> ./classes/IoModel.cpp:2866
[13] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[14] ??? Error using ==> ./classes/IoModel.cpp:2866
[14] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file

[15] ??? Error using ==> ./classes/IoModel.cpp:2866
[15] SetFilePointerToData error message: could not find data with name "md.hydrologyshreve.spcwatercolumn" in binary file