> TURB_NCAR, transfer coefficients with NCAR 3.0/3.5

Use turb_ncar() of module mod_blk_ncar (mod_blk_ncar.f90).
Example of a call:

          PROGRAM TEST_COEFF
              USE mod_const
              USE mod_blk_ncar
              ...
              jpi = Ni ! x-shape of the 2D domain
              jpj = Nj ! y-shape of the 2D domain
              ...
              CALL TURB_NCAR( zt, zu, sst, t_zt, ssq, q_zt, U_zu, &
              &               Cd, Ch, Ce, t_zu, q_zu, U_blk       &
              &               [ , xz0=z0, xu_star=u_s, xL=L ] )
              ...
          END PROGRAM TEST_COEFF
      

INPUT ARGUMENTS:

• zt : (Sc,real) height for air temperature and humidity [m]
• zu : (Sc,real) height for wind speed (generally 10m) [m]
• sst : (2D,real) SST [K]
• t_zt : (2D,real) potential air temperature at zt [K]
• ssq : (2D,real) surface saturation spec. humidity [kg/kg]
• q_zt : (2D,real) air spec. humidity of at zt [kg/kg]
• U_zu : (2D,real) scalar wind speed at zu [m/s]

OUTPUT ARGUMENTS:

• Cd : (2D,real) drag coefficient
• Ch : (2D,real) sensible heat transfer coefficient
• Ce : (2D,real) moisture transfer (evaporation) coefficient
• t_zu : (2D,real) air pot. temperature adjusted at zu [K]
• q_zu : (2D,real) air spec. humidity adjusted at zu [kg/kg]
• Ublk : (2D,real) bulk wind speed at 10m [m/s]

[ OPTIONAL OUTPUT ARGUMENTS: ]

• z0 : (2D,real) roughness length of the sea surface [m]
• u_s : (2D,real) friction velocity [m/s]
• L : (2D,real) Monin-Obukhov length [m]