Changeset 24870 for issm/trunk-jpl/test/NightlyRun/test244.py

Timestamp:

05/17/20 17:15:00 (5 years ago)

Author:

jdquinn

Message:

CHG: Translated MATLAB changes to Python that allow for distributing variable over partition.

File:

• issm/trunk-jpl/test/NightlyRun/test244.py (modified) (2 diffs)

issm/trunk-jpl/test/NightlyRun/test244.py ¶

@@ -73 +73 @@
 
 #variables
-md.qmu.variables.surface_mass_balanceC = normal_uncertain.normal_uncertain('scaled_SmbC', 1, 0.5)
+md.qmu.variables.surface_mass_balance = normal_uncertain.normal_uncertain(
+    'descriptor', 'scaled_SmbC',
+    'mean', np.ones(md.qmu.numberofpartitions),
+    'stddev', .5 * np.ones(md.qmu.numberofpartitions),
+    'partition', md.qmu.epartition
+    )
 Tmin = 273.
 telms = np.atleast_2d(np.min(md.smb.Ta[0:-1, :], 1))
@@ -85 +90 @@
 
 mint_on_partition[np.where(np.isnan(mint_on_partition))] = 10**-10
-md.qmu.variables.surface_mass_balanceTa = uniform_uncertain.uniform_uncertain('scaled_SmbTa', 1, 0.05)
-md.qmu.variables.surface_mass_balanceTa[0].lower = 0.95
-md.qmu.variables.surface_mass_balanceTa[0].upper = np.maximum(np.minimum(np.maximum(1.05, mint_on_partition), 0.9999), 0.0001)
+md.qmu.variables.surface_mass_balanceTa = uniform_uncertain.uniform_uncertain(
+    'descriptor', 'scaled_SmbTa',
+    'lower', .95 * np.ones(md.qmu.numberofpartitions),
+    'upper', np.maximum(np.minimum(np.maximum(1.05, mint_on_partition), 0.9999), 0.0001),
+    'partition', md.qmu.epartition
+    )
 
 #responses
-md.qmu.responses.IceVolume = response_function.response_function('IceVolume', [], [0.0001, 0.001, 0.01, 0.25, 0.5, 0.75, 0.99, 0.999, 0.9999])
-md.qmu.responses.IceMass = response_function.response_function('IceMass', [], [0.0001, 0.001, 0.01, 0.25, 0.5, 0.75, 0.99, 0.999, 0.9999])
-md.qmu.responses.TotalSmb = response_function.response_function('TotalSmb', [], [0.0001, 0.001, 0.01, 0.25, 0.5, 0.75, 0.99, 0.999, 0.9999])
+md.qmu.responses.IceVolume = response_function.response_function('descriptor','IceVolume')
+md.qmu.responses.IceMass = response_function.response_function('descriptor','IceMass')
+md.qmu.responses.TotalSmb = response_function.response_function('descriptor','TotalSmb')
 
 #  nond_sampling study