Merge pull request #9 from natgeo-wong/main

Debugging fixed wind speed flux specification

Author: natgeo-wong

surface.f90

@@ -44,10 +44,10 @@ subroutine surface()
             q_s = salt_factor*qsatw(sstxy(i,j)+t00,presi(1))
 
             if (dobulksfcflx) then
-              call oceflx(rho(1), u_h, v_h, ta_h, q_h, t_h, z(1)-zi(1), t_s, q_s, &
+              call oceflx(rho(1), bulksfcflx_u, 0, ta_h, q_h, t_h, z(1)-zi(1), t_s, q_s, &
                 fluxt0, fluxq0, taux0, tauy0)
             else
-              call oceflx(rho(1), bulksfcflx_u, 0, ta_h, q_h, t_h, z(1)-zi(1), t_s, q_s, &
+              call oceflx(rho(1), u_h, v_h, ta_h, q_h, t_h, z(1)-zi(1), t_s, q_s, &
                 fluxt0, fluxq0, taux0, tauy0)
             endif