Changeset 18201
- Timestamp:
- 06/30/14 16:42:56 (11 years ago)
- Location:
- issm/trunk-jpl/examples/Jakobshavn
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/examples/Jakobshavn/Jks.par
r17799 r18201 6 6 %Load SeaRISE dataset 7 7 disp(' Loading SeaRISE data from NetCDF'); 8 ncdata = '. /Greenland_5km_v1.1.nc';8 ncdata = '../Data/Greenland_5km_dev1.2.nc'; 9 9 x1 = ncread(ncdata,'x1'); 10 10 y1 = ncread(ncdata,'y1'); … … 58 58 disp(' Set other boundary conditions'); 59 59 md=SetMarineIceSheetBC(md,'./Front.exp'); 60 md.basalforcings. melting_rate=zeros(md.mesh.numberofvertices,1);60 md.basalforcings.floatingice_melting_rate=zeros(md.mesh.numberofvertices,1); -
issm/trunk-jpl/examples/Jakobshavn/runme.m
r17799 r18201 1 steps=[ 1:3];1 steps=[4]; 2 2 3 3 if any(steps==1) … … 6 6 7 7 %Get observed velocity field on mesh nodes 8 ncdata=' Greenland_5km_v1.1.nc';8 ncdata='../Data/Greenland_5km_dev1.2.nc'; 9 9 if ~exist(ncdata,'file'), 10 error(' Download Greenland5km_v1.1.nc first on http://websrv.cs.umt.edu/isis/index.php/Present_Day_Greenland');10 error('File Greenland_5km_dev1.2.nc not downloaded in Data Directory. Please cd ../Data and run the download script.'); 11 11 end 12 12 x1 = ncread(ncdata,'x1'); … … 84 84 'colorbar#1','off','colorbar#2','on','colorbartitle#2','[m/yr]',... 85 85 'caxis#1-2',[0,7000],... 86 'data',md.geometry.b ed,'title','Bedelevation',...86 'data',md.geometry.base,'title','Base elevation',... 87 87 'data',md.results.StressbalanceSolution.FrictionCoefficient,... 88 88 'title','Friction Coefficient',...
Note:
See TracChangeset
for help on using the changeset viewer.