00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040 #include "plumeinput.h"
00041 #include "plumevars.h"
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051
00052 int plumejump( Plume_river *river )
00053 {
00054 double froude, d1, u1;
00055
00056 froude = river->u0/sqrt(grv*river->d0);
00057
00058 if( froude > 1.0 )
00059 {
00060 d1 = -(river->d0/2)
00061 + sqrt( sq(river->d0/2) + ((2*river->d0*sq(river->u0))/grv) );
00062 u1 = river->u0*river->d0/d1;
00063
00064 fprintf(stderr, " PlumeJump: Hydraulic Jump detected and fixed \n");
00065 fprintf(fidlog, " PlumeJump: Hydraulic Jump detected and fixed \n");
00066 fprintf(fidlog, " Fr = %g \n", froude);
00067 fprintf(fidlog, " Uold = %g \n", river->u0);
00068 fprintf(fidlog, " Unew = %g \n", u1);
00069 fprintf(fidlog, " Dold = %g \n", river->d0);
00070 fprintf(fidlog, " Dnew = %g \n\n", d1);
00071
00072 river->u0 = u1;
00073 river->d0 = d1;
00074 }
00075
00076 #ifdef DBG
00077 fprintf(fidlog, " PlumeJump: \n");
00078 fprintf(fidlog, " Fr = %g \n", froude);
00079 fprintf(fidlog, " U = %g \n", river->u0);
00080 fprintf(fidlog, " D = %g \n", river->d0);
00081 #endif
00082
00083 return 0;
00084
00085 }
00086