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

```asm
<calc_psi_zeta>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  mov    %rdi,-0x18(%rbp)
  mov    %rsi,-0x20(%rbp)
  mov    %rdx,-0x28(%rbp)
  mov    %rcx,-0x30(%rbp)
  mov    %r8,-0x38(%rbp)
  mov    %r9d,-0x3c(%rbp)
  movss  %xmm0,-0x40(%rbp)
  movss  %xmm1,-0x44(%rbp)
  movl   $0x1,-0x8(%rbp)
  jmp    L75
L135:
  movl   $0x1,-0x4(%rbp)
  jmp    L42
L118:
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4(%rbp),%eax
  cltq
  add    %rdx,%rax
  movzbl (%rax),%eax
  movsbl %al,%eax
  and    $0x10,%eax
  test   %eax,%eax
  je     L129
  mov    -0x8(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4(%rbp),%eax
  cltq
  add    %rdx,%rax
  movzbl (%rax),%eax
  movsbl %al,%eax
  and    $0x10,%eax
  test   %eax,%eax
  je     L129
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4(%rbp),%eax
  cltq
  add    $0x1,%rax
  add    %rdx,%rax
  movzbl (%rax),%eax
  movsbl %al,%eax
  and    $0x10,%eax
  test   %eax,%eax
  je     L129
  mov    -0x8(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4(%rbp),%eax
  cltq
  add    $0x1,%rax
  add    %rdx,%rax
  movzbl (%rax),%eax
  movsbl %al,%eax
  and    $0x10,%eax
  test   %eax,%eax
  je     L129
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4(%rbp),%eax
  cltq
  add    $0x1,%rax
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  subss  %xmm1,%xmm0
  divss  -0x44(%rbp),%xmm0
  mov    -0x8(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm2
  subss  %xmm2,%xmm1
  divss  -0x40(%rbp),%xmm1
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  subss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  jmp    L93
L129:
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
L93:
  addl   $0x1,-0x4(%rbp)
L42:
  mov    0x10(%rbp),%eax
  cmp    -0x4(%rbp),%eax
  jg     L118
  addl   $0x1,-0x8(%rbp)
L75:
  mov    -0x3c(%rbp),%eax
  cmp    -0x8(%rbp),%eax
  jg     L135
  movl   $0x0,-0x8(%rbp)
  jmp    L145
L27:
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  movl   $0x1,-0x4(%rbp)
  jmp    L49
L74:
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4(%rbp),%eax
  cltq
  shl    $0x2,%rax
  sub    $0x4,%rax
  add    %rax,%rdx
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x28(%rbp),%rax
  add    %rcx,%rax
  mov    (%rax),%rcx
  mov    -0x4(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rcx,%rax
  movss  (%rdx),%xmm0
  movss  %xmm0,(%rax)
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4(%rbp),%eax
  cltq
  add    %rdx,%rax
  movzbl (%rax),%eax
  movsbl %al,%eax
  and    $0x10,%eax
  test   %eax,%eax
  jne    L10
  mov    -0x8(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4(%rbp),%eax
  cltq
  add    %rdx,%rax
  movzbl (%rax),%eax
  movsbl %al,%eax
  and    $0x10,%eax
  test   %eax,%eax
  je     L19
L10:
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mulss  -0x44(%rbp),%xmm0
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  addss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
L19:
  addl   $0x1,-0x4(%rbp)
L49:
  mov    -0x4(%rbp),%eax
  cmp    0x10(%rbp),%eax
  jle    L74
  addl   $0x1,-0x8(%rbp)
L145:
  mov    -0x8(%rbp),%eax
  cmp    -0x3c(%rbp),%eax
  jle    L27
  nop
  nop
  pop    %rbp
  ret
```

assistant: ```c
void calc_psi_zeta(float **u, float **v, float **psi, float **zeta, char **flag, int imax, int jmax, float delx, float dely) {
  int i, j;
  for (i = 1; i <= imax - 1; i++) {
    for (j = 1; j <= jmax - 1; j++) {
      if ((flag[i][j] & 0x0010) && (flag[i + 1][j] & 0x0010) && (flag[i][j + 1] & 0x0010) && (flag[i + 1][j + 1] & 0x0010)) {
        zeta[i][j] = (u[i][j + 1] - u[i][j]) / dely - (v[i + 1][j] - v[i][j]) / delx;
      } else {
        zeta[i][j] = 0.0;
      }
    }
  }
  for (i = 0; i <= imax; i++) {
    psi[i][0] = 0.0;
    for (j = 1; j <= jmax; j++) {
      psi[i][j] = psi[i][j - 1];
      if ((flag[i][j] & 0x0010) || (flag[i + 1][j] & 0x0010)) {
        psi[i][j] += u[i][j] * dely;
      }
    }
  }
}
```