Index: sm/trunk/src/m/shared/GaussPenta.m =================================================================== --- /issm/trunk/src/m/shared/GaussPenta.m (revision 6297) +++ (revision ) @@ -1,13 +1,0 @@ -function [num_area_gauss,first_area_gauss_coord,second_area_gauss_coord,third_area_gauss_coord,area_gauss_weights, vert_gauss_coord,vert_gauss_weights]=GaussPenta(area_order,num_vert_gauss) -%GAUSSPENTA - get Gauss point for penta elements -% -% This routine calls the routines GAUSSTRIA and GAUSSSEGMENT -% -% Usage: -% [num_area_gauss,first_area_gauss_coord,second_area_gauss_coord,third_area_gauss_coord,area_gauss_weights, vert_gauss_coord,vert_gauss_weights]=GaussPenta(area_order,num_vert_gauss) - -%Get Tria area gauss points -[num_area_gauss,first_area_gauss_coord,second_area_gauss_coord,third_area_gauss_coord,area_gauss_weights]=GaussTria(area_order); - -%Get vertical segment gauss points -[vert_gauss_coord vert_gauss_weights]=GaussSegment(num_vert_gauss); Index: sm/trunk/src/m/shared/GaussQuad.m =================================================================== --- /issm/trunk/src/m/shared/GaussQuad.m (revision 6297) +++ (revision ) @@ -1,36 +1,0 @@ -function [num_gauss,first_gauss_coord,second_gauss_coord,gauss_weights]=GaussQuad(order) -%GAUSSQUAD - get Gauss point for quad elements -% -% This routine computes gaussian points on a reference qud elements -% The number of gaussian points returned depends on the order of integration ('order'). The order of -% integration can be computed from the polynomial degree p that needs to be integrated. The formula is: -% order = (p+1) /2 -% order=1, num_gauss=1. Can integrate polynomials of degree 0 to 1 -% -% Usage: -% [num_gauss,first_gauss_coord,second_gauss_coord,gauss_weights]=GaussQuad(order) - - -if order==1, - num_gauss=1; - gauss_weights=4; - first_gauss_coord= 0; - second_gauss_coord= 0; - -%order=2, num_gauss=2. Can integrate polynomials of degree 0 to 3 -elseif order==2, - num_gauss=4; - gauss_weights=[ 1, 1, 1, 1]; - first_gauss_coord=[ -0.57735026918962, 0.57735026918962, -0.57735026918962, 0.57735026918962]; - second_gauss_coord=[0.57735026918962, 0.57735026918962, -0.57735026918962, -0.57735026918962]; - -%order=3, num_gauss=3. Can integrate polynomials of degree 0 to 5 -elseif order==3, - num_gauss=9; - gauss_weights=[ 0.30864197530864, 0.49382716049383, 0.30864197530864, 0.49382716049383, 0.79012345679012, 0.49382716049383, 0.30864197530864, 0.49382716049383, 0.30864197530864]; - first_gauss_coord=[ -0.77459666924148, 0, 0.77459666924148, -0.77459666924148, 0, 0.77459666924148, -0.77459666924148, 0, 0.77459666924148]; - second_gauss_coord=[ -0.77459666924148, -0.77459666924148, -0.77459666924148, 0, 0, 0, 0.77459666924148, 0.77459666924148, 0.77459666924148]; - -else - error('GaussQuad error message: order not supported yet'); -end Index: sm/trunk/src/m/shared/GaussSegment.m =================================================================== --- /issm/trunk/src/m/shared/GaussSegment.m (revision 6297) +++ (revision ) @@ -1,41 +1,0 @@ -function [segment_coord,gauss_weights]=GaussSegment(numgauss) -%GAUSSSEGMENT - get Gauss point for segment elements -% -% This routine computes numgauss gaussian points on a reference segment from -1 to 1. -% The number of gaussian points returned depends on the order of integration ('order'). The order of -% integration can be computed from the polynomial degree p that needs to be integrated. The formula is: -% order = (p+1) /2 -% order=1, num_gauss=1. Can integrate polynomials of degree 0 to 1 -% -% Usage: -% [segment_coord,gauss_weights]=GaussSegment(numgauss) - -if numgauss==1, - - segment_coord=0; - gauss_weights=2; - -%p=3, numgauss=2. Can integrate polynomials of degree 0 to 3 -elseif numgauss==2, - gauss_weights=[ 1.0, 1.0]; - segment_coord=[-0.577350269189626, 0.577350269189626]; - -%p=5, numgauss=3. Can integrate polynomials of degree 0 to 5 -elseif numgauss==3, - gauss_weights=[ 0.555555555555555556, 0.88888888888889, 0.555555555555555556]; - segment_coord=[ -0.774596669241483, 0.000000000000000, 0.774596669241483]; - -%p=7, numgauss=4. Can integrate polynomials of degree 0 to 7 -elseif numgauss==4, - gauss_weights=[ 0.347854845137454, 0.652145154862546, 0.652145154862546, 0.347854845137454]; - segment_coord=[ -0.861136311594053, -0.339981043584856, 0.339981043584856, 0.861136311594053]; - -%p=9, numgauss=5. Can integrate polynomials of degree 0 to 9 -elseif numgauss==5, - gauss_weights=[ 0.236926885056189, 0.478628670499366, 0.568888888888889, 0.478628670499366, 0.236926885056189]; - segment_coord=[ -0.906179845938664, -0.538469310105683, 0, 0.538469310105683, 0.906179845938664]; - - -else - error('GaussSegment error message: order not supported yet'); -end Index: sm/trunk/src/m/shared/GaussTetra.m =================================================================== --- /issm/trunk/src/m/shared/GaussTetra.m (revision 6297) +++ (revision ) @@ -1,128 +1,0 @@ -function [num_gauss,first_gauss_vol_coord,second_gauss_vol_coord,third_gauss_vol_coord,fourth_gauss_vol_coord,gauss_weights]=GaussTetra(order) -%GAUSSTETRA - get Gauss point for tetra elements -% -% This routine computes gaussian points on a reference tetra, in terms of their volume coordinates. -% The number of gaussian points returned depends on the order of integration ('order'). The order of -% integration can be computed from the polynomial degree p that needs to be integrated. The formula is: -% order = (p+1) /2 -% order=1, num_gauss=1. Can integrate polynomials of degree 0 to 1 -% -% Usage: -% [num_gauss,first_gauss_vol_coord,second_gauss_vol_coord,third_gauss_vol_coord,fourth_gauss_vol_coord,gauss_weights]=GaussTetra(order) - -if order==1, - num_gauss=1; - gauss_weights=2*sqrt(2)/3*1; - first_gauss_vol_coord=0.250000000000000; - second_gauss_vol_coord=0.250000000000000; - third_gauss_vol_coord=0.250000000000000; - fourth_gauss_vol_coord=0.250000000000000; - -elseif order==2, - - num_gauss=4; - gauss_weights=2*sqrt(2)/3*[0.250000000000000, 0.250000000000000, 0.250000000000000, 0.250000000000000]; - first_gauss_vol_coord=[ 0.585410196624969, 0.138196601125011, 0.138196601125011, 0.138196601125011]; - second_gauss_vol_coord=[ 0.138196601125011, 0.585410196624969, 0.138196601125011, 0.138196601125011]; - third_gauss_vol_coord=[ 0.138196601125011, 0.138196601125011, 0.585410196624969, 0.138196601125011]; - fourth_gauss_vol_coord=[0.138196601125011, 0.138196601125011, 0.138196601125011, 0.585410196624969]; - -elseif order==3, - - num_gauss=5; - gauss_weights=2*sqrt(2)/3*[ -0.800000000000000, 0.450000000000000, 0.450000000000000, 0.450000000000000, 0.450000000000000]; - first_gauss_vol_coord=[ 0.250000000000000, 0.500000000000000, 0.166666666666667,0.166666666666667, 0.166666666666667]; - second_gauss_vol_coord=[ 0.250000000000000, 0.166666666666667, 0.500000000000000,0.166666666666667, 0.166666666666667]; - third_gauss_vol_coord=[ 0.250000000000000, 0.166666666666667, 0.166666666666667,0.500000000000000, 0.166666666666667]; - fourth_gauss_vol_coord=[0.250000000000000, 0.166666666666667, 0.166666666666667,0.166666666666667, 0.500000000000000]; - -elseif order==4, - - num_gauss=11; - gauss_weights=2*sqrt(2)/3*[-0.013155555555556, 0.007622222222222, 0.007622222222222,0.007622222222222, 0.007622222222222, 0.024888888888889, 0.024888888888889, 0.024888888888889, 0.024888888888889, 0.024888888888889, 0.024888888888889]; - first_gauss_vol_coord=[ 0.250000000000000, 0.785714285714286, 0.071428571428571,0.071428571428571, 0.071428571428571, 0.399403576166799,0.399403576166799, 0.399403576166799, 0.100596423833201, 0.100596423833201, 0.100596423833201]; - second_gauss_vol_coord=[0.250000000000000, 0.071428571428571, 0.785714285714286,0.071428571428571, 0.071428571428571, 0.399403576166799,0.100596423833201, 0.100596423833201, 0.399403576166799, 0.399403576166799, 0.100596423833201]; - third_gauss_vol_coord=[ 0.250000000000000, 0.071428571428571, 0.071428571428571,0.785714285714286, 0.071428571428571, 0.100596423833201,0.399403576166799, 0.100596423833201, 0.399403576166799, 0.100596423833201, 0.399403576166799]; - fourth_gauss_vol_coord=[0.250000000000000, 0.071428571428571, 0.071428571428571,0.071428571428571, 0.785714285714286, 0.100596423833201,0.100596423833201, 0.399403576166799, 0.100596423833201, 0.399403576166799, 0.399403576166799]; - - -elseif order==5, - - num_gauss=15; - gauss_weights=2*sqrt(2)/3*[0.030283678097089, 0.006026785714286, 0.006026785714286,... - 0.006026785714286, 0.006026785714286, 0.011645249086029,... - 0.011645249086029, 0.011645249086029, 0.011645249086029,... - 0.010949141561386, 0.010949141561386, 0.010949141561386,... - 0.010949141561386, 0.010949141561386, 0.010949141561386]; - first_gauss_vol_coord=[ 0.250000000000000, 0.000000000000000, 0.333333333333333,... - 0.333333333333333, 0.333333333333333, 0.727272727272727,... - 0.090909090909091, 0.090909090909091, 0.090909090909091,... - 0.066550153573664, 0.066550153573664, 0.066550153573664,... - 0.433449846426336, 0.433449846426336, 0.433449846426336]; - second_gauss_vol_coord=[0.250000000000000, 0.333333333333333, 0.000000000000000,... - 0.333333333333333, 0.333333333333333, 0.090909090909091,... - 0.727272727272727, 0.090909090909091, 0.090909090909091,... - 0.066550153573664, 0.433449846426336, 0.433449846426336,... - 0.066550153573664, 0.066550153573664, 0.433449846426336]; - third_gauss_vol_coord=[0.250000000000000, 0.333333333333333, 0.333333333333333,... - 0.000000000000000, 0.333333333333333, 0.090909090909091,... - 0.090909090909091, 0.727272727272727, 0.090909090909091,... - 0.433449846426336, 0.066550153573664, 0.433449846426336,... - 0.066550153573664, 0.433449846426336, 0.066550153573664]; - fourth_gauss_vol_coord=[ 0.250000000000000, 0.333333333333333, 0.333333333333333,... - 0.333333333333333, 0.000000000000000, 0.090909090909091,... - 0.090909090909091, 0.090909090909091, 0.727272727272727,... - 0.433449846426336, 0.433449846426336, 0.066550153573664,... - 0.433449846426336, 0.066550153573664, 0.066550153573664]; - - -elseif order==6, - - num_gauss=24; - gauss_weights=2*sqrt(2)/3*[ 0.006653791709695, 0.006653791709695, 0.006653791709695,... - 0.006653791709695, 0.001679535175887, 0.001679535175887,... - 0.001679535175887, 0.001679535175887, 0.009226196923942,... - 0.009226196923942, 0.009226196923942, 0.009226196923942,... - 0.008035714285714, 0.008035714285714, 0.008035714285714,... - 0.008035714285714, 0.008035714285714, 0.008035714285714,... - 0.008035714285714, 0.008035714285714, 0.008035714285714,... - 0.008035714285714, 0.008035714285714, 0.008035714285714]; - first_gauss_vol_coord=[ 0.356191386222545, 0.214602871259152, 0.214602871259152,... - 0.214602871259152, 0.877978124396166, 0.040673958534611,... - 0.040673958534611, 0.040673958534611, 0.032986329573173,... - 0.322337890142276, 0.322337890142276, 0.322337890142276,... - - 0.063661001875018, 0.063661001875018, 0.063661001875018,... - 0.063661001875018, 0.063661001875018, 0.063661001875018,... - 0.269672331458316, 0.603005664791649, 0.269672331458316,... - 0.603005664791649, 0.269672331458316, 0.603005664791649]; - second_gauss_vol_coord=[0.214602871259152, 0.356191386222545, 0.214602871259152, ... - 0.214602871259152, 0.040673958534611, 0.877978124396166, ... - 0.040673958534611, 0.040673958534611, 0.322337890142276, ... - 0.032986329573173, 0.322337890142276, 0.322337890142276,... - - 0.063661001875018, 0.063661001875018, 0.269672331458316, ... - 0.603005664791649, 0.269672331458316, 0.603005664791649, ... - 0.063661001875018, 0.063661001875018, 0.063661001875018, ... - 0.063661001875018, 0.603005664791649, 0.269672331458316]; - third_gauss_vol_coord=[ 0.214602871259152, 0.214602871259152, 0.356191386222545, ... - 0.214602871259152, 0.040673958534611, 0.040673958534611, ... - 0.877978124396166, 0.040673958534611, 0.322337890142276, ... - 0.322337890142276, 0.032986329573173, 0.322337890142276,... - - 0.269672331458316, 0.603005664791649, 0.063661001875018, ... - 0.063661001875018, 0.603005664791649, 0.269672331458316,... - 0.063661001875018, 0.063661001875018, 0.603005664791649,... - 0.269672331458316, 0.063661001875018, 0.063661001875018]; - fourth_gauss_vol_coord=[0.214602871259152, 0.214602871259152, 0.214602871259152,... - 0.356191386222545, 0.040673958534611, 0.040673958534611,... - 0.040673958534611, 0.877978124396166, 0.322337890142276,... - 0.322337890142276, 0.322337890142276, 0.032986329573173,... - - 0.603005664791649, 0.269672331458316, 0.603005664791649,... - 0.269672331458316, 0.063661001875018, 0.063661001875018,... - 0.603005664791649, 0.269672331458316, 0.063661001875018,... - 0.063661001875018, 0.063661001875018, 0.063661001875018]; -else - error('GaussTria error message: order not supported yet'); -end Index: sm/trunk/src/m/shared/GaussTria.m =================================================================== --- /issm/trunk/src/m/shared/GaussTria.m (revision 6297) +++ (revision ) @@ -1,67 +1,0 @@ -function [num_gauss,first_gauss_area_coord,second_gauss_area_coord,third_gauss_area_coord,gauss_weights]=GaussTria(order) -%GAUSSTRIA - get Gauss point for triangle elements -% -% This routine computes gaussian points on a reference triangle, in terms of their area coordinates. -% The number of gaussian points returned depends on the order of integration ('order'). The order of -% integration can be computed from the polynomial degree p that needs to be integrated. The formula is: -% order = (p+1) /2 -% order=1, num_gauss=1. Can integrate polynomials of degree 0 to 1 -% -% Usage: -% [num_gauss,first_gauss_area_coord,second_gauss_area_coord,third_gauss_area_coord,gauss_weights]=GaussTria(order) - -if order==1, - num_gauss=1; - gauss_weights=1.732050807568877; - first_gauss_area_coord= 0.333333333333333; - second_gauss_area_coord= 0.333333333333333; - third_gauss_area_coord= 0.333333333333333; - -%order=2, num_gauss=3. Can integrate polynomials of degree 0 to 3 -elseif order==2, - - num_gauss=3; - gauss_weights=[ 0.577350269189625, 0.577350269189625, 0.577350269189625]; - first_gauss_area_coord=[ 0.666666666666667, 0.166666666666667, 0.166666666666667]; - second_gauss_area_coord=[ 0.166666666666667, 0.666666666666667, 0.166666666666667]; - third_gauss_area_coord=[ 0.166666666666667, 0.166666666666667, 0.666666666666667]; - -%order=3, num_gauss=4. Can integrate polynomials of degree 0 to 5 -elseif order==3, - num_gauss=4; - gauss_weights=[ -0.974278579257493, 0.902109795608790, 0.902109795608790, 0.902109795608790]; - first_gauss_area_coord=[ 0.333333333333333, 0.600000000000000, 0.200000000000000, 0.200000000000000]; - second_gauss_area_coord=[ 0.333333333333333, 0.200000000000000, 0.600000000000000, 0.200000000000000]; - third_gauss_area_coord=[ 0.333333333333333, 0.200000000000000, 0.200000000000000, 0.600000000000000]; - -%order=4, num_gauss=6. Can integrate polynomials of degree 0 to 7 -elseif order==4, - num_gauss=6; - gauss_weights=[ 0.386908262797819, 0.386908262797819, 0.386908262797819, 0.190442006391807, 0.190442006391807, 0.190442006391807]; - first_gauss_area_coord=[ 0.108103018168070, 0.445948490915965, 0.445948490915965, 0.816847572980459, 0.091576213509771, 0.091576213509771]; - second_gauss_area_coord=[ 0.445948490915965, 0.108103018168070, 0.445948490915965, 0.091576213509771, 0.816847572980459, 0.091576213509771]; - third_gauss_area_coord=[ 0.445948490915965, 0.445948490915965, 0.108103018168070, 0.091576213509771, 0.091576213509771, 0.816847572980459]; - -%order=5, num_gauss=7. Can integrate polynomials of degree 0 to 9 -elseif order==5, - num_gauss=7; - - gauss_weights=[ 0.389711431702997, 0.229313399254729, 0.229313399254729, 0.229313399254729, 0.218133059367230, 0.218133059367230, 0.218133059367230]; - first_gauss_area_coord=[ 0.333333333333333, 0.059715871789770, 0.470142064105115, 0.470142064105115, 0.797426985353087, 0.101286507323456, 0.101286507323456]; - second_gauss_area_coord=[ 0.333333333333333, 0.470142064105115, 0.059715871789770, 0.470142064105115, 0.101286507323456, 0.797426985353087, 0.101286507323456]; - third_gauss_area_coord=[ 0.333333333333333, 0.470142064105115, 0.470142064105115, 0.059715871789770, 0.101286507323456, 0.101286507323456, 0.797426985353087]; - -%order=6, num_gauss=12. Can integrate polynomials of degree 0 to 11 -elseif order==6, - num_gauss=12; - gauss_weights=[ 0.202279763184836, 0.202279763184836, 0.202279763184836, 0.088065961139281, 0.088065961139281, 0.088065961139281, ... - 0.143502272432755, 0.143502272432755, 0.143502272432755, 0.143502272432755, 0.143502272432755, 0.143502272432755]; - first_gauss_area_coord=[ 0.501426509658179, 0.249286745170910, 0.249286745170910, 0.873821971016996, 0.063089014491502, 0.063089014491502, ... - 0.053145049844817, 0.053145049844817, 0.310352451033784, 0.636502499121399, 0.310352451033784, 0.636502499121399]; - second_gauss_area_coord=[ 0.249286745170910, 0.501426509658179, 0.249286745170910, 0.063089014491502, 0.873821971016996, 0.063089014491502, ... - 0.310352451033784, 0.636502499121399, 0.053145049844817, 0.053145049844817, 0.636502499121399, 0.310352451033784]; - third_gauss_area_coord=[ 0.249286745170910, 0.249286745170910, 0.501426509658179, 0.063089014491502, 0.063089014491502, 0.873821971016996, ... - 0.636502499121399, 0.310352451033784, 0.636502499121399, 0.310352451033784, 0.053145049844817, 0.053145049844817]; -else - error('GaussTria error message: order not supported yet'); -end Index: /issm/trunk/src/m/shared/GetVerbosityLevel.m =================================================================== --- /issm/trunk/src/m/shared/GetVerbosityLevel.m (revision 6298) +++ /issm/trunk/src/m/shared/GetVerbosityLevel.m (revision 6298) @@ -0,0 +1,18 @@ +function level=GetVerbosityLevel(level_in) +%GETVERBOSITYLEVEL - get verbosity level +% +% Usage: +% level=GetVerbosityLevel(); + +persistent verbositylevel; + +if nargin==1 & nargout==0, + verbositylevel=level_in; +elseif nargin==0 & nargout==1, + if isempty(verbositylevel), error('verbositylevel has not been set. Call SetVerbositylevel first!'); end + level=verbositylevel; + end +else + help GetVerbosityLevel + error('Bad usage'); +end Index: sm/trunk/src/m/shared/Getalpha2.m =================================================================== --- /issm/trunk/src/m/shared/Getalpha2.m (revision 6297) +++ (revision ) @@ -1,49 +1,0 @@ -function alpha2=Getalpha2(parameters) -%GETALPHA2 - compute the friction parameter alpha2 -% -% This routine calculates the basal friction coefficient -% o alpha2= drag^2 * Neff ^r * vel ^s -% with -% o Neff=rho_ice*g*thickness+rho_ice*g*bed, -% o r=q/p and s=1/p -% -% Usage: -% alpha2=Getalpha2(parameters) - -%recover parameters -gravity=parameters.g; -rho_ice=parameters.rho_ice; -rho_water=parameters.rho_water; -K=parameters.k; -bed=parameters.b; -thickness=parameters.h; -velocity_x=parameters.velocities(:,1); -velocity_y=parameters.velocities(:,2); -if size(parameters.velocities,2)==3, %we have to consider vz if Stokes model - velocity_z=parameters.velocities(:,3); -end - -%From bed and thickness, compute effective pressure when drag is viscous: -Neff=gravity*(rho_ice*thickness+rho_water*bed); - -%If effective pressure becomes negative, sliding becomes unstable (Paterson 4th edition p 148). This is because -%the water pressure is so high, the ice sheet elevates over its ice bumps and slides. But the limit behaviour -%for this should be an ice shelf sliding (no basal drag). Therefore, for any effective pressure Neff < 0, we should -%replace it by Neff=0 (ie, equival it to an ice shelf)*/ -pos=find(Neff<0); -Neff(pos)=0; - -%recover p and q in the drag law u = k * sigma ^ p * Neff ^ q : -pcoeff=parameters.p; -qcoeff=parameters.q; -rcoeff=qcoeff./pcoeff; -scoeff=1./pcoeff; - -%We need the velocity magnitude to evaluate the basal stress: -if size(parameters.velocities,2)==2, - velocity_mag=sqrt(velocity_x.^2+velocity_y.^2); -elseif size(parameters.velocities,2)==3, - velocity_mag=sqrt(velocity_x.^2+velocity_y.^2+velocity_z.^2); -end - -alpha2=(K.^2).*(Neff.^rcoeff).*velocity_mag.^(scoeff-1); Index: sm/trunk/src/m/shared/MatlabSolver.m =================================================================== --- /issm/trunk/src/m/shared/MatlabSolver.m (revision 6297) +++ (revision ) @@ -1,22 +1,0 @@ -function X=Solver(A,B,solver_type); -%ICESOLVER - solve the matrix equation AX=B -% -% Solver AX=B, with Chol, Lu or general solvers from matlab. -% We can use either the LU or the Cholesky decomposition, but the -% Cholesky decomposition is twice as efficient as LU for symmetric -% definite positive matrix -% -% Usage: -% X=IceSolver(A,B,solver_type); - -if strcmpi(solver_type,'lu'), - % Solve by LU decomposition. - [L,U] = lu(A); - X = U\(L\B); -elseif strcmpi(solver_type,'cholesky'), - % Solve by Choleski decomposition. - L = chol(A); X = L\(L'\B); -else - % use matlab's generic solver - X = A\B; -end Index: /issm/trunk/src/m/shared/SetVerbosityLevel.m =================================================================== --- /issm/trunk/src/m/shared/SetVerbosityLevel.m (revision 6298) +++ /issm/trunk/src/m/shared/SetVerbosityLevel.m (revision 6298) @@ -0,0 +1,10 @@ +function SetVerbosityLevel(level) +%SETVERBOSITYLEVEL - Set verbosity level +% +% Usage: +% SetVerbosityLevel(level) + + %Set persistent variable in GetVerbositylevel + GetVerbosityLevel(level); + +end Index: sm/trunk/src/m/shared/recover_input.m =================================================================== --- /issm/trunk/src/m/shared/recover_input.m (revision 6297) +++ (revision ) @@ -1,22 +1,0 @@ -function [field,fieldispresent]=recover_input(inputs,fieldname) -%RECOVER_INPUT - recover the field of a structure -% -% From a certain structure, recover the field that corresponds to a certain -% field name. For example: field=recover_input(input,'bof') where input=struct('bof',4,'baf',5) -% will return the value 4, and 1 for fieldispresent. If 'bof' field did not exist, field would be -% NaN and fieldispresent would be 0. -% -% Usage: -% [field,fieldispresent]=recover_input(inputs,fieldname) - -if isfield(inputs,fieldname), - field=getfield(inputs,fieldname); - if ~isempty(field), - fieldispresent=1; - else - fieldispresent=0; - end -else - field=NaN; - fieldispresent=0; -end