Changeset 19043 for issm/trunk-jpl/examples/Pig/Pig.par

Timestamp:

01/27/15 05:10:08 (10 years ago)

Author:

Mathieu Morlighem

Message:

CHG: removed all folds

File:

• issm/trunk-jpl/examples/Pig/Pig.par (modified) (7 diffs)

issm/trunk-jpl/examples/Pig/Pig.par ¶

@@ -7 +7 @@
 md.mesh.hemisphere='s';
 
-% {{{ NetCdf Loading
+%NetCdf Loading
 disp('   Loading SeaRISE data from NetCDF');
 ncdata='../Data/Antarctica_5km_withshelves_v0.75.nc';
@@ -47 +47 @@
 y2=(ymax-ny*spacing)+(0:1:ny)'*spacing; 
 y2=double(y2);
-% }}}
-% {{{ Geometry
+
+%Geometry
 disp('   Interpolating surface and ice base');
 md.geometry.base    = InterpFromGridToMesh(x1,y1,topg,md.mesh.x,md.mesh.y,0);
@@ -73 +73 @@
 md.geometry.thickness(pos0)=1;
 md.geometry.surface=md.geometry.thickness+md.geometry.base;
-% }}} 
-% {{{ Initialization parameters
+
+%Initialization parameters
 disp('   Interpolating temperatures');
 md.initialization.temperature=InterpFromGridToMesh(x1,y1,temp,md.mesh.x,md.mesh.y,0)+273.15+Temp_change;
@@ -93 +93 @@
 md.initialization.pressure=md.materials.rho_ice*md.constants.g*md.geometry.thickness;
 
-% }}}
-
 disp('   Construct ice rheological properties');
 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
 md.materials.rheology_B=paterson(md.initialization.temperature);
 
-% {{{ Forcings
+%Forcings
 disp('   Interpolating surface mass balance');
 md.surfaceforcings.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0);
@@ -108 +106 @@
 md.basalforcings.geothermalflux=InterpFromGridToMesh(x1,y1,gflux,md.mesh.x,md.mesh.y,0);
 clear gflux;
-% }}}
 
-% {{{ Friction and inversion set up
+%Friction and inversion set up
 disp('   Construct basal friction parameters');
 md.friction.coefficient=friction_coefficient*ones(md.mesh.numberofvertices,1);
@@ -124 +121 @@
 md.inversion.vy_obs=vy_obs;
 md.inversion.vel_obs=vel_obs;
-% }}}
-
 
 disp('   Set boundary conditions');
@@ -133 +128 @@
 md.thermal.spctemperature     = [md.initialization.temperature;1]; %impose observed temperature on surface
 md.masstransport.spcthickness    = NaN*ones(md.mesh.numberofvertices,1);
-