Changeset 22398

Timestamp:

02/06/18 14:03:37 (7 years ago)

Author:

kruegern

Message:

BUG: fixed test2085.py and paresresultsfromdisk.py

Location:

issm/trunk-jpl

Files:

• src/m/solve/parseresultsfromdisk.py (modified) (1 diff)
• test/NightlyRun/test2085.py (added)

issm/trunk-jpl/src/m/solve/parseresultsfromdisk.py

@@ -221 +221 @@
                 elif fieldname=='CalvingCalvingrate':
                         field = field*yts
+                elif fieldname=='LoveKernelsReal' or fieldname=='LoveKernelsImag':
+                        nlayer = md.materials.numlayers
+                        degmax = md.love.sh_nmax
+                        nfreq  = md.love.nfreq
+                        temp_field = np.full((degmax+1,nfreq,nlayer+1,6),0.0)
+                        for ii in range(degmax+1):
+                                for jj in range(nfreq):
+                                        for kk in range(nlayer+1):
+                                                ll = ii*(nlayer+1)*6 + (kk*6+1)
+                                                for mm in range(6):
+                                                        temp_field[ii,jj,kk,mm] = field[ll+mm-1,jj]
+                        field=temp_field
 
                 if time !=-9999: