source: issm/oecreview/Archive/24684-25833/ISSM-25708-25709.diff@ 25834

Last change on this file since 25834 was 25834, checked in by Mathieu Morlighem, 4 years ago

CHG: added 24684-25833
File size: 21.0 KB

  • ../trunk-jpl/test/NightlyRun/test243.py

     
    6363#Run transient
    6464md = solve(md, 'Transient')
    6565
    66 nlayers = np.size(md.results.TransientSolution[0].SmbT, 1)
     66nlayers = md.results.TransientSolution[0].SmbT.shape[1]
    6767for i in range(1, len(md.results.TransientSolution)):
    68     nlayers = np.minimum(np.size(md.results.TransientSolution[i].SmbT, 1), nlayers)
     68    nlayers=np.minimum(md.results.TransientSolution[i].SmbT.shape[1], nlayers)
    6969
    7070#Fields and tolerances to track changes
    7171field_names = ['Layers', 'SmbDz', 'SmbT', 'SmbD', 'SmbRe', 'SmbGdn', 'SmbGsp', 'SmbA', 'SmbEC', 'SmbMassBalance', 'SmbMAdd', 'SmbDzAdd', 'SmbFAC', 'SmbMeanSHF', 'SmbMeanLHF', 'SmbMeanULW', 'SmbNetLW', 'SmbNetSW']
    7272field_tolerances = [1e-12, 2e-12, 1e-12, 1e-11, 1e-11, 2e-11, 1e-11, 1e-12, 1e-11, 1e-12, 1e-12, 1e-12, 1e-11, 2e-11, 2e-11, 1e-11, 9e-10, 2e-11]
    73 #shape is different in python solution (fixed using reshape) which can cause test failure:
     73# Shape is different in python solution (fixed using reshape) which can cause test failure
    7474field_values = [
    7575    nlayers,
    76     md.results.TransientSolution[-1].SmbDz[0, 0:nlayers],
    77     md.results.TransientSolution[-1].SmbT[0, 0:nlayers],
    78     md.results.TransientSolution[-1].SmbD[0, 0:nlayers],
    79     md.results.TransientSolution[-1].SmbRe[0, 0:nlayers],
    80     md.results.TransientSolution[-1].SmbGdn[0, 0:nlayers],
    81     md.results.TransientSolution[-1].SmbGsp[0, 0:nlayers],
    82     md.results.TransientSolution[-1].SmbA[0, 0:nlayers],
     76    md.results.TransientSolution[-1].SmbDz[0, 0:nlayers].reshape(1, -1),
     77    md.results.TransientSolution[-1].SmbT[0, 0:nlayers].reshape(1, -1),
     78    md.results.TransientSolution[-1].SmbD[0, 0:nlayers].reshape(1, -1),
     79    md.results.TransientSolution[-1].SmbRe[0, 0:nlayers].reshape(1, -1),
     80    md.results.TransientSolution[-1].SmbGdn[0, 0:nlayers].reshape(1, -1),
     81    md.results.TransientSolution[-1].SmbGsp[0, 0:nlayers].reshape(1, -1),
     82    md.results.TransientSolution[-1].SmbA[0, 0:nlayers].reshape(1, -1),
    8383    md.results.TransientSolution[-1].SmbEC[0],
    8484    md.results.TransientSolution[-1].SmbMassBalance[0],
    8585    md.results.TransientSolution[-1].SmbMAdd[0],

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

     
    6969#Run transient
    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]

  • ../trunk-jpl/test/NightlyRun/test244.py

     
    9494    'lower', .95 * np.ones((npart, 1)),
    9595    'upper', upper,
    9696    'partition', partition
    97     )
     97)
    9898
    9999#responses
    100100md.qmu.responses.IceVolume = response_function.response_function('descriptor','IceVolume')

  • ../trunk-jpl/test/NightlyRun/test253.py

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