01.12.2012 Views

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 ...

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!