for(int i = 0; i < 576 ; i++)//The input data is complete source data

{

VpwdPrjProcess()

{

Thread_Count=strtol(argv[1],MULL,10);//Get the number of threads

for（thread=0；thread

{

//Create a thread, the inversion calculation into Thread_Count sub-tasks, by Thread_Count threads in parallel

pthread_create(VpwdProcess);

}

for（thread=0；thread

{

//Blocking the call thread, until all the implementation of the end of the sub-thread, the main thread to perform the remaining operations.

pthread_join(thread_handles[thread]);

}

……Other non-parallel operations

}

}

WriteFile();//Generate precipitation products

VpwdDelete();//Release memory