Not sure if this has been reported before, but the scale factor passed from write_sta.F to extract.F (roms-2.1), or extract_sta.F (roms-2.2), does not appear to be applied. In my application, history and station files contain predictably different values (O(10^3)) for bottom stress in the same approximate location.
Presumably, the fix is best applied in extract*.F. I believe a correction is:
Apos(np)=Ascl*Apos(np)
On lines 165, 219, 273, 474, 589, 704, and 811
Though this seems a bit redundant. I welcome alternative input as to the proper correction.
Scale factor in extract*.F (roms-2.1 and 2.2)
-
arango
- Site Admin
- Posts: 1367
- Joined: Wed Feb 26, 2003 4:41 pm
- Location: DMCS, Rutgers University
- Contact:
Yes, good catch 
We need to multiply by the scale factor Ascl in the interpolation routines extract_sta2d and extrac_sta3d . Fortunately, this only affected the scaling of surface heat fluxes and surface and bottom momentum stresses since this scale is unity for other output station data. The corrected routine can be downloaded from the corrections tar.
Thank you
We need to multiply by the scale factor Ascl in the interpolation routines extract_sta2d and extrac_sta3d . Fortunately, this only affected the scaling of surface heat fluxes and surface and bottom momentum stresses since this scale is unity for other output station data. The corrected routine can be downloaded from the corrections tar.Thank you