Changeset 27035 for issm/trunk/test/NightlyRun/test257.py
- Timestamp:
- 06/01/22 05:01:48 (3 years ago)
- Location:
- issm/trunk
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk
- Property svn:mergeinfo changed
/issm/trunk-jpl merged: 26745-26955,26957-27031
- Property svn:mergeinfo changed
-
issm/trunk/test
- Property svn:mergeinfo changed
-
issm/trunk/test/NightlyRun/test257.py
r26744 r27035 43 43 md.smb = SMBautoregression() 44 44 md.smb.num_basins = 3 # number of basins 45 md.smb.basin_id = idbasin - 1 # prescribe basin ID number to elements; # NOTE: offset because of 1-based vertex indexing45 md.smb.basin_id = idbasin # prescribe basin ID number to elements; 46 46 md.smb.beta0 = np.array([[0.5, 1.2, 1.5]]) # intercept values of SMB in basins [m ice eq./yr] 47 47 md.smb.beta1 = np.array([[0.0, 0.01, -0.01]]) # trend values of SMB in basins [m ice eq./yr^2] … … 50 50 md.smb.ar_timestep = 2.0 #timestep of the autoregressive model [yr] 51 51 md.smb.phi = np.array([[0.2, 0.1, 0.05, 0.01], [0.4, 0.2, -0.2, 0.1], [0.4, -0.4, 0.1, -0.1]]) 52 md.smb.lapserates = np.array([[0.01,0.0],[0.01,-0.01],[0.0,-0.01]]) 53 md.smb.elevationbins = np.array([100,150,100]).reshape(md.smb.num_basins,1) 52 54 53 55 # Stochastic forcing … … 56 58 md.stochasticforcing.covariance = np.array([[0.15, 0.08, -0.02], [0.08, 0.12, -0.05], [-0.02, -0.05, 0.1]]) # global covariance among- and between-fields 57 59 md.stochasticforcing.randomflag = 0 # fixed random seeds 58 60 md.stochasticforcing.stochastictimestep = 1.0 59 61 60 62 md = solve(md, 'Transient')
Note:
See TracChangeset
for help on using the changeset viewer.