Scale factor in extract*.F (roms-2.1 and 2.2)

Bug reports, work arounds and fixes

Moderators: arango, robertson

Post Reply
Message
Author
User avatar
cvl
Posts: 18
Joined: Tue Jun 03, 2003 7:39 pm

Scale factor in extract*.F (roms-2.1 and 2.2)

#1 Unread post by cvl »

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.

User avatar
arango
Site Admin
Posts: 1351
Joined: Wed Feb 26, 2003 4:41 pm
Location: DMCS, Rutgers University
Contact:

#2 Unread post by arango »

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 :)

Post Reply