sedi.c File Reference

Calculate sedimentation velocity. More...

#include "mptrac.h"

Go to the source code of this file.

Functions
int	main (int argc, char *argv[])

Detailed Description

Calculate sedimentation velocity.

Definition in file sedi.c.

Function Documentation

main()

int main	(	int	argc,
		char *	argv[]
	)

Definition at line 27 of file sedi.c.

                {
 
  /* Check arguments... */
  if (argc < 5)
    ERRMSG("Give parameters: <p> <T> <r_p> <rho_p>");
 
  /* Read arguments... */
  double p = atof(argv[1]);
  double T = atof(argv[2]);
  double r_p = atof(argv[3]);
  double rho_p = atof(argv[4]);
 
  /* Calculate sedimentation velocity... */
  double vs = sedi(p, T, r_p, rho_p);
 
  /* Density of dry air [kg / m^3]... */
  double rho = 100. * p / (RA * T);
 
  /* Dynamic viscosity of air [kg / (m s)]... */
  double eta = 1.8325e-5 * (416.16 / (T + 120.)) * pow(T / 296.16, 1.5);
 
  /* Particle Reynolds number... */
  double Re = 2e-6 * r_p * vs * rho / eta;
 
  /* Write output... */
  printf("    p= %g hPa\n", p);
  printf("    T= %g K\n", T);
  printf("  r_p= %g microns\n", r_p);
  printf("rho_p= %g kg/m^3\n", rho_p);
  printf("rho_a= %g kg/m^3\n", RHO(p, T));
  printf("  v_s= %g m/s\n", vs);
  printf("   Re= %g\n", Re);
 
  return EXIT_SUCCESS;
}

Here is the call graph for this function: