[cig-commits] [commit] master: adding min and max radius/colat to netcdf file (5c46c43)

[cig_noreply at geodynamics.org]

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

On branch  : master
Link       : https://github.com/geodynamics/axisem/compare/f5897a60c20dba957e3501960ebe1c1a1234a7f3...5c46c43fff66f9350665d71c46395c8cc2550df8

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

commit 5c46c43fff66f9350665d71c46395c8cc2550df8
Author: martinvandriel <martin at vandriel.de>
Date:   Thu Sep 4 14:58:41 2014 +0200

    adding min and max radius/colat to netcdf file


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

5c46c43fff66f9350665d71c46395c8cc2550df8
 SOLVER/parameters.F90 | 6 ++++++
 1 file changed, 6 insertions(+)

diff --git a/SOLVER/parameters.F90 b/SOLVER/parameters.F90
index c1c9c9d..2306904 100644
--- a/SOLVER/parameters.F90
+++ b/SOLVER/parameters.F90
@@ -1450,6 +1450,12 @@ subroutine write_parameters
            call nc_write_att_int(nstrain,              'number of strain dumps')
            call nc_write_att_dble(deltat_coarse,       'strain dump sampling rate in sec')
            call nc_write_att_char(trim(dump_type),     'dump type (displ_only, displ_velo, fullfields)')
+           if (dump_type == 'displ_only') then
+              call nc_write_att_dble(kwf_rmin / 1000,  'kernel wavefield rmin')
+              call nc_write_att_dble(min(kwf_rmax, router) / 1000,  'kernel wavefield rmax')
+              call nc_write_att_dble(kwf_thetamin / pi * 180.,  'kernel wavefield colatmin')
+              call nc_write_att_dble(kwf_thetamax / pi * 180.,  'kernel wavefield colatmax')
+           endif
         else
            call nc_write_att_int(0,                    'number of strain dumps')       
            call nc_write_att_dble(0.d0,                'strain dump sampling rate in sec' )