Alternative small scale meteorology input to a chemical transport ...
Alternative small scale meteorology input to a chemical transport ...
Alternative small scale meteorology input to a chemical transport ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
DO 10 I=1,NWORDS<br />
XARR(I)=X<br />
10 CONTINUE<br />
C<br />
RETURN<br />
END<br />
c----------------------------------------------------------------------<br />
FUNCTION JULDAY(iyear,imonth,iiday)<br />
c----------------------------------------------------------------------<br />
c --- PURPOSE: Compute the Julian day number from the Gregorian date<br />
c (month, day)<br />
c<br />
c --- INPUTS:<br />
c IYEAR - integer - Year (YYYY)<br />
c IMONTH - integer - Month (MM)<br />
c IDAY - integer - Day (DD)<br />
c<br />
c --- OUTPUT:<br />
c JULDAY - integer - Julian day (DDD)<br />
c<br />
c --- JULDAY called by: MAIN<br />
c --- JULDAY calls: none<br />
c---------------------------------------------------------------------<br />
integer iyear ! Year (YYYY)<br />
integer imonth ! Month (MM)<br />
integer iiday ! Day (DD)<br />
c<br />
integer kday(12)<br />
data kday/0,31,59,90,120,151,181,212,243,273,304,334/<br />
if(imonth.lt.1.or.imonth.gt.12.or.iiday.gt.31)then<br />
write(*,10)imonth<br />
10 format(//2x,'ERROR in SUBR. JULDAY -- invalid date - imonth =<br />
',<br />
1 i8,5x,'iiday = ',i8)<br />
s<strong>to</strong>p<br />
endif<br />
c<br />
JULDAY=kday(imonth)+iiday<br />
if(imonth.le.2) return<br />
if(mod(iyear,4).EQ.0) JULDAY=JULDAY+1<br />
c<br />
return<br />
end function JULDAY<br />
90