
p r o g r a m  I S O _ 1 9 9 0 5 _ 1   w32 version 00.11     
------------------------------------   at 07:04:48 on 19-JUL-2024
jackup leg penetration assessment during preload according to ISO 19905-1 (2016)
(c) S.Kay Consultant, Leiderdorp, The Netherlands.
Please read the ISO_19905_1 Users Manual [ref SKA/MAN/001]
        and the ISO_19905_1 FAQ File     [ref SKA/FAQ/001].


 pidata: < --- BEGIN
 pidata: input  text file: "JOB.I01"
 pidata: output text file: "JOB.O01"
 pidata: output plot file: "Q23-AS-PWD_SOFF_CPT002_LB_ENSCO.O01.PDF"
 pidata: program name    : "ISO_19905_1"

 pidata: rtitle          : "Q23-AS-PWD_SOFF_CPT002_LB_ENSCO"
       : ptitle          : "QAT-0689"
       : projno          : "Q-NOC-QAT-0689"
       : locatn          : "Q23-AS-PWD_SOFF_CPT002_LB_ENSC"

 pidata: END --- >

 gmdata: < --- BEGIN
 numsl
   10
  zdep      bden     su_top    su_bot     phi  
   [m]   [kN/m^3]    [kPa]      [kPa]    [deg] 
  0.00      2.40     0.800      1.30      0.00    
 0.500      7.65      0.00      0.00      20.0    
  3.60      7.00      90.0      90.0      0.00    
  4.80      7.00      0.00      0.00      15.0    
  6.30      7.25      48.0      48.0      0.00    
  13.0      9.75      0.00      0.00      25.0    
  14.4      7.00      71.0      71.0      0.00    
  16.4      8.50      0.00      0.00      30.0    
  17.3      7.00      78.0      78.0      0.00    
  17.9      9.00      0.00      0.00      17.5    
  18.4    
 gmdata: END --- >

 spfdat: < --- BEGIN
 spname
 ENSCO110
 spgeom
 CIRCULAR
 v_spud_total [m^3]
  421.    
 p_spud [kN]
 0.669E+05
 numsp_BH
    3
  sp_B      sp_H 
   [m]       [m] 
  0.00      0.00    
  3.00      1.22    
  15.3      3.10    
 spfdat: END --- >

 ssidat: < --- BEGIN
 sp_alpha
 0.600    
 i_bf i sqc
    1    1
 ssidat: END --- >

 resdat: < --- BEGIN
 alpha_b_tan_phi   alpha_b_su
 0.250     0.500    
 a_sq   b_sq   r_sq
  5.00     0.330      1.50    
 num_n_s
    1
 v_n_s
  4.00    
 resdat: END --- >

 s__out: < --- BEGIN
 ioutd
    2
 z_tip_verify_min   z_tip_verify_max
  7.90      7.90    
 s__out: END --- >

spudcan_mat_X_sub: < --- BEGIN
   V_D_maxB                            [m^3]:       145.    
   V_D_missing                         [m^3]:       276.    
   sp_H_missing                          [m]:       4.51    
   sp_H                                  [m]:       7.61    
 > sp_H              sp_B
 > [m]               [m]
       0.00           0.00    
       3.00           1.22    
       15.3           3.10    
       0.00           7.61    
spudcan_mat_X_sub: END --- >


D_B_eq_sub: < --- BEGIN
 < z_tip                                 [m]:       7.90    
   sp_H_max                              [m]:       3.10    
   sp_B_max                              [m]:       15.3    
 > D                                     [m]:       4.80    
 > B_eq                                  [m]:       15.3    
D_B_eq_sub: END --- >


sp_beta_fun: < --- BEGIN
 < z_tip                                 [m]:       7.90    
   sp_H_max                              [m]:       3.10    
   z_tip_dummy                           [m]:       3.10    
   V_eq                                [m^3]:       145.    
   B_eq                                  [m]:       15.3    
   A_eq                                [m^2]:       184.    
   H_eq                                  [m]:       2.37    
   sp_H_max                              [m]:       2.54    
   z_tip_dummy                           [m]:       146.    
   V_eq                                [m^3]:
 > sp_beta                             [rad]:       2.54    
 > sp_beta/deg                         [deg]:       146.    
sp_beta_fun: END --- >

   gamma_sub_av                     [kN/m^3]:       6.94    
   A                                   [m^2]:       184.    
   H_cav                                 [m]:       4.30    
   V_spud                              [m^3]:       227.    
   V_D                                 [m^3]:       145.    
   W_BF = gamma_sub_av [ A (D - H_cav) - (V_spud - V_D ) ]
 > W_BF                                 [kN]:   68.10682    
W_BF_fun: END --- >


Qv_Dpt_sub: < --- BEGIN
 < D                                     [m]:       4.80    
 < B_eq                                  [m]:       15.3    

Q_ub_sub: < --- BEGIN     
 < D                                     [m]:       4.80    
 < D_pt                                  [m]:       4.80    
 < B_pt                                  [m]:       15.3    

Qv_GSS_fun: < --- BEGIN
 < D                                     [m]:       4.80    
 < B_eq                                  [m]:       15.3    

phi_av_fun: < --- BEGIN
 < z1                                    [m]:       4.80    
 < z2                                    [m]:       8.62    
   is_D_clay                             [-]:         F
   lay_num_z1                            [-]:         4
   lay_num_z2                            [-]:         5
   tan_phi_av_sum                    [deg-m]:      0.402    
   T_sum                                 [m]:       1.50    
 > phi_av                              [rad]:      0.262    
 > phi_av/deg                          [deg]:       15.0    
phi_av_fun : END --- >

   is_D_clay                           [t/f]:         F
   gamma_sub_av                     [kN/m^3]:       7.15    
   phi_rad_av                          [rad]:      0.262    
   phi_rad_av/deg                      [deg]:       15.0    
   N_gamma                               [-]:       1.80    
   N_q                                   [-]:       7.20    
   d_gamma                               [-]:       1.00    
   d_q                                   [-]:       1.09    
   p0                                  [kPa]:       33.3    
   gamma_sub_av*D_gamma*N_gamma*pi*B_eq*B_eq*B_eq/8
                                        [kN]:      0.181E+05
   p0*d_q*N_q*pi*B_eq*B_eq/4            [kN]:      0.480E+05
   Qv_GSS = gamma_sub_av*D_gamma*N_gamma*pi*B_eq*B_eq*B_eq/8 + p0*d_q*N_q*pi*B_eq*B_eq/4
 > Qv_GSS                               [kN]:   66154.04    
Qv_GSS_fun: END --- >

   phi_lay_deg1                        [deg]:       15.0    
   phi_lay_deg2                        [deg]:       15.0    
 > Q_ub                                 [kN]:   66154.04    
 > Qv_fmode_pt                              :GSS
Q_ub_sub: END --- >


Rs_pt_fun: < --- BEGIN
 < D                                     [m]:       4.80    
 < B_pt                                  [m]:       15.3    
 < T_pt                                  [m]:       0.00    
   D + T_pt                              [m]:       4.80    
   su_av                                kPa]:       0.00    
   P_pt                                  [m]:       0.00    
   T_pt                                  [m]:       0.00    
   Rs_pt = ( 3/4 ) * su_av * P_pt * T_pt;    
 > Rs_pt                                [kN]:   0.000000    
Rs_pt_fun: END --- >


dW_pt_fun: < --- BEGIN
 < D                                     [m]:       4.80    
 < B_pt                                  [m]:       15.3    
 < T_pt                                  [m]:       0.00    
   D + T_pt                              [m]:       4.80    
   A_pt                                [m^2]:       184.    
   gamma_sub_av                     [kN/m^3]:       7.00    
   dW_pt = A_pt * T_pt * gamma_sub_av;       
 > dW_pt = A_pt * T_pt * gamma_sub_av   [kN]:   0.000000    
dW_pt_fun: END --- >


Q_ub_sub: < --- BEGIN     
 < D                                     [m]:       4.80    
 < D_pt                                  [m]:       6.30    
 < B_pt                                  [m]:       16.1    

Qv_SQC_fun: < --- BEGIN
 < D                                     [m]:       6.30    
 < B_eq                                  [m]:       16.1    

ZTOP_sq_fun: < --- BEGIN
 < D                                     [m]:       6.30    
 < B_eq                                  [m]:
   z1                                [m bsf]:       6.30    
   z0                                [m bsf]:       0.00    
   lay_num_z1                            [-]:         5
   lay_num_z0                            [-]:         1
   r_sq                                  [-]:       1.50    
   su0_bot                             [kPa]:      0.100E-02
   su1_top                             [kPa]:      0.100E-02
   su1_top/su0_bot                       [-]:       1.00    
 > ZTOP_sq_fun                       [m bsf]:       6.30    
ZTOP_sq_fun: END --- >


T_sq_fun: < --- BEGIN
 < D                                     [m]:       6.30    
 < B_eq                                  [m]:       16.1    
   z1                                    [m]:       6.30    
   z2                                    [m]:       14.3    
   lay_num_z1                            [-]:         5
   lay_num_z2                            [-]:         6
   r_sq                                  [-]:       1.50    
   su1_bot                             [kPa]:       48.0    
   su2_top                             [kPa]:       48.0    
   su2_top/su1_bot                       [-]:       1.00    
 > T_sq                                  [m]:       6.70    
T_sq_fun: END --- >

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:       6.70    
   b_sq                                  [-]:      0.330    
   B_eq                                  [m]:       16.1    
   T_sq/b_sq                             [m]:       20.3    
 * B_eq > T_sq/b_sq                    [t/f]:         F
   ztop_sq                           [m bsf]:       6.30    
   B_eq                                  [m]:       16.1    
   ztop_sq/B_eq (i.e. 'D/B')             [-]:      0.393    
 * 'D/B' <= 2.5                        [t/f]:         T
   Dtip (ie D corresponding to tip depth)[m]:       4.80    
   Dtip/B                                [-]:      0.299    
 * is_SQC_occurring at Dtip/B=2.5      [t/f]:         F
   [B_eq > T_sq/b_sq] & ['D/B' <= 2.5]       
 * is_SQC_occurring                    [t/f]:         F

 > Qv_SQC                               [kN]:   1000000.    
Qv_SQC_fun: END --- >


Qv_GSC_fun: < --- BEGIN
 < D                                 [m bsf]:       6.30    
 < B_eq                                  [m]:       16.1    
 < sp_beta                             [rad]:       3.14    
 < sp_beta/deg                         [deg]:       180.    

rho_av_fun: < --- BEGIN     
 < D                                     [m]:       6.30    
 < B                                     [m]:       16.1    
   is_D_clay                           [t/f]:         T
   z1                                    [m]:       6.30    
   z2                                    [m]:       14.3    
   lay_num_z1                            [-]:         5
   lay_num_z2                            [-]:         6
   rho_av_sum                          [kPa]:       0.00    
   T_sum                                 [m]:       6.70    
 > rho_av                            [kPa/m]:       0.00    
rho_av_fun: END --- >

   is_D_clay                           [t/f]:         T
   rho_av                            [kPa/m]:       0.00    
   su0                                 [kPa]:       48.0    
   D                                     [m]:       6.30    
   z_top_first_CLAY_layer_aboveD     [m bsf]:       6.30    
   sum = su0 - rho_av*D                [kPa]:       48.0    
   sum_modified = max(sum,0)           [kPa]:       48.0    
   kappa_m = rho_av*B_eq/sum_modified    [-]:       0.00    
   kappa_m_modified = max(kappa_m,5)     [-]:       0.00    
   B_eq                                  [m]:       16.1    
   lambda = D/B_eq                       [-]:      0.393    
   kappa_0=kappa_mod/(lambda*kappa_mod+1)[-]:       0.00    
   sp_beta/deg                         [deg]:       180.    
   sp_alpha                              [-]:      0.600    
   Nc0                                   [-]:       6.87    
   su0 * Nc0                           [kPa]:       330.    
   p0                                  [kPa]:       43.8    
   su0*Nc0 + p0                        [kPa]:       373.    
   pi * B_eq * B_eq /4                 [m^2]:       202.    
   Qv_GSC = ( su0*Nc0 + p0 ) * pi * B_eq * B_eq / 4 
 > Qv_GSC                               [kN]:   75561.31    
Qv_GSC_fun: END --- >

   phi_lay_deg1                        [deg]:       15.0    
   phi_lay_deg2                        [deg]:       0.00    
 > Q_ub                                 [kN]:   75561.31    
 > Qv_fmode_pt                              :GSC
Q_ub_sub: END --- >


Rs_pt_fun: < --- BEGIN
 < D                                     [m]:       4.80    
 < B_pt                                  [m]:       16.1    
 < T_pt                                  [m]:       1.50    
   D + T_pt                              [m]:       6.30    
   su_av                                kPa]:       0.00    
   P_pt                                  [m]:       0.00    
   T_pt                                  [m]:       1.50    
   Rs_pt = ( 3/4 ) * su_av * P_pt * T_pt;    
 > Rs_pt                                [kN]:   0.000000    
Rs_pt_fun: END --- >


dW_pt_fun: < --- BEGIN
 < D                                     [m]:       4.80    
 < B_pt                                  [m]:       16.1    
 < T_pt                                  [m]:       1.50    
   D + T_pt                              [m]:       6.30    
   A_pt                                [m^2]:       202.    
   gamma_sub_av                     [kN/m^3]:       7.00    
   dW_pt = A_pt * T_pt * gamma_sub_av;       
 > dW_pt = A_pt * T_pt * gamma_sub_av   [kN]:   2124.563    
dW_pt_fun: END --- >


Q_ub_sub: < --- BEGIN     
 < D                                     [m]:       6.30    
 < D_pt                                  [m]:       13.0    
 < B_pt                                  [m]:       16.1    

Qv_GSS_fun: < --- BEGIN
 < D                                     [m]:       13.0    
 < B_eq                                  [m]:       16.1    

phi_av_fun: < --- BEGIN
 < z1                                    [m]:       13.0    
 < z2                                    [m]:       17.0    
   is_D_clay                             [-]:         F
   lay_num_z1                            [-]:         6
   lay_num_z2                            [-]:         8
   tan_phi_av_sum                    [deg-m]:       1.01    
   T_sum                                 [m]:       2.01    
 > phi_av                              [rad]:      0.464    
 > phi_av/deg                          [deg]:       26.6    
phi_av_fun : END --- >

   is_D_clay                           [t/f]:         F
   gamma_sub_av                     [kN/m^3]:       8.19    
   phi_rad_av                          [rad]:      0.464    
   phi_rad_av/deg                      [deg]:       26.6    
   N_gamma                               [-]:       8.15    
   N_q                                   [-]:       22.9    
   d_gamma                               [-]:       1.00    
   d_q                                   [-]:       1.21    
   p0                                  [kPa]:       92.4    
   gamma_sub_av*D_gamma*N_gamma*pi*B_eq*B_eq*B_eq/8
                                        [kN]:      0.108E+06
   p0*d_q*N_q*pi*B_eq*B_eq/4            [kN]:      0.516E+06
   Qv_GSS = gamma_sub_av*D_gamma*N_gamma*pi*B_eq*B_eq*B_eq/8 + p0*d_q*N_q*pi*B_eq*B_eq/4
 > Qv_GSS                               [kN]:   624740.1    
Qv_GSS_fun: END --- >

   phi_lay_deg1                        [deg]:       0.00    
   phi_lay_deg2                        [deg]:       25.0    
 > Q_ub                                 [kN]:   624740.1    
 > Qv_fmode_pt                              :GSS
Q_ub_sub: END --- >


Rs_pt_fun: < --- BEGIN
 < D                                     [m]:       6.30    
 < B_pt                                  [m]:       16.1    
 < T_pt                                  [m]:       6.70    
   D + T_pt                              [m]:       13.0    
   su_av                                kPa]:       699.    
   P_pt                                  [m]:       50.4    
   T_pt                                  [m]:       6.70    
   Rs_pt = ( 3/4 ) * su_av * P_pt * T_pt;    
 > Rs_pt                                [kN]:   177110.6    
Rs_pt_fun: END --- >


dW_pt_fun: < --- BEGIN
 < D                                     [m]:       6.30    
 < B_pt                                  [m]:       16.1    
 < T_pt                                  [m]:       6.70    
   D + T_pt                              [m]:       13.0    
   A_pt                                [m^2]:       202.    
   gamma_sub_av                     [kN/m^3]:       7.25    
   dW_pt = A_pt * T_pt * gamma_sub_av;       
 > dW_pt = A_pt * T_pt * gamma_sub_av   [kN]:   9829.926    
dW_pt_fun: END --- >


Q_ub_sub: < --- BEGIN     
 < D                                     [m]:       13.0    
 < D_pt                                  [m]:       14.4    
 < B_pt                                  [m]:       16.8    

Qv_SQC_fun: < --- BEGIN
 < D                                     [m]:       14.4    
 < B_eq                                  [m]:       16.8    

ZTOP_sq_fun: < --- BEGIN
 < D                                     [m]:       14.4    
 < B_eq                                  [m]:
   z1                                [m bsf]:       14.4    
   z0                                [m bsf]:       0.00    
   lay_num_z1                            [-]:         7
   lay_num_z0                            [-]:         1
   r_sq                                  [-]:       1.50    
   su0_bot                             [kPa]:      0.100E-02
   su1_top                             [kPa]:      0.100E-02
   su1_top/su0_bot                       [-]:       1.00    
 > ZTOP_sq_fun                       [m bsf]:       14.4    
ZTOP_sq_fun: END --- >


T_sq_fun: < --- BEGIN
 < D                                     [m]:       14.4    
 < B_eq                                  [m]:       16.8    
   z1                                    [m]:       14.4    
   z2                                    [m]:       18.4    
   lay_num_z1                            [-]:         7
   lay_num_z2                            [-]:        10
   r_sq                                  [-]:       1.50    
   su1_bot                             [kPa]:       71.0    
   su2_top                             [kPa]:       71.0    
   su2_top/su1_bot                       [-]:       1.00    
 > T_sq                                  [m]:       2.00    
T_sq_fun: END --- >

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:       2.00    
   b_sq                                  [-]:      0.330    
   B_eq                                  [m]:       16.8    
   T_sq/b_sq                             [m]:       6.06    
 * B_eq > T_sq/b_sq                    [t/f]:         T
   ztop_sq                           [m bsf]:       14.4    
   B_eq                                  [m]:       16.8    
   ztop_sq/B_eq (i.e. 'D/B')             [-]:      0.860    
 * 'D/B' <= 2.5                        [t/f]:         T
   Dtip (ie D corresponding to tip depth)[m]:       4.80    
   Dtip/B                                [-]:      0.287    
 * is_SQC_occurring at Dtip/B=2.5      [t/f]:         F
   [B_eq > T_sq/b_sq] & ['D/B' <= 2.5]       
 * is_SQC_occurring                    [t/f]:         T

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:       2.00    
   D + T_sq                              [m]:       16.4    
   a_sq                                  [-]:       5.00    
   b_sq                                  [-]:      0.330    
   a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq    [-]:       8.80    
   su_av                               [kPa]:       71.0    
   p'0                                 [kPa]:       106.    
   (a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su0 [kPa]:       624.    
   pi*B_eq*B_eq/4                      [m^2]:       220.    
   Qv_SQC = (( a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su_av + p'0 ) *pi*B_eq*B_eq/4
 > Qv_SQC                               [kN]:   160970.9    
Qv_SQC_fun: END --- >


Qv_GSC_fun: < --- BEGIN
 < D                                 [m bsf]:       14.4    
 < B_eq                                  [m]:       16.8    
 < sp_beta                             [rad]:       3.14    
 < sp_beta/deg                         [deg]:       180.    

rho_av_fun: < --- BEGIN     
 < D                                     [m]:       14.4    
 < B                                     [m]:       16.8    
   is_D_clay                           [t/f]:         T
   z1                                    [m]:       14.4    
   z2                                    [m]:       18.4    
   lay_num_z1                            [-]:         7
   lay_num_z2                            [-]:        10
   rho_av_sum                          [kPa]:       0.00    
   T_sum                                 [m]:       2.60    
 > rho_av                            [kPa/m]:       0.00    
rho_av_fun: END --- >

   is_D_clay                           [t/f]:         T
   rho_av                            [kPa/m]:       0.00    
   su0                                 [kPa]:       71.0    
   D                                     [m]:       14.4    
   z_top_first_CLAY_layer_aboveD     [m bsf]:       14.4    
   sum = su0 - rho_av*D                [kPa]:       71.0    
   sum_modified = max(sum,0)           [kPa]:       71.0    
   kappa_m = rho_av*B_eq/sum_modified    [-]:       0.00    
   kappa_m_modified = max(kappa_m,5)     [-]:       0.00    
   B_eq                                  [m]:       16.8    
   lambda = D/B_eq                       [-]:      0.860    
   kappa_0=kappa_mod/(lambda*kappa_mod+1)[-]:       0.00    
   sp_beta/deg                         [deg]:       180.    
   sp_alpha                              [-]:      0.600    
   Nc0                                   [-]:       7.52    
   su0 * Nc0                           [kPa]:       534.    
   p0                                  [kPa]:       106.    
   su0*Nc0 + p0                        [kPa]:       640.    
   pi * B_eq * B_eq /4                 [m^2]:       220.    
   Qv_GSC = ( su0*Nc0 + p0 ) * pi * B_eq * B_eq / 4 
 > Qv_GSC                               [kN]:   140951.0    
Qv_GSC_fun: END --- >

   phi_lay_deg1                        [deg]:       25.0    
   phi_lay_deg2                        [deg]:       0.00    
 > Q_ub                                 [kN]:   160970.9    
 > Qv_fmode_pt                              :SQC
Q_ub_sub: END --- >


Rs_pt_fun: < --- BEGIN
 < D                                     [m]:       13.0    
 < B_pt                                  [m]:       16.8    
 < T_pt                                  [m]:       1.40    
   D + T_pt                              [m]:       14.4    
   su_av                                kPa]:       0.00    
   P_pt                                  [m]:       0.00    
   T_pt                                  [m]:       1.40    
   Rs_pt = ( 3/4 ) * su_av * P_pt * T_pt;    
 > Rs_pt                                [kN]:   0.000000    
Rs_pt_fun: END --- >


dW_pt_fun: < --- BEGIN
 < D                                     [m]:       13.0    
 < B_pt                                  [m]:       16.8    
 < T_pt                                  [m]:       1.40    
   D + T_pt                              [m]:       14.4    
   A_pt                                [m^2]:       220.    
   gamma_sub_av                     [kN/m^3]:       9.75    
   dW_pt = A_pt * T_pt * gamma_sub_av;       
 > dW_pt = A_pt * T_pt * gamma_sub_av   [kN]:   3007.270    
dW_pt_fun: END --- >


Qv_SQC_fun: < --- BEGIN
 < D                                     [m]:       14.4    
 < B_eq                                  [m]:       16.8    

ZTOP_sq_fun: < --- BEGIN
 < D                                     [m]:       14.4    
 < B_eq                                  [m]:
   z1                                [m bsf]:       14.4    
   z0                                [m bsf]:       0.00    
   lay_num_z1                            [-]:         7
   lay_num_z0                            [-]:         1
   r_sq                                  [-]:       1.50    
   su0_bot                             [kPa]:      0.100E-02
   su1_top                             [kPa]:      0.100E-02
   su1_top/su0_bot                       [-]:       1.00    
 > ZTOP_sq_fun                       [m bsf]:       14.4    
ZTOP_sq_fun: END --- >


T_sq_fun: < --- BEGIN
 < D                                     [m]:       14.4    
 < B_eq                                  [m]:       16.8    
   z1                                    [m]:       14.4    
   z2                                    [m]:       18.4    
   lay_num_z1                            [-]:         7
   lay_num_z2                            [-]:        10
   r_sq                                  [-]:       1.50    
   su1_bot                             [kPa]:       71.0    
   su2_top                             [kPa]:       71.0    
   su2_top/su1_bot                       [-]:       1.00    
 > T_sq                                  [m]:       2.00    
T_sq_fun: END --- >

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:       2.00    
   b_sq                                  [-]:      0.330    
   B_eq                                  [m]:       16.8    
   T_sq/b_sq                             [m]:       6.06    
 * B_eq > T_sq/b_sq                    [t/f]:         T
   ztop_sq                           [m bsf]:       14.4    
   B_eq                                  [m]:       16.8    
   ztop_sq/B_eq (i.e. 'D/B')             [-]:      0.860    
 * 'D/B' <= 2.5                        [t/f]:         T
   Dtip (ie D corresponding to tip depth)[m]:       4.80    
   Dtip/B                                [-]:      0.287    
 * is_SQC_occurring at Dtip/B=2.5      [t/f]:         F
   [B_eq > T_sq/b_sq] & ['D/B' <= 2.5]       
 * is_SQC_occurring                    [t/f]:         T

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:       2.00    
   D + T_sq                              [m]:       16.4    
   a_sq                                  [-]:       5.00    
   b_sq                                  [-]:      0.330    
   a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq    [-]:       8.80    
   su_av                               [kPa]:       71.0    
   p'0                                 [kPa]:       106.    
   (a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su0 [kPa]:       624.    
   pi*B_eq*B_eq/4                      [m^2]:       220.    
   Qv_SQC = (( a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su_av + p'0 ) *pi*B_eq*B_eq/4
 > Qv_SQC                               [kN]:   160970.9    
Qv_SQC_fun: END --- >


Qv_SQC_fun: < --- BEGIN
 < D                                     [m]:       14.4    
 < B_eq                                  [m]:       15.3    

ZTOP_sq_fun: < --- BEGIN
 < D                                     [m]:       14.4    
 < B_eq                                  [m]:
   z1                                [m bsf]:       14.4    
   z0                                [m bsf]:       0.00    
   lay_num_z1                            [-]:         7
   lay_num_z0                            [-]:         1
   r_sq                                  [-]:       1.50    
   su0_bot                             [kPa]:      0.100E-02
   su1_top                             [kPa]:      0.100E-02
   su1_top/su0_bot                       [-]:       1.00    
 > ZTOP_sq_fun                       [m bsf]:       14.4    
ZTOP_sq_fun: END --- >


T_sq_fun: < --- BEGIN
 < D                                     [m]:       14.4    
 < B_eq                                  [m]:       15.3    
   z1                                    [m]:       14.4    
   z2                                    [m]:       18.4    
   lay_num_z1                            [-]:         7
   lay_num_z2                            [-]:        10
   r_sq                                  [-]:       1.50    
   su1_bot                             [kPa]:       71.0    
   su2_top                             [kPa]:       71.0    
   su2_top/su1_bot                       [-]:       1.00    
 > T_sq                                  [m]:       2.00    
T_sq_fun: END --- >

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:       2.00    
   b_sq                                  [-]:      0.330    
   B_eq                                  [m]:       15.3    
   T_sq/b_sq                             [m]:       6.06    
 * B_eq > T_sq/b_sq                    [t/f]:         T
   ztop_sq                           [m bsf]:       14.4    
   B_eq                                  [m]:       15.3    
   ztop_sq/B_eq (i.e. 'D/B')             [-]:      0.941    
 * 'D/B' <= 2.5                        [t/f]:         T
   Dtip (ie D corresponding to tip depth)[m]:       4.80    
   Dtip/B                                [-]:      0.314    
 * is_SQC_occurring at Dtip/B=2.5      [t/f]:         F
   [B_eq > T_sq/b_sq] & ['D/B' <= 2.5]       
 * is_SQC_occurring                    [t/f]:         T

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:       2.00    
   D + T_sq                              [m]:       16.4    
   a_sq                                  [-]:       5.00    
   b_sq                                  [-]:      0.330    
   a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq    [-]:       8.65    
   su_av                               [kPa]:       71.0    
   p'0                                 [kPa]:       106.    
   (a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su0 [kPa]:       614.    
   pi*B_eq*B_eq/4                      [m^2]:       184.    
   Qv_SQC = (( a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su_av + p'0 ) *pi*B_eq*B_eq/4
 > Qv_SQC                               [kN]:   132460.8    
Qv_SQC_fun: END --- >


Rs_pt_fun: < --- BEGIN
 < D                                     [m]:       13.0    
 < B_pt                                  [m]:       15.3    
 < T_pt                                  [m]:       1.40    
   D + T_pt                              [m]:       14.4    
   su_av                                kPa]:       0.00    
   P_pt                                  [m]:       0.00    
   T_pt                                  [m]:       1.40    
   Rs_pt = ( 3/4 ) * su_av * P_pt * T_pt;    
 > Rs_pt                                [kN]:   0.000000    
Rs_pt_fun: END --- >


dW_pt_fun: < --- BEGIN
 < D                                     [m]:       13.0    
 < B_pt                                  [m]:       15.3    
 < T_pt                                  [m]:       1.40    
   D + T_pt                              [m]:       14.4    
   A_pt                                [m^2]:       184.    
   gamma_sub_av                     [kN/m^3]:       9.75    
   dW_pt = A_pt * T_pt * gamma_sub_av;       
 > dW_pt = A_pt * T_pt * gamma_sub_av   [kN]:   2509.144    
dW_pt_fun: END --- >


   D, D_pt                               [m]:       13.0           14.4    
   D_pt layer #                          [-]:    4
   Qv_Fmode_PT                           [-]:PTS>SQC
   decreasing dW for punch-through layer#[-]:    3
   using B_eq (was B_pt), B_eq , B_pt    [m]:       15.3           16.8    
   (B_eq/B_pt)^2                         [-]:      0.834    
   old dW,     new dW                   [kN]:      0.301E+04      0.251E+04
   old Qv_SQC, new Qv_SQC (ie Q_ub)     [kN]:      0.161E+06      0.132E+06

Q_ub_sub: < --- BEGIN     
 < D                                     [m]:       14.4    
 < D_pt                                  [m]:       16.4    
 < B_pt                                  [m]:       16.8    

Qv_GSS_fun: < --- BEGIN
 < D                                     [m]:       16.4    
 < B_eq                                  [m]:       16.8    

phi_av_fun: < --- BEGIN
 < z1                                    [m]:       16.4    
 < z2                                    [m]:       18.4    
   is_D_clay                             [-]:         F
   lay_num_z1                            [-]:         8
   lay_num_z2                            [-]:        10
   tan_phi_av_sum                    [deg-m]:      0.677    
   T_sum                                 [m]:       1.40    
 > phi_av                              [rad]:      0.451    
 > phi_av/deg                          [deg]:       25.8    
phi_av_fun : END --- >

   is_D_clay                           [t/f]:         F
   gamma_sub_av                     [kN/m^3]:       8.18    
   phi_rad_av                          [rad]:      0.451    
   phi_rad_av/deg                      [deg]:       25.8    
   N_gamma                               [-]:       7.08    
   N_q                                   [-]:       20.6    
   d_gamma                               [-]:       1.00    
   d_q                                   [-]:       1.24    
   p0                                  [kPa]:       120.    
   gamma_sub_av*D_gamma*N_gamma*pi*B_eq*B_eq*B_eq/8
                                        [kN]:      0.107E+06
   p0*d_q*N_q*pi*B_eq*B_eq/4            [kN]:      0.675E+06
   Qv_GSS = gamma_sub_av*D_gamma*N_gamma*pi*B_eq*B_eq*B_eq/8 + p0*d_q*N_q*pi*B_eq*B_eq/4
 > Qv_GSS                               [kN]:   782035.7    
Qv_GSS_fun: END --- >

   phi_lay_deg1                        [deg]:       0.00    
   phi_lay_deg2                        [deg]:       30.0    
 > Q_ub                                 [kN]:   782035.7    
 > Qv_fmode_pt                              :GSS
Q_ub_sub: END --- >


Rs_pt_fun: < --- BEGIN
 < D                                     [m]:       14.4    
 < B_pt                                  [m]:       16.8    
 < T_pt                                  [m]:       2.00    
   D + T_pt                              [m]:       16.4    
   su_av                                kPa]:       722.    
   P_pt                                  [m]:       52.6    
   T_pt                                  [m]:       2.00    
   Rs_pt = ( 3/4 ) * su_av * P_pt * T_pt;    
 > Rs_pt                                [kN]:   56988.92    
Rs_pt_fun: END --- >


dW_pt_fun: < --- BEGIN
 < D                                     [m]:       14.4    
 < B_pt                                  [m]:       16.8    
 < T_pt                                  [m]:       2.00    
   D + T_pt                              [m]:       16.4    
   A_pt                                [m^2]:       220.    
   gamma_sub_av                     [kN/m^3]:       7.00    
   dW_pt = A_pt * T_pt * gamma_sub_av;       
 > dW_pt = A_pt * T_pt * gamma_sub_av   [kN]:   3085.376    
dW_pt_fun: END --- >


Q_ub_sub: < --- BEGIN     
 < D                                     [m]:       16.4    
 < D_pt                                  [m]:       17.3    
 < B_pt                                  [m]:       17.2    

Qv_SQC_fun: < --- BEGIN
 < D                                     [m]:       17.3    
 < B_eq                                  [m]:       17.2    

ZTOP_sq_fun: < --- BEGIN
 < D                                     [m]:       17.3    
 < B_eq                                  [m]:
   z1                                [m bsf]:       17.3    
   z0                                [m bsf]:       0.00    
   lay_num_z1                            [-]:         9
   lay_num_z0                            [-]:         1
   r_sq                                  [-]:       1.50    
   su0_bot                             [kPa]:      0.100E-02
   su1_top                             [kPa]:      0.100E-02
   su1_top/su0_bot                       [-]:       1.00    
 > ZTOP_sq_fun                       [m bsf]:       17.3    
ZTOP_sq_fun: END --- >


T_sq_fun: < --- BEGIN
 < D                                     [m]:       17.3    
 < B_eq                                  [m]:       17.2    
   z1                                    [m]:       17.3    
   z2                                    [m]:       18.4    
   lay_num_z1                            [-]:         9
   lay_num_z2                            [-]:        10
   r_sq                                  [-]:       1.50    
   su1_bot                             [kPa]:       78.0    
   su2_top                             [kPa]:       78.0    
   su2_top/su1_bot                       [-]:       1.00    
 > T_sq                                  [m]:      0.600    
T_sq_fun: END --- >

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:      0.600    
   b_sq                                  [-]:      0.330    
   B_eq                                  [m]:       17.2    
   T_sq/b_sq                             [m]:       1.82    
 * B_eq > T_sq/b_sq                    [t/f]:         T
   ztop_sq                           [m bsf]:       17.3    
   B_eq                                  [m]:       17.2    
   ztop_sq/B_eq (i.e. 'D/B')             [-]:       1.01    
 * 'D/B' <= 2.5                        [t/f]:         T
   Dtip (ie D corresponding to tip depth)[m]:       4.80    
   Dtip/B                                [-]:      0.279    
 * is_SQC_occurring at Dtip/B=2.5      [t/f]:         F
   [B_eq > T_sq/b_sq] & ['D/B' <= 2.5]       
 * is_SQC_occurring                    [t/f]:         T

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:      0.600    
   D + T_sq                              [m]:       17.9    
   a_sq                                  [-]:       5.00    
   b_sq                                  [-]:      0.330    
   a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq    [-]:       15.7    
   su_av                               [kPa]:       78.0    
   p'0                                 [kPa]:       128.    
   (a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su0 [kPa]:      0.122E+04
   pi*B_eq*B_eq/4                      [m^2]:       232.    
   Qv_SQC = (( a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su_av + p'0 ) *pi*B_eq*B_eq/4
 > Qv_SQC                               [kN]:   313609.1    
Qv_SQC_fun: END --- >


Qv_GSC_fun: < --- BEGIN
 < D                                 [m bsf]:       17.3    
 < B_eq                                  [m]:       17.2    
 < sp_beta                             [rad]:       3.14    
 < sp_beta/deg                         [deg]:       180.    

rho_av_fun: < --- BEGIN     
 < D                                     [m]:       17.3    
 < B                                     [m]:       17.2    
   is_D_clay                           [t/f]:         T
   z1                                    [m]:       17.3    
   z2                                    [m]:       18.4    
   lay_num_z1                            [-]:         9
   lay_num_z2                            [-]:        10
   rho_av_sum                          [kPa]:       0.00    
   T_sum                                 [m]:      0.600    
 > rho_av                            [kPa/m]:       0.00    
rho_av_fun: END --- >

   is_D_clay                           [t/f]:         T
   rho_av                            [kPa/m]:       0.00    
   su0                                 [kPa]:       78.0    
   D                                     [m]:       17.3    
   z_top_first_CLAY_layer_aboveD     [m bsf]:       17.3    
   sum = su0 - rho_av*D                [kPa]:       78.0    
   sum_modified = max(sum,0)           [kPa]:       78.0    
   kappa_m = rho_av*B_eq/sum_modified    [-]:       0.00    
   kappa_m_modified = max(kappa_m,5)     [-]:       0.00    
   B_eq                                  [m]:       17.2    
   lambda = D/B_eq                       [-]:       1.01    
   kappa_0=kappa_mod/(lambda*kappa_mod+1)[-]:       0.00    
   sp_beta/deg                         [deg]:       180.    
   sp_alpha                              [-]:      0.600    
   Nc0                                   [-]:       7.70    
   su0 * Nc0                           [kPa]:       600.    
   p0                                  [kPa]:       128.    
   su0*Nc0 + p0                        [kPa]:       728.    
   pi * B_eq * B_eq /4                 [m^2]:       232.    
   Qv_GSC = ( su0*Nc0 + p0 ) * pi * B_eq * B_eq / 4 
 > Qv_GSC                               [kN]:   169178.6    
Qv_GSC_fun: END --- >

   phi_lay_deg1                        [deg]:       30.0    
   phi_lay_deg2                        [deg]:       0.00    
 > Q_ub                                 [kN]:   313609.1    
 > Qv_fmode_pt                              :SQC
Q_ub_sub: END --- >


Rs_pt_fun: < --- BEGIN
 < D                                     [m]:       16.4    
 < B_pt                                  [m]:       17.2    
 < T_pt                                  [m]:      0.900    
   D + T_pt                              [m]:       17.3    
   su_av                                kPa]:       0.00    
   P_pt                                  [m]:       0.00    
   T_pt                                  [m]:      0.900    
   Rs_pt = ( 3/4 ) * su_av * P_pt * T_pt;    
 > Rs_pt                                [kN]:   0.000000    
Rs_pt_fun: END --- >


dW_pt_fun: < --- BEGIN
 < D                                     [m]:       16.4    
 < B_pt                                  [m]:       17.2    
 < T_pt                                  [m]:      0.900    
   D + T_pt                              [m]:       17.3    
   A_pt                                [m^2]:       232.    
   gamma_sub_av                     [kN/m^3]:       8.50    
   dW_pt = A_pt * T_pt * gamma_sub_av;       
 > dW_pt = A_pt * T_pt * gamma_sub_av   [kN]:   1777.290    
dW_pt_fun: END --- >


Qv_SQC_fun: < --- BEGIN
 < D                                     [m]:       17.3    
 < B_eq                                  [m]:       17.2    

ZTOP_sq_fun: < --- BEGIN
 < D                                     [m]:       17.3    
 < B_eq                                  [m]:
   z1                                [m bsf]:       17.3    
   z0                                [m bsf]:       0.00    
   lay_num_z1                            [-]:         9
   lay_num_z0                            [-]:         1
   r_sq                                  [-]:       1.50    
   su0_bot                             [kPa]:      0.100E-02
   su1_top                             [kPa]:      0.100E-02
   su1_top/su0_bot                       [-]:       1.00    
 > ZTOP_sq_fun                       [m bsf]:       17.3    
ZTOP_sq_fun: END --- >


T_sq_fun: < --- BEGIN
 < D                                     [m]:       17.3    
 < B_eq                                  [m]:       17.2    
   z1                                    [m]:       17.3    
   z2                                    [m]:       18.4    
   lay_num_z1                            [-]:         9
   lay_num_z2                            [-]:        10
   r_sq                                  [-]:       1.50    
   su1_bot                             [kPa]:       78.0    
   su2_top                             [kPa]:       78.0    
   su2_top/su1_bot                       [-]:       1.00    
 > T_sq                                  [m]:      0.600    
T_sq_fun: END --- >

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:      0.600    
   b_sq                                  [-]:      0.330    
   B_eq                                  [m]:       17.2    
   T_sq/b_sq                             [m]:       1.82    
 * B_eq > T_sq/b_sq                    [t/f]:         T
   ztop_sq                           [m bsf]:       17.3    
   B_eq                                  [m]:       17.2    
   ztop_sq/B_eq (i.e. 'D/B')             [-]:       1.01    
 * 'D/B' <= 2.5                        [t/f]:         T
   Dtip (ie D corresponding to tip depth)[m]:       4.80    
   Dtip/B                                [-]:      0.279    
 * is_SQC_occurring at Dtip/B=2.5      [t/f]:         F
   [B_eq > T_sq/b_sq] & ['D/B' <= 2.5]       
 * is_SQC_occurring                    [t/f]:         T

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:      0.600    
   D + T_sq                              [m]:       17.9    
   a_sq                                  [-]:       5.00    
   b_sq                                  [-]:      0.330    
   a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq    [-]:       15.7    
   su_av                               [kPa]:       78.0    
   p'0                                 [kPa]:       128.    
   (a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su0 [kPa]:      0.122E+04
   pi*B_eq*B_eq/4                      [m^2]:       232.    
   Qv_SQC = (( a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su_av + p'0 ) *pi*B_eq*B_eq/4
 > Qv_SQC                               [kN]:   313609.1    
Qv_SQC_fun: END --- >


Qv_SQC_fun: < --- BEGIN
 < D                                     [m]:       17.3    
 < B_eq                                  [m]:       15.3    

ZTOP_sq_fun: < --- BEGIN
 < D                                     [m]:       17.3    
 < B_eq                                  [m]:
   z1                                [m bsf]:       17.3    
   z0                                [m bsf]:       0.00    
   lay_num_z1                            [-]:         9
   lay_num_z0                            [-]:         1
   r_sq                                  [-]:       1.50    
   su0_bot                             [kPa]:      0.100E-02
   su1_top                             [kPa]:      0.100E-02
   su1_top/su0_bot                       [-]:       1.00    
 > ZTOP_sq_fun                       [m bsf]:       17.3    
ZTOP_sq_fun: END --- >


T_sq_fun: < --- BEGIN
 < D                                     [m]:       17.3    
 < B_eq                                  [m]:       15.3    
   z1                                    [m]:       17.3    
   z2                                    [m]:       18.4    
   lay_num_z1                            [-]:         9
   lay_num_z2                            [-]:        10
   r_sq                                  [-]:       1.50    
   su1_bot                             [kPa]:       78.0    
   su2_top                             [kPa]:       78.0    
   su2_top/su1_bot                       [-]:       1.00    
 > T_sq                                  [m]:      0.600    
T_sq_fun: END --- >

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:      0.600    
   b_sq                                  [-]:      0.330    
   B_eq                                  [m]:       15.3    
   T_sq/b_sq                             [m]:       1.82    
 * B_eq > T_sq/b_sq                    [t/f]:         T
   ztop_sq                           [m bsf]:       17.3    
   B_eq                                  [m]:       15.3    
   ztop_sq/B_eq (i.e. 'D/B')             [-]:       1.13    
 * 'D/B' <= 2.5                        [t/f]:         T
   Dtip (ie D corresponding to tip depth)[m]:       4.80    
   Dtip/B                                [-]:      0.314    
 * is_SQC_occurring at Dtip/B=2.5      [t/f]:         F
   [B_eq > T_sq/b_sq] & ['D/B' <= 2.5]       
 * is_SQC_occurring                    [t/f]:         T

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:      0.600    
   D + T_sq                              [m]:       17.9    
   a_sq                                  [-]:       5.00    
   b_sq                                  [-]:      0.330    
   a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq    [-]:       14.8    
   su_av                               [kPa]:       78.0    
   p'0                                 [kPa]:       128.    
   (a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su0 [kPa]:      0.115E+04
   pi*B_eq*B_eq/4                      [m^2]:       184.    
   Qv_SQC = (( a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su_av + p'0 ) *pi*B_eq*B_eq/4
 > Qv_SQC                               [kN]:   235313.7    
Qv_SQC_fun: END --- >


Rs_pt_fun: < --- BEGIN
 < D                                     [m]:       16.4    
 < B_pt                                  [m]:       15.3    
 < T_pt                                  [m]:      0.900    
   D + T_pt                              [m]:       17.3    
   su_av                                kPa]:       0.00    
   P_pt                                  [m]:       0.00    
   T_pt                                  [m]:      0.900    
   Rs_pt = ( 3/4 ) * su_av * P_pt * T_pt;    
 > Rs_pt                                [kN]:   0.000000    
Rs_pt_fun: END --- >


dW_pt_fun: < --- BEGIN
 < D                                     [m]:       16.4    
 < B_pt                                  [m]:       15.3    
 < T_pt                                  [m]:      0.900    
   D + T_pt                              [m]:       17.3    
   A_pt                                [m^2]:       184.    
   gamma_sub_av                     [kN/m^3]:       8.50    
   dW_pt = A_pt * T_pt * gamma_sub_av;       
 > dW_pt = A_pt * T_pt * gamma_sub_av   [kN]:   1406.320    
dW_pt_fun: END --- >


   D, D_pt                               [m]:       16.4           17.3    
   D_pt layer #                          [-]:    6
   Qv_Fmode_PT                           [-]:PTS>SQC
   decreasing dW for punch-through layer#[-]:    5
   using B_eq (was B_pt), B_eq , B_pt    [m]:       15.3           17.2    
   (B_eq/B_pt)^2                         [-]:      0.791    
   old dW,     new dW                   [kN]:      0.178E+04      0.141E+04
   old Qv_SQC, new Qv_SQC (ie Q_ub)     [kN]:      0.314E+06      0.235E+06

Q_ub_sub: < --- BEGIN     
 < D                                     [m]:       17.3    
 < D_pt                                  [m]:       17.9    
 < B_pt                                  [m]:       17.2    

Qv_GSS_fun: < --- BEGIN
 < D                                     [m]:       17.9    
 < B_eq                                  [m]:       17.2    

phi_av_fun: < --- BEGIN
 < z1                                    [m]:       17.9    
 < z2                                    [m]:       18.4    
   is_D_clay                             [-]:         F
   lay_num_z1                            [-]:        10
   lay_num_z2                            [-]:        10
   tan_phi_av_sum                    [deg-m]:      0.158    
   T_sum                                 [m]:      0.500    
 > phi_av                              [rad]:      0.305    
 > phi_av/deg                          [deg]:       17.5    
phi_av_fun : END --- >

   is_D_clay                           [t/f]:         F
   gamma_sub_av                     [kN/m^3]:       9.00    
   phi_rad_av                          [rad]:      0.305    
   phi_rad_av/deg                      [deg]:       17.5    
   N_gamma                               [-]:       2.10    
   N_q                                   [-]:       8.40    
   d_gamma                               [-]:       1.00    
   d_q                                   [-]:       1.25    
   p0                                  [kPa]:       132.    
   gamma_sub_av*D_gamma*N_gamma*pi*B_eq*B_eq*B_eq/8
                                        [kN]:      0.378E+05
   p0*d_q*N_q*pi*B_eq*B_eq/4            [kN]:      0.321E+06
   Qv_GSS = gamma_sub_av*D_gamma*N_gamma*pi*B_eq*B_eq*B_eq/8 + p0*d_q*N_q*pi*B_eq*B_eq/4
 > Qv_GSS                               [kN]:   359111.6    
Qv_GSS_fun: END --- >

   phi_lay_deg1                        [deg]:       0.00    
   phi_lay_deg2                        [deg]:       17.5    
 > Q_ub                                 [kN]:   359111.6    
 > Qv_fmode_pt                              :GSS
Q_ub_sub: END --- >


Rs_pt_fun: < --- BEGIN
 < D                                     [m]:       17.3    
 < B_pt                                  [m]:       17.2    
 < T_pt                                  [m]:      0.600    
   D + T_pt                              [m]:       17.9    
   su_av                                kPa]:       729.    
   P_pt                                  [m]:       54.0    
   T_pt                                  [m]:      0.600    
   Rs_pt = ( 3/4 ) * su_av * P_pt * T_pt;    
 > Rs_pt                                [kN]:   17726.09    
Rs_pt_fun: END --- >


dW_pt_fun: < --- BEGIN
 < D                                     [m]:       17.3    
 < B_pt                                  [m]:       17.2    
 < T_pt                                  [m]:      0.600    
   D + T_pt                              [m]:       17.9    
   A_pt                                [m^2]:       232.    
   gamma_sub_av                     [kN/m^3]:       7.00    
   dW_pt = A_pt * T_pt * gamma_sub_av;       
 > dW_pt = A_pt * T_pt * gamma_sub_av   [kN]:   976.0607    
dW_pt_fun: END --- >

   D        D_pt      B_pt      T_pt        Q_ub          Rs_pt          dW_pt     lay_D_pt  sum_W_pt       sum_Rs_pt     Qv_D         .   Qv_fmode
  [m]        [m]       [m]      [m]         [kN]          [kN]           [kN]           [-]  [kN]           [kN]          [kN]
  4.80      4.80      15.3      0.00       66154.04       0.000000       0.000000        4   0.000000       0.000000       66154.04    *  1 GSS
  4.80      6.30      16.1      1.50       75561.31       0.000000       2124.563        5   2124.563       0.000000       73436.75       2 PTS>GSC
  6.30      13.0      16.1      6.70       624740.1       177110.6       9829.926        6   11954.49       177110.6       789896.3       6 PTS>PTC>GSS
  13.0      14.4      15.3      1.40       132460.8       0.000000       2509.144        7   13372.50       177110.6       296198.9       3 PTS>PTC>PTS>SQC
  14.4      16.4      16.8      2.00       782035.7       56988.92       3085.376        8   19112.63       234099.5       997022.6       7 PTS>PTC>PTS>PTC>GSS
  16.4      17.3      15.3     0.900       235313.7       0.000000       1406.320        9   17353.13       234099.5       452060.1       4 PTS>PTC>PTS>PTC>PTS>SQC
  17.3      17.9      17.2     0.600       359111.6       17726.09       976.0607       10   22906.73       251825.6       588030.5       5 PTS>PTC>PTS>PTC>PTS>PTC>GSS
 > Qv_Dpt                               [kN]:   66154.04    
Qv_Dpt_sub: END --- >


ioutd_z_tip_verify: < --- BEGIN
 < z_tip_verify [m]:  7.90    

Qv, GROSS resistances at foundation depth D ------:
--------------------------------------------------:

Qv_Dpt_sub: < --- BEGIN
 < D                                     [m]:       4.80    
 < B_eq                                  [m]:       15.3    

Q_ub_sub: < --- BEGIN     
 < D                                     [m]:       4.80    
 < D_pt                                  [m]:       4.80    
 < B_pt                                  [m]:       15.3    

Qv_GSS_fun: < --- BEGIN
 < D                                     [m]:       4.80    
 < B_eq                                  [m]:       15.3    

phi_av_fun: < --- BEGIN
 < z1                                    [m]:       4.80    
 < z2                                    [m]:       8.62    
   is_D_clay                             [-]:         F
   lay_num_z1                            [-]:         4
   lay_num_z2                            [-]:         5
   tan_phi_av_sum                    [deg-m]:      0.402    
   T_sum                                 [m]:       1.50    
 > phi_av                              [rad]:      0.262    
 > phi_av/deg                          [deg]:       15.0    
phi_av_fun : END --- >

   is_D_clay                           [t/f]:         F
   gamma_sub_av                     [kN/m^3]:       7.15    
   phi_rad_av                          [rad]:      0.262    
   phi_rad_av/deg                      [deg]:       15.0    
   N_gamma                               [-]:       1.80    
   N_q                                   [-]:       7.20    
   d_gamma                               [-]:       1.00    
   d_q                                   [-]:       1.09    
   p0                                  [kPa]:       33.3    
   gamma_sub_av*D_gamma*N_gamma*pi*B_eq*B_eq*B_eq/8
                                        [kN]:      0.181E+05
   p0*d_q*N_q*pi*B_eq*B_eq/4            [kN]:      0.480E+05
   Qv_GSS = gamma_sub_av*D_gamma*N_gamma*pi*B_eq*B_eq*B_eq/8 + p0*d_q*N_q*pi*B_eq*B_eq/4
 > Qv_GSS                               [kN]:   66154.04    
Qv_GSS_fun: END --- >

   phi_lay_deg1                        [deg]:       15.0    
   phi_lay_deg2                        [deg]:       15.0    
 > Q_ub                                 [kN]:   66154.04    
 > Qv_fmode_pt                              :GSS
Q_ub_sub: END --- >


Rs_pt_fun: < --- BEGIN
 < D                                     [m]:       4.80    
 < B_pt                                  [m]:       15.3    
 < T_pt                                  [m]:       0.00    
   D + T_pt                              [m]:       4.80    
   su_av                                kPa]:       0.00    
   P_pt                                  [m]:       0.00    
   T_pt                                  [m]:       0.00    
   Rs_pt = ( 3/4 ) * su_av * P_pt * T_pt;    
 > Rs_pt                                [kN]:   0.000000    
Rs_pt_fun: END --- >


dW_pt_fun: < --- BEGIN
 < D                                     [m]:       4.80    
 < B_pt                                  [m]:       15.3    
 < T_pt                                  [m]:       0.00    
   D + T_pt                              [m]:       4.80    
   A_pt                                [m^2]:       184.    
   gamma_sub_av                     [kN/m^3]:       7.00    
   dW_pt = A_pt * T_pt * gamma_sub_av;       
 > dW_pt = A_pt * T_pt * gamma_sub_av   [kN]:   0.000000    
dW_pt_fun: END --- >


Q_ub_sub: < --- BEGIN     
 < D                                     [m]:       4.80    
 < D_pt                                  [m]:       6.30    
 < B_pt                                  [m]:       16.1    

Qv_SQC_fun: < --- BEGIN
 < D                                     [m]:       6.30    
 < B_eq                                  [m]:       16.1    

ZTOP_sq_fun: < --- BEGIN
 < D                                     [m]:       6.30    
 < B_eq                                  [m]:
   z1                                [m bsf]:       6.30    
   z0                                [m bsf]:       0.00    
   lay_num_z1                            [-]:         5
   lay_num_z0                            [-]:         1
   r_sq                                  [-]:       1.50    
   su0_bot                             [kPa]:      0.100E-02
   su1_top                             [kPa]:      0.100E-02
   su1_top/su0_bot                       [-]:       1.00    
 > ZTOP_sq_fun                       [m bsf]:       6.30    
ZTOP_sq_fun: END --- >


T_sq_fun: < --- BEGIN
 < D                                     [m]:       6.30    
 < B_eq                                  [m]:       16.1    
   z1                                    [m]:       6.30    
   z2                                    [m]:       14.3    
   lay_num_z1                            [-]:         5
   lay_num_z2                            [-]:         6
   r_sq                                  [-]:       1.50    
   su1_bot                             [kPa]:       48.0    
   su2_top                             [kPa]:       48.0    
   su2_top/su1_bot                       [-]:       1.00    
 > T_sq                                  [m]:       6.70    
T_sq_fun: END --- >

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:       6.70    
   b_sq                                  [-]:      0.330    
   B_eq                                  [m]:       16.1    
   T_sq/b_sq                             [m]:       20.3    
 * B_eq > T_sq/b_sq                    [t/f]:         F
   ztop_sq                           [m bsf]:       6.30    
   B_eq                                  [m]:       16.1    
   ztop_sq/B_eq (i.e. 'D/B')             [-]:      0.393    
 * 'D/B' <= 2.5                        [t/f]:         T
   Dtip (ie D corresponding to tip depth)[m]:       4.80    
   Dtip/B                                [-]:      0.299    
 * is_SQC_occurring at Dtip/B=2.5      [t/f]:         F
   [B_eq > T_sq/b_sq] & ['D/B' <= 2.5]       
 * is_SQC_occurring                    [t/f]:         F

 > Qv_SQC                               [kN]:   1000000.    
Qv_SQC_fun: END --- >


Qv_GSC_fun: < --- BEGIN
 < D                                 [m bsf]:       6.30    
 < B_eq                                  [m]:       16.1    
 < sp_beta                             [rad]:       3.14    
 < sp_beta/deg                         [deg]:       180.    

rho_av_fun: < --- BEGIN     
 < D                                     [m]:       6.30    
 < B                                     [m]:       16.1    
   is_D_clay                           [t/f]:         T
   z1                                    [m]:       6.30    
   z2                                    [m]:       14.3    
   lay_num_z1                            [-]:         5
   lay_num_z2                            [-]:         6
   rho_av_sum                          [kPa]:       0.00    
   T_sum                                 [m]:       6.70    
 > rho_av                            [kPa/m]:       0.00    
rho_av_fun: END --- >

   is_D_clay                           [t/f]:         T
   rho_av                            [kPa/m]:       0.00    
   su0                                 [kPa]:       48.0    
   D                                     [m]:       6.30    
   z_top_first_CLAY_layer_aboveD     [m bsf]:       6.30    
   sum = su0 - rho_av*D                [kPa]:       48.0    
   sum_modified = max(sum,0)           [kPa]:       48.0    
   kappa_m = rho_av*B_eq/sum_modified    [-]:       0.00    
   kappa_m_modified = max(kappa_m,5)     [-]:       0.00    
   B_eq                                  [m]:       16.1    
   lambda = D/B_eq                       [-]:      0.393    
   kappa_0=kappa_mod/(lambda*kappa_mod+1)[-]:       0.00    
   sp_beta/deg                         [deg]:       180.    
   sp_alpha                              [-]:      0.600    
   Nc0                                   [-]:       6.87    
   su0 * Nc0                           [kPa]:       330.    
   p0                                  [kPa]:       43.8    
   su0*Nc0 + p0                        [kPa]:       373.    
   pi * B_eq * B_eq /4                 [m^2]:       202.    
   Qv_GSC = ( su0*Nc0 + p0 ) * pi * B_eq * B_eq / 4 
 > Qv_GSC                               [kN]:   75561.31    
Qv_GSC_fun: END --- >

   phi_lay_deg1                        [deg]:       15.0    
   phi_lay_deg2                        [deg]:       0.00    
 > Q_ub                                 [kN]:   75561.31    
 > Qv_fmode_pt                              :GSC
Q_ub_sub: END --- >


Rs_pt_fun: < --- BEGIN
 < D                                     [m]:       4.80    
 < B_pt                                  [m]:       16.1    
 < T_pt                                  [m]:       1.50    
   D + T_pt                              [m]:       6.30    
   su_av                                kPa]:       0.00    
   P_pt                                  [m]:       0.00    
   T_pt                                  [m]:       1.50    
   Rs_pt = ( 3/4 ) * su_av * P_pt * T_pt;    
 > Rs_pt                                [kN]:   0.000000    
Rs_pt_fun: END --- >


dW_pt_fun: < --- BEGIN
 < D                                     [m]:       4.80    
 < B_pt                                  [m]:       16.1    
 < T_pt                                  [m]:       1.50    
   D + T_pt                              [m]:       6.30    
   A_pt                                [m^2]:       202.    
   gamma_sub_av                     [kN/m^3]:       7.00    
   dW_pt = A_pt * T_pt * gamma_sub_av;       
 > dW_pt = A_pt * T_pt * gamma_sub_av   [kN]:   2124.563    
dW_pt_fun: END --- >


Q_ub_sub: < --- BEGIN     
 < D                                     [m]:       6.30    
 < D_pt                                  [m]:       13.0    
 < B_pt                                  [m]:       16.1    

Qv_GSS_fun: < --- BEGIN
 < D                                     [m]:       13.0    
 < B_eq                                  [m]:       16.1    

phi_av_fun: < --- BEGIN
 < z1                                    [m]:       13.0    
 < z2                                    [m]:       17.0    
   is_D_clay                             [-]:         F
   lay_num_z1                            [-]:         6
   lay_num_z2                            [-]:         8
   tan_phi_av_sum                    [deg-m]:       1.01    
   T_sum                                 [m]:       2.01    
 > phi_av                              [rad]:      0.464    
 > phi_av/deg                          [deg]:       26.6    
phi_av_fun : END --- >

   is_D_clay                           [t/f]:         F
   gamma_sub_av                     [kN/m^3]:       8.19    
   phi_rad_av                          [rad]:      0.464    
   phi_rad_av/deg                      [deg]:       26.6    
   N_gamma                               [-]:       8.15    
   N_q                                   [-]:       22.9    
   d_gamma                               [-]:       1.00    
   d_q                                   [-]:       1.21    
   p0                                  [kPa]:       92.4    
   gamma_sub_av*D_gamma*N_gamma*pi*B_eq*B_eq*B_eq/8
                                        [kN]:      0.108E+06
   p0*d_q*N_q*pi*B_eq*B_eq/4            [kN]:      0.516E+06
   Qv_GSS = gamma_sub_av*D_gamma*N_gamma*pi*B_eq*B_eq*B_eq/8 + p0*d_q*N_q*pi*B_eq*B_eq/4
 > Qv_GSS                               [kN]:   624740.1    
Qv_GSS_fun: END --- >

   phi_lay_deg1                        [deg]:       0.00    
   phi_lay_deg2                        [deg]:       25.0    
 > Q_ub                                 [kN]:   624740.1    
 > Qv_fmode_pt                              :GSS
Q_ub_sub: END --- >


Rs_pt_fun: < --- BEGIN
 < D                                     [m]:       6.30    
 < B_pt                                  [m]:       16.1    
 < T_pt                                  [m]:       6.70    
   D + T_pt                              [m]:       13.0    
   su_av                                kPa]:       699.    
   P_pt                                  [m]:       50.4    
   T_pt                                  [m]:       6.70    
   Rs_pt = ( 3/4 ) * su_av * P_pt * T_pt;    
 > Rs_pt                                [kN]:   177110.6    
Rs_pt_fun: END --- >


dW_pt_fun: < --- BEGIN
 < D                                     [m]:       6.30    
 < B_pt                                  [m]:       16.1    
 < T_pt                                  [m]:       6.70    
   D + T_pt                              [m]:       13.0    
   A_pt                                [m^2]:       202.    
   gamma_sub_av                     [kN/m^3]:       7.25    
   dW_pt = A_pt * T_pt * gamma_sub_av;       
 > dW_pt = A_pt * T_pt * gamma_sub_av   [kN]:   9829.926    
dW_pt_fun: END --- >


Q_ub_sub: < --- BEGIN     
 < D                                     [m]:       13.0    
 < D_pt                                  [m]:       14.4    
 < B_pt                                  [m]:       16.8    

Qv_SQC_fun: < --- BEGIN
 < D                                     [m]:       14.4    
 < B_eq                                  [m]:       16.8    

ZTOP_sq_fun: < --- BEGIN
 < D                                     [m]:       14.4    
 < B_eq                                  [m]:
   z1                                [m bsf]:       14.4    
   z0                                [m bsf]:       0.00    
   lay_num_z1                            [-]:         7
   lay_num_z0                            [-]:         1
   r_sq                                  [-]:       1.50    
   su0_bot                             [kPa]:      0.100E-02
   su1_top                             [kPa]:      0.100E-02
   su1_top/su0_bot                       [-]:       1.00    
 > ZTOP_sq_fun                       [m bsf]:       14.4    
ZTOP_sq_fun: END --- >


T_sq_fun: < --- BEGIN
 < D                                     [m]:       14.4    
 < B_eq                                  [m]:       16.8    
   z1                                    [m]:       14.4    
   z2                                    [m]:       18.4    
   lay_num_z1                            [-]:         7
   lay_num_z2                            [-]:        10
   r_sq                                  [-]:       1.50    
   su1_bot                             [kPa]:       71.0    
   su2_top                             [kPa]:       71.0    
   su2_top/su1_bot                       [-]:       1.00    
 > T_sq                                  [m]:       2.00    
T_sq_fun: END --- >

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:       2.00    
   b_sq                                  [-]:      0.330    
   B_eq                                  [m]:       16.8    
   T_sq/b_sq                             [m]:       6.06    
 * B_eq > T_sq/b_sq                    [t/f]:         T
   ztop_sq                           [m bsf]:       14.4    
   B_eq                                  [m]:       16.8    
   ztop_sq/B_eq (i.e. 'D/B')             [-]:      0.860    
 * 'D/B' <= 2.5                        [t/f]:         T
   Dtip (ie D corresponding to tip depth)[m]:       4.80    
   Dtip/B                                [-]:      0.287    
 * is_SQC_occurring at Dtip/B=2.5      [t/f]:         F
   [B_eq > T_sq/b_sq] & ['D/B' <= 2.5]       
 * is_SQC_occurring                    [t/f]:         T

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:       2.00    
   D + T_sq                              [m]:       16.4    
   a_sq                                  [-]:       5.00    
   b_sq                                  [-]:      0.330    
   a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq    [-]:       8.80    
   su_av                               [kPa]:       71.0    
   p'0                                 [kPa]:       106.    
   (a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su0 [kPa]:       624.    
   pi*B_eq*B_eq/4                      [m^2]:       220.    
   Qv_SQC = (( a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su_av + p'0 ) *pi*B_eq*B_eq/4
 > Qv_SQC                               [kN]:   160970.9    
Qv_SQC_fun: END --- >


Qv_GSC_fun: < --- BEGIN
 < D                                 [m bsf]:       14.4    
 < B_eq                                  [m]:       16.8    
 < sp_beta                             [rad]:       3.14    
 < sp_beta/deg                         [deg]:       180.    

rho_av_fun: < --- BEGIN     
 < D                                     [m]:       14.4    
 < B                                     [m]:       16.8    
   is_D_clay                           [t/f]:         T
   z1                                    [m]:       14.4    
   z2                                    [m]:       18.4    
   lay_num_z1                            [-]:         7
   lay_num_z2                            [-]:        10
   rho_av_sum                          [kPa]:       0.00    
   T_sum                                 [m]:       2.60    
 > rho_av                            [kPa/m]:       0.00    
rho_av_fun: END --- >

   is_D_clay                           [t/f]:         T
   rho_av                            [kPa/m]:       0.00    
   su0                                 [kPa]:       71.0    
   D                                     [m]:       14.4    
   z_top_first_CLAY_layer_aboveD     [m bsf]:       14.4    
   sum = su0 - rho_av*D                [kPa]:       71.0    
   sum_modified = max(sum,0)           [kPa]:       71.0    
   kappa_m = rho_av*B_eq/sum_modified    [-]:       0.00    
   kappa_m_modified = max(kappa_m,5)     [-]:       0.00    
   B_eq                                  [m]:       16.8    
   lambda = D/B_eq                       [-]:      0.860    
   kappa_0=kappa_mod/(lambda*kappa_mod+1)[-]:       0.00    
   sp_beta/deg                         [deg]:       180.    
   sp_alpha                              [-]:      0.600    
   Nc0                                   [-]:       7.52    
   su0 * Nc0                           [kPa]:       534.    
   p0                                  [kPa]:       106.    
   su0*Nc0 + p0                        [kPa]:       640.    
   pi * B_eq * B_eq /4                 [m^2]:       220.    
   Qv_GSC = ( su0*Nc0 + p0 ) * pi * B_eq * B_eq / 4 
 > Qv_GSC                               [kN]:   140951.0    
Qv_GSC_fun: END --- >

   phi_lay_deg1                        [deg]:       25.0    
   phi_lay_deg2                        [deg]:       0.00    
 > Q_ub                                 [kN]:   160970.9    
 > Qv_fmode_pt                              :SQC
Q_ub_sub: END --- >


Rs_pt_fun: < --- BEGIN
 < D                                     [m]:       13.0    
 < B_pt                                  [m]:       16.8    
 < T_pt                                  [m]:       1.40    
   D + T_pt                              [m]:       14.4    
   su_av                                kPa]:       0.00    
   P_pt                                  [m]:       0.00    
   T_pt                                  [m]:       1.40    
   Rs_pt = ( 3/4 ) * su_av * P_pt * T_pt;    
 > Rs_pt                                [kN]:   0.000000    
Rs_pt_fun: END --- >


dW_pt_fun: < --- BEGIN
 < D                                     [m]:       13.0    
 < B_pt                                  [m]:       16.8    
 < T_pt                                  [m]:       1.40    
   D + T_pt                              [m]:       14.4    
   A_pt                                [m^2]:       220.    
   gamma_sub_av                     [kN/m^3]:       9.75    
   dW_pt = A_pt * T_pt * gamma_sub_av;       
 > dW_pt = A_pt * T_pt * gamma_sub_av   [kN]:   3007.270    
dW_pt_fun: END --- >


Qv_SQC_fun: < --- BEGIN
 < D                                     [m]:       14.4    
 < B_eq                                  [m]:       16.8    

ZTOP_sq_fun: < --- BEGIN
 < D                                     [m]:       14.4    
 < B_eq                                  [m]:
   z1                                [m bsf]:       14.4    
   z0                                [m bsf]:       0.00    
   lay_num_z1                            [-]:         7
   lay_num_z0                            [-]:         1
   r_sq                                  [-]:       1.50    
   su0_bot                             [kPa]:      0.100E-02
   su1_top                             [kPa]:      0.100E-02
   su1_top/su0_bot                       [-]:       1.00    
 > ZTOP_sq_fun                       [m bsf]:       14.4    
ZTOP_sq_fun: END --- >


T_sq_fun: < --- BEGIN
 < D                                     [m]:       14.4    
 < B_eq                                  [m]:       16.8    
   z1                                    [m]:       14.4    
   z2                                    [m]:       18.4    
   lay_num_z1                            [-]:         7
   lay_num_z2                            [-]:        10
   r_sq                                  [-]:       1.50    
   su1_bot                             [kPa]:       71.0    
   su2_top                             [kPa]:       71.0    
   su2_top/su1_bot                       [-]:       1.00    
 > T_sq                                  [m]:       2.00    
T_sq_fun: END --- >

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:       2.00    
   b_sq                                  [-]:      0.330    
   B_eq                                  [m]:       16.8    
   T_sq/b_sq                             [m]:       6.06    
 * B_eq > T_sq/b_sq                    [t/f]:         T
   ztop_sq                           [m bsf]:       14.4    
   B_eq                                  [m]:       16.8    
   ztop_sq/B_eq (i.e. 'D/B')             [-]:      0.860    
 * 'D/B' <= 2.5                        [t/f]:         T
   Dtip (ie D corresponding to tip depth)[m]:       4.80    
   Dtip/B                                [-]:      0.287    
 * is_SQC_occurring at Dtip/B=2.5      [t/f]:         F
   [B_eq > T_sq/b_sq] & ['D/B' <= 2.5]       
 * is_SQC_occurring                    [t/f]:         T

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:       2.00    
   D + T_sq                              [m]:       16.4    
   a_sq                                  [-]:       5.00    
   b_sq                                  [-]:      0.330    
   a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq    [-]:       8.80    
   su_av                               [kPa]:       71.0    
   p'0                                 [kPa]:       106.    
   (a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su0 [kPa]:       624.    
   pi*B_eq*B_eq/4                      [m^2]:       220.    
   Qv_SQC = (( a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su_av + p'0 ) *pi*B_eq*B_eq/4
 > Qv_SQC                               [kN]:   160970.9    
Qv_SQC_fun: END --- >


Qv_SQC_fun: < --- BEGIN
 < D                                     [m]:       14.4    
 < B_eq                                  [m]:       15.3    

ZTOP_sq_fun: < --- BEGIN
 < D                                     [m]:       14.4    
 < B_eq                                  [m]:
   z1                                [m bsf]:       14.4    
   z0                                [m bsf]:       0.00    
   lay_num_z1                            [-]:         7
   lay_num_z0                            [-]:         1
   r_sq                                  [-]:       1.50    
   su0_bot                             [kPa]:      0.100E-02
   su1_top                             [kPa]:      0.100E-02
   su1_top/su0_bot                       [-]:       1.00    
 > ZTOP_sq_fun                       [m bsf]:       14.4    
ZTOP_sq_fun: END --- >


T_sq_fun: < --- BEGIN
 < D                                     [m]:       14.4    
 < B_eq                                  [m]:       15.3    
   z1                                    [m]:       14.4    
   z2                                    [m]:       18.4    
   lay_num_z1                            [-]:         7
   lay_num_z2                            [-]:        10
   r_sq                                  [-]:       1.50    
   su1_bot                             [kPa]:       71.0    
   su2_top                             [kPa]:       71.0    
   su2_top/su1_bot                       [-]:       1.00    
 > T_sq                                  [m]:       2.00    
T_sq_fun: END --- >

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:       2.00    
   b_sq                                  [-]:      0.330    
   B_eq                                  [m]:       15.3    
   T_sq/b_sq                             [m]:       6.06    
 * B_eq > T_sq/b_sq                    [t/f]:         T
   ztop_sq                           [m bsf]:       14.4    
   B_eq                                  [m]:       15.3    
   ztop_sq/B_eq (i.e. 'D/B')             [-]:      0.941    
 * 'D/B' <= 2.5                        [t/f]:         T
   Dtip (ie D corresponding to tip depth)[m]:       4.80    
   Dtip/B                                [-]:      0.314    
 * is_SQC_occurring at Dtip/B=2.5      [t/f]:         F
   [B_eq > T_sq/b_sq] & ['D/B' <= 2.5]       
 * is_SQC_occurring                    [t/f]:         T

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:       2.00    
   D + T_sq                              [m]:       16.4    
   a_sq                                  [-]:       5.00    
   b_sq                                  [-]:      0.330    
   a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq    [-]:       8.65    
   su_av                               [kPa]:       71.0    
   p'0                                 [kPa]:       106.    
   (a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su0 [kPa]:       614.    
   pi*B_eq*B_eq/4                      [m^2]:       184.    
   Qv_SQC = (( a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su_av + p'0 ) *pi*B_eq*B_eq/4
 > Qv_SQC                               [kN]:   132460.8    
Qv_SQC_fun: END --- >


Rs_pt_fun: < --- BEGIN
 < D                                     [m]:       13.0    
 < B_pt                                  [m]:       15.3    
 < T_pt                                  [m]:       1.40    
   D + T_pt                              [m]:       14.4    
   su_av                                kPa]:       0.00    
   P_pt                                  [m]:       0.00    
   T_pt                                  [m]:       1.40    
   Rs_pt = ( 3/4 ) * su_av * P_pt * T_pt;    
 > Rs_pt                                [kN]:   0.000000    
Rs_pt_fun: END --- >


dW_pt_fun: < --- BEGIN
 < D                                     [m]:       13.0    
 < B_pt                                  [m]:       15.3    
 < T_pt                                  [m]:       1.40    
   D + T_pt                              [m]:       14.4    
   A_pt                                [m^2]:       184.    
   gamma_sub_av                     [kN/m^3]:       9.75    
   dW_pt = A_pt * T_pt * gamma_sub_av;       
 > dW_pt = A_pt * T_pt * gamma_sub_av   [kN]:   2509.144    
dW_pt_fun: END --- >


   D, D_pt                               [m]:       13.0           14.4    
   D_pt layer #                          [-]:    4
   Qv_Fmode_PT                           [-]:PTS>SQC
   decreasing dW for punch-through layer#[-]:    3
   using B_eq (was B_pt), B_eq , B_pt    [m]:       15.3           16.8    
   (B_eq/B_pt)^2                         [-]:      0.834    
   old dW,     new dW                   [kN]:      0.301E+04      0.251E+04
   old Qv_SQC, new Qv_SQC (ie Q_ub)     [kN]:      0.161E+06      0.132E+06

Q_ub_sub: < --- BEGIN     
 < D                                     [m]:       14.4    
 < D_pt                                  [m]:       16.4    
 < B_pt                                  [m]:       16.8    

Qv_GSS_fun: < --- BEGIN
 < D                                     [m]:       16.4    
 < B_eq                                  [m]:       16.8    

phi_av_fun: < --- BEGIN
 < z1                                    [m]:       16.4    
 < z2                                    [m]:       18.4    
   is_D_clay                             [-]:         F
   lay_num_z1                            [-]:         8
   lay_num_z2                            [-]:        10
   tan_phi_av_sum                    [deg-m]:      0.677    
   T_sum                                 [m]:       1.40    
 > phi_av                              [rad]:      0.451    
 > phi_av/deg                          [deg]:       25.8    
phi_av_fun : END --- >

   is_D_clay                           [t/f]:         F
   gamma_sub_av                     [kN/m^3]:       8.18    
   phi_rad_av                          [rad]:      0.451    
   phi_rad_av/deg                      [deg]:       25.8    
   N_gamma                               [-]:       7.08    
   N_q                                   [-]:       20.6    
   d_gamma                               [-]:       1.00    
   d_q                                   [-]:       1.24    
   p0                                  [kPa]:       120.    
   gamma_sub_av*D_gamma*N_gamma*pi*B_eq*B_eq*B_eq/8
                                        [kN]:      0.107E+06
   p0*d_q*N_q*pi*B_eq*B_eq/4            [kN]:      0.675E+06
   Qv_GSS = gamma_sub_av*D_gamma*N_gamma*pi*B_eq*B_eq*B_eq/8 + p0*d_q*N_q*pi*B_eq*B_eq/4
 > Qv_GSS                               [kN]:   782035.7    
Qv_GSS_fun: END --- >

   phi_lay_deg1                        [deg]:       0.00    
   phi_lay_deg2                        [deg]:       30.0    
 > Q_ub                                 [kN]:   782035.7    
 > Qv_fmode_pt                              :GSS
Q_ub_sub: END --- >


Rs_pt_fun: < --- BEGIN
 < D                                     [m]:       14.4    
 < B_pt                                  [m]:       16.8    
 < T_pt                                  [m]:       2.00    
   D + T_pt                              [m]:       16.4    
   su_av                                kPa]:       722.    
   P_pt                                  [m]:       52.6    
   T_pt                                  [m]:       2.00    
   Rs_pt = ( 3/4 ) * su_av * P_pt * T_pt;    
 > Rs_pt                                [kN]:   56988.92    
Rs_pt_fun: END --- >


dW_pt_fun: < --- BEGIN
 < D                                     [m]:       14.4    
 < B_pt                                  [m]:       16.8    
 < T_pt                                  [m]:       2.00    
   D + T_pt                              [m]:       16.4    
   A_pt                                [m^2]:       220.    
   gamma_sub_av                     [kN/m^3]:       7.00    
   dW_pt = A_pt * T_pt * gamma_sub_av;       
 > dW_pt = A_pt * T_pt * gamma_sub_av   [kN]:   3085.376    
dW_pt_fun: END --- >


Q_ub_sub: < --- BEGIN     
 < D                                     [m]:       16.4    
 < D_pt                                  [m]:       17.3    
 < B_pt                                  [m]:       17.2    

Qv_SQC_fun: < --- BEGIN
 < D                                     [m]:       17.3    
 < B_eq                                  [m]:       17.2    

ZTOP_sq_fun: < --- BEGIN
 < D                                     [m]:       17.3    
 < B_eq                                  [m]:
   z1                                [m bsf]:       17.3    
   z0                                [m bsf]:       0.00    
   lay_num_z1                            [-]:         9
   lay_num_z0                            [-]:         1
   r_sq                                  [-]:       1.50    
   su0_bot                             [kPa]:      0.100E-02
   su1_top                             [kPa]:      0.100E-02
   su1_top/su0_bot                       [-]:       1.00    
 > ZTOP_sq_fun                       [m bsf]:       17.3    
ZTOP_sq_fun: END --- >


T_sq_fun: < --- BEGIN
 < D                                     [m]:       17.3    
 < B_eq                                  [m]:       17.2    
   z1                                    [m]:       17.3    
   z2                                    [m]:       18.4    
   lay_num_z1                            [-]:         9
   lay_num_z2                            [-]:        10
   r_sq                                  [-]:       1.50    
   su1_bot                             [kPa]:       78.0    
   su2_top                             [kPa]:       78.0    
   su2_top/su1_bot                       [-]:       1.00    
 > T_sq                                  [m]:      0.600    
T_sq_fun: END --- >

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:      0.600    
   b_sq                                  [-]:      0.330    
   B_eq                                  [m]:       17.2    
   T_sq/b_sq                             [m]:       1.82    
 * B_eq > T_sq/b_sq                    [t/f]:         T
   ztop_sq                           [m bsf]:       17.3    
   B_eq                                  [m]:       17.2    
   ztop_sq/B_eq (i.e. 'D/B')             [-]:       1.01    
 * 'D/B' <= 2.5                        [t/f]:         T
   Dtip (ie D corresponding to tip depth)[m]:       4.80    
   Dtip/B                                [-]:      0.279    
 * is_SQC_occurring at Dtip/B=2.5      [t/f]:         F
   [B_eq > T_sq/b_sq] & ['D/B' <= 2.5]       
 * is_SQC_occurring                    [t/f]:         T

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:      0.600    
   D + T_sq                              [m]:       17.9    
   a_sq                                  [-]:       5.00    
   b_sq                                  [-]:      0.330    
   a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq    [-]:       15.7    
   su_av                               [kPa]:       78.0    
   p'0                                 [kPa]:       128.    
   (a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su0 [kPa]:      0.122E+04
   pi*B_eq*B_eq/4                      [m^2]:       232.    
   Qv_SQC = (( a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su_av + p'0 ) *pi*B_eq*B_eq/4
 > Qv_SQC                               [kN]:   313609.1    
Qv_SQC_fun: END --- >


Qv_GSC_fun: < --- BEGIN
 < D                                 [m bsf]:       17.3    
 < B_eq                                  [m]:       17.2    
 < sp_beta                             [rad]:       3.14    
 < sp_beta/deg                         [deg]:       180.    

rho_av_fun: < --- BEGIN     
 < D                                     [m]:       17.3    
 < B                                     [m]:       17.2    
   is_D_clay                           [t/f]:         T
   z1                                    [m]:       17.3    
   z2                                    [m]:       18.4    
   lay_num_z1                            [-]:         9
   lay_num_z2                            [-]:        10
   rho_av_sum                          [kPa]:       0.00    
   T_sum                                 [m]:      0.600    
 > rho_av                            [kPa/m]:       0.00    
rho_av_fun: END --- >

   is_D_clay                           [t/f]:         T
   rho_av                            [kPa/m]:       0.00    
   su0                                 [kPa]:       78.0    
   D                                     [m]:       17.3    
   z_top_first_CLAY_layer_aboveD     [m bsf]:       17.3    
   sum = su0 - rho_av*D                [kPa]:       78.0    
   sum_modified = max(sum,0)           [kPa]:       78.0    
   kappa_m = rho_av*B_eq/sum_modified    [-]:       0.00    
   kappa_m_modified = max(kappa_m,5)     [-]:       0.00    
   B_eq                                  [m]:       17.2    
   lambda = D/B_eq                       [-]:       1.01    
   kappa_0=kappa_mod/(lambda*kappa_mod+1)[-]:       0.00    
   sp_beta/deg                         [deg]:       180.    
   sp_alpha                              [-]:      0.600    
   Nc0                                   [-]:       7.70    
   su0 * Nc0                           [kPa]:       600.    
   p0                                  [kPa]:       128.    
   su0*Nc0 + p0                        [kPa]:       728.    
   pi * B_eq * B_eq /4                 [m^2]:       232.    
   Qv_GSC = ( su0*Nc0 + p0 ) * pi * B_eq * B_eq / 4 
 > Qv_GSC                               [kN]:   169178.6    
Qv_GSC_fun: END --- >

   phi_lay_deg1                        [deg]:       30.0    
   phi_lay_deg2                        [deg]:       0.00    
 > Q_ub                                 [kN]:   313609.1    
 > Qv_fmode_pt                              :SQC
Q_ub_sub: END --- >


Rs_pt_fun: < --- BEGIN
 < D                                     [m]:       16.4    
 < B_pt                                  [m]:       17.2    
 < T_pt                                  [m]:      0.900    
   D + T_pt                              [m]:       17.3    
   su_av                                kPa]:       0.00    
   P_pt                                  [m]:       0.00    
   T_pt                                  [m]:      0.900    
   Rs_pt = ( 3/4 ) * su_av * P_pt * T_pt;    
 > Rs_pt                                [kN]:   0.000000    
Rs_pt_fun: END --- >


dW_pt_fun: < --- BEGIN
 < D                                     [m]:       16.4    
 < B_pt                                  [m]:       17.2    
 < T_pt                                  [m]:      0.900    
   D + T_pt                              [m]:       17.3    
   A_pt                                [m^2]:       232.    
   gamma_sub_av                     [kN/m^3]:       8.50    
   dW_pt = A_pt * T_pt * gamma_sub_av;       
 > dW_pt = A_pt * T_pt * gamma_sub_av   [kN]:   1777.290    
dW_pt_fun: END --- >


Qv_SQC_fun: < --- BEGIN
 < D                                     [m]:       17.3    
 < B_eq                                  [m]:       17.2    

ZTOP_sq_fun: < --- BEGIN
 < D                                     [m]:       17.3    
 < B_eq                                  [m]:
   z1                                [m bsf]:       17.3    
   z0                                [m bsf]:       0.00    
   lay_num_z1                            [-]:         9
   lay_num_z0                            [-]:         1
   r_sq                                  [-]:       1.50    
   su0_bot                             [kPa]:      0.100E-02
   su1_top                             [kPa]:      0.100E-02
   su1_top/su0_bot                       [-]:       1.00    
 > ZTOP_sq_fun                       [m bsf]:       17.3    
ZTOP_sq_fun: END --- >


T_sq_fun: < --- BEGIN
 < D                                     [m]:       17.3    
 < B_eq                                  [m]:       17.2    
   z1                                    [m]:       17.3    
   z2                                    [m]:       18.4    
   lay_num_z1                            [-]:         9
   lay_num_z2                            [-]:        10
   r_sq                                  [-]:       1.50    
   su1_bot                             [kPa]:       78.0    
   su2_top                             [kPa]:       78.0    
   su2_top/su1_bot                       [-]:       1.00    
 > T_sq                                  [m]:      0.600    
T_sq_fun: END --- >

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:      0.600    
   b_sq                                  [-]:      0.330    
   B_eq                                  [m]:       17.2    
   T_sq/b_sq                             [m]:       1.82    
 * B_eq > T_sq/b_sq                    [t/f]:         T
   ztop_sq                           [m bsf]:       17.3    
   B_eq                                  [m]:       17.2    
   ztop_sq/B_eq (i.e. 'D/B')             [-]:       1.01    
 * 'D/B' <= 2.5                        [t/f]:         T
   Dtip (ie D corresponding to tip depth)[m]:       4.80    
   Dtip/B                                [-]:      0.279    
 * is_SQC_occurring at Dtip/B=2.5      [t/f]:         F
   [B_eq > T_sq/b_sq] & ['D/B' <= 2.5]       
 * is_SQC_occurring                    [t/f]:         T

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:      0.600    
   D + T_sq                              [m]:       17.9    
   a_sq                                  [-]:       5.00    
   b_sq                                  [-]:      0.330    
   a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq    [-]:       15.7    
   su_av                               [kPa]:       78.0    
   p'0                                 [kPa]:       128.    
   (a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su0 [kPa]:      0.122E+04
   pi*B_eq*B_eq/4                      [m^2]:       232.    
   Qv_SQC = (( a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su_av + p'0 ) *pi*B_eq*B_eq/4
 > Qv_SQC                               [kN]:   313609.1    
Qv_SQC_fun: END --- >


Qv_SQC_fun: < --- BEGIN
 < D                                     [m]:       17.3    
 < B_eq                                  [m]:       15.3    

ZTOP_sq_fun: < --- BEGIN
 < D                                     [m]:       17.3    
 < B_eq                                  [m]:
   z1                                [m bsf]:       17.3    
   z0                                [m bsf]:       0.00    
   lay_num_z1                            [-]:         9
   lay_num_z0                            [-]:         1
   r_sq                                  [-]:       1.50    
   su0_bot                             [kPa]:      0.100E-02
   su1_top                             [kPa]:      0.100E-02
   su1_top/su0_bot                       [-]:       1.00    
 > ZTOP_sq_fun                       [m bsf]:       17.3    
ZTOP_sq_fun: END --- >


T_sq_fun: < --- BEGIN
 < D                                     [m]:       17.3    
 < B_eq                                  [m]:       15.3    
   z1                                    [m]:       17.3    
   z2                                    [m]:       18.4    
   lay_num_z1                            [-]:         9
   lay_num_z2                            [-]:        10
   r_sq                                  [-]:       1.50    
   su1_bot                             [kPa]:       78.0    
   su2_top                             [kPa]:       78.0    
   su2_top/su1_bot                       [-]:       1.00    
 > T_sq                                  [m]:      0.600    
T_sq_fun: END --- >

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:      0.600    
   b_sq                                  [-]:      0.330    
   B_eq                                  [m]:       15.3    
   T_sq/b_sq                             [m]:       1.82    
 * B_eq > T_sq/b_sq                    [t/f]:         T
   ztop_sq                           [m bsf]:       17.3    
   B_eq                                  [m]:       15.3    
   ztop_sq/B_eq (i.e. 'D/B')             [-]:       1.13    
 * 'D/B' <= 2.5                        [t/f]:         T
   Dtip (ie D corresponding to tip depth)[m]:       4.80    
   Dtip/B                                [-]:      0.314    
 * is_SQC_occurring at Dtip/B=2.5      [t/f]:         F
   [B_eq > T_sq/b_sq] & ['D/B' <= 2.5]       
 * is_SQC_occurring                    [t/f]:         T

   is_D_clay                           [t/f]:         T
   T_sq                                  [m]:      0.600    
   D + T_sq                              [m]:       17.9    
   a_sq                                  [-]:       5.00    
   b_sq                                  [-]:      0.330    
   a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq    [-]:       14.8    
   su_av                               [kPa]:       78.0    
   p'0                                 [kPa]:       128.    
   (a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su0 [kPa]:      0.115E+04
   pi*B_eq*B_eq/4                      [m^2]:       184.    
   Qv_SQC = (( a_sq + b_sq*B_eq/T_sq + 1.2*D/B_eq)*su_av + p'0 ) *pi*B_eq*B_eq/4
 > Qv_SQC                               [kN]:   235313.7    
Qv_SQC_fun: END --- >


Rs_pt_fun: < --- BEGIN
 < D                                     [m]:       16.4    
 < B_pt                                  [m]:       15.3    
 < T_pt                                  [m]:      0.900    
   D + T_pt                              [m]:       17.3    
   su_av                                kPa]:       0.00    
   P_pt                                  [m]:       0.00    
   T_pt                                  [m]:      0.900    
   Rs_pt = ( 3/4 ) * su_av * P_pt * T_pt;    
 > Rs_pt                                [kN]:   0.000000    
Rs_pt_fun: END --- >


dW_pt_fun: < --- BEGIN
 < D                                     [m]:       16.4    
 < B_pt                                  [m]:       15.3    
 < T_pt                                  [m]:      0.900    
   D + T_pt                              [m]:       17.3    
   A_pt                                [m^2]:       184.    
   gamma_sub_av                     [kN/m^3]:       8.50    
   dW_pt = A_pt * T_pt * gamma_sub_av;       
 > dW_pt = A_pt * T_pt * gamma_sub_av   [kN]:   1406.320    
dW_pt_fun: END --- >


   D, D_pt                               [m]:       16.4           17.3    
   D_pt layer #                          [-]:    6
   Qv_Fmode_PT                           [-]:PTS>SQC
   decreasing dW for punch-through layer#[-]:    5
   using B_eq (was B_pt), B_eq , B_pt    [m]:       15.3           17.2    
   (B_eq/B_pt)^2                         [-]:      0.791    
   old dW,     new dW                   [kN]:      0.178E+04      0.141E+04
   old Qv_SQC, new Qv_SQC (ie Q_ub)     [kN]:      0.314E+06      0.235E+06

Q_ub_sub: < --- BEGIN     
 < D                                     [m]:       17.3    
 < D_pt                                  [m]:       17.9    
 < B_pt                                  [m]:       17.2    

Qv_GSS_fun: < --- BEGIN
 < D                                     [m]:       17.9    
 < B_eq                                  [m]:       17.2    

phi_av_fun: < --- BEGIN
 < z1                                    [m]:       17.9    
 < z2                                    [m]:       18.4    
   is_D_clay                             [-]:         F
   lay_num_z1                            [-]:        10
   lay_num_z2                            [-]:        10
   tan_phi_av_sum                    [deg-m]:      0.158    
   T_sum                                 [m]:      0.500    
 > phi_av                              [rad]:      0.305    
 > phi_av/deg                          [deg]:       17.5    
phi_av_fun : END --- >

   is_D_clay                           [t/f]:         F
   gamma_sub_av                     [kN/m^3]:       9.00    
   phi_rad_av                          [rad]:      0.305    
   phi_rad_av/deg                      [deg]:       17.5    
   N_gamma                               [-]:       2.10    
   N_q                                   [-]:       8.40    
   d_gamma                               [-]:       1.00    
   d_q                                   [-]:       1.25    
   p0                                  [kPa]:       132.    
   gamma_sub_av*D_gamma*N_gamma*pi*B_eq*B_eq*B_eq/8
                                        [kN]:      0.378E+05
   p0*d_q*N_q*pi*B_eq*B_eq/4            [kN]:      0.321E+06
   Qv_GSS = gamma_sub_av*D_gamma*N_gamma*pi*B_eq*B_eq*B_eq/8 + p0*d_q*N_q*pi*B_eq*B_eq/4
 > Qv_GSS                               [kN]:   359111.6    
Qv_GSS_fun: END --- >

   phi_lay_deg1                        [deg]:       0.00    
   phi_lay_deg2                        [deg]:       17.5    
 > Q_ub                                 [kN]:   359111.6    
 > Qv_fmode_pt                              :GSS
Q_ub_sub: END --- >


Rs_pt_fun: < --- BEGIN
 < D                                     [m]:       17.3    
 < B_pt                                  [m]:       17.2    
 < T_pt                                  [m]:      0.600    
   D + T_pt                              [m]:       17.9    
   su_av                                kPa]:       729.    
   P_pt                                  [m]:       54.0    
   T_pt                                  [m]:      0.600    
   Rs_pt = ( 3/4 ) * su_av * P_pt * T_pt;    
 > Rs_pt                                [kN]:   17726.09    
Rs_pt_fun: END --- >


dW_pt_fun: < --- BEGIN
 < D                                     [m]:       17.3    
 < B_pt                                  [m]:       17.2    
 < T_pt                                  [m]:      0.600    
   D + T_pt                              [m]:       17.9    
   A_pt                                [m^2]:       232.    
   gamma_sub_av                     [kN/m^3]:       7.00    
   dW_pt = A_pt * T_pt * gamma_sub_av;       
 > dW_pt = A_pt * T_pt * gamma_sub_av   [kN]:   976.0607    
dW_pt_fun: END --- >

 < z_tip_verify                          [m]:       7.90    
   D                                     [m]:       4.80    
   B_eq                                  [m]:       15.3    
   lay_num (D)                           [-]:         4
   is_D_clay                           [t/f]:         F

   Qv_GSS, ISO 19905-1 Equation A.9.3-8,
      GSS general shear resistance in SAND [kN]:   66154.04    *
   Qv_GSC, ISO 19905-1 Equation A.9.3-7,
      GSC general shear resistance in CLAY [kN]:   1000000.     
   Qv_SQC, ISO 19905-1 Equation A.9.3-10,
      SQC squeezing     resistance in CLAY [kN]:   1000000.     
   Qv_D_clay, GSC or SQC           in CLAY [kN]:   1000000.    
   Qub_D, GROSS resistance (GSS, GSC/SQC)  [kN]:   66154.04    
                       Qv_fmode, failure mode  :   'GSS'


Qv_GSS_fun: < --- BEGIN
 < D                                     [m]:       4.80    
 < B_eq                                  [m]:       15.3    

phi_av_fun: < --- BEGIN
 < z1                                    [m]:       4.80    
 < z2                                    [m]:       8.62    
   is_D_clay                             [-]:         F
   lay_num_z1                            [-]:         4
   lay_num_z2                            [-]:         5
   tan_phi_av_sum                    [deg-m]:      0.402    
   T_sum                                 [m]:       1.50    
 > phi_av                              [rad]:      0.262    
 > phi_av/deg                          [deg]:       15.0    
phi_av_fun : END --- >

   is_D_clay                           [t/f]:         F
   gamma_sub_av                     [kN/m^3]:       7.15    
   phi_rad_av                          [rad]:      0.262    
   phi_rad_av/deg                      [deg]:       15.0    
   N_gamma                               [-]:       1.80    
   N_q                                   [-]:       7.20    
   d_gamma                               [-]:       1.00    
   d_q                                   [-]:       1.09    
   p0                                  [kPa]:       33.3    
   gamma_sub_av*D_gamma*N_gamma*pi*B_eq*B_eq*B_eq/8
                                        [kN]:      0.181E+05
   p0*d_q*N_q*pi*B_eq*B_eq/4            [kN]:      0.480E+05
   Qv_GSS = gamma_sub_av*D_gamma*N_gamma*pi*B_eq*B_eq*B_eq/8 + p0*d_q*N_q*pi*B_eq*B_eq/4
 > Qv_GSS                               [kN]:   66154.04    
Qv_GSS_fun: END --- >


V_L, spudcan GROSS ultimate bearing resistance at foundation depth D ------:

D_B_eq_sub: < --- BEGIN
 < z_tip                                 [m]:       7.90    
   sp_H_max                              [m]:       3.10    
   sp_B_max                              [m]:       15.3    
 > D                                     [m]:       4.80    
 > B_eq                                  [m]:       15.3    
D_B_eq_sub: END --- >


H_cav_fun: < --- BEGIN
 < D                                     [m]:       4.80    
 < B_eq                                  [m]:       15.3    
 < I_BF                                  [-]:         1

 < D_top_first_SAND_layer                [m]:      0.500    
 < D_top_first_CLAY_layer                [m]:       0.00    
   all_CLAY_above_D                    [t/f]:         F
 > H_cav, ISO 19905-1 Eq. 9.3-3          [m]:       4.30    
H_cav_fun: END --- >


B_S_fun: < --- BEGIN
 < z_tip                                 [m]:       7.90    

D_B_eq_sub: < --- BEGIN
 < z_tip                                 [m]:       7.90    
   sp_H_max                              [m]:       3.10    
   sp_B_max                              [m]:       15.3    
 > D                                     [m]:       4.80    
 > B_eq                                  [m]:       15.3    
D_B_eq_sub: END --- >

   z1                                    [m]:       7.90    
   z2                                    [m]:       4.80    
   gamma_sub_av                     [kN/m^3]:       7.13    
   V_D                                  [kN]:       145.    
 > B_S = gamma_sub_av*V_D               [kN]:   1034.508    
B_S_fun: END --- >


W_BF_fun: < --- BEGIN
 < z_tip                                 [m]:       7.90    
 < I_BF                                  [-]:         1

D_B_eq_sub: < --- BEGIN
 < z_tip                                 [m]:       7.90    
   sp_H_max                              [m]:       3.10    
   sp_B_max                              [m]:       15.3    
 > D                                     [m]:       4.80    
 > B_eq                                  [m]:       15.3    
D_B_eq_sub: END --- >

   gamma_sub_av                     [kN/m^3]:       6.94    
   A                                   [m^2]:       184.    
   H_cav                                 [m]:       4.30    
   V_spud                              [m^3]:       227.    
   V_D                                 [m^3]:       145.    
   W_BF = gamma_sub_av [ A (D - H_cav) - (V_spud - V_D ) ]
 > W_BF                                 [kN]:   68.10682    
W_BF_fun: END --- >

 < z_tip_verify                          [m]:       7.90    
   D                                     [m]:       4.80    
   B_eq                                  [m]:       15.3    
   lay_num (D)                           [-]:         4
   is_D_clay                           [t/f]:         F

GROSS spudcan resistance at depth D --- >
-----------------------------------------
   H_CAV, ISO 19905-1 Eq.A.9.3-3,
     cavity depth                                        [m bsf]:       4.30    
   V_SPUD, ISO 19905-1 Fig.A.9.3-6, row1,
     spudcan total volume BENEATH BACKFILL                 [m^3]:       227.    
   V_D, ISO 19905-1 Fig.A.9.3-6, row2, 
     spudcan total volume BELOW MAXIMUM BEARING AREA       [m^3]:       145.    
   W_BF, ISO 19905-1 Eq.A.9.3-2,
     backfill weight due to backflow,
   W_BF = gamma [A (D - H_CAV) - (V_SPUD - V_D)]            [kN]:   68.10682    

   QV_DPT, spudcan NETT resistance at depth D               [kN]:   66154.04    
   B_S, soil buoyancy of spudcan BELOW MAXIMUM BEARING AREA [kN]:   1034.508    
        i.e. soil submerged weight displaced by spudcan below D [D is MAX AREA].
   W_BF, backfill weight due to backflow                    [kN]:   68.10682    
 > V_L, spudcan GROSS ultimate bearing resistance
        ie V_L = QV_DPT - W_BF + B_S                        [kN]:   67120.44    


ioutd_z_tip_verify: END --- >


soil_mat_adjust: < --- BEGIN
 < z_tip                                 [m]:  7.90    
 < i_sqc                                 [-]:  1
 < Qv_fmode                              [-]: GSS
 -------------------------------------------------------------------------------------------------------------
 output  file: "JOB.O01"
 load spread factor #=    1; load spread factor =  4.00     (#rows V_L results=  184)
 -------------------------------------------------------------------------------------------------------------
 z_tip        D      B_eq     H_cav    V_spud       V_D      B_S           W_BF         Qv_Dpt           V_L  
  [m]        [m]      [m]       [m]      [m^3]     [m^3]      [kN]         [kN]          [kN]            [kN] 
 -------------------------------------------------------------------------------------------------------------
  0.00      0.00     0.100E-02  0.00      0.00     0.100E-02 0.240E-02  0.2400000E-02  0.1000000E-02  0.1000000E-02
 0.100      0.00     0.246      0.00     0.158E-02 0.158E-02 0.380E-02  0.1000000E-02  0.2338996      0.2366989    
 0.200      0.00     0.492      0.00     0.127E-01 0.127E-01 0.304E-01  0.1000000E-02  0.9776733       1.007068    
 0.300      0.00     0.738      0.00     0.427E-01 0.427E-01 0.103      0.1000000E-02   2.294434       2.396015    
 0.400      0.00     0.984      0.00     0.101     0.101     0.243      0.1000000E-02   4.247293       4.489448    
 0.500      0.00      1.23      0.00     0.198     0.198     0.476      0.1000000E-02   6.899364       7.373952    
 0.600      0.00      1.48      0.00     0.342     0.342      1.12      0.1000000E-02   10.31376       11.43377    
 0.700      0.00      1.72      0.00     0.543     0.543      2.12      0.1000000E-02   14.98585       17.10274    
 0.800      0.00      1.97      0.00     0.811     0.811      3.54      0.1000000E-02   20.09105       23.63376    
 0.900      0.00      2.21      0.00      1.15      1.15      5.46      0.1000000E-02   26.08295       31.54570    
  1.00      0.00      2.46      0.00      1.58      1.58      7.96      0.1000000E-02   33.01009       40.96927    
  1.10      0.00      2.70      0.00      2.11      2.11      11.1      0.1000000E-02   41.02954       52.12766    
  1.20      0.00      2.95      0.00      2.74      2.74      14.9      0.1000000E-02   50.34324       65.27551    
  1.30      0.00      3.52      0.00      3.54      3.54      20.0      0.1000000E-02   74.96031       94.93951    
  1.40      0.00      4.18      0.00      4.71      4.71      27.2      0.1000000E-02   111.5955       138.8306    
  1.50      0.00      4.83      0.00      6.31      6.31      37.2      0.1000000E-02   157.5960       194.8368    
  1.60      0.00      5.49      0.00      8.40      8.40      50.5      0.1000000E-02   213.8760       264.3927    
  1.70      0.00      6.14      0.00      11.1      11.1      67.6      0.1000000E-02   281.3497       348.9006    
  1.80      0.00      6.79      0.00      14.3      14.3      88.8      0.1000000E-02   360.9312       449.7424    
  1.90      0.00      7.45      0.00      18.3      18.3      115.      0.1000000E-02   453.5346       568.5285    
  2.00      0.00      8.10      0.00      23.1      23.1      146.      0.1000000E-02   560.0740       706.4524    
  2.10      0.00      8.76      0.00      28.7      28.7      184.      0.1000000E-02   681.4636       865.0268    
  2.20      0.00      9.41      0.00      35.2      35.2      227.      0.1000000E-02   818.6175       1045.606    
  2.30      0.00      10.1      0.00      42.6      42.6      277.      0.1000000E-02   972.4499       1249.602    
  2.40      0.00      10.7      0.00      51.1      51.1      335.      0.1000000E-02   1143.875       1479.029    
  2.50      0.00      11.4      0.00      60.7      60.7      401.      0.1000000E-02   1333.807       1734.341    
  2.60      0.00      12.0      0.00      71.4      71.4      474.      0.1000000E-02   1543.159       2017.626    
  2.70      0.00      12.7      0.00      83.4      83.4      558.      0.1000000E-02   1772.847       2330.392    
  2.80      0.00      13.3      0.00      96.7      96.7      649.      0.1000000E-02   2023.783       2672.837    
  2.90      0.00      14.0      0.00      111.      111.      751.      0.1000000E-02   2296.883       3048.158    
  3.00      0.00      14.6      0.00      128.      128.      863.      0.1000000E-02   2593.060       3456.438    
  3.10      0.00      15.3      0.00      145.      145.      988.      0.1000000E-02   2913.229       3900.798    
  3.20     0.100      15.3      0.00      163.      145.     0.101E+04  0.9669725       3608.051       4618.461    
  3.30     0.200      15.3      0.00      180.      145.     0.104E+04   3.870591       4705.117       5737.422    
  3.40     0.300      15.3      0.00      197.      145.     0.106E+04   8.609188       6807.575       7858.948    
  3.50     0.400      15.3      0.00      212.      145.     0.109E+04   15.19531       12931.17       14001.75    
  3.60     0.500      15.3      0.00      227.      145.     0.111E+04   23.59727       27544.27       28630.69    
  3.70     0.600      15.3     0.100      227.      145.     0.111E+04   32.18821       28861.25       29936.13    
  3.80     0.700      15.3     0.200      227.      145.     0.110E+04   38.29583       30180.92       31246.50    
  3.90     0.800      15.3     0.300      227.      145.     0.110E+04   42.92688       31503.28       32561.28    
  4.00     0.900      15.3     0.400      227.      145.     0.110E+04   46.48406       32836.89       33888.26    
  4.10      1.00      15.3     0.500      227.      145.     0.109E+04   49.37008       34094.06       35139.60    
  4.20      1.10      15.3     0.600      227.      145.     0.109E+04   51.71916       35100.51       36140.50    
  4.30      1.20      15.3     0.700      227.      145.     0.109E+04   53.59843       36092.21       37127.38    
  4.40      1.30      15.3     0.800      227.      145.     0.109E+04   55.34346       37057.64       38087.99    
  4.50      1.40      15.3     0.900      227.      145.     0.108E+04   56.75292       38004.91       39030.90    
  4.60      1.50      15.3      1.00      227.      145.     0.108E+04   57.96102       38933.65       39955.25    
  4.70      1.60      15.3      1.10      227.      145.     0.108E+04   59.03488       39835.44       40853.02    
  4.80      1.70      15.3      1.20      227.      145.     0.107E+04   59.97452       40718.13       41731.70    
  4.90      1.80      15.3      1.30      227.      145.     0.107E+04   60.77992       41576.08       42585.89    
  5.00      1.90      15.3      1.40      227.      145.     0.107E+04   61.58532       42416.90       43422.72    
  5.10      2.00      15.3      1.50      227.      145.     0.106E+04   62.25648       43232.40       44234.59    
  5.20      2.10      15.3      1.60      227.      145.     0.106E+04   62.86053       44022.35       45020.87    
  5.30      2.20      15.3      1.70      227.      145.     0.106E+04   63.39747       44794.12       45789.16    
  5.40      2.30      15.3      1.80      227.      145.     0.106E+04   63.86728       45539.75       46531.12    
  5.50      2.40      15.3      1.90      227.      145.     0.105E+04   64.40422       46259.02       47246.91    
  5.60      2.50      15.3      2.00      227.      145.     0.105E+04   64.80692       47020.09       48004.51    
  5.70      2.60      15.3      2.10      227.      145.     0.105E+04   65.20962       47924.12       48905.19    
  5.80      2.70      15.3      2.20      227.      145.     0.104E+04   65.61232       48793.69       49771.16    
  5.90      2.80      15.3      2.30      227.      145.     0.104E+04   65.88079       49611.89       50586.15    
  6.00      2.90      15.3      2.40      227.      145.     0.104E+04   66.21637       50376.92       51347.77    
  6.10      3.00      15.3      2.50      227.      145.     0.103E+04   66.48484       51076.47       52044.10    
  6.20      3.10      15.3      2.60      227.      145.     0.103E+04   66.82042       51695.34       52659.44    
  6.30      3.20      15.3      2.70      227.      145.     0.103E+04   67.08889       52236.42       53197.33    
  6.40      3.30      15.3      2.80      227.      145.     0.103E+04   67.29024       52672.75       53631.52    
  6.50      3.40      15.3      2.90      227.      145.     0.102E+04   67.49159       52980.18       53937.03    
  6.60      3.50      15.3      3.00      227.      145.     0.102E+04   67.69294       53150.94       54105.53    
  6.70      3.60      15.3      3.10      227.      145.     0.102E+04   67.89013       96664.40       97616.89    
  6.80      3.70      15.3      3.20      227.      145.     0.102E+04   67.92529       94121.87       95075.46    
  6.90      3.80      15.3      3.30      227.      145.     0.102E+04   67.96875       91579.34       92534.02    
  7.00      3.90      15.3      3.40      227.      145.     0.102E+04   68.00391       89036.81       89992.59    
  7.10      4.00      15.3      3.50      227.      145.     0.103E+04   68.05569       86494.28       87451.36    
  7.20      4.10      15.3      3.60      227.      145.     0.103E+04   68.05696       83951.75       84909.96    
  7.30      4.20      15.3      3.70      227.      145.     0.103E+04   68.11705       81409.22       82368.60    
  7.40      4.30      15.3      3.80      227.      145.     0.103E+04   68.05952       78866.69       79827.26    
  7.50      4.40      15.3      3.90      227.      145.     0.103E+04   68.06911       76324.16       77285.97    
  7.60      4.50      15.3      4.00      227.      145.     0.103E+04   68.16243       73781.63       74744.48    
  7.70      4.60      15.3      4.10      227.      145.     0.103E+04   68.11321       71239.10       72203.13    
  7.80      4.70      15.3      4.20      227.      145.     0.103E+04   68.14773       68696.57       69661.70    
  7.90      4.80      15.3      4.30      227.      145.     0.103E+04   68.10682       66154.04       67120.44    
  8.00      4.90      15.3      4.40      227.      145.     0.104E+04   68.14965       67258.39       68225.89    
  8.10      5.00      15.3      4.50      227.      145.     0.104E+04   68.11705       68365.68       69334.34    
  8.20      5.10      15.3      4.60      227.      145.     0.104E+04   68.16819       69476.29       70446.03    
  8.30      5.20      15.3      4.70      227.      145.     0.104E+04   68.14390       70593.09       71564.11    
  8.40      5.30      15.3      4.80      227.      145.     0.104E+04   68.19503       71707.41       72679.50    
  8.50      5.40      15.3      4.90      227.      145.     0.104E+04   68.18736       71641.31       72614.64    
  8.60      5.50      15.3      5.00      227.      145.     0.104E+04   68.24681       71340.58       72314.99    
  8.70      5.60      15.3      5.10      227.      145.     0.104E+04   68.23914       71039.45       72015.21    
  8.80      5.70      15.3      5.20      227.      145.     0.105E+04   68.29858       70737.92       71714.76    
  8.90      5.80      15.3      5.30      227.      145.     0.105E+04   68.22380       70436.01       71414.05    
  9.00      5.90      15.3      5.40      227.      145.     0.105E+04   68.22444       70133.71       71112.89    
  9.10      6.00      15.3      5.50      227.      145.     0.105E+04   68.30050       69831.05       70811.29    
  9.20      6.10      15.3      5.60      227.      145.     0.105E+04   68.30114       69528.03       70509.40    
  9.30      6.20      15.3      5.70      227.      145.     0.105E+04   68.24297       69224.65       70207.21    
  9.40      6.30      15.3      5.80      227.      145.     0.105E+04   68.31233       66780.57       67764.29    
  9.50      6.40      15.3      5.90      227.      145.     0.105E+04   68.36794       66995.12       67978.78    
  9.60      6.50      15.3      6.00      227.      145.     0.105E+04   68.34493       67209.44       68193.13    
  9.70      6.60      15.3      6.10      227.      145.     0.105E+04   68.46765       67423.55       68407.11    
  9.80      6.70      15.3      6.20      227.      145.     0.105E+04   68.44976       67637.44       68621.02    
  9.90      6.80      15.3      6.30      227.      145.     0.105E+04   68.50217       67851.11       68834.64    
  10.0      6.90      15.3      6.40      227.      145.     0.105E+04   68.48427       68064.57       69048.12    
  10.1      7.00      15.3      6.50      227.      145.     0.105E+04   68.60381       68277.81       69261.24    
  10.2      7.10      15.3      6.60      227.      145.     0.105E+04   68.58271       68490.84       69474.29    
  10.3      7.20      15.3      6.70      227.      145.     0.105E+04   68.63193       68703.66       69687.06    
  10.4      7.30      15.3      6.80      227.      145.     0.105E+04   68.61915       68916.27       69899.68    
  10.5      7.40      15.3      6.90      227.      145.     0.105E+04   68.74891       69128.67       70111.95    
  10.6      7.50      15.3      7.00      227.      145.     0.105E+04   68.72462       69340.86       70324.17    
  10.7      7.60      15.3      7.10      227.      145.     0.105E+04   68.77895       69552.85       70536.11    
  10.8      7.70      15.3      7.20      227.      145.     0.105E+04   68.77447       69764.64       70747.89    
  10.9      7.80      15.3      7.30      227.      145.     0.105E+04   68.88953       69976.22       70959.36    
  11.0      7.90      15.3      7.40      227.      145.     0.105E+04   68.88186       70187.60       71170.75    
  11.1      8.00      15.3      7.50      227.      145.     0.105E+04   68.92980       70398.78       71381.88    
  11.2      8.10      15.3      7.60      227.      145.     0.105E+04   68.92213       70609.76       71592.87    
  11.3      8.20      15.3      7.70      227.      145.     0.105E+04   68.90295       70820.54       71803.67    
  11.4      8.30      15.3      7.80      227.      145.     0.105E+04   69.03783       71031.12       72014.12    
  11.5      8.40      15.3      7.90      227.      145.     0.105E+04   69.02696       71241.51       72224.52    
  11.6      8.50      15.3      8.00      227.      145.     0.105E+04   69.00459       71451.71       72434.74    
  11.7      8.60      15.3      8.10      227.      145.     0.105E+04   69.06595       71661.71       72644.68    
  11.8      8.70      15.3      8.20      227.      145.     0.105E+04   69.05828       71871.52       72854.50    
  11.9      8.80      15.3      8.30      227.      145.     0.105E+04   69.05253       72212.70       73195.68    
  12.0      8.90      15.3      8.40      227.      145.     0.105E+04   69.18101       72673.96       73656.81    
  12.1      9.00      15.3      8.50      227.      145.     0.105E+04   69.15864       73148.16       74131.03    
  12.2      9.10      15.3      8.60      227.      145.     0.105E+04   69.15288       73636.29       74619.17    
  12.3      9.20      15.3      8.70      227.      145.     0.105E+04   69.21425       74139.46       75122.28    
  12.4      9.30      15.3      8.80      227.      145.     0.105E+04   69.20849       74658.88       75641.70    
  12.5      9.40      15.3      8.90      227.      145.     0.105E+04   69.19635       75195.90       76178.74    
  12.6      9.50      15.3      9.00      227.      145.     0.105E+04   69.32483       75752.04       76734.75    
  12.7      9.60      15.3      9.10      227.      145.     0.105E+04   69.31077       76328.98       77311.71    
  12.8      9.70      15.3      9.20      227.      145.     0.105E+04   69.31013       76928.62       77911.34    
  12.9      9.80      15.3      9.30      227.      145.     0.105E+04   69.29798       77553.07       78535.81    
  13.0      9.90      15.3      9.40      227.      145.     0.105E+04   69.35935       78204.75       79187.66    
  13.1      10.0      15.3      9.50      227.      145.     0.106E+04   69.35040       78886.37       79880.62    
  13.2      10.1      15.3      9.60      227.      145.     0.108E+04   69.33825       79601.03       80607.58    
  13.3      10.2      15.3      9.70      227.      145.     0.109E+04   69.32930       80352.28       81370.17    
  13.4      10.3      15.3      9.80      227.      145.     0.110E+04   69.45139       81144.17       82173.75    
  13.5      10.4      15.3      9.90      227.      145.     0.111E+04   69.45906       81981.40       83022.31    
  13.6      10.5      15.3      10.0      227.      145.     0.112E+04   69.44692       82869.41       83922.61    
  13.7      10.6      15.3      10.1      227.      145.     0.113E+04   69.43797       83814.54       84879.09    
  13.8      10.7      15.3      10.2      227.      145.     0.115E+04   69.50125       84824.25       85900.54    
  13.9      10.8      15.3      10.3      227.      145.     0.116E+04   69.49230       85907.34       86994.98    
  14.0      10.9      15.3      10.4      227.      145.     0.117E+04   69.49678       87074.29       88174.21    
  14.1      11.0      15.3      10.5      227.      145.     0.118E+04   69.48144       88337.69       89448.96    
  14.2      11.1      15.3      10.6      227.      145.     0.119E+04   69.48080       89712.76       90835.84    
  14.3      11.2      15.3      10.7      227.      145.     0.120E+04   69.60289       91218.12       92352.41    
  14.4      11.3      15.3      10.8      227.      145.     0.122E+04   69.60416       92876.75       94023.06    
  14.5      11.4      15.3      10.9      227.      145.     0.121E+04   69.59521       94717.40       95862.59    
  14.6      11.5      15.3      11.0      227.      145.     0.121E+04   69.58818       96776.46       97920.97    
  14.7      11.6      15.3      11.1      227.      145.     0.121E+04   69.57284       99100.75       100244.1    
  14.8      11.7      15.3      11.2      227.      145.     0.121E+04   69.64443       101751.5       102893.1    
  14.9      11.8      15.3      11.3      227.      145.     0.121E+04   69.63740       104810.2       105950.7    
  15.0      11.9      15.3      11.4      227.      145.     0.121E+04   69.62206       108388.3       109527.1    
  15.1      12.0      15.3      11.5      227.      145.     0.121E+04   69.62973       112641.5       113779.1    
  15.2      12.1      15.3      11.6      227.      145.     0.121E+04   69.62270       117794.8       118931.7    
  15.3      12.2      15.3      11.7      227.      145.     0.121E+04   69.75821       124185.8       125321.5    
  15.4      12.3      15.3      11.8      227.      145.     0.120E+04   69.74287       132345.0       133479.0    
  15.5      12.4      15.3      11.9      227.      145.     0.120E+04   69.73584       143156.4       144289.3    
  15.6      12.5      15.3      12.0      227.      145.     0.120E+04   69.73712       141884.6       143016.7    
  15.7      12.6      15.3      12.1      227.      145.     0.120E+04   69.90651       143850.1       144982.1    
  15.8      12.7      15.3      12.2      227.      145.     0.120E+04   70.04202       146008.3       147140.1    
  15.9      12.8      15.3      12.3      227.      145.     0.120E+04   70.14685       148393.2       149525.0    
  16.0      12.9      15.3      12.4      227.      145.     0.120E+04   70.27405       151047.4       152179.0    
  16.1      13.0      15.3      12.5      227.      145.     0.120E+04   70.40956       154024.6       155156.1    
  16.2      13.1      15.3      12.6      227.      145.     0.120E+04   70.67100       157394.2       158525.4    
  16.3      13.2      15.3      12.7      227.      145.     0.120E+04   70.78989       159720.9       160852.0    
  16.4      13.3      15.3      12.8      227.      145.     0.120E+04   70.62242       158859.9       159990.2    
  16.5      13.4      15.3      12.9      227.      145.     0.121E+04   70.81674       163315.9       164452.8    
  16.6      13.5      15.3      13.0      227.      145.     0.122E+04   70.93563       168540.0       169684.1    
  16.7      13.6      15.3      13.1      227.      145.     0.122E+04   71.12739       174762.7       175913.4    
  16.8      13.7      15.3      13.2      227.      145.     0.123E+04   71.31340       182316.8       183474.5    
  16.9      13.8      15.3      13.3      227.      145.     0.124E+04   71.65378       183490.4       184654.5    
  17.0      13.9      15.3      13.4      227.      145.     0.124E+04   71.65378       189315.2       190484.9    
  17.1      14.0      15.3      13.5      227.      145.     0.125E+04   71.51220       201315.5       202492.2    
  17.2      14.1      15.3      13.6      227.      145.     0.126E+04   71.59657       217238.9       218422.6    
  17.3      14.2      15.3      13.7      227.      145.     0.126E+04   71.71546       221615.4       222807.2    
  17.4      14.3      15.3      13.8      227.      145.     0.126E+04   71.69117       237052.6       238243.3    
  17.5      14.4      15.3      13.9      227.      145.     0.126E+04   71.78194       246544.0       247734.6    
  17.6      14.5      15.3      14.0      227.      145.     0.126E+04   72.15492       270468.7       271658.9    
  17.7      14.6      15.3      14.1      227.      145.     0.126E+04   72.15332       268327.2       269517.5    
  17.8      14.7      15.3      14.2      227.      145.     0.126E+04   72.17377       296968.3       298158.6    
  17.9      14.8      15.3      14.3      227.      145.     0.126E+04   72.29043       302384.7       303572.4    
  18.0      14.9      15.3      14.4      227.      145.     0.127E+04   72.23386       323698.4       324895.3    
  18.1      15.0      15.3      14.5      227.      145.     0.128E+04   72.40900       339601.4       340807.2    
  18.2      15.1      15.3      14.6      227.      145.     0.129E+04   72.42626       337235.6       338450.6    
  18.3      15.2      15.3      14.7      227.      145.     0.128E+04   72.46909       335446.0       336657.1    
 -------------------------------------------------------------------------------------------------------------
 z_tip        D      B_eq     H_cav    V_spud       V_D      B_S           W_BF         Qv_Dpt           V_L  
  [m]        [m]      [m]       [m]      [m^3]     [m^3]      [kN]         [kN]          [kN]            [kN] 
 -------------------------------------------------------------------------------------------------------------
 -------------------------------------------------------------------
 output  file: "JOB.O01"
 load spread factor #=    1; load spread factor =  4.00     (#rows QV_FMODE results=  184)
 -------------------------------------------------------------------
  z_tip        D      B_eq     Qv_Dpt           V_L        Qv_fmode 
   [m]        [m]      [m]      [kN]            [kN]                
 -------------------------------------------------------------------
  0.00      0.00     0.100E-02  0.1000000E-02  0.1000000E-02 n/r
 0.100      0.00     0.246      0.2338996      0.2366989     GSC
 0.200      0.00     0.492      0.9776733       1.007068     GSC
 0.300      0.00     0.738       2.294434       2.396015     GSC
 0.400      0.00     0.984       4.247293       4.489448     GSC
 0.500      0.00      1.23       6.899364       7.373952     GSC
 0.600      0.00      1.48       10.31376       11.43377     GSC
 0.700      0.00      1.72       14.98585       17.10274     SQC
 0.800      0.00      1.97       20.09105       23.63376     SQC
 0.900      0.00      2.21       26.08295       31.54570     SQC
  1.00      0.00      2.46       33.01009       40.96927     SQC
  1.10      0.00      2.70       41.02954       52.12766     GSC
  1.20      0.00      2.95       50.34324       65.27551     GSC
  1.30      0.00      3.52       74.96031       94.93951     SQC
  1.40      0.00      4.18       111.5955       138.8306     SQC
  1.50      0.00      4.83       157.5960       194.8368     SQC
  1.60      0.00      5.49       213.8760       264.3927     SQC
  1.70      0.00      6.14       281.3497       348.9006     SQC
  1.80      0.00      6.79       360.9312       449.7424     SQC
  1.90      0.00      7.45       453.5346       568.5285     SQC
  2.00      0.00      8.10       560.0740       706.4524     SQC
  2.10      0.00      8.76       681.4636       865.0268     SQC
  2.20      0.00      9.41       818.6175       1045.606     SQC
  2.30      0.00      10.1       972.4499       1249.602     SQC
  2.40      0.00      10.7       1143.875       1479.029     SQC
  2.50      0.00      11.4       1333.807       1734.341     SQC
  2.60      0.00      12.0       1543.159       2017.626     SQC
  2.70      0.00      12.7       1772.847       2330.392     SQC
  2.80      0.00      13.3       2023.783       2672.837     SQC
  2.90      0.00      14.0       2296.883       3048.158     SQC
  3.00      0.00      14.6       2593.060       3456.438     SQC
  3.10      0.00      15.3       2913.229       3900.798     SQC
  3.20     0.100      15.3       3608.051       4618.461     SQC
  3.30     0.200      15.3       4705.117       5737.422     SQC
  3.40     0.300      15.3       6807.575       7858.948     SQC
  3.50     0.400      15.3       12931.17       14001.75     SQC
  3.60     0.500      15.3       27544.27       28630.69     GSS
  3.70     0.600      15.3       28861.25       29936.13     GSS
  3.80     0.700      15.3       30180.92       31246.50     GSS
  3.90     0.800      15.3       31503.28       32561.28     GSS
  4.00     0.900      15.3       32836.89       33888.26     GSS
  4.10      1.00      15.3       34094.06       35139.60     GSS
  4.20      1.10      15.3       35100.51       36140.50     GSS
  4.30      1.20      15.3       36092.21       37127.38     GSS
  4.40      1.30      15.3       37057.64       38087.99     GSS
  4.50      1.40      15.3       38004.91       39030.90     GSS
  4.60      1.50      15.3       38933.65       39955.25     GSS
  4.70      1.60      15.3       39835.44       40853.02     GSS
  4.80      1.70      15.3       40718.13       41731.70     GSS
  4.90      1.80      15.3       41576.08       42585.89     GSS
  5.00      1.90      15.3       42416.90       43422.72     GSS
  5.10      2.00      15.3       43232.40       44234.59     GSS
  5.20      2.10      15.3       44022.35       45020.87     GSS
  5.30      2.20      15.3       44794.12       45789.16     GSS
  5.40      2.30      15.3       45539.75       46531.12     GSS
  5.50      2.40      15.3       46259.02       47246.91     GSS
  5.60      2.50      15.3       47020.09       48004.51     GSS
  5.70      2.60      15.3       47924.12       48905.19     GSS
  5.80      2.70      15.3       48793.69       49771.16     GSS
  5.90      2.80      15.3       49611.89       50586.15     GSS
  6.00      2.90      15.3       50376.92       51347.77     GSS
  6.10      3.00      15.3       51076.47       52044.10     GSS
  6.20      3.10      15.3       51695.34       52659.44     GSS
  6.30      3.20      15.3       52236.42       53197.33     GSS
  6.40      3.30      15.3       52672.75       53631.52     GSS
  6.50      3.40      15.3       52980.18       53937.03     GSS
  6.60      3.50      15.3       53150.94       54105.53     GSS
  6.70      3.60      15.3       96664.40       97616.89     PTC>GSS
  6.80      3.70      15.3       94121.87       95075.46     PTC>GSS
  6.90      3.80      15.3       91579.34       92534.02     PTC>GSS
  7.00      3.90      15.3       89036.81       89992.59     PTC>GSS
  7.10      4.00      15.3       86494.28       87451.36     PTC>GSS
  7.20      4.10      15.3       83951.75       84909.96     PTC>GSS
  7.30      4.20      15.3       81409.22       82368.60     PTC>GSS
  7.40      4.30      15.3       78866.69       79827.26     PTC>GSS
  7.50      4.40      15.3       76324.16       77285.97     PTC>GSS
  7.60      4.50      15.3       73781.63       74744.48     PTC>GSS
  7.70      4.60      15.3       71239.10       72203.13     PTC>GSS
  7.80      4.70      15.3       68696.57       69661.70     PTC>GSS
  7.90      4.80      15.3       66154.04       67120.44     GSS
  8.00      4.90      15.3       67258.39       68225.89     GSS
  8.10      5.00      15.3       68365.68       69334.34     GSS
  8.20      5.10      15.3       69476.29       70446.03     GSS
  8.30      5.20      15.3       70593.09       71564.11     GSS
  8.40      5.30      15.3       71707.41       72679.50     GSS
  8.50      5.40      15.3       71641.31       72614.64     PTS>GSC
  8.60      5.50      15.3       71340.58       72314.99     PTS>GSC
  8.70      5.60      15.3       71039.45       72015.21     PTS>GSC
  8.80      5.70      15.3       70737.92       71714.76     PTS>GSC
  8.90      5.80      15.3       70436.01       71414.05     PTS>GSC
  9.00      5.90      15.3       70133.71       71112.89     PTS>GSC
  9.10      6.00      15.3       69831.05       70811.29     PTS>GSC
  9.20      6.10      15.3       69528.03       70509.40     PTS>GSC
  9.30      6.20      15.3       69224.65       70207.21     PTS>GSC
  9.40      6.30      15.3       66780.57       67764.29     GSC
  9.50      6.40      15.3       66995.12       67978.78     GSC
  9.60      6.50      15.3       67209.44       68193.13     GSC
  9.70      6.60      15.3       67423.55       68407.11     GSC
  9.80      6.70      15.3       67637.44       68621.02     GSC
  9.90      6.80      15.3       67851.11       68834.64     GSC
  10.0      6.90      15.3       68064.57       69048.12     GSC
  10.1      7.00      15.3       68277.81       69261.24     GSC
  10.2      7.10      15.3       68490.84       69474.29     GSC
  10.3      7.20      15.3       68703.66       69687.06     GSC
  10.4      7.30      15.3       68916.27       69899.68     GSC
  10.5      7.40      15.3       69128.67       70111.95     GSC
  10.6      7.50      15.3       69340.86       70324.17     GSC
  10.7      7.60      15.3       69552.85       70536.11     GSC
  10.8      7.70      15.3       69764.64       70747.89     GSC
  10.9      7.80      15.3       69976.22       70959.36     GSC
  11.0      7.90      15.3       70187.60       71170.75     GSC
  11.1      8.00      15.3       70398.78       71381.88     GSC
  11.2      8.10      15.3       70609.76       71592.87     GSC
  11.3      8.20      15.3       70820.54       71803.67     GSC
  11.4      8.30      15.3       71031.12       72014.12     GSC
  11.5      8.40      15.3       71241.51       72224.52     GSC
  11.6      8.50      15.3       71451.71       72434.74     GSC
  11.7      8.60      15.3       71661.71       72644.68     GSC
  11.8      8.70      15.3       71871.52       72854.50     GSC
  11.9      8.80      15.3       72212.70       73195.68     SQC
  12.0      8.90      15.3       72673.96       73656.81     SQC
  12.1      9.00      15.3       73148.16       74131.03     SQC
  12.2      9.10      15.3       73636.29       74619.17     SQC
  12.3      9.20      15.3       74139.46       75122.28     SQC
  12.4      9.30      15.3       74658.88       75641.70     SQC
  12.5      9.40      15.3       75195.90       76178.74     SQC
  12.6      9.50      15.3       75752.04       76734.75     SQC
  12.7      9.60      15.3       76328.98       77311.71     SQC
  12.8      9.70      15.3       76928.62       77911.34     SQC
  12.9      9.80      15.3       77553.07       78535.81     SQC
  13.0      9.90      15.3       78204.75       79187.66     SQC
  13.1      10.0      15.3       78886.37       79880.62     SQC
  13.2      10.1      15.3       79601.03       80607.58     SQC
  13.3      10.2      15.3       80352.28       81370.17     SQC
  13.4      10.3      15.3       81144.17       82173.75     SQC
  13.5      10.4      15.3       81981.40       83022.31     SQC
  13.6      10.5      15.3       82869.41       83922.61     SQC
  13.7      10.6      15.3       83814.54       84879.09     SQC
  13.8      10.7      15.3       84824.25       85900.54     SQC
  13.9      10.8      15.3       85907.34       86994.98     SQC
  14.0      10.9      15.3       87074.29       88174.21     SQC
  14.1      11.0      15.3       88337.69       89448.96     SQC
  14.2      11.1      15.3       89712.76       90835.84     SQC
  14.3      11.2      15.3       91218.12       92352.41     SQC
  14.4      11.3      15.3       92876.75       94023.06     SQC
  14.5      11.4      15.3       94717.40       95862.59     SQC
  14.6      11.5      15.3       96776.46       97920.97     SQC
  14.7      11.6      15.3       99100.75       100244.1     SQC
  14.8      11.7      15.3       101751.5       102893.1     SQC
  14.9      11.8      15.3       104810.2       105950.7     SQC
  15.0      11.9      15.3       108388.3       109527.1     SQC
  15.1      12.0      15.3       112641.5       113779.1     SQC
  15.2      12.1      15.3       117794.8       118931.7     SQC
  15.3      12.2      15.3       124185.8       125321.5     SQC
  15.4      12.3      15.3       132345.0       133479.0     SQC
  15.5      12.4      15.3       143156.4       144289.3     SQC
  15.6      12.5      15.3       141884.6       143016.7     PTC>PTS>SQC
  15.7      12.6      15.3       143850.1       144982.1     PTC>PTS>SQC
  15.8      12.7      15.3       146008.3       147140.1     PTC>PTS>SQC
  15.9      12.8      15.3       148393.2       149525.0     PTC>PTS>SQC
  16.0      12.9      15.3       151047.4       152179.0     PTC>PTS>SQC
  16.1      13.0      15.3       154024.6       155156.1     PTC>PTS>SQC
  16.2      13.1      15.3       157394.2       158525.4     PTC>PTS>SQC
  16.3      13.2      15.3       159720.9       160852.0     SQC
  16.4      13.3      15.3       158859.9       159990.2     PTC>PTS>SQC
  16.5      13.4      15.3       163315.9       164452.8     PTC>PTS>SQC
  16.6      13.5      15.3       168540.0       169684.1     PTC>PTS>SQC
  16.7      13.6      15.3       174762.7       175913.4     PTC>PTS>SQC
  16.8      13.7      15.3       182316.8       183474.5     PTC>PTS>SQC
  16.9      13.8      15.3       183490.4       184654.5     SQC
  17.0      13.9      15.3       189315.2       190484.9     PTC>PTS>SQC
  17.1      14.0      15.3       201315.5       202492.2     PTC>PTS>SQC
  17.2      14.1      15.3       217238.9       218422.6     PTC>PTS>SQC
  17.3      14.2      15.3       221615.4       222807.2     SQC
  17.4      14.3      15.3       237052.6       238243.3     PTC>PTS>SQC
  17.5      14.4      15.3       246544.0       247734.6     PTC>PTS>PTC>PTS>SQC
  17.6      14.5      15.3       270468.7       271658.9     PTC>PTS>SQC
  17.7      14.6      15.3       268327.2       269517.5     PTC>PTS>PTC>PTS>SQC
  17.8      14.7      15.3       296968.3       298158.6     SQC
  17.9      14.8      15.3       302384.7       303572.4     PTC>PTS>PTC>PTS>SQC
  18.0      14.9      15.3       323698.4       324895.3     PTC>PTS>SQC
  18.1      15.0      15.3       339601.4       340807.2     PTC>PTS>PTC>PTS>PTC>GSS
  18.2      15.1      15.3       337235.6       338450.6     PTS>PTC>PTS>PTC>GSS
  18.3      15.2      15.3       335446.0       336657.1     PTS>PTC>PTS>PTC>GSS
 -------------------------------------------------------------------
  z_tip        D      B_eq     Qv_Dpt           V_L        Qv_fmode 
   [m]        [m]      [m]      [kN]            [kN]                
 -------------------------------------------------------------------

out_sp_pen_sub: < --- BEGIN
 -------------------------------------------------------------------
 load spread factor #=    1; load spread factor =  4.00    
 (#rows QV_FMODE summary results=   25)
 -------------------------------------------------------------------
  from      to        from      to     thickness             
  z_tip     z_tip      D         D         T        Qv_fmode 
   [m]       [m]      [m]       [m]       [m]        [-]     
 -------------------------------------------------------------------
 0.100     0.700      0.00      0.00     0.600       GSC
 0.700      1.10      0.00      0.00     0.400       SQC
  1.10      1.30      0.00      0.00     0.200       GSC
  1.30      3.60      0.00     0.500      2.30       SQC
  3.60      6.70     0.500      3.60      3.10       GSS
  6.70      7.90      3.60      4.80      1.20       PTC>GSS
  7.90      8.50      4.80      5.40     0.600       GSS
  8.50      9.40      5.40      6.30     0.900       PTS>GSC
  9.40      11.9      6.30      8.80      2.50       GSC
  11.9      15.6      8.80      12.5      3.70       SQC
  15.6      16.3      12.5      13.2     0.700       PTC>PTS>SQC
  16.3      16.4      13.2      13.3     0.100       SQC
  16.4      16.9      13.3      13.8     0.500       PTC>PTS>SQC
  16.9      17.0      13.8      13.9     0.100       SQC
  17.0      17.3      13.9      14.2     0.300       PTC>PTS>SQC
  17.3      17.4      14.2      14.3     0.100       SQC
  17.4      17.5      14.3      14.4     0.100       PTC>PTS>SQC
  17.5      17.6      14.4      14.5     0.100       PTC>PTS>PTC>PTS>SQC
  17.6      17.7      14.5      14.6     0.100       PTC>PTS>SQC
  17.7      17.8      14.6      14.7     0.100       PTC>PTS>PTC>PTS>SQC
  17.8      17.9      14.7      14.8     0.100       SQC
  17.9      18.0      14.8      14.9     0.100       PTC>PTS>PTC>PTS>SQC
  18.0      18.1      14.9      15.0     0.100       PTC>PTS>SQC
  18.1      18.2      15.0      15.1     0.100       PTC>PTS>PTC>PTS>PTC>GSS
  18.2      18.3      15.1      15.2     0.100       PTS>PTC>PTS>PTC>GSS
 -------------------------------------------------------------------
  from      to        from      to     thickness             
  z_tip     z_tip      D         D         T        Qv_fmode 
   [m]       [m]      [m]       [m]       [m]        [-]     
 -------------------------------------------------------------------

iso_19905_1: nlines_errors  :    0
iso_19905_1: nlines_warnings:    5

WARNING. gmdata: soil layer # 3; BDEN=7.00 kN/m^3 is inconsistent with SU_TOP=90.00 kPa.
WARNING. gmdata: soil layer # 4; PHI < 20 degrees [ISO 19905-1 Table A.9.3-3].
WARNING. gmdata: soil layer # 7; BDEN=7.00 kN/m^3 is inconsistent with SU_TOP=71.00 kPa.
WARNING. gmdata: soil layer # 9; BDEN=7.00 kN/m^3 is inconsistent with SU_TOP=78.00 kPa.
WARNING. gmdata: soil layer # 10; PHI < 20 degrees [ISO 19905-1 Table A.9.3-3].

 < P_spud                [kN]:   66910.00    
   load spread factor #=    1; load spread factor =  4.00    
   is_PT_occurring      [t/f]:         F
 > z_tip_penetration      [m]:       6.63    
out_sp_pen_sub: END --- >


input   file: "JOB.I01"
output  file: "JOB.O01"
iso_19905_1: normal program end
