Skip to content

Commit

Permalink
plot for thick inversion data added
Browse files Browse the repository at this point in the history
  • Loading branch information
GHilmarG committed Oct 22, 2023
1 parent dde6364 commit a60c470
Showing 1 changed file with 53 additions and 0 deletions.
53 changes: 53 additions & 0 deletions PlotDataForThicknessInversion.m
Original file line number Diff line number Diff line change
@@ -0,0 +1,53 @@












function PlotDataForThicknessInversion(UserVar,CtrlVar,MUA,F,BCs,Priors,Meas,htrue)


narginchk(8,8)



cbar=UaPlots(CtrlVar,MUA,F,Meas.h,FigureTitle="h measured") ;
xlabel("$x$ (km) ") ; ylabel("$y$ (km) ") ; title(cbar,"(m a.s.l.)")
title("Measured Ice thickness",Interpreter="latex")
hold on ; PlotGroundingLines();

cbar=UaPlots(CtrlVar,MUA,F,"-uv-",FigureTitle="measured velocities");
xlabel("$x$ (km) ") ; ylabel("$y$ (km) ")
title("Measured surface velocities",Interpreter="latex")
hold on ; PlotGroundingLines();

cbar=UaPlots(CtrlVar,MUA,F,F.as,FigureTitle="measured mass balance");
xlabel("$x$ (km) ") ; ylabel("$y$ (km) ") ; title(cbar,"(m/yr)",interpreter="latex")
title("Surface mass balance",Interpreter="latex")
hold on ; PlotGroundingLines();

cbar=UaPlots(CtrlVar,MUA,F,F.dsdt,FigureTitle="surface rate of elevation changes");
title("Surface elevation changes",Interpreter="latex")
title(cbar,"(m/yr)",interpreter="latex")
xlabel("$x$ (km) ") ; ylabel("$y$ (km) ")
hold on ; PlotGroundingLines();













end

0 comments on commit a60c470

Please sign in to comment.