Changeset 5569 for issm/trunk

Timestamp:

08/25/10 09:20:45 (15 years ago)

Author:

seroussi

Message:

minor improvements in validation .par

Location:

issm/trunk/test/Par

Files:

• ISMIPC.par (modified) (1 diff)
• ISMIPE.par (modified) (1 diff)

issm/trunk/test/Par/ISMIPC.par

@@ -9 +9 @@
 disp('      creating drag');
 md.drag_type=2; %0 none 1 plastic 2 viscous
-md.drag_coefficient=sqrt(md.yts.*(1000+1000*sin(md.x*2*pi/max(md.x/2)).*sin(md.y*2*pi/max(md.x/2)))./(md.g*(md.rho_ice*md.thickness+md.rho_water*md.bed)));
-%md.drag_coefficient=sqrt(md.yts.*(1000+1000*sin(md.x*2*pi/max(md.x)).*sin(md.y*2*pi/max(md.x)))./(md.g*(md.rho_ice*md.thickness+md.rho_water*md.bed)));
+%md.drag_coefficient=sqrt(md.yts.*(1000+1000*sin(md.x*2*pi/max(md.x/2)).*sin(md.y*2*pi/max(md.x/2)))./(md.g*(md.rho_ice*md.thickness+md.rho_water*md.bed)));
+md.drag_coefficient=sqrt(md.yts.*(1000+1000*sin(md.x*2*pi/max(md.x)).*sin(md.y*2*pi/max(md.x)))./(md.g*(md.rho_ice*md.thickness+md.rho_water*md.bed)));
 %Take care of iceshelves: no basal drag
 pos=find(md.elementoniceshelf);

issm/trunk/test/Par/ISMIPE.par

@@ -22 +22 @@
 md.drag_type=2; %0 none 1 plastic 2 viscous
 md.drag_coefficient=zeros(md.numberofgrids,1);
-%Take care of iceshelves: no basal drag
-pos=find(md.elementoniceshelf);
-md.drag_coefficient(md.elements(pos,:))=0;
 md.drag_p=ones(md.numberofelements,1);
 md.drag_q=ones(md.numberofelements,1);