Alternative small scale meteorology input to a chemical transport ...

C met. given simplified inputs such as USTAR, M-O Length, PBL height
C
C INPUTS:
C >
C ISESN ! season index for table
C SDATE ! starting date, format YYYYDDD
C STIME ! starting time, format HHMMSS
C ISESN ! season code
C GMT !G.M. TIME in hours (e.g., 10.75 = 10:45 GMT)
C
C >
C NDEP
C INPUTS:
C >
C XLUSE (NCOLS_X,NROWS_X,LURADM)!landuse fractions
C X3HTF (NCOLS_X,NROWS_X,0:MAXK)!height at model layerface
C X3HTM (NCOLS_X,NROWS_X, MAXK) ! height at layer middle
C XTEMPM (NCOLS_X,NROWS_X, MAXK) ! temperature profiles
C XDENSAM (NCOLS_X,NROWS_X, MAXK) ! air density (kg/m^3)
C XDENSWM (NCOLS_X,NROWS_X, MAXK) ! vapor density (kg/m^3)
C XTEMPG (NCOLS_X,NROWS_X) ! temperature at ground
C XRAINN (NCOLS_X,NROWS_X) ! accumulated nonc. rain
C XRAINC (NCOLS_X,NROWS_X) ! accumulated conv. rain
C XUSTAR (NCOLS_X,NROWS_X,0) ! ustar (cell avg) (m/s)
C XMOL (NCOLS_X,NROWS_X) ! MONIN-OBUKHOV length
C XPBL (NCOLS_X,NROWS_X) ! PBL height (m)
C
C OUTPUT:
C Compute following
C XUSTAR (NCOLS_X,NROWS_X,0:LURADM)! 0 = avg,
C others = LUSE specific
C XRIB (NCOLS_X,NROWS_X)
C XHFX (NCOLS_X,NROWS_X)
C XQFX (NCOLS_X,NROWS_X)
C XWSTAR (NCOLS_X,NROWS_X,0) ! ustar (cell avg) (m/s)
C XRADYN (NCOLS_X,NROWS_X,0:LURADM)! 0 = avg,
C others = LUSE specific
C XRBNDY (NCOLS_X,NROWS_X,0:LURADM) ! 0 = avg,
C others = LUSE specific
C XRSTOM (NCOLS_X,NROWS_X,0:LURADM)
C XTEMP10 (NCOLS_X,NROWS_X)
C useful to compute following if possible
C ??? XSURF1 (NCOLS_X,NROWS_X)
C ??? XSURF2 (NCOLS_X,NROWS_X)
C ??? XFSOIL (NCOLS_X,NROWS_X)
C Now, from solar.F routine XALBEDO (NCOLS_X,NROWS_X)
C Now, from solar.F routine XGLW (NCOLS_X,NROWS_X)
C Now, from solar.F routine XGSW (NCOLS_X,NROWS_X)
C Now, from solar.F routine XRGRND (NCOLS_X,NROWS_X)
C Now, from solar.F routine XRNET (NCOLS_X,NROWS_X)
C ISTAT ! error status for sanity check
98