Using rain as LwSrc

Suggest improvements/optimizations to the ROMS code.

Moderators: arango, robertson

Post Reply
Message
Author
jcwarner
Posts: 1172
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Using rain as LwSrc

#1 Unread post by jcwarner »

We are using ROMS to simulate rainfall and urban flooding during hurricanes. Currently the 'rain' input from the surface forcing is used for several physics such as a BC for salinity. This will allow the water to get fresher. But the volume of that rain is not accounted for. The LwSrc point source allows for the inflow to be incorporated into the model Volume, but it would be cumbersome to have Lm x Mm point sources.
So can we modify the code, maybe something like this in Omega:

IF (LwSrc(ng)) THEN

#ifdef RAIN_FLUX (or some other name)
DO ii=Istr,Iend
DO jj=Jstr,Jend
#else
DO is=1,Nsrc(ng)
IF (INT(SOURCES(ng)%Dsrc(is)).eq.2) THEN
ii=SOURCES(ng)%Isrc(is)
jj=SOURCES(ng)%Jsrc(is)
#endif

....

#ifdef RAIN_FLUX (or some other name)
& FORCES(ng)%rain(ii,jj)/N(ng) !not sure the units here, *om_r*on_r?
#else
& SOURCES(ng)%Qsrc(is,k)
#endif


this might be a simple way to use rain as a source over the whole grid.
-j

Post Reply