Ignore:
Timestamp:
10/22/20 12:31:59 (4 years ago)
Author:
jdquinn
Message:

BUG: Addressing of TransientSolution results structure; clean up

File:
1 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/test/NightlyRun/test252.py

    r25398 r25709  
    7070md = solve(md, 'Transient')
    7171
    72 nlayers=np.size(md.results.TransientSolution[0].SmbT,1)
    73 for i in range(1, np.size(md.results.TransientSolution,0)):
    74     nlayers=np.minimum(np.size(md.results.TransientSolution[i].SmbT,1),nlayers)
     72nlayers = md.results.TransientSolution[0].SmbT.shape[1]
     73for i in range(1, len(md.results.TransientSolution)):
     74    nlayers=np.minimum(md.results.TransientSolution[i].SmbT.shape[1], nlayers)
    7575
    7676#Fields and tolerances to track changes
    77 
    7877field_names = ['Layers','SmbDz1', 'SmbT1', 'SmbD1', 'SmbRe1', 'SmbGdn1', 'SmbGsp1', 'SmbA1', 'SmbEC1', 'SmbMassBalance1', 'SmbMAdd1', 'SmbDzAdd1', 'SmbFAC1',
    7978               'SmbDz2', 'SmbT2', 'SmbD2', 'SmbRe2', 'SmbGdn2', 'SmbGsp2', 'SmbA2', 'SmbEC2', 'SmbMassBalance2', 'SmbMAdd2', 'SmbDzAdd2', 'SmbFAC2',
    8079               'SmbDz3', 'SmbT3', 'SmbD3', 'SmbRe3', 'SmbGdn3', 'SmbGsp3', 'SmbA3', 'SmbEC3', 'SmbMassBalance3', 'SmbMAdd3', 'SmbDzAdd3', 'SmbFAC3',
    8180               'SmbDz4', 'SmbT4', 'SmbD4', 'SmbRe4', 'SmbGdn4', 'SmbGsp4', 'SmbA4', 'SmbEC4', 'SmbMassBalance4', 'SmbMAdd4', 'SmbDzAdd4', 'SmbFAC4']
    82 field_tolerances = [1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,
    83                    1e-12,1e-12,1e-11,1e-10,2e-11,1e-11,1e-12,1e-11,1e-12,1e-12,1e-12,1e-11,
    84                    1e-12,1e-12,2e-12,2e-11,4e-11,1e-11,1e-12,1e-11,1e-11,1e-12,1e-12,1e-11,
    85                    1e-11,1e-11,4e-11,4e-11,1e-12,4e-11,1e-12,1e-12,1e-10,1e-12,1e-12,2e-11]
    86 
    87 
    88 #shape is different in python solution (fixed using reshape) which can cause test failure:
    89 field_values = [nlayers,md.results.TransientSolution[0].SmbDz[0, 0:nlayers].reshape(1, -1),
    90                 md.results.TransientSolution[0].SmbT[0, 0:nlayers].reshape(1, -1),
    91                 md.results.TransientSolution[0].SmbD[0, 0:nlayers].reshape(1, -1),
    92                 md.results.TransientSolution[0].SmbRe[0, 0:nlayers].reshape(1, -1),
    93                 md.results.TransientSolution[0].SmbGdn[0, 0:nlayers].reshape(1, -1),
    94                 md.results.TransientSolution[0].SmbGsp[0, 0:nlayers].reshape(1, -1),
    95                 md.results.TransientSolution[0].SmbA[0, 0:nlayers].reshape(1, -1),
    96                 md.results.TransientSolution[0].SmbEC[0],
    97                 md.results.TransientSolution[0].SmbMassBalance[0],
    98                 md.results.TransientSolution[0].SmbMAdd[0],
    99                 md.results.TransientSolution[0].SmbDzAdd[0],
    100                 md.results.TransientSolution[0].SmbFAC[0],
    101                 md.results.TransientSolution[145].SmbDz[0, 0:nlayers].reshape(1, -1),
    102                 md.results.TransientSolution[145].SmbT[0, 0:nlayers].reshape(1, -1),
    103                 md.results.TransientSolution[145].SmbD[0, 0:nlayers].reshape(1, -1),
    104                 md.results.TransientSolution[145].SmbRe[0, 0:nlayers].reshape(1, -1),
    105                 md.results.TransientSolution[145].SmbGdn[0, 0:nlayers].reshape(1, -1),
    106                 md.results.TransientSolution[145].SmbGsp[0, 0:nlayers].reshape(1, -1),
    107                 md.results.TransientSolution[145].SmbA[0, 0:nlayers].reshape(1, -1),
    108                 md.results.TransientSolution[145].SmbEC[0],
    109                 md.results.TransientSolution[145].SmbMassBalance[0],
    110                 md.results.TransientSolution[145].SmbMAdd[0],
    111                 md.results.TransientSolution[145].SmbDzAdd[0],
    112                 md.results.TransientSolution[145].SmbFAC[0],
    113                 md.results.TransientSolution[146].SmbDz[0, 0:nlayers].reshape(1, -1),
    114                 md.results.TransientSolution[146].SmbT[0, 0:nlayers].reshape(1, -1),
    115                 md.results.TransientSolution[146].SmbD[0, 0:nlayers].reshape(1, -1),
    116                 md.results.TransientSolution[146].SmbRe[0, 0:nlayers].reshape(1, -1),
    117                 md.results.TransientSolution[146].SmbGdn[0, 0:nlayers].reshape(1, -1),
    118                 md.results.TransientSolution[146].SmbGsp[0, 0:nlayers].reshape(1, -1),
    119                 md.results.TransientSolution[146].SmbA[0, 0:nlayers].reshape(1, -1),
    120                 md.results.TransientSolution[146].SmbEC[0],
    121                 md.results.TransientSolution[146].SmbMassBalance[0],
    122                 md.results.TransientSolution[146].SmbMAdd[0],
    123                 md.results.TransientSolution[146].SmbDzAdd[0],
    124                 md.results.TransientSolution[146].SmbFAC[0],
    125                 md.results.TransientSolution[-1].SmbDz[0, 0:nlayers].reshape(1, -1),
    126                 md.results.TransientSolution[-1].SmbT[0, 0:nlayers].reshape(1, -1),
    127                 md.results.TransientSolution[-1].SmbD[0, 0:nlayers].reshape(1, -1),
    128                 md.results.TransientSolution[-1].SmbRe[0, 0:nlayers].reshape(1, -1),
    129                 md.results.TransientSolution[-1].SmbGdn[0, 0:nlayers].reshape(1, -1),
    130                 md.results.TransientSolution[-1].SmbGsp[0, 0:nlayers].reshape(1, -1),
    131                 md.results.TransientSolution[-1].SmbA[0, 0:nlayers].reshape(1, -1),
    132                 md.results.TransientSolution[-1].SmbEC[0],
    133                 md.results.TransientSolution[-1].SmbMassBalance[0],
    134                 md.results.TransientSolution[-1].SmbMAdd[0],
    135                 md.results.TransientSolution[-1].SmbDzAdd[0],
    136                 md.results.TransientSolution[-1].SmbFAC[0]]
     81field_tolerances = [
     82    1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,
     83    1e-12,1e-12,1e-11,1e-10,2e-11,1e-11,1e-12,1e-11,1e-12,1e-12,1e-12,1e-11,
     84    1e-12,1e-12,2e-12,2e-11,4e-11,1e-11,1e-12,1e-11,1e-11,1e-12,1e-12,1e-11,
     85    1e-11,1e-11,4e-11,4e-11,1e-12,4e-11,1e-12,1e-12,1e-10,1e-12,1e-12,2e-11
     86]
     87# Shape is different in python solution (fixed using reshape) which can cause test failure
     88field_values = [
     89    nlayers,
     90    md.results.TransientSolution[0].SmbDz[0, 0:nlayers].reshape(1, -1),
     91    md.results.TransientSolution[0].SmbT[0, 0:nlayers].reshape(1, -1),
     92    md.results.TransientSolution[0].SmbD[0, 0:nlayers].reshape(1, -1),
     93    md.results.TransientSolution[0].SmbRe[0, 0:nlayers].reshape(1, -1),
     94    md.results.TransientSolution[0].SmbGdn[0, 0:nlayers].reshape(1, -1),
     95    md.results.TransientSolution[0].SmbGsp[0, 0:nlayers].reshape(1, -1),
     96    md.results.TransientSolution[0].SmbA[0, 0:nlayers].reshape(1, -1),
     97    md.results.TransientSolution[0].SmbEC[0],
     98    md.results.TransientSolution[0].SmbMassBalance[0],
     99    md.results.TransientSolution[0].SmbMAdd[0],
     100    md.results.TransientSolution[0].SmbDzAdd[0],
     101    md.results.TransientSolution[0].SmbFAC[0],
     102    md.results.TransientSolution[145].SmbDz[0, 0:nlayers].reshape(1, -1),
     103    md.results.TransientSolution[145].SmbT[0, 0:nlayers].reshape(1, -1),
     104    md.results.TransientSolution[145].SmbD[0, 0:nlayers].reshape(1, -1),
     105    md.results.TransientSolution[145].SmbRe[0, 0:nlayers].reshape(1, -1),
     106    md.results.TransientSolution[145].SmbGdn[0, 0:nlayers].reshape(1, -1),
     107    md.results.TransientSolution[145].SmbGsp[0, 0:nlayers].reshape(1, -1),
     108    md.results.TransientSolution[145].SmbA[0, 0:nlayers].reshape(1, -1),
     109    md.results.TransientSolution[145].SmbEC[0],
     110    md.results.TransientSolution[145].SmbMassBalance[0],
     111    md.results.TransientSolution[145].SmbMAdd[0],
     112    md.results.TransientSolution[145].SmbDzAdd[0],
     113    md.results.TransientSolution[145].SmbFAC[0],
     114    md.results.TransientSolution[146].SmbDz[0, 0:nlayers].reshape(1, -1),
     115    md.results.TransientSolution[146].SmbT[0, 0:nlayers].reshape(1, -1),
     116    md.results.TransientSolution[146].SmbD[0, 0:nlayers].reshape(1, -1),
     117    md.results.TransientSolution[146].SmbRe[0, 0:nlayers].reshape(1, -1),
     118    md.results.TransientSolution[146].SmbGdn[0, 0:nlayers].reshape(1, -1),
     119    md.results.TransientSolution[146].SmbGsp[0, 0:nlayers].reshape(1, -1),
     120    md.results.TransientSolution[146].SmbA[0, 0:nlayers].reshape(1, -1),
     121    md.results.TransientSolution[146].SmbEC[0],
     122    md.results.TransientSolution[146].SmbMassBalance[0],
     123    md.results.TransientSolution[146].SmbMAdd[0],
     124    md.results.TransientSolution[146].SmbDzAdd[0],
     125    md.results.TransientSolution[146].SmbFAC[0],
     126    md.results.TransientSolution[-1].SmbDz[0, 0:nlayers].reshape(1, -1),
     127    md.results.TransientSolution[-1].SmbT[0, 0:nlayers].reshape(1, -1),
     128    md.results.TransientSolution[-1].SmbD[0, 0:nlayers].reshape(1, -1),
     129    md.results.TransientSolution[-1].SmbRe[0, 0:nlayers].reshape(1, -1),
     130    md.results.TransientSolution[-1].SmbGdn[0, 0:nlayers].reshape(1, -1),
     131    md.results.TransientSolution[-1].SmbGsp[0, 0:nlayers].reshape(1, -1),
     132    md.results.TransientSolution[-1].SmbA[0, 0:nlayers].reshape(1, -1),
     133    md.results.TransientSolution[-1].SmbEC[0],
     134    md.results.TransientSolution[-1].SmbMassBalance[0],
     135    md.results.TransientSolution[-1].SmbMAdd[0],
     136    md.results.TransientSolution[-1].SmbDzAdd[0],
     137    md.results.TransientSolution[-1].SmbFAC[0]
     138]
Note: See TracChangeset for help on using the changeset viewer.