A_Kowaleski
New member
Hello,
I am trying to modify module_sf_sfclayrev.F to eliminate surface fluxes based on distance from the tropical cyclone center. To find the TC center, I seek to use the minimum pressure in the domain (simple, but it should work for stronger storms, which I am simulating). I also want to find the location of maximum 10-meter wind speed, so that I can eliminate fluxes as a function of the radius of maximum winds.
How do I go about computing the min pressure/max wind over the entire domain? Although my modified module_sf_sfclayrev.F compiles correctly, I think that I am currently having difficulty because it computes min pressure/max wind over tiles, rather than the entire domain. For example, within SFCLAYREV:
PMIN_LOC
)=MINLOC(PSFC)
WSPD_MAX=MAXVAL(WSPD10)
Is computing the min/max over the entire domain possible within module_sf_sfclayrev.F? If so, how? If not, what do you suggest?
Thank you
I am trying to modify module_sf_sfclayrev.F to eliminate surface fluxes based on distance from the tropical cyclone center. To find the TC center, I seek to use the minimum pressure in the domain (simple, but it should work for stronger storms, which I am simulating). I also want to find the location of maximum 10-meter wind speed, so that I can eliminate fluxes as a function of the radius of maximum winds.
How do I go about computing the min pressure/max wind over the entire domain? Although my modified module_sf_sfclayrev.F compiles correctly, I think that I am currently having difficulty because it computes min pressure/max wind over tiles, rather than the entire domain. For example, within SFCLAYREV:
PMIN_LOC
WSPD_MAX=MAXVAL(WSPD10)
Is computing the min/max over the entire domain possible within module_sf_sfclayrev.F? If so, how? If not, what do you suggest?
Thank you