Physical Climatology Programming Exercise 1

You are required to do only two of the following four questions. Feel free to use any programming language. Please turn in your results or outputs plus the code for each question.

1. Modify my sample FORTRAN code (austin_solar2.pdf) that computes daily mean top-of-the-atmosphere (TOA) solar radiation. My code uses the hourly interval. Please modify the code to use a series of different time intervals and then find at what interval the discrete solution is within 0.01 Wm-2 of the analytical solution. You may want to program the code so that it searches the desired solution automatically.

2. Modify the above code to calculate the length of sunlight during the day for the first day of each month from January through December given an input of city name and its latitude. For example, the daylight length at 30ºN (Austin) on October 1 (DOY=274) is 11 hours and 41 minutes. You need to have your output similar to the content in DayLengthAustin.pdf

3.Modify the code in #2 to calculate the length of sunlight during the day for the first day of each season (spring, summer, fall, and winter). You need to have your output similar to the content in DayLength_LatSeason.pdf

4. Modify the code (austin_solar5.pdf) to compute the daily average insolation at the top of the atmosphere as a function of latitude and season (i.e. the mid-month of each month from January to December). Your output should have a format similar to that in TOA_Insolation_LatMon.pdf. Compare your results with Figure 2.6 in Hartmann.

[Useful links:

FORTRAN 77 tutorial

Websites with overview of unix and vi commands (Courtesy of Lindsey Gulden):

]