user: What is the c source code of the assembly code below?

```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  D78(%rip),%xmm0
  lea    D6(%rip),%rdi
  mov    $0x1,%eax
  call   <printf@plt>
  lea    0x50(%rsp),%rsi
  lea    D103(%rip),%rdi
  call   <scanf@plt>
  add    $0x10,%rsp
  cmp    $0x1,%eax
  je     L21
  lea    D122(%rip),%rdi
  call   <puts@plt>
  mov    $0x1,%ebx
L35:
  mov    0x78(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L48
  mov    %ebx,%eax
  add    $0x88,%rsp
  pop    %rbx
  pop    %rbp
  ret
L21:
  mov    %eax,%ebx
  movsd  D78(%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    L59
  movsd  0x40(%rsp),%xmm0
  movsd  0x10(%rsp),%xmm1
  comisd %xmm1,%xmm0
  ja     L10
L0:
  lea    D51(%rip),%rdi
  call   <puts@plt>
  movsd  D123(%rip),%xmm2
  movapd %xmm2,%xmm3
  mulsd  0x60(%rsp),%xmm3
  movsd  0x70(%rsp),%xmm5
  movsd  0x68(%rsp),%xmm4
  mulsd  D76(%rip),%xmm3
  mulsd  0x58(%rsp),%xmm2
  movsd  0x50(%rsp),%xmm1
  movsd  0x48(%rsp),%xmm0
  lea    D50(%rip),%rdi
  mov    $0x6,%eax
  call   <printf@plt>
  mov    %ebp,%ebx
  jmp    L35
L59:
  movsd  0x8(%rsp),%xmm0
  divsd  D78(%rip),%xmm0
  lea    D71(%rip),%rdi
  mov    $0x1,%eax
  call   <printf@plt>
  jmp    L35
L10:
  movsd  D78(%rip),%xmm2
  divsd  %xmm2,%xmm0
  divsd  %xmm2,%xmm1
  lea    D82(%rip),%rdi
  mov    $0x2,%eax
  call   <printf@plt>
  jmp    L0
L48:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: read_data
arguments: data_tag: D78
data_type: f64
name: read_data
arguments: data_tag: D6
data_type: string
name: read_data
arguments: data_tag: D103
data_type: string
name: read_data
arguments: data_tag: D122
data_type: string
name: read_data
arguments: data_tag: D51
data_type: string
name: read_data
arguments: data_tag: D123
data_type: f64
name: read_data
arguments: data_tag: D76
data_type: f64
name: read_data
arguments: data_tag: D50
data_type: string
name: read_data
arguments: data_tag: D71
data_type: string
name: read_data
arguments: data_tag: D82
data_type: string

tool: D78: 1000

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

tool: D103: %lf

tool: D122: ... input error; exiting

tool: D51: Results from Test N:

tool: D123: 31556926

tool: D76: 100

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


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


tool: D82: 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/year)\n"
         "     M    = surface mass balance = %12.5f (cm/year)\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;
}
```