user: 请显示与下面列出的汇编代码相匹配的C语言实现：

```asm
<main>:
  endbr64
  push   %rbp
  push   %rbx
  sub    $0x88,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0x78(%rsp)
  xor    %eax,%eax
  lea    0x18(%rsp),%rcx
  lea    0x10(%rsp),%rdx
  lea    0x8(%rsp),%rsi
  mov    %rsp,%rdi
  lea    0x38(%rsp),%rax
  push   %rax
  lea    0x38(%rsp),%rax
  push   %rax
  lea    0x38(%rsp),%r9
  lea    0x30(%rsp),%r8
  call   <params_exactN@plt>
  movsd  0x18(%rsp),%xmm0
  divsd  D23(%rip),%xmm0
  lea    D93(%rip),%rdi
  mov    $0x1,%eax
  call   <printf@plt>
  lea    0x50(%rsp),%rsi
  lea    D56(%rip),%rdi
  call   <scanf@plt>
  add    $0x10,%rsp
  cmp    $0x1,%eax
  je     L35
  lea    D66(%rip),%rdi
  call   <puts@plt>
  mov    $0x1,%ebx
L40:
  mov    0x78(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L30
  mov    %ebx,%eax
  add    $0x88,%rsp
  pop    %rbx
  pop    %rbp
  ret
L35:
  mov    %eax,%ebx
  movsd  D23(%rip),%xmm0
  mulsd  0x40(%rsp),%xmm0
  movsd  %xmm0,0x40(%rsp)
  lea    0x60(%rsp),%rcx
  lea    0x58(%rsp),%rdx
  lea    0x50(%rsp),%rsi
  lea    0x48(%rsp),%rdi
  lea    0x70(%rsp),%r9
  lea    0x68(%rsp),%r8
  call   <exactN@plt>
  mov    %eax,%ebp
  test   %eax,%eax
  jne    L39
  movsd  0x40(%rsp),%xmm0
  movsd  0x10(%rsp),%xmm1
  comisd %xmm1,%xmm0
  ja     L49
L28:
  lea    D129(%rip),%rdi
  call   <puts@plt>
  movsd  D95(%rip),%xmm2
  movapd %xmm2,%xmm3
  mulsd  0x60(%rsp),%xmm3
  movsd  0x70(%rsp),%xmm5
  movsd  0x68(%rsp),%xmm4
  mulsd  D68(%rip),%xmm3
  mulsd  0x58(%rsp),%xmm2
  movsd  0x50(%rsp),%xmm1
  movsd  0x48(%rsp),%xmm0
  lea    D40(%rip),%rdi
  mov    $0x6,%eax
  call   <printf@plt>
  mov    %ebp,%ebx
  jmp    L40
L39:
  movsd  0x8(%rsp),%xmm0
  divsd  D23(%rip),%xmm0
  lea    D80(%rip),%rdi
  mov    $0x1,%eax
  call   <printf@plt>
  jmp    L40
L49:
  movsd  D23(%rip),%xmm2
  divsd  %xmm2,%xmm0
  divsd  %xmm2,%xmm1
  lea    D13(%rip),%rdi
  mov    $0x2,%eax
  call   <printf@plt>
  jmp    L28
L30:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: read_data
arguments: label: D23
type: f64
name: read_data
arguments: label: D93
type: string
name: read_data
arguments: label: D56
type: string
name: read_data
arguments: label: D66
type: string
name: read_data
arguments: label: D129
type: string
name: read_data
arguments: label: D95
type: f64
name: read_data
arguments: label: D68
type: f64
name: read_data
arguments: label: D40
type: string
name: read_data
arguments: label: D80
type: string
name: read_data
arguments: label: D13
type: string

tool: D23: 1000

tool: D93: Enter  x  (in km; 0.0 <= x <= %5.1f):   

tool: D56: %lf

tool: D66: ... input error; exiting

tool: D129: Results from Test N:

tool: D95: 31556926

tool: D68: 100

tool: D40:      H    = ice thickness        = %12.5f (m)
     h_x  = surface slope        = %12.5e
     u    = ice velocity         = %12.5f (m a-1)
     M    = surface mass balance = %12.5f (cm a-1)
     B    = ice hardness         = %12.5e (Pa s^(1/3))
     beta = ice hardness         = %12.5e (Pa s m-1)


tool: D80: SIMPLEN ERROR:  x  out of allowed domain  0.0 <= x <= %5.1f km
   ... ending ...


tool: D13: WARNING:  x = %5.1f km is past the calving front at xc = %5.1f km


assistant: ```c
int main() {
  double H0, L0, xc, a, Hela, k, Hc, Tc, x, H, hx, u, M, B, beta;
  const double secpera = 31556926.0;
  int scanret, retvalN;
  params_exactN(&H0, &L0, &xc, &a, &Hela, &k, &Hc, &Tc);
  printf("Enter  x  (in km; 0.0 <= x <= %5.1f):   ", L0 / 1000.0);
  scanret = scanf("%lf", &x);
  if (scanret != 1) {
    printf("... input error; exiting\n");
    return 1;
  }
  x = x * 1000.0;
  retvalN = exactN(x, &H, &hx, &u, &M, &B, &beta);
  if (retvalN) {
    printf("SIMPLEN ERROR:  x  out of allowed domain  0.0 <= x <= %5.1f km\n"
           "   ... ending ...\n",
           L0 / 1000.0);
    return 1;
  }
  if (x > xc) {
    printf("WARNING:  x = %5.1f km is past the calving front at xc = %5.1f km\n", x / 1000.0, xc / 1000.0);
  }
  printf("Results from Test N:\n");
  printf("     H    = ice thickness        = %12.5f (m)\n"
         "     h_x  = surface slope        = %12.5e\n"
         "     u    = ice velocity         = %12.5f (m a-1)\n"
         "     M    = surface mass balance = %12.5f (cm a-1)\n"
         "     B    = ice hardness         = %12.5e (Pa s^(1/3))\n"
         "     beta = ice hardness         = %12.5e (Pa s m-1)\n",
         H, hx, u * secpera, M * secpera * 100.0, B, beta);
  return 0;
}
```