Changeset 27411

Timestamp:

11/22/22 07:43:43 (2 years ago)

Author:

vverjans

Message:

CHG: adding NightlyRun files as examples of ARMA and Stochastic capabilities

Location:

issm/trunk-jpl/test

Files:

• Archives/Archive543.arch (modified) ( previous)
• Archives/Archive545.arch (added)
• NightlyRun/test543.m (modified) (8 diffs)
• NightlyRun/test543.py (modified) (9 diffs)
• NightlyRun/test545.m (added)
• NightlyRun/test545.py (added)

issm/trunk-jpl/test/NightlyRun/test543.m

@@ -1 @@
-%Test Name: PigTranStochasticforcing
+%Test Name: PigTranRignotarma
 md=triangle(model(),'../Exp/Pig.exp',10000.);
 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp');
@@ -21 +21 @@
     end
 end
-% Basin separation default
-idb_df = zeros(md.mesh.numberofelements,1);
-iid1   = find(md.mesh.x<=-1.62e6);
-for ii=1:md.mesh.numberofelements
-    for vertex=1:3
-        if any(iid1==md.mesh.elements(ii,vertex)) %one vertex in basin 1
-            idb_df(ii) = 1;
-        end
-    end
-    if idb_df(ii)==0 %no vertex was found in basin 1
-        idb_df(ii) = 2;
-    end
-end
 % Dimensionalities
 nb_tf    = 2;
-nb_clv   = 2;
-nb_flmlt = 2;
 
 %Calving parameters
 md.mask.ice_levelset = 1e4*(md.mask.ice_levelset + 0.5);
-md.calving.calvingrate = 0.3*ones(md.mesh.numberofvertices,1);
+md.calving.calvingrate = 0*ones(md.mesh.numberofvertices,1);
 md.levelset.spclevelset = NaN(md.mesh.numberofvertices,1);
 md.levelset.migration_max = 10.0; %avoid fast advance/retreat of the front
@@ -115 +100 @@
 end
 % Floating Ice Melt parameters
-md.basalforcings.floatingice_melting_rate = 0.1*ones(md.mesh.numberofvertices,1);
+md.basalforcings.floatingice_melting_rate = 0*ones(md.mesh.numberofvertices,1);
 
 
 % Covariance matrix
 covtf       = 1e-4*eye(nb_tf);
-covclv      = 1e-1*eye(nb_clv);
-covclv(1,1) = 1/10*covclv(1,1);
-covflmlt    = 0.05*eye(nb_flmlt);
 covsd       = 1e3*eye(nb_tf);
-covglob     = blkdiag(covtf,covclv,covflmlt,covsd);
+covglob     = blkdiag(covtf,covsd);
 
 %Stochastic forcing
 md.stochasticforcing.isstochasticforcing = 1;
-md.stochasticforcing.fields              = [{'FrontalForcingsRignotarma'},{'DefaultCalving'},{'FloatingMeltRate'},{'FrontalForcingsSubglacialDischargearma'}];
+md.stochasticforcing.fields              = [{'FrontalForcingsRignotarma'},{'FrontalForcingsSubglacialDischargearma'}];
 md.stochasticforcing.defaultdimension    = 2;
-md.stochasticforcing.default_id          = idb_df;
+md.stochasticforcing.default_id          = idb_tf;
 md.stochasticforcing.covariance          = covglob; %global covariance among- and between-fields
 md.stochasticforcing.randomflag          = 0; %determines true/false randomness
@@ -137 +119 @@
 md.transient.isgroundingline = 1;
 
-md.transient.requested_outputs = {'default','CalvingCalvingrate','CalvingMeltingrate','BasalforcingsFloatingiceMeltingRate'};
+md.transient.requested_outputs = {'default','CalvingMeltingrate'};
 md.cluster=generic('name',oshostname(),'np',2);
 md=solve(md,'Transient');
@@ -143 +125 @@
 %Fields and tolerances to track changes
 field_names ={...
-   'Vx1' ,'Vy1' ,'Vel1' ,'Thickness1' ,'MaskIceLevelset1' ,'CalvingCalvingrate1' ,'CalvingMeltingrate1' ,'BasalforcingsFloatingiceMeltingRate1',...
-   'Vx2' ,'Vy2' ,'Vel2' ,'Thickness2' ,'MaskIceLevelset2' ,'CalvingCalvingrate2' ,'CalvingMeltingrate2' ,'BasalforcingsFloatingiceMeltingRate2',...
-   'Vx10','Vy10','Vel10','Thickness10','MaskIceLevelset10','CalvingCalvingrate10','CalvingMeltingrate10','BasalforcingsFloatingiceMeltingRate10',...
+   'Vx1' ,'Vy1' ,'Vel1' ,'Thickness1' ,'MaskIceLevelset1' ,'CalvingMeltingrate1' ,...
+   'Vx2' ,'Vy2' ,'Vel2' ,'Thickness2' ,'MaskIceLevelset2' ,'CalvingMeltingrate2' ,...
+   'Vx10','Vy10','Vel10','Thickness10','MaskIceLevelset10','CalvingMeltingrate10',...
    };
 field_tolerances={...
-   1e-11,2e-11,2e-11,1e-11,1e-9,1e-10,1e-10,1e-10,...
-   2e-11,1e-11,1e-11,9e-11,2e-9,1e-10,1e-10,1e-10,...
-   2e-6,1e-6,1e-6,1e-6,5e-6,1e-6,1e-6,1e-6,...
+   1e-11,2e-11,2e-11,1e-11,1e-9,1e-10,...
+   2e-11,1e-11,1e-11,9e-11,2e-9,1e-10,...
+   2e-6,1e-6,1e-6,1e-6,5e-6,1e-6,...
    };
 field_values={...
@@ -158 +140 @@
    (md.results.TransientSolution(1).Thickness),...
    (md.results.TransientSolution(1).MaskIceLevelset),...
-   (md.results.TransientSolution(1).CalvingCalvingrate),...
    (md.results.TransientSolution(1).CalvingMeltingrate),...
-   (md.results.TransientSolution(1).BasalforcingsFloatingiceMeltingRate),...
    (md.results.TransientSolution(20).Vx),...
    (md.results.TransientSolution(20).Vy),...
@@ -166 +146 @@
    (md.results.TransientSolution(20).Thickness),...
    (md.results.TransientSolution(20).MaskIceLevelset),...
-   (md.results.TransientSolution(20).CalvingCalvingrate),...
    (md.results.TransientSolution(20).CalvingMeltingrate),...
-   (md.results.TransientSolution(20).BasalforcingsFloatingiceMeltingRate),...
         (md.results.TransientSolution(40).Vx),...
         (md.results.TransientSolution(40).Vy),...
@@ -174 +152 @@
         (md.results.TransientSolution(40).Thickness),...
         (md.results.TransientSolution(40).MaskIceLevelset),...
-        (md.results.TransientSolution(40).CalvingCalvingrate),...
         (md.results.TransientSolution(40).CalvingMeltingrate),...
-        (md.results.TransientSolution(40).BasalforcingsFloatingiceMeltingRate),...
         };

issm/trunk-jpl/test/NightlyRun/test543.py

@@ -1 @@
-#Test Name: PigTranStochasticforcing
+#Test Name: PigTranRignotarma 
 import numpy as np
 from frontalforcingsrignotarma import frontalforcingsrignotarma
@@ -29 +29 @@
         for vertex in range(3):
             idb_tf[ii] = 2
-# Basin separation default
-idb_df = np.zeros((md.mesh.numberofelements,))
-iid1 = np.where(md.mesh.x <= -1.62e6)[0]
-for ii in range(md.mesh.numberofelements):
-    for vertex in range(3):
-        if md.mesh.elements[ii][vertex] - 1 in iid1:  # one vertex in basin 1; NOTE: offset because of 1-based vertex indexing
-            idb_df[ii] = 1
-    if idb_df[ii] == 0:  # no vertex was found in basin 1
-        for vertex in range(3):
-            idb_df[ii] = 2
 #Dimensionalities
 nb_tf = 2
-nb_clv  = 2
-nb_flmlt = 2
 
 # Calving parameters
 md.mask.ice_levelset = 1e4 * (md.mask.ice_levelset + 0.5)
-md.calving.calvingrate = 0.3 * np.ones((md.mesh.numberofvertices,))
+md.calving.calvingrate = 0 * np.ones((md.mesh.numberofvertices,))
 md.levelset.spclevelset = np.full((md.mesh.numberofvertices,), np.nan)
 md.levelset.migration_max = 10.0
@@ -119 +107 @@
     md.frontalforcings.sd_polynomialparams   = sd_polyparam
 #Floating Ice Melt parameters
-md.basalforcings.floatingice_melting_rate = 0.1 * np.ones((md.mesh.numberofvertices,))
+md.basalforcings.floatingice_melting_rate = 0.0 * np.ones((md.mesh.numberofvertices,))
 
 #Covariance matrix
 covtf = 1e-4 * np.identity(nb_tf)
-covclv = 1e-1 * np.identity(nb_clv)
-covclv[0, 0] = 1 / 10 * covclv[0, 0]
-covflmlt = 0.05 * np.identity(nb_flmlt)
+covsd = 1e3 * np.identity(nb_tf)
 #covglob          = np.zeros([6,6])
 #covglob[0:2,0:2] = covtf
@@ -132 +118 @@
 
 #Hard-coding covariance matrix because python is complaining
-covglob = np.array([[1e-4, 0., 0., 0., 0., 0.,0.,0.],
-                    [0., 1e-4, 0., 0., 0., 0.,0.,0.],
-                    [0., 0., 1e-2, 0., 0., 0.,0.,0.],
-                    [0., 0., 0., 1e-1, 0., 0.,0.,0.],
-                    [0., 0., 0., 0., 0.05, 0.,0.,0.],
-                    [0., 0., 0., 0., 0., 0.05,0.,0.],
-                    [0., 0., 0., 0., 0., 0., 1e3,0.],
-                    [0., 0., 0., 0., 0., 0.,0., 1e3]])
+covglob = np.array([[1e-4, 0., 0., 0.],
+                    [0., 1e-4, 0., 0.],
+                    [0., 0., 1e3, 0.],
+                    [0., 0., 0., 1e3]])
 #testchol = np.linalg.cholesky(covglob)
 #print(testchol)
@@ -145 +127 @@
 # Stochastic forcing
 md.stochasticforcing.isstochasticforcing = 1
-md.stochasticforcing.fields = ['FrontalForcingsRignotarma', 'DefaultCalving', 'FloatingMeltRate','FrontalForcingsSubglacialDischargearma']
+md.stochasticforcing.fields = ['FrontalForcingsRignotarma','FrontalForcingsSubglacialDischargearma']
 md.stochasticforcing.defaultdimension = 2
-md.stochasticforcing.default_id = idb_df
+md.stochasticforcing.default_id = idb_tf
 md.stochasticforcing.covariance = covglob  # global covariance among- and between-fields
 md.stochasticforcing.randomflag = 0  # determines true/false randomness
@@ -160 +142 @@
 # Fields and tolerances to track changes
 field_names = [
-    'Vx1', 'Vy1', 'Vel1', 'Thickness1', 'MaskIceLevelset1', 'CalvingCalvingrate1', 'CalvingMeltingrate1', 'BasalforcingsFloatingiceMeltingRate1',
-    'Vx2', 'Vy2', 'Vel2', 'Thickness2', 'MaskIceLevelset2', 'CalvingCalvingrate2', 'CalvingMeltingrate2', 'BasalforcingsFloatingiceMeltingRate2',
-    'Vx10', 'Vy10', 'Vel10', 'Thickness10', 'MaskIceLevelset10', 'CalvingCalvingrate10', 'CalvingMeltingrate10', 'BasalforcingsFloatingiceMeltingRate10']
+    'Vx1', 'Vy1', 'Vel1', 'Thickness1', 'MaskIceLevelset1', 'CalvingMeltingrate1',
+    'Vx2', 'Vy2', 'Vel2', 'Thickness2', 'MaskIceLevelset2', 'CalvingMeltingrate2',
+    'Vx10', 'Vy10', 'Vel10', 'Thickness10', 'MaskIceLevelset10', 'CalvingMeltingrate10']
 
 field_tolerances = [
-    1e-11, 2e-11, 2e-11, 1e-11, 1e-9, 1e-10, 1e-10, 1e-10,
-    2e-11, 1e-11, 1e-11, 9e-11, 2e-9, 1e-10, 1e-10, 1e-10,
-    2e-6, 1e-6, 1e-6, 1e-6, 5e-6, 1e-6, 1e-6, 1e-6]
+        1e-11,2e-11,2e-11,1e-11,1e-9,1e-10,
+        2e-11,1e-11,1e-11,9e-11,2e-9,1e-10,
+        2e-6,1e-6,1e-6,1e-6,5e-6,1e-6]
 field_values = [
     md.results.TransientSolution[0].Vx,
@@ -174 +156 @@
     md.results.TransientSolution[0].Thickness,
     md.results.TransientSolution[0].MaskIceLevelset,
-    md.results.TransientSolution[0].CalvingCalvingrate,
     md.results.TransientSolution[0].CalvingMeltingrate,
-    md.results.TransientSolution[0].BasalforcingsFloatingiceMeltingRate,
     md.results.TransientSolution[19].Vx,
     md.results.TransientSolution[19].Vy,
@@ -182 +162 @@
     md.results.TransientSolution[19].Thickness,
     md.results.TransientSolution[19].MaskIceLevelset,
-    md.results.TransientSolution[19].CalvingCalvingrate,
     md.results.TransientSolution[19].CalvingMeltingrate,
-    md.results.TransientSolution[19].BasalforcingsFloatingiceMeltingRate,
     md.results.TransientSolution[39].Vx,
     md.results.TransientSolution[39].Vy,
@@ -190 +168 @@
     md.results.TransientSolution[39].Thickness,
     md.results.TransientSolution[39].MaskIceLevelset,
-    md.results.TransientSolution[39].CalvingCalvingrate,
-    md.results.TransientSolution[39].CalvingMeltingrate,
-    md.results.TransientSolution[39].BasalforcingsFloatingiceMeltingRate]
+    md.results.TransientSolution[39].CalvingMeltingrate]