kiwi_info

1
2
3
4
5
6
7
8
delta=min(100.max(shortest_doi,0.5)/2.,4000.)   ! delta_max=4km     ; delta_min= 1000.25=25m
nxf=floor(2.0bord_semi_x/delta)                 ! 2.02265/50=90    ;
nyf=floor(2.0bord_semi_y/delta)                 ! 2.03150/50=126
if(nxf.le.1) nxf=2 ; if(bord_semi_x.eq.0) nxf=1  ! line source
if(nyf.le.1) nyf=2 ; if(bord_semi_y.eq.0) nyf=1  ! line source
!print ,‘–> nxf,nyf’ ,nxf,nyf  ! 90,126
delta_x = 2.0bord_semi_x/nxf  !
delta_y = 2.0*bord_semi_y/nyf  !

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
set_inversion_method sa1 -C0/0/0 -N10/100 -P0.25 -R0.5/0.25/0.5/0.5 -S0.1 -Ts0.03 -Tr0.95 -V
  [-C<iselect>/<ioptimize>/<imisfit>/[<ireduction>[/<imodel>]]]  :
      iselect:   -1: one params 0: one fault 1: all params :
      ioptimize: 0: GSA   1: ASA   2: MCMC :
      imisfit:   0: MSE   1: BAYESIAN :
      irecduction: 0:ASA  1: EXP
      imodel: 0:UNIFORM  1: RANDOM 3: OLD
  [-L]  : Layer model instead of halfspace
  [-Mi<0|1|2>]  : initial model: 0:uniform, 1:random, 2:old value,  have not implement.
  [-N<nsteps>/<niterations>/[<init_ntiral>]]   :
  [-O]  : output result each step to sa_out.hd5
  [-P<val>] : poisson ratio
  [-R<moment>[/<rake>[/<fronttime>[/<risetime>]]]]:factor set range (upper and lower boundary)
  [-Q<val>] : quench factor -Q1.0
  [-S] : smooth_factor [0.0 ~ +inf)
  [-Ts<val>] : Ts: set temperature scale factor
  [-Tr<val>] : Tr: set reduction factor, for -C <ireduction>=0
  [-T<i1>/<i2>/<i3>]  : set terminal condition, have not implement
  [-T<n>/<m>]  : tune facor -T<100.0>/<100.0>  n=kf, m=tb/t
  [-V]  : Verbose mode will plot result after each circle
  [-W<mcmc_ratio>]  : Makov Chain Width
  [-Z]  : Layer model instead of halfspace

1
2
3
4
5
6
7
set_inversion_method slipgen  [-K<kvalue>] [-P<factor>] [-R] [-W<c|e>] [-G] [-S<value>]
  [-K<1.0>]  : k^-2 parameter: (0.0~1.0)
  [-P<0.0>]  : rake perturbation factor: (0.0~1.0)
  [-R]       : texture risetime, auto-determined by moment distribution
  [-W<c>]    : save in which grid, c:cgrid, e:egrid
  [-S<0.5>]  : spike model, value spike_location
  [-G]       : random generator layout of derterminstic part 5X3 array

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
          

         current_node=>list_head ; i_sign=1
         do while(associated(current_node))
             psm=>current_node%psm
             current_node=>current_node%next


         ....

         call psm_transfer_params(psm%sourcetype,psm%params,i_sign)  
         ! i.e. copy psm%params =>  psm_params_cache_eikonal_modify

         grid%psm_sourcetype=psm%sourcetype 
         ! save a copy to grid%psm_params
         call psm_transfer_params(psm%sourcetype,grid%psm_params,-i_sign)

         ...

         call psm_get(psm,rake=rake0,risetime=risetime )

         ...

         call texture_moment()   # dummy function

         call texture_slip_slipgen(ctrl_slipgen%sdin,ctrl_slipgen%kvalue)
         call texture_rake(rake0)
         call texture_risetime(risetime)! rake, risetime,moment.
     end do<span class="w">

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
          

        current_node=>list_head ; i_sign=1
         do while(associated(current_node))
             psm=>current_node%psm
             current_node=>current_node%next
             grid=>psm%cgrid


         call psm_transfer_params(psm%sourcetype,psm%params,i_sign) 
         ! update specific params structure

         grid%psm_sourcetype=psm%sourcetype 
         ! save a copy to grid%psm_params
         call psm_transfer_params(psm%sourcetype,grid%psm_params,-i_sign)

         call texture_moment()   # dummy function 

     end do<span class="w">

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
          

       ! count  Number of discrete-source
        current_node=>list_head ; is=0; i_sign=1
        do while(associated(current_node))
           …
       is=is+1
           …
        end do  ! while


   ! count Number of obs data.
    do ireceiver=1, nreceivers.
       nobs = nobs + span_ref/delta 
    end do 

   ! allocate memory  Ax=B 
      allocate(A,B,X, ...) 

   ! calculate green A.
      make_seismogram(sub_tdsm, ... )

   ! calculate obs B
      B(i:j)=strip%data(i′:j′)

   ! set_constrains
      .... 

   ! invoke
     call nnls(...)

   ! deallocate memory.
     ....

   ! save result.
     ....<span class="w">

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
          

      ! count  Number of discrete-source
        current_node=>list_head ; is=0; i_sign=1
        do while(associated(current_node))
           …
       is=is+1
           …
        end do  ! while


   ! count N-obs
     nobs=0; ndataset=size(obs_datasets)
     do loop=1,ndataset
        nobs=nobs+size((obs_datasets(loop)%obs))
     end do

   ! allocate memory
     allocate(AA,BB,XX,DD ...) 

   ! construct AA ... DD ...
     .... 

   ! call SVD
     ....

   ! deallocate momery
     ....

   ! save result
     ....<span class="w">

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
          

subroutine calculate_and_scale_seismogram()
   delta:
      db%dt,


   receiver::
      type(t_receiver), dimension(:), allocatable :: receivers {   !   receiver.f90
         type(t_strip),dimension(:) allocatable :: displacement {    ! sparse_trace.f90
              real(kind=real_kind), dimension(:), allocatable :: data   

         }
      }
   !—————————————————————————————————————–
   !             real, dimension(:) :: t_receiver :: receivers%t_strip displacement%data
   !   should update_to:

   !             real, dimension(:) :: t_receiver :: receivers%t_trace :: traces%t_strip :: stripsdata
   !—————————————————————————————————————–

   ! source ::
     type(t_list_node), pointer :: current_node = list_head  {
         type(t_list_node), pointer :: next=>NULL()
         type(t_psm) :: psm  {
              type(t_eikonal_grid) :: egrid {
                  grid%moments(nsd,nx,ny,nw)
              }
              type(t_eikonal_grid) :: cgrid
         }
         logical :: xxx_dirty,xxx_inited
     }
   !—————————————————————————————
   if(ctrl_inversion%imethod(2)==slip_inversion_method_slipgen) then
       if(.not.ctrl_slipgen%cgrid) grid=>psm%egrid   ! slipgen may save result in find grid        

   end if
   !—————————————————————————————
   grid=>cgrid or  grid=>egrid

! info:  step_1: if cgrid%duration > effective_dt.  then discretize each subfault to several part.
! info:  step_2: apply risetime
! info:  step_3:    

{
   do iy=1:ny,ix=1:nx

  ! step_1:
    duration=grid%durations(ix,iy) ; egrid%duration=0;  cgrid%duration
    call discretize_subfault_time( duration, 0.0,  effective_dt, t_weights, t_offsets,nsrc)    ! rectangular
    allocate(sub_tdsm%centroids(nsrc))

  ! step_2:
    if(ctrl_inversion%imethod(2)==slip_inversion_method_nnls) then  ! n triangular stf
                    call discretize_subfault_time(risetime,  deltat, tweights, toffsets,nt) ! for strip_fold
    else  ! one rectangular stf( source time function )
                    call discretize_subfault_time(0.0,risetime,  deltat, tweights, toffsets,nt) ! for strip_fold 
    end if

  ! step_3: 
         forall(it=1:nsrc) ! same geocoordinate.
                    sub_tdsm%centroids(it)%north = grid%points(1,ix,iy)
                    sub_tdsm%centroids(it)%east  = grid%points(2,ix,iy)
                    sub_tdsm%centroids(it)%depth = grid%points(3,ix,iy)
         end forall


      do iw=1,nw ! asa,slipgen,uniform:  nw=svd =1,   mtw: nw >= 1
                    !print *, 'debug minimizer_engine.f90 1566 iw',iw
          FORALL(it=1:nsrc) sub_tdsm%centroids(it)%time  =  &
                    grid%times(ix,iy) + &  ! times reference, relative time
                    origin_time-grid%center_time  + & ! due to origin%(time,x,y,z)  is centeroid location
                    t_offsets(it)       + & ! offset by nsrc        
                    (iw-1)*0.5*risetime + & ! offset by nw
                    real(psm%ref_time)   ! new ref scheme
      end do


   end do
}

end subroutine