Changeset 18974


Ignore:
Timestamp:
01/05/15 10:29:41 (10 years ago)
Author:
bdef
Message:

CHG: reverting nightly 3300

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/test/NightlyRun/test3300.m

    r18973 r18974  
    2222md.hydrology.mask_eplactive_node=0*ones(md.mesh.numberofvertices,1);
    2323
    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;
     24md.hydrology.epl_conductivity=1.5e-2;
    2925md.hydrology.epl_initial_thickness=1.0;
    30 md.hydrology.epl_colapse_thickness=1.0e-5;
     26md.hydrology.epl_colapse_thickness=1.0e-6;
    3127md.hydrology.epl_max_thickness=5.0;
    3228
    3329md.hydrology.transfer_flag=1.0;
    34 md.hydrology.leakage_factor=1.0e-1;
     30md.hydrology.leakage_factor=3.9e-12;
    3531
    36 times=0:0.2:30.0;
     32times=0:0.2:8.0;
    3733md.basalforcings.groundedice_melting_rate=ones(md.mesh.numberofvertices+1,length(times));
    3834
    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 %       else
    45 %               md.basalforcings.groundedice_melting_rate(:,i)=0.0;
    46 % end
    47 % end   
    48 md.basalforcings.groundedice_melting_rate(:,:)=1.0;
     35for 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
     43end     
     44
    4945md.basalforcings.groundedice_melting_rate(end,:)=times;
    5046
    5147
    5248md.timestepping.time_step=0.2;
    53 md.timestepping.final_time=30.0;
     49md.timestepping.final_time=8.0;
    5450
    5551md=solve(md,TransientSolutionEnum());
Note: See TracChangeset for help on using the changeset viewer.