sedi.c File Reference

Calculate sedimentation velocity. More...

#include "mptrac.h"

Go to the source code of this file.

Functions
void	usage (void)
	Print command-line help. More...
int	main (int argc, char *argv[])

Detailed Description

Calculate sedimentation velocity.

Definition in file sedi.c.

Function Documentation

usage()

void usage

(

void

)

Print command-line help.

Definition at line 84 of file sedi.c.

        {
 
  printf("\nMPTRAC sedi tool.\n\n");
  printf
    ("Calculate sedimentation velocity for aerosol or cloud particles.\n");
  printf("\n");
  printf("Usage:\n");
  printf("  sedi <p> <T> <r_p> <rho_p>\n");
  printf("\n");
  printf("Arguments:\n");
  printf("  <p>      Pressure [hPa].\n");
  printf("  <T>      Temperature [K].\n");
  printf("  <r_p>    Particle radius [um].\n");
  printf("  <rho_p>  Particle density [kg/m3].\n");
  printf("\nFurther information:\n");
  printf("  Manual: https://slcs-jsc.github.io/mptrac/\n");
}

main()

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

Definition at line 39 of file sedi.c.

                {
 
  /* Print usage information... */
  USAGE;
 
  /* Check arguments... */
  if (argc < 5)
    ERRMSG("Missing or invalid command-line arguments.\n\n"
           "Usage: sedi <p> <T> <r_p> <rho_p>\n\n" "Use -h for full help.");
 
  /* Read arguments... */
  const double p = atof(argv[1]);
  const double T = atof(argv[2]);
  const double r_p = atof(argv[3]);
  const double rho_p = atof(argv[4]);
 
  /* Calculate sedimentation velocity... */
  const double vs = sedi(p, T, r_p, rho_p);
 
  /* Density of dry air [kg / m^3]... */
  const double rho = 100. * p / (RA * T);
 
  /* Dynamic viscosity of air [kg / (m s)]... */
  const double eta = 1.8325e-5 * (416.16 / (T + 120.)) * pow(T / 296.16, 1.5);
 
  /* Particle Reynolds number... */
  const 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: