Difference between revisions of "River Runoff"
| Line 19: | Line 19: | ||
In ana_psource.h, set the values for | In ana_psource.h, set the values for | ||
Nsrc: Number of sources | Nsrc: Number of sources. | ||
Dsrc(Nsrc): direction. (0 for along xi-axis, 1 for eta-axis) | Dsrc(Nsrc): direction. (0 for along xi-axis, 1 for eta-axis). | ||
Isrc(Nsrc): location of i index (u or v point) | Isrc(Nsrc): location of i index (u or v point) | ||
Jsrc(Nsrc): location of j index (u or v point) | Jsrc(Nsrc): location of j index (u or v point) | ||
Revision as of 02:29, 9 April 2008
River Input
First you need to set the appropriate CPP options. Typically, this would be:
TS_PSOURCE for salinity and/or temperature sources/sinks
and/or
UV_PSOURCE for mass sources/sinks due to river runoff.
Then you need to decide how to provide the information of river flow, vertical distribution of flow, direction, location, and tracer (temp, salt, sed, etc) values. You can provide the data analytically or by providing the data in a NetCDF forcing file.
To provide this information analytically, you need to activate ANA_PSOURCE. In ana_psource.h, set the values for
Nsrc: Number of sources. Dsrc(Nsrc): direction. (0 for along xi-axis, 1 for eta-axis). Isrc(Nsrc): location of i index (u or v point) Jsrc(Nsrc): location of j index (u or v point) Lsrc(Nsrc, itemp): T or F if temp information will be provided. If not, then the model uses the value at the inflowing rho point. Lsrc(Nsrc, isalt): T or F if salt information will be provided. If not, then the model uses the value at the inflowing rho point. Lsrc(Nsrc, ised): T or F if sediment information will be provided. If not, then the model uses the value at the inflowing rho point. Qshape(Nsrc,N): vertical distribution of the flow. Qbar(Nsrc): depth-integrated mass flow. Tsrc(Nsrc,N,itemp): actual vertical distribution of the temperature values. Tsrc(Nsrc,N,isalt): actual vertical distribution of the salinity values. Tsrc(Nsrc,N,ised): actual vertical distribution of the sediment values.
To provide this information in a netcdf file, do NOT define the ANA_PSOURCE. Instead, create the netcdf file and provide the file name in ocean.in. You can find examples on how to write such a river forcing file in the "matlab" package at https://www.myroms.org/index.php?page=Processing you may also be interested in this https://www.myroms.org/forum/viewtopic.php?t=63
Here are some forum discussions about the vertical profile of the flow. https://www.myroms.org/forum/viewtopic.php?t=60 and https://www.myroms.org/forum/viewtopic.php?t=79
