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

```asm
<latsens>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  movq   %xmm0,%r15
  movq   %xmm1,%r14
  movsd  %xmm2,0x8(%rsp)
  movq   %xmm3,%rbx
  movsd  %xmm4,0x10(%rsp)
  movq   %xmm5,%rbp
  movsd  %xmm6,0x20(%rsp)
  movsd  %xmm7,0x28(%rsp)
  mov    %rdi,%r13
  mov    %rsi,%r12
  movsd  0x8(%rsp),%xmm7
  comisd D70(%rip),%xmm7
  jbe    L43
  movapd %xmm4,%xmm2
  movsd  0x70(%rsp),%xmm4
  movsd  D62(%rip),%xmm3
  movq   %rbx,%xmm1
  addsd  D162(%rip),%xmm1
  call   <lkdrag@plt>
  movss  %xmm0,0x1c(%rsp)
L39:
  pxor   %xmm0,%xmm0
  comisd 0x8(%rsp),%xmm0
  jb     L3
  movq   %r15,%xmm7
  movq   %r14,%xmm3
  comisd %xmm3,%xmm7
  ja     L9
L3:
  movq   %r15,%xmm6
  movq   %r15,%xmm0
  subsd  D162(%rip),%xmm0
  mulsd  D127(%rip),%xmm0
  subsd  D67(%rip),%xmm6
  divsd  %xmm6,%xmm0
  call   <exp@plt>
  mulsd  D32(%rip),%xmm0
L45:
  movsd  D110(%rip),%xmm2
  movsd  0x20(%rsp),%xmm5
  movapd %xmm5,%xmm1
  mulsd  %xmm2,%xmm1
  movq   %rbp,%xmm3
  subsd  %xmm1,%xmm3
  divsd  %xmm3,%xmm5
  movapd %xmm5,%xmm1
  movsd  D177(%rip),%xmm4
  mulsd  %xmm4,%xmm1
  mulsd  %xmm0,%xmm2
  movapd %xmm2,%xmm3
  movq   %rbp,%xmm2
  subsd  %xmm3,%xmm2
  divsd  %xmm2,%xmm0
  mulsd  %xmm4,%xmm0
  subsd  %xmm0,%xmm1
  movss  0x1c(%rsp),%xmm7
  movaps %xmm7,%xmm0
  xorps  D133(%rip),%xmm0
  cvtss2sd %xmm0,%xmm0
  movsd  0x10(%rsp),%xmm4
  mulsd  %xmm4,%xmm0
  movsd  0x28(%rsp),%xmm6
  mulsd  %xmm6,%xmm0
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,0x0(%r13)
  pxor   %xmm0,%xmm0
  cvtss2sd %xmm7,%xmm0
  mulsd  %xmm4,%xmm0
  mulsd  %xmm6,%xmm0
  mulsd  D159(%rip),%xmm0
  movq   %rbx,%xmm1
  addsd  D162(%rip),%xmm1
  movq   %r15,%xmm5
  subsd  %xmm5,%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,(%r12)
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L43:
  movsd  0x70(%rsp),%xmm4
  movsd  D25(%rip),%xmm3
  movsd  0x10(%rsp),%xmm2
  movq   %rbx,%xmm1
  addsd  D162(%rip),%xmm1
  call   <lkdrag@plt>
  movss  %xmm0,0x1c(%rsp)
  jmp    L39
L9:
  movsd  D162(%rip),%xmm2
  movq   %r15,%xmm5
  movapd %xmm7,%xmm0
  subsd  %xmm2,%xmm0
  mulsd  D131(%rip),%xmm0
  addsd  D134(%rip),%xmm5
  movapd %xmm5,%xmm1
  subsd  %xmm2,%xmm1
  divsd  %xmm1,%xmm0
  call   <exp@plt>
  mulsd  D32(%rip),%xmm0
  jmp    L45
```

tool_calls: name: unpack_data
arguments: tag: D70
type_guessed: f64
name: unpack_data
arguments: tag: D62
type_guessed: f64
name: unpack_data
arguments: tag: D162
type_guessed: f64
name: unpack_data
arguments: tag: D127
type_guessed: f64
name: unpack_data
arguments: tag: D67
type_guessed: f64
name: unpack_data
arguments: tag: D32
type_guessed: byte[8]
name: unpack_data
arguments: tag: D110
type_guessed: f64
name: unpack_data
arguments: tag: D177
type_guessed: f64
name: unpack_data
arguments: tag: D133
type_guessed: byte[16]
name: unpack_data
arguments: tag: D159
type_guessed: f64
name: unpack_data
arguments: tag: D25
type_guessed: f64
name: unpack_data
arguments: tag: D131
type_guessed: f64
name: unpack_data
arguments: tag: D134
type_guessed: f64

tool: D70: 0

tool: D62: 0.005

tool: D162: 273.15

tool: D127: 21.874

tool: D67: 7.66

tool: {"D32": [0x5a,0x64,0x3b,0xdf,0x4f,0x8d,0xe3,0x3f]}

tool: D110: 0.378

tool: D177: 0.622

tool: {"D133": [0x00,0x00,0x00,0x80,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00]}

tool: D159: 1013

tool: D25: 0.0045

tool: D131: 17.269

tool: D134: 237.3

assistant: ```c
void latsens(double Tsurf, double Tcutk, double hice, double tair, double wind, double pressure, double vp, double air_density, double *evap, double *qsen, double wind_h) {
  float dragcoeff;
  double eog, delT;
  double qair, qlake;
  double delq;
  if (hice > 0.)
    dragcoeff = lkdrag(Tsurf, tair + 273.15, wind, 0.005, wind_h);
  else
    dragcoeff = lkdrag(Tsurf, tair + 273.15, wind, 0.0045, wind_h);
  if ((hice <= 0.) && (Tsurf > Tcutk)) {
    eog = .611 * exp(17.269 * (Tsurf - 273.15) / (Tsurf + 237.3 - 273.15));
  } else {
    eog = .611 * exp(21.874 * (Tsurf - 273.15) / (Tsurf - 7.66));
  }
  qlake = 0.622 * (eog / (pressure - 0.378 * eog));
  qair = 0.622 * (vp / (pressure - 0.378 * vp));
  delq = qair - qlake;
  *evap = -1 * dragcoeff * wind * air_density * delq;
  delT = tair + 273.15 - Tsurf;
  *qsen = dragcoeff * wind * air_density * 1013.0;
  *qsen = *qsen * delT;
}
```