[cig-commits] [commit] master: adding strain_only mode to fieldtransform (d23dbce)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/axisem

On branch  : master
Link       : https://github.com/geodynamics/axisem/compare/ca4b3002407f17f95bbf7013f7622a43ceb4db35...7a0ed4115e0445ccb0219c894eedf4bf157d8689

>---------------------------------------------------------------

commit d23dbce69af409a8c69dec8b681687a8e17f93f7
Author: martinvandriel <martin at vandriel.de>
Date:   Mon Oct 13 16:49:12 2014 +0200

    adding strain_only mode to fieldtransform


>---------------------------------------------------------------

d23dbce69af409a8c69dec8b681687a8e17f93f7
 SOLVER/UTILS/field_transform.F90 | 13 +++++++++++++
 1 file changed, 13 insertions(+)

diff --git a/SOLVER/UTILS/field_transform.F90 b/SOLVER/UTILS/field_transform.F90
index 778b40f..a743ff7 100644
--- a/SOLVER/UTILS/field_transform.F90
+++ b/SOLVER/UTILS/field_transform.F90
@@ -163,6 +163,19 @@ program field_transformation
            varnamelist = ['disp_s     ', 'disp_p     ', 'disp_z     ']
        end if
 
+    elseif (trim(dump_type) == 'strain_only') then
+       if (sourcetype=='monopole')  then
+           nvar = 4
+           allocate(varnamelist(nvar))
+           varnamelist = (/'strain_dsus', 'strain_dsuz', 'strain_dpup', &
+                           'straintrace'/)
+       else
+           nvar = 6
+           allocate(varnamelist(nvar))
+           varnamelist = (/'strain_dsus', 'strain_dsuz', 'strain_dpup', &
+                           'strain_dsup', 'strain_dzup', 'straintrace'/)
+       end if
+
     elseif (trim(dump_type) == 'fullfields') then
        if (sourcetype=='monopole')  then
            nvar = 6