Changeset 25268

Timestamp:

07/10/20 23:08:28 (5 years ago)

Author:

Eric.Larour

Message:

CHG: diverse

Location:

issm/trunk-jpl/test/SandBox

Files:

• Makefile (modified) (1 diff)
• hier2002.m (modified) (6 diffs)
• partmme2002.m (modified) (4 diffs)
• scale2002_element.m (added)
• scale2002_vertices.m (added)

issm/trunk-jpl/test/SandBox/Makefile ¶

@@ -1 +1 @@
 clean: 
-        rm -rf *.tar.gz *.errlog *.outlog *.bin *.qmu.in *.qmu.out *qmu.err *.queue *.toolkits *.out
+        rm -rf *.tar.gz *.errlog *.outlog *.bin *.qmu.in *.qmu.out *qmu.err *.queue *.toolkits *.out *.outbin.*

issm/trunk-jpl/test/SandBox/hier2002.m ¶

@@ -3 +3 @@
 %mesh earth:
 md=model;
-md.mesh=gmshplanet('radius',6.371012*10^3,'resolution',1000.); %500 km resolution mesh
-md.cluster.np=16;
+md.mesh=gmshplanet('radius',6.371012*10^3,'resolution',700.); %500 km resolution mesh
+md.cluster.np=2;
 
 %parameterize slr solution:
 %slr loading:  {{{
-nt=100;
-md.solidearth.surfaceload.icethicknesschange=zeros(md.mesh.numberofelements+1,nt);
+nsteps=100;
+md.solidearth.surfaceload.icethicknesschange=zeros(md.mesh.numberofelements+1,nsteps);
 md.solidearth.surfaceload.icethicknesschange(end,:)=2000.5:1:2099.5;
 md.solidearth.sealevel=zeros(md.mesh.numberofvertices,1);
@@ -24 +24 @@
 
 %greenland
-posgre=find(late > 70 &  late < 80 & longe>-60 & longe<-30);
+posgre=find(late > 70 &  late < 80 & longe>-50 & longe<-30);
 md.solidearth.surfaceload.icethicknesschange(posgre,:)=-100;
 
 %alaska : 
-posala=find(late > 62 &  late < 68 & longe>-152 & longe<-147);
+posala=find(late > 55 &  late < 68 & longe>-162 & longe<-140);
 md.solidearth.surfaceload.icethicknesschange(posala,:)=-150;
 
-%himalaya : 
-poshim=find(late > 25 &  late < 32 & longe>81 & longe<86);
-md.solidearth.surfaceload.icethicknesschange(poshim,:)=-200;
-
+
+%hawaii : 
+poshaw=find(late > 15 &  late < 25 & longe>-170 & longe<-130);
+md.solidearth.surfaceload.icethicknesschange(poshaw,:)=10;
 
 %elastic loading from love numbers:
@@ -78 +78 @@
 %initialize GIA: 
 md.gia=giamme();
-md.gia.Ngia=ones(md.mesh.numberofvertices,1);
-md.gia.Ugia=ones(md.mesh.numberofvertices,1);
+md.gia.Ngia=zeros(md.mesh.numberofvertices,1);
+md.gia.Ugia=zeros(md.mesh.numberofvertices,1);
 md.gia.modelid=1;
 
-%Uncertainty Quantification%{{{
+%Uncertainty Quantification
 md.qmu.variables=struct();;
 
-%partition vector: 
-npart=4;
+%partitioning scaling
+npart=1;
 partition=-ones(md.mesh.numberofelements,1);
-partition(posgre)=0;
-partition(posant)=1;
+partition(posala)=0;
+
+%variable scaling
+md.qmu.variables.surfaceload0=normal_uncertain(...
+        'descriptor','scaled_SurfaceloadIceThicknessChange',...
+        'mean',ones(npart,nsteps),...
+        'stddev',.1*ones(npart,nsteps),...
+        'partition',partition,'nsteps',nsteps);
+
+%partitioning distributed
+npart=3;
+partition=-ones(md.mesh.numberofelements,1);
+
+partition(posant)=0;
+partition(posgre)=1;
 partition(posala)=2;
-partition(poshim)=3;
+
+vposant=unique(md.mesh.elements(posant,:));
+vposgre=unique(md.mesh.elements(posgre,:));
+vposala=unique(md.mesh.elements(posala,:));
+
 
 %prepare arrays: 
-background=md.solidearth.surfaceload.icethicknesschange;
-antarray={}; for i=1:10, antarray{end+1}= (-50+i)*ones(length(posant),nsteps); end
-grearray={}; for i=1:20, grearray{end+1}= (-100+i)*ones(length(posgre),nsteps); end
-alaarray={}; for i=1:30, alaarray{end+1}= (-150+i)*ones(length(posala),nsteps); end
-himarray={-200*ones(length(poshim),nsteps)};
-arrays={antarray,grearray,alaarray,himarray};
-
-partarray={antarray,grearray,alaarray,himarray};
-histarray={equiprobable_histogram_uncertain('ant',10), equiprobable_histogram_uncertain('gre',20), equiprobable_histogram_uncertain('ala',30), equiprobable_histogram_uncertain('him',1)};
-pdfant=normal_uncertain('descriptor','pdfant','mean',1*ones(1,nt),'stddev',.1*ones(1,nt));
-pdfgre=uniform_uncertain('descriptor','pdfgre','lower',.8*ones(1,nt),'upper',1.2*ones(1,nt));
-pdfala=dirac_uncertain('descriptor','pdfala',2);
-pdfhim=dirac_uncertain('descriptor','pdfhim',1);
+nant=10; ngre=20; nala=30; nmax=max([nant,ngre,nala]);
+
+%make ocean levelset temporal:
+md.mask.ocean_levelset=[md.mask.ocean_levelset;2000];
+
+%background arrays: 
+bkgd=md.solidearth.surfaceload.icethicknesschange; 
+bkgdnan=bkgd; bkgdnan(1:end-1,:)=NaN;
+
+bkgdocean=md.mask.ocean_levelset;
+bkgdoceannan=bkgdocean; bkgdoceannan(1:end-1,1)=NaN;
+
+%mme arrays:
+md.solidearth.surfaceload.icethicknesschange=cell(30,1);
+md.mask.ocean_levelset=cell(30,1);
+for i=1:nmax,
+        if i==1,
+                md.solidearth.surfaceload.icethicknesschange{i}=bkgd;
+                md.mask.ocean_levelset{i}=bkgdocean;
+        else
+                bi=bkgdnan;
+                oi=bkgdoceannan;
+                if i<=10, 
+                        bi(posant,:)=(-50+i)*ones(length(posant),nsteps); 
+                        %oi(vposant,:)=bkgdocean(vposant,:)-50+i;
+                        oi(vposant,:)=i;
+                end
+                if i<=20, 
+                        bi(posgre,:)=(-100+i)*ones(length(posgre),nsteps); 
+                        %oi(vposgre,:)=bkgdocean(vposgre,:)-100+i;
+                        oi(vposgre,:)=2*i;
+                end
+                if i<=30, 
+                        bi(posala,:)=(-150+i)*ones(length(posala),nsteps); 
+                        %oi(vposala,:)=bkgdocean(vposala,:)-150+i;
+                        oi(vposala,:)=3*i;
+                end
+                md.solidearth.surfaceload.icethicknesschange{i}=bi;
+                ois=repmat(oi,1,nsteps); ois(end,:)=bi(end,:);
+                md.mask.ocean_levelset{i}=ois;
+        end
+        md.mask.ocean_levelset{i}=bkgdocean;
+end
+%md.mask.ice_levelset=md.mask.ocean_levelset;
+%md.solidearth.surfaceload.icethicknesschange=bkgd;
+%md.mask.ocean_levelset=bkgdocean;
+
+%create distributed histograms for each partition:
+pairs_per_variable=zeros(npart,1);
+abscissas=cell(npart,1);
+counts=cell(npart,1);
+[abscissas{1} counts{1} pairs_per_variable(1)]=equiprobable_histogram_uncertain(10); 
+[abscissas{2} counts{2} pairs_per_variable(2)]=equiprobable_histogram_uncertain(20); 
+[abscissas{3} counts{3} pairs_per_variable(3)]=equiprobable_histogram_uncertain(30); 
 
 %variables: 
-md.qmu.variables.deltathickness=normal_uncertain('descriptor','SurfaceloadIceThicknessChange',...
-        'partition',partition,'nsteps',nt,...
-        'background',md.solidearth.surfaceload.icethicknesschange,...
-        'histogram_bin_uncertain_perpartition',histarray,...
-        'matarray_perpartition',partarray,...
-
-        );
-
-md.qmu.correlation_matrix=zeros(npart*nt,npart*nt);
-for i=1:npart,
-        for j=1:nt,
-                indi=(i-1)*nt+j;
-                for k=1:npart,
-                        for l=1:nt,
-                                indj=(k-1)*nt+l;
-                                if i~=k,
-                                        md.qmu.correlation_matrix(indi,indj)=0;
-                                else
-                                        %same partition:
-                                        if j==l, 
-                                                md.qmu.correlation_matrix(indi,indj)=1;
-                                        else
-                                                md.qmu.correlation_matrix(indi,indj)=.2;
-                                        end
-                                end
-                        end
-                end
-        end
-end
+md.qmu.variables.surfaceload1=histogram_bin_uncertain(...
+        'descriptor','distributed_SurfaceloadModelid',...
+        'pairs_per_variable',pairs_per_variable,'abscissas',abscissas,'counts',counts,'partition',partition);
+
+
 md.qmu.correlation_matrix=[];
-
-
 
 %responses 
@@ -179 +211 @@
 md.qmu.params.tabular_graphics_data=true;
 md.qmu.output=1;
-%}}}
 
 %transient: 
 md.timestepping.start_time=2000;
-md.timestepping.interp_forcings=0;
+md.timestepping.interp_forcings=1;
 md.timestepping.final_time=2002;
 md.transient.issmb=0;
@@ -191 +222 @@
 md.transient.isslr=1;
 md.transient.isgia=1;
-md.slr.requested_outputs= {'default',...
-                'SealevelriseDeltathickness','Sealevel','SealevelRSLRate','SealevelriseCumDeltathickness',...
+md.solidearth.requested_outputs= {'default',...
+                'SurfaceloadIceThicknessChange','Sealevel','SealevelRSLRate',...
                 'SealevelNEsaRate', 'SealevelUEsaRate', 'NGiaRate', 'UGiaRate',...
-                'SealevelEustaticMask','SealevelEustaticOceanMask'};
+                'SealevelEustaticMask','SealevelEustaticOceanMask','MaskOceanLevelset','MaskIceLevelset'};
 
 %hack: 
@@ -203 +234 @@
 %eustatic + rigid + elastic + rotation run:
 md.verbose=verbose('11111111111');
-%md.verbose.qmu=1;
+%md.verbose=verbose(0);
+md.verbose.qmu=1;
 md.solidearth.settings.rigid=1; md.solidearth.settings.elastic=1;md.solidearth.settings.rotation=1;
 md.qmu.isdakota=1;
-md=solve(md,'slr');
-
+md=solve(md,'tr');
+

issm/trunk-jpl/test/SandBox/partmme2002.m ¶

@@ -3 +3 @@
 %mesh earth:
 md=model;
-md.mesh=gmshplanet('radius',6.371012*10^3,'resolution',400.); %500 km resolution mesh
-md.cluster.np=16;
+md.mesh=gmshplanet('radius',6.371012*10^3,'resolution',500.); %500 km resolution mesh
+md.cluster.np=2;
 
 %parameterize slr solution:
@@ -140 +140 @@
                 md.mask.ocean_levelset{i}=ois;
         end
+        md.mask.ocean_levelset{i}=bkgdocean;
 end
+%md.mask.ice_levelset=md.mask.ocean_levelset;
+%md.solidearth.surfaceload.icethicknesschange=bkgd;
+%md.mask.ocean_levelset=bkgdocean;
 
 %create distributed histograms for each partition:
@@ -190 +194 @@
 md.qmu.params.analysis_driver='matlab';
 md.qmu.params.evaluation_scheduling='master';
-md.qmu.params.processors_per_evaluation=5;
+md.qmu.params.processors_per_evaluation=1;
 md.qmu.params.tabular_graphics_data=true;
 md.qmu.output=1;
@@ -216 +220 @@
 
 %eustatic + rigid + elastic + rotation run:
-%md.verbose=verbose('11111111111');
-md.verbose=verbose(0);
+md.verbose=verbose('11111111111');
+%md.verbose=verbose(0);
 md.verbose.qmu=1;
 md.solidearth.settings.rigid=1; md.solidearth.settings.elastic=1;md.solidearth.settings.rotation=1;