/*
  To compile c function:
  Example #1
  cc -c -KPIC call_bvw99.c
  ld -shared -o call_bvw99.so call_bvw99.o
  example #2
  cc -I/usr/local/rsi/idl/external -c -O -32 -KPIC call_bvw99.c
  ld -shared -32 -o call_bvw99.so call_bvw99.o
*/

#include "bvw99.h"
#include "bvw99.c"

float call_bvw99(argc, argv)
  int argc;
  void *argv[];
{
  int    *dyn_in_prm;
  float  *dyn_in_val, *rel_wind_ang, *CON_P, *conv_crit, *press_sfc;
  int    *air_moist_prm;
  float  *air_moist_val;
  int    *sfc_moist_prm;
  float  *sfc_moist_val, *salinity;
  int    *ss_prm;
  float  *ss_val, *ta, *tskin, *z_wind_ref, *z_temp_ref;
  int    *astab, *warn;
  int    bvw_flag;
  float  *dom_phs_spd, *h_sig, *lhf, *q_star, *shf, *tau, *t_star,
         *u_star, *wave_age, *ww_stab, *z_over_L, *zo_m;
 
  dyn_in_prm     = (int *)   argv[0];
  dyn_in_val     = (float *) argv[1];
  rel_wind_ang   = (float *) argv[2]; 
  CON_P          = (float *) argv[3];
  conv_crit      = (float *) argv[4]; 
  press_sfc      = (float *) argv[5];
  air_moist_prm  = (int *)   argv[6];
  air_moist_val  = (float *) argv[7]; 
  sfc_moist_prm  = (int *)   argv[8];
  sfc_moist_val  = (float *) argv[9];
  salinity       = (float *) argv[10];
  ss_prm         = (int *)   argv[11];
  ss_val         = (float *) argv[12];
  ta             = (float *) argv[13];
  tskin          = (float *) argv[14];
  z_wind_ref     = (float *) argv[15];
  z_temp_ref     = (float *) argv[16];
  astab          = (int *)   argv[17]; 
  warn           = (int *)   argv[18];
  shf            = (float *) argv[19];  
  lhf            = (float *) argv[20]; 
  tau            = (float *) argv[21]; 
  u_star         = (float *) argv[22]; 
  t_star         = (float *) argv[23]; 
  q_star         = (float *) argv[24]; 
  z_over_L       = (float *) argv[25]; 
  wave_age       = (float *) argv[26]; 
  dom_phs_spd    = (float *) argv[27]; 
  h_sig          = (float *) argv[28]; 
  ww_stab        = (float *) argv[29]; 
  zo_m           = (float *) argv[30];

/* 
  printf( "argc = %i\n", argc ); 
  printf( "call_bvw99.c: %i %f %f %f %f %f %i %f %i %f %f %i %f %f %f %f %f %i %i\n", 
   *dyn_in_prm, *dyn_in_val, *rel_wind_ang,
   *CON_P, *conv_crit, *press_sfc, *air_moist_prm, *air_moist_val, 
   *sfc_moist_prm, *sfc_moist_val, *salinity, *ss_prm, *ss_val, *ta, 
   *tskin, *z_wind_ref, *z_temp_ref, *astab, *warn ); 
*/ 

  bvw_flag = pmix_(*dyn_in_prm, *dyn_in_val, *rel_wind_ang, *CON_P, 
   *conv_crit, *press_sfc, *air_moist_prm, *air_moist_val, 
   *sfc_moist_prm, *sfc_moist_val, *salinity, *ss_prm, *ss_val, *ta, 
   *tskin, *z_wind_ref, *z_temp_ref, *astab, *warn, shf, lhf, tau, 
   u_star, t_star, q_star, z_over_L, wave_age, dom_phs_spd, h_sig, 
   ww_stab, zo_m);
/*
  printf( "call_bvw99.c: %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f\n",
   *shf, *lhf, tau[0], tau[1], u_star[0], u_star[1], *t_star, *q_star, 
   *z_over_L, *wave_age, *dom_phs_spd, *h_sig, *ww_stab, zo_m[0], zo_m[1] );  
*/
  return( bvw_flag );
}