user: Please show me the C implementation that matches the assembly code listed below:

```asm
<g_0>:
  endbr64
  push   %rbx
  movapd %xmm0,%xmm5
  mov    %edi,%ebx
  sub    $0x20,%rsp
  movsd  D40(%rip),%xmm0
  movsd  %xmm3,0x8(%rsp)
  comisd %xmm4,%xmm0
  jae    L51
  comisd D78(%rip),%xmm4
  jbe    L35
L51:
  movapd %xmm5,%xmm3
  divsd  D72(%rip),%xmm3
L74:
  movsd  D25(%rip),%xmm0
  mov    $0x4,%edi
  movsd  %xmm2,0x18(%rsp)
  movsd  %xmm5,0x10(%rsp)
  mulsd  %xmm3,%xmm0
  mulsd  %xmm3,%xmm0
  mulsd  D9(%rip),%xmm3
  addsd  %xmm3,%xmm0
  movsd  %xmm0,(%rsp)
  movapd %xmm1,%xmm0
  call   <pow@plt>
  movsd  0x18(%rsp),%xmm2
  subsd  D5(%rip),%xmm2
  movapd %xmm0,%xmm1
  movsd  0x8(%rsp),%xmm0
  movsd  0x10(%rsp),%xmm5
  mulsd  D13(%rip),%xmm1
  mulsd  %xmm2,%xmm0
  movapd %xmm1,%xmm2
  movsd  D79(%rip),%xmm1
  divsd  %xmm5,%xmm0
  subsd  %xmm2,%xmm1
  mulsd  (%rsp),%xmm0
  mulsd  %xmm1,%xmm0
  test   %ebx,%ebx
  je     L90
  mulsd  D146(%rip),%xmm0
  subsd  D8(%rip),%xmm0
L90:
  add    $0x20,%rsp
  pop    %rbx
  ret
  nopl   0x0(%rax)
L35:
  comisd %xmm0,%xmm4
  movsd  D112(%rip),%xmm0
  ja     L46
L59:
  comisd %xmm0,%xmm4
  movsd  D91(%rip),%xmm0
  jbe    L48
  comisd %xmm4,%xmm0
  jae    L16
L48:
  comisd %xmm0,%xmm4
  jbe    L20
  movsd  D78(%rip),%xmm6
  movapd %xmm5,%xmm3
  comisd %xmm4,%xmm6
  jae    L74
L20:
  movapd %xmm5,%xmm3
  divsd  D121(%rip),%xmm3
  jmp    L74
  nop
L46:
  comisd %xmm4,%xmm0
  movapd %xmm5,%xmm3
  jae    L74
  jmp    L59
L16:
  movapd %xmm5,%xmm3
  divsd  D58(%rip),%xmm3
  jmp    L74
```

tool_calls: name: extract_symbolic_data
arguments: data_label: D40
type: f64
name: extract_symbolic_data
arguments: data_label: D78
type: f64[2]
name: extract_symbolic_data
arguments: data_label: D72
type: f64
name: extract_symbolic_data
arguments: data_label: D25
type: f64
name: extract_symbolic_data
arguments: data_label: D9
type: f64
name: extract_symbolic_data
arguments: data_label: D5
type: f64
name: extract_symbolic_data
arguments: data_label: D13
type: f64
name: extract_symbolic_data
arguments: data_label: D79
type: f64
name: extract_symbolic_data
arguments: data_label: D146
type: f64
name: extract_symbolic_data
arguments: data_label: D8
type: f64
name: extract_symbolic_data
arguments: data_label: D112
type: f64
name: extract_symbolic_data
arguments: data_label: D91
type: f64
name: extract_symbolic_data
arguments: data_label: D121
type: byte[8]
name: extract_symbolic_data
arguments: data_label: D58
type: f64

tool: D40: 9

tool: {"D78": [15.0,11.0]}

tool: D72: 1.1

tool: D25: 0.0062

tool: D9: 0.0032

tool: D5: 273.15

tool: D13: 0.978

tool: D79: 1

tool: D146: 1.43

tool: D8: 0.0845

tool: D112: 11

tool: D91: 13

tool: {"D121": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00]}

tool: D58: 0.9

assistant: ```c
double g_0(double bbalb, double ndvi, double tempk, double rnet, double time, int roerink) {
  double a, b, result;
  double r0_coef;
  if (time <= 9.0 || time > 15.0) {
    r0_coef = 1.1;
  } else if (time > 9.0 && time <= 11.0) {
    r0_coef = 1.0;
  } else if (time > 11.0 && time <= 13.0) {
    r0_coef = 0.9;
  } else if (time > 13.0 && time <= 15.0) {
    r0_coef = 1.0;
  }
  a = (0.0032 * (bbalb / r0_coef) + 0.0062 * (bbalb / r0_coef) * (bbalb / r0_coef));
  b = (1 - 0.978 * pow(ndvi, 4));
  result = (rnet * (tempk - 273.15) / bbalb) * a * b;
  if (roerink) {
    result = result * 1.430 - 0.0845;
  }
  return result;
}
```