Changeset 921


Ignore:
Timestamp:
03/21/13 17:40:34 (6 years ago)
Author:
marcus
Message:
  • local magnitude formula parameters are now configurable
  • preparing next version

(internal: r423 | walther | 2013-03-21 17:25:05 +0100 (Thu, 21 Mar 2013) | 3 lines)

Location:
SH_SHM/trunk
Files:
8 edited

Legend:

Unmodified
Added
Removed
  • SH_SHM/trunk/inputs/shm-config.txt

    r898 r921  
    6565!                            More information: http://www.seismic-handler.org/wiki/HowTo/LocalMagnitudeFormula 
    6666local_magnitude_formula      classic 
     67 
     68!                            The following values define the default parameters for 
     69!                            the new formula. You are asked to alter theses values to 
     70!                            match your local conditions. 
     71local_magnitude_param_n      1.11 
     72local_magnitude_param_k      0.00189 
     73local_magnitude_param_offset 2.09 
    6774 
    6875!                            Stations in read dialog box for MiniSEED data. 
  • SH_SHM/trunk/source/globalparams.c

    r898 r921  
    33 *      ============== 
    44 * 
    5  * $Revision: 406 $, $Date: 2013-02-07 16:03:35 +0100 (Do, 07. Feb 2013) $ 
    6  * 
    75 * Management of global parameters 
    8  * K. Stammler, 9-May-2006 
     6 * K. Stammler 
    97 */ 
    108 
     
    212210        "axis_label_rounding", 
    213211        "autofilter_sepdist", 
     212        "local_magnitude_param_n", 
     213        "local_magnitude_param_k", 
     214        "local_magnitude_param_offset", 
    214215        "f_last_parameter" 
    215216}; 
     
    277278        0.01,      /* axis_label_rounding */ 
    278279        5.0,       /* autofilter_sepdist */ 
     280        1.11,      /* local magnitude 1st parameter */ 
     281        0.00189,   /* local magnitude 2nd parameter */ 
     282        2.09,      /* local magnitude 3rd parameter */ 
    279283        0.0        /* last parameter */ 
    280284}; 
  • SH_SHM/trunk/source/globalparams.h

    r898 r921  
    33 *      ============== 
    44 * 
    5  * $Revision: 406 $, $Date: 2013-02-07 16:03:35 +0100 (Do, 07. Feb 2013) $ 
    6  * 
    75 * Management of global parameters 
    8  * K. Stammler, 9-May-2006 
     6 * K. Stammler 
    97 */ 
    108 
     
    160158        cGpF_axis_label_rounding, 
    161159        cGpF_autofilter_sepdist, 
     160        cGpF_local_magnitude_param_n, 
     161        cGpF_local_magnitude_param_k, 
     162        cGpF_local_magnitude_param_offset, 
    162163        cGpF_last 
    163164} TGpFloatIndex; 
  • SH_SHM/trunk/source/motif/cblib2.c

    r900 r921  
    33 *      ======= 
    44 * 
    5  * version 42,9-Nov-2006 
    6  * 
    75 * library of utility routines for callbacks 
    8  * K. Stammler, 28-Oct-93 
     6 * K. Stammler 
    97 */ 
    108 
     
    399397        formula_mode = GpGetString(cGpS_local_magnitude_formula); 
    400398        if (strcmp(formula_mode, "iaspei") == 0) 
    401             ml = si_magn_ml_analytic( ml_ampl, distance_km, status ); 
     399            ml = si_magn_ml_analytic( 
     400                    ml_ampl, distance_km,  
     401                    GpGetFloat(cGpF_local_magnitude_param_n), 
     402                    GpGetFloat(cGpF_local_magnitude_param_k), 
     403                    GpGetFloat(cGpF_local_magnitude_param_offset), status 
     404            ); 
    402405        else 
    403406            ml = si_magn_ml( ml_ampl, distance_km, status ); 
  • SH_SHM/trunk/source/motif/seismics.c

    r902 r921  
    33 *      ========== 
    44 * 
    5  * version 34, 16-Jun-2006 
    6  * 
    75 * seismic utilities 
    8  * K. Stammler, 25-Mar-93 
     6 * K. Stammler 
    97 */ 
    108 
     
    499497 
    500498 
    501 float si_magn_ml_analytic( float ampl, float distance, STATUS *status ) 
     499float si_magn_ml_analytic( float ampl, float distance, float n, float k, float offset, STATUS *status ) 
    502500 
    503501/* Returns magnitude ml using analytical formula.  Amplitude and period must 
     
    507505 * float      ampl;        input; amplitude of WOOD-ANDERSON simulation 
    508506 * float      distance;    input; epicentral distance in km 
     507 * float      n;           input; 1st parameter for distance correction 
     508 * float      k;           input; 2nd parameter for distance correction 
     509 * float      offset;      input; offset parameter 
    509510 * STATUS     *status;     output; return status 
    510511 *                         returns value of ml 
     
    530531        ampl *= 1000; /* nm */ 
    531532        magn = log10((double)ampl); 
    532         magn += 1.11 * log10((double)distance); 
    533         magn += 0.00189 * (double)distance; 
    534         magn -= 2.09; 
     533        magn += n * log10((double)distance); 
     534        magn += k * (double)distance; 
     535        magn -= offset; 
    535536 
    536537        if  (GpGetInt(cGpI_debug_level) > 0) 
  • SH_SHM/trunk/source/motif/seismics.h

    r898 r921  
    33 *      ========== 
    44 * 
    5  * version 11, 22-May-2005 
    6  * 
    75 * header file of seismics.c (seismic utilities) 
    8  * K. Stammler, 25-Mar-93 
     6 * K. Stammler 
    97 */ 
    108 
     
    148146 
    149147 
    150 float si_magn_ml_analytic( float ampl, float distance, STATUS *status ); 
     148float si_magn_ml_analytic( float ampl, float distance, float n, float k, float offset, STATUS *status ); 
    151149 
    152150/* Returns magnitude ml using analytical formula.  Amplitude and period must 
     
    156154 * float      ampl;        input; amplitude of WOOD-ANDERSON simulation 
    157155 * float      distance;    input; epicentral distance in km 
     156 * float      n;           input; 1st parameter for distance correction 
     157 * float      k;           input; 2nd parameter for distance correction 
     158 * float      offset;      input; offset parameter 
    158159 * STATUS     *status;     output; return status 
    159160 *                         returns value of ml 
  • SH_SHM/trunk/source/motif/shm.c

    r898 r921  
    33 *      ===== 
    44 * 
    5  * $Revision: 406 $, $Date: 2013-02-07 16:03:35 +0100 (Do, 07. Feb 2013) $ 
    6  * 
    75 * main module of shm 
    8  * K. Stammler, 15-Feb-93 
     6 * K. Stammler 
    97 */ 
    108 
     
    683681 
    684682        printf( "\n" ); 
    685         printf( "SeismicHandler Motif version 2013a, Copyright (C) 1990-2013\n" ); 
     683        printf( "SeismicHandler Motif version 2013b, Copyright (C) 1990-2013\n" ); 
    686684        printf( "Klaus Stammler, Seismological Observatory Graefenberg (SZGRF)\n" ); 
    687685        printf( "of the Federal Institute for Geosciences and Natural Resources (BGR).\n" ); 
  • SH_SHM/trunk/source/shconst.h

    r898 r921  
    33 *      ========= 
    44 * 
    5  * $Revision: 406 $, $Date: 2013-02-07 16:03:35 +0100 (Do, 07. Feb 2013) $ 
    6  * 
    7  * constants of seismhandler program 
    8  * K. Stammler, 22-MAY-91 
     5 * constants of Seismic Handler program 
    96 */ 
    107 
     
    3633#define __SHCONST 
    3734 
    38 #define SHC_VERSION "2013a" 
     35#define SHC_VERSION "2013b" 
    3936        /* version number */ 
    4037#define TW 1 
Note: See TracChangeset for help on using the changeset viewer.