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

```asm
<main>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x90,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  movsd  D2(%rip),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  lea    -0x60(%rbp),%r9
  lea    -0x68(%rbp),%r8
  lea    -0x70(%rbp),%rcx
  lea    -0x78(%rbp),%rdx
  lea    -0x80(%rbp),%rsi
  lea    -0x88(%rbp),%rax
  lea    -0x50(%rbp),%rdi
  push   %rdi
  lea    -0x58(%rbp),%rdi
  push   %rdi
  mov    %rax,%rdi
  call   <params_exactN@plt>
  add    $0x10,%rsp
  movsd  -0x80(%rbp),%xmm0
  movsd  D27(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  lea    D109(%rip),%rax
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <printf@plt>
  lea    -0x48(%rbp),%rax
  mov    %rax,%rsi
  lea    D0(%rip),%rax
  mov    %rax,%rdi
  call   <scanf@plt>
  mov    %eax,-0x90(%rbp)
  cmpl   $0x1,-0x90(%rbp)
  je     L11
  lea    D92(%rip),%rax
  mov    %rax,%rdi
  call   <puts@plt>
  mov    $0x1,%eax
  jmp    L30
L11:
  movsd  -0x48(%rbp),%xmm1
  movsd  D27(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  mov    -0x48(%rbp),%rax
  lea    -0x18(%rbp),%r9
  lea    -0x20(%rbp),%r8
  lea    -0x28(%rbp),%rcx
  lea    -0x30(%rbp),%rdx
  lea    -0x38(%rbp),%rsi
  lea    -0x40(%rbp),%rdi
  movq   %rax,%xmm0
  call   <exactN@plt>
  mov    %eax,-0x8c(%rbp)
  cmpl   $0x0,-0x8c(%rbp)
  je     L27
  movsd  -0x80(%rbp),%xmm0
  movsd  D27(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  lea    D7(%rip),%rax
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <printf@plt>
  mov    $0x1,%eax
  jmp    L30
L27:
  movsd  -0x48(%rbp),%xmm0
  movsd  -0x78(%rbp),%xmm1
  comisd %xmm1,%xmm0
  jbe    L24
  movsd  -0x78(%rbp),%xmm0
  movsd  D27(%rip),%xmm1
  movapd %xmm0,%xmm2
  divsd  %xmm1,%xmm2
  movsd  -0x48(%rbp),%xmm0
  movsd  D27(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movapd %xmm2,%xmm1
  movq   %rax,%xmm0
  lea    D69(%rip),%rax
  mov    %rax,%rdi
  mov    $0x2,%eax
  call   <printf@plt>
L24:
  lea    D97(%rip),%rax
  mov    %rax,%rdi
  call   <puts@plt>
  movsd  -0x18(%rbp),%xmm4
  movsd  -0x20(%rbp),%xmm3
  movsd  -0x28(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x10(%rbp),%xmm1
  movsd  D139(%rip),%xmm0
  movapd %xmm1,%xmm2
  mulsd  %xmm0,%xmm2
  movsd  -0x30(%rbp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x10(%rbp),%xmm1
  movsd  -0x38(%rbp),%xmm0
  mov    -0x40(%rbp),%rax
  movapd %xmm4,%xmm5
  movapd %xmm3,%xmm4
  movapd %xmm2,%xmm3
  movapd %xmm1,%xmm2
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  lea    D79(%rip),%rax
  mov    %rax,%rdi
  mov    $0x6,%eax
  call   <printf@plt>
  mov    $0x0,%eax
L30:
  mov    -0x8(%rbp),%rdx
  sub    %fs:0x28,%rdx
  je     L19
  call   <__stack_chk_fail@plt>
L19:
  leave
  ret
```

tool_calls: name: retrieve_data
arguments: tag: D2
type_infered: f64
name: retrieve_data
arguments: tag: D27
type_infered: f64
name: retrieve_data
arguments: tag: D109
type_infered: string
name: retrieve_data
arguments: tag: D0
type_infered: string
name: retrieve_data
arguments: tag: D92
type_infered: string
name: retrieve_data
arguments: tag: D7
type_infered: string
name: retrieve_data
arguments: tag: D69
type_infered: string
name: retrieve_data
arguments: tag: D97
type_infered: string
name: retrieve_data
arguments: tag: D139
type_infered: f64
name: retrieve_data
arguments: tag: D79
type_infered: string

tool: D2: 31556926

tool: D27: 1000

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

tool: D0: %lf

tool: D92: ... input error; exiting

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


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


tool: D97: Results from Test N:

tool: D139: 100

tool: D79:      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)


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;
}
```