Changeset 18974
- Timestamp:
- 01/05/15 10:29:41 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test3300.m
r18973 r18974 22 22 md.hydrology.mask_eplactive_node=0*ones(md.mesh.numberofvertices,1); 23 23 24 pos=find(md.mesh.y==0 & md.mesh.x<6.0e5 & md.mesh.x>4.0e5); 25 md.hydrology.spcepl_head(pos)=0.0; 26 27 28 md.hydrology.epl_conductivity=1000; 24 md.hydrology.epl_conductivity=1.5e-2; 29 25 md.hydrology.epl_initial_thickness=1.0; 30 md.hydrology.epl_colapse_thickness=1.0e- 5;26 md.hydrology.epl_colapse_thickness=1.0e-6; 31 27 md.hydrology.epl_max_thickness=5.0; 32 28 33 29 md.hydrology.transfer_flag=1.0; 34 md.hydrology.leakage_factor= 1.0e-1;30 md.hydrology.leakage_factor=3.9e-12; 35 31 36 times=0:0.2: 30.0;32 times=0:0.2:8.0; 37 33 md.basalforcings.groundedice_melting_rate=ones(md.mesh.numberofvertices+1,length(times)); 38 34 39 %for i=1:length(times)40 %if(times(i)<=1.0)41 %md.basalforcings.groundedice_melting_rate(:,i)=1.0;42 %elseif(times(i)<=6.0)43 %md.basalforcings.groundedice_melting_rate(:,i)=-0.2;44 %else45 %md.basalforcings.groundedice_melting_rate(:,i)=0.0;46 %end47 %end48 md.basalforcings.groundedice_melting_rate(:,:)=1.0; 35 for i=1:length(times) 36 if(times(i)<=1.0) 37 md.basalforcings.groundedice_melting_rate(:,i)=1.0; 38 elseif(times(i)<=6.0) 39 md.basalforcings.groundedice_melting_rate(:,i)=-0.2; 40 else 41 md.basalforcings.groundedice_melting_rate(:,i)=0.0; 42 end 43 end 44 49 45 md.basalforcings.groundedice_melting_rate(end,:)=times; 50 46 51 47 52 48 md.timestepping.time_step=0.2; 53 md.timestepping.final_time= 30.0;49 md.timestepping.final_time=8.0; 54 50 55 51 md=solve(md,TransientSolutionEnum());
Note:
See TracChangeset
for help on using the changeset viewer.