Opened 17 months ago

Closed 17 months ago

Last modified 17 months ago

#809 closed bug (Fixed)

IMPORTANT: Corrected bug in ana_srflux.h

Reported by: arango Owned by:
Priority: major Milestone: Release ROMS/TOMS 3.7
Component: Nonlinear Version: 3.7
Keywords: Cc:

Description

Finally, I have time to browse the literature, Googling, and write a Matlab script to check a bug in ana_srflx.h as discussed in the Forum Indeed, the 15.0 factor used to compute the zenith is incorrect. We need to have instead:

        zenith=cff1+cff2*COS(Hangle-lonr(i,j)*deg2rad)

In my Matlab script, I compared against a WRF-ROMS coupled solution for the US West Coast at a particular point in the grid:

http://www.myroms.org/trac/swrad.png

Here, the X-axis is the time since initialization and it is in GMT time zone. We need to subtract 8 hours for local Pacific Standard Time (PST) to analyze the results. As you can see, the WRF-ROMS coupled solution values agree well with the red curve that computes zenith with updated above equation. I am now satisfied that I proved to myself that it is a bug.

Many thanks to all that brought this bug to my attention.

Change History (2)

comment:1 Changed 17 months ago by arango

  • Resolution set to Fixed
  • Status changed from new to closed

comment:2 Changed 17 months ago by arango

Forgot to correct how Dangle is computed. We need to have:

      Dangle=23.44_dp*COS((172.0_dp-yday)*2.0_dp*pi/365.2425_dp)

The length of the year is 648 seconds less to a value of 365.2425 days.

Note: See TracTickets for help on using tickets.