user: Translate the below assembly code into C programming language:

```asm
<solar_azimuth>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x48,%rsp
  mov    %edi,%r12d
  mov    %esi,%ebx
  movsd  %xmm0,0x8(%rsp)
  movsd  %xmm1,0x10(%rsp)
  movsd  %xmm2,0x18(%rsp)
  lea    0x38(%rsp),%r14
  lea    0x30(%rsp),%r15
  movapd %xmm0,%xmm3
  mulsd  D72(%rip),%xmm3
  movapd %xmm3,%xmm0
  divsd  D13(%rip),%xmm0
  mov    %r15,%rsi
  mov    %r14,%rdi
  call   <sincos@plt>
  movsd  0x30(%rsp),%xmm4
  movsd  %xmm4,0x28(%rsp)
  mov    0x38(%rsp),%r13
  movsd  0x18(%rsp),%xmm1
  movsd  0x10(%rsp),%xmm0
  mov    %ebx,%esi
  mov    %r12d,%edi
  call   <LAT@plt>
  mov    %eax,%ebp
  mov    %ebx,%edi
  call   <declination@plt>
  movsd  %xmm0,0x20(%rsp)
  mov    %ebp,%edi
  call   <hour_angle@plt>
  movsd  0x18(%rsp),%xmm2
  movsd  0x10(%rsp),%xmm1
  movsd  0x8(%rsp),%xmm0
  mov    %ebx,%esi
  mov    %r12d,%edi
  call   <solar_zenith@plt>
  mulsd  D107(%rip),%xmm0
  mov    %r15,%rsi
  mov    %r14,%rdi
  call   <sincos@plt>
  movsd  D107(%rip),%xmm0
  mulsd  0x20(%rsp),%xmm0
  call   <sin@plt>
  movapd %xmm0,%xmm1
  movq   %r13,%xmm0
  mulsd  0x30(%rsp),%xmm0
  subsd  %xmm1,%xmm0
  divsd  0x38(%rsp),%xmm0
  divsd  0x28(%rsp),%xmm0
  comisd D90(%rip),%xmm0
  jbe    L53
  movsd  D45(%rip),%xmm1
  comisd %xmm0,%xmm1
  ja     L24
L53:
  movsd  D68(%rip),%xmm1
  comisd %xmm0,%xmm1
  jbe    L3
  comisd D9(%rip),%xmm0
  jbe    L3
  movapd %xmm1,%xmm0
  jmp    L3
L24:
  movsd  D90(%rip),%xmm0
L3:
  call   <acos@plt>
  cmp    $0xa8c0,%ebp
  ja     L9
  xorpd  D66(%rip),%xmm0
L9:
  mulsd  D63(%rip),%xmm0
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: disassemble_data
arguments: label: D72
type_guessed: f64
name: disassemble_data
arguments: label: D13
type_guessed: f64
name: disassemble_data
arguments: label: D107
type_guessed: byte[8]
name: disassemble_data
arguments: label: D90
type_guessed: f64[2]
name: disassemble_data
arguments: label: D45
type_guessed: byte[8]
name: disassemble_data
arguments: label: D68
type_guessed: f32[2]
name: disassemble_data
arguments: label: D9
type_guessed: byte[16]
name: disassemble_data
arguments: label: D66
type_guessed: f32[4]
name: disassemble_data
arguments: label: D63
type_guessed: byte[8]

tool: D72: 3.1415926

tool: D13: 180

tool: {"D107": [0x18,0x39,0x35,0x9d,0x46,0xdf,0x91,0x3f]}

tool: {"D90": [1.0,3.1415926]}

tool: {"D45": [0x0b,0x7a,0x6f,0x0c,0x01,0x00,0xf0,0x3f]}

tool: {"D68": [0.0,-1.875]}

tool: {"D9": [0x0b,0x7a,0x6f,0x0c,0x01,0x00,0xf0,0xbf,0xaa,0x9f,0x96,0x22,0xdc,0xa5,0x4c,0x40]}

tool: {"D66": [0.0,-0.0,0.0,0.0]}

tool: {"D63": [0xaa,0x9f,0x96,0x22,0xdc,0xa5,0x4c,0x40]}

assistant: ```c
double solar_azimuth(int time, int day, double latitude, double longitude, double long_std) {
  double theta = 0, omega = 0, delta = 0, phi = latitude * 3.1415926 / 180.0, psi = 0, cospsi = 0;
  int lat = LAT(time, day, longitude, long_std);
  delta = 3.1415926 / 180.0 * declination(day);
  omega = 3.1415926 / 180.0 * hour_angle(lat);
  theta = 3.1415926 / 180.0 * solar_zenith(time, day, latitude, longitude, long_std);
  cospsi = (cos(theta) * sin(phi) - sin(delta)) / sin(theta) / cos(phi);
  if (cospsi > 1.0 && cospsi < 1.0 + 1e-6)
    cospsi = 1.0;
  if (cospsi < -1.0 && cospsi > -1.0 - 1e-6)
    cospsi = -1.0;
  psi = -acos(cospsi);
  if (lat > 43200 || lat < 0)
    psi = -psi;
  psi *= (180.0 / 3.1415926);
  return psi;
}
```