Changeset 21346


Ignore:
Timestamp:
11/04/16 17:45:26 (8 years ago)
Author:
adhikari
Message:

CHG: nightly updated for rotational feedback

Location:
issm/trunk-jpl/test/NightlyRun
Files:
2 edited

Legend:

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

    r21333 r21346  
    1212late=sum(md.mesh.lat(md.mesh.elements),2)/3;
    1313longe=sum(md.mesh.long(md.mesh.elements),2)/3;
    14 pos=find(late <-80);
    15 md.slr.deltathickness(pos)=-100;
     14pos=find(late <-75 & longe <0);
     15md.slr.deltathickness(pos)=-1;
    1616
    1717%elastic loading from love numbers:
     
    6464
    6565%eustatic + rigid + elastic run:
    66 md.slr.rigid=1; md.slr.elastic=1; md.slr.rotation=0;
     66md.slr.rigid=1; md.slr.elastic=1; md.slr.rotation=0;
     67md.cluster=generic('name',oshostname(),'np',1);
     68%md.verbose=verbose('111111111');
    6769md=solve(md,'Sealevelrise');
    6870SnoRotation=md.results.SealevelriseSolution.Sealevel;
     
    7072%eustatic + rigid + elastic + rotation run:
    7173md.slr.rigid=1; md.slr.elastic=1; md.slr.rotation=1;
     74md.cluster=generic('name',oshostname(),'np',1);
     75%md.verbose=verbose('111111111');
    7276md=solve(md,'Sealevelrise');
    73 Srotation=md.results.SealevelriseSolution.Sealevel;
     77SRotation=md.results.SealevelriseSolution.Sealevel;
    7478
    7579%Fields and tolerances to track changes
    7680field_names     ={'noRotation','Rotation'};
    7781field_tolerances={1e-13,1e-13};
    78 field_values={SnoRotation,Srotation};
     82field_values={SnoRotation,SRotation};
  • issm/trunk-jpl/test/NightlyRun/test2110.m

    r21332 r21346  
    1717%love numbers:
    1818nlov=10001;     % horizontal displacements do not work for low degree truncation, e.g., 101
    19 %md.esa.love_h = love_numbers('h','CM'); md.esa.love_h(nlov+1:end)=[];
    20 %md.esa.love_l = love_numbers('l','CM'); md.esa.love_l(nlov+1:end)=[];
    2119md.esa.love_h = love_numbers('h'); md.esa.love_h(nlov+1:end)=[];
    2220md.esa.love_l = love_numbers('l'); md.esa.love_l(nlov+1:end)=[];
Note: See TracChangeset for help on using the changeset viewer.