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

#include "bvw04x2.c"

int call_ht_adj04x2(argc, argv)
  int argc;
  void *argv[];
{
  int    *dyn_in_prm;
  float  *dyn_in_val, *dyn_in_val2, *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; 
  int    *sst_prm;
  float  *tskin, *z_wind_ref, *z_temp_ref;
  int    *astab; 
  float  *Qnet;
  int    *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;
  int    *eqv_neut;
  float  *z_wanted, *u_at_z, *t_at_z, *q_at_z;
  int    *flux_model;

  float  phi_u;
 
  dyn_in_prm     = (int *)   argv[0];
  dyn_in_val     = (float *) argv[1];
  dyn_in_val2    = (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];
  sst_prm        = (int *)   argv[14];
  tskin          = (float *) argv[15];
  z_wind_ref     = (float *) argv[16];
  z_temp_ref     = (float *) argv[17];
  astab          = (int *)   argv[18]; 
  Qnet           = (float *) argv[19];
  warn           = (int *)   argv[20];
  shf            = (float *) argv[21];  
  lhf            = (float *) argv[22]; 
  tau            = (float *) argv[23]; 
  u_star         = (float *) argv[24]; 
  t_star         = (float *) argv[25]; 
  q_star         = (float *) argv[26]; 
  z_over_L       = (float *) argv[27]; 
  wave_age       = (float *) argv[28]; 
  dom_phs_spd    = (float *) argv[29]; 
  h_sig          = (float *) argv[30]; 
  ww_stab        = (float *) argv[31]; 
  zo_m           = (float *) argv[32];
  eqv_neut       = (int *)   argv[33];
  z_wanted       = (float *) argv[34];
  u_at_z         = (float *) argv[35];
  t_at_z         = (float *) argv[36];
  q_at_z         = (float *) argv[37];
  flux_model     = (int *)   argv[38];

/* 
  printf( "argc = %i\n", argc ); 
  printf( "call_ht_adj03x2.c: %f %f %f %f %f %i %f %i %f %f %i %f %f %i %f %f %f %i %f %i %i %f\n", 
   *dyn_in_val, *dyn-in_val2,
   *CON_P, *conv_crit, *press_sfc, *air_moist_prm, *air_moist_val, 
   *sfc_moist_prm, *sfc_moist_val, *salinity, *ss_prm, *ss_val, *ta, 
   *sst_prm, *tskin, *z_wind_ref, *z_temp_ref, *astab, *Qnet, *warn, *eqv_neut, *z_wanted); 
*/

  bvw_flag = ht_adj_( *dyn_in_prm, *dyn_in_val, *dyn_in_val2, *CON_P, 
              *conv_crit, *press_sfc,
              *air_moist_prm, *air_moist_val, *sfc_moist_prm, *sfc_moist_val,
              *salinity, *ss_prm, *ss_val, *ta, *sst_prm, *tskin, *z_wind_ref,
              *z_temp_ref, *astab, *Qnet, *warn, shf, lhf, tau, u_star, t_star,
              q_star, z_over_L, wave_age, dom_phs_spd, h_sig, ww_stab, zo_m, 
              *eqv_neut, *z_wanted, u_at_z, t_at_z, q_at_z, *flux_model );

/*
  printf( "call_ht_adj04x2.c: %f %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, *u_at_z,
   *t_at_z, *q_at_z, *flux_model );  

  printf( "call_ht_adj04x2.c %i\n", bvw_flag );
*/

  return( bvw_flag );
}