!
! ProTeX sample FORTRAN source: Style 2
!
! NOTE: Styles 1 and 2 differ only on the order of the !DESCRIPTION
! keyword.
!
!-------------------------------------------------------------------------
! NASA/GSFC, Data Assimilation Office, Code 910.3, GEOS/DAS !
!-------------------------------------------------------------------------
!BOP
!
! !ROUTINE: ATESS_RunMeans() --- Update current estimates of O-F grids
!
! !DESCRIPTION: This routine reads previous estimates of O-F grids from
! disk and updates them according to the new estimates
! given on input, which are generally computed with a much smaller sample
! (1 or 2 days, say). Denote by $\overline{\delta}_{old}$ the current
! estimate of the time-mean innovation, and by $\delta$ a small sample
! estimate (sample size $N_{syn}$). The updated time-mean innovation
! $\overline{\delta}_{new}$ is given by
! $$
! \overline{\delta}_{new} := (1-\eta)\overline{\delta}_{old}
! + \eta \delta
! $$
! where $\eta = N_{syn}/N_* $, and $N_*$ is the number of synoptic
! times in the ``running-mean'' window. Typically, $N_* = 30\times4=120$
! corresponding to a 1 month long ``running-mean''. The parameter
! $N_*$ is determined from the resource file.
!
! !INTERFACE:
!
subroutine ATESS_RunMeans ( delgrid_ifname,
. bias_raob, rtms_raob, nobs_raob,
. bias_tovs, rtms_tovs, nobs_tovs,
. nsyn_raob, nsyn_tovs )
! !USES:
Use config ! System wide configuration
Use m_inpak ! Resource files
! !INPUT PARAMETERS:
!
Implicit NONE
include 'atess_dim.h' ! Dimension of grids
character*(255) delgrid_ifname ! Name of the input file with
! current estimate of the
! O-F grids. This is a GrADS
! control (ctl) file.
! Radiosonde height O-F grids:
real bias_raob(idm,jdm,kdm,ldm) ! Time mean
real rtms_raob(idm,jdm,kdm,ldm) ! Root time mean square
real nobs_raob(idm,jdm,kdm,ldm) ! No. obs per gridbox
! NESDIS TOVS A/B height O-F Grids
real bias_tovs(idm,jdm,kdm,ldm) ! Time mean
real rtms_tovs(idm,jdm,kdm,ldm) ! Root time mean square
real nobs_tovs(idm,jdm,kdm,ldm) ! No. obs per gridbox
integer nsyn_raob ! Number of synoptic times used
! in the computation of
! radiosonde grids.
integer nsyn_tovs ! Number of synoptic times used
! in the computation of
! TOVS grids.
! NOTE: The last dimension of the grids above (ldm) specify whether or
! not radiosonde-TOVS co-location has ben used to produce
! the grids:
!
! ldm = 1: all data (non-colocated)
! ldm = 2: co-located RAOB/TOVS observations
!
! !OUTPUT PARAMETERS:
!
! The O-F grids bias_raob, rtms_raob, nobs_raob, bias_tovs, rtms_tovs,
! and nobs_tovs are re-written with the updated values.
!
! !FILES USED:
!
! The O-F grid files on input are opened, read from and closed
! on return. Routine i90_lua() is used to determine the temporary
! logical unit numbers.
!
! Internal configuration parameters are obtained from the resource
! file through the I90 package which is assumed to have been properly
! initialized.
!
! !REVISION HISTORY:
!
! 09Nov96 da Silva Specification and prologue.
!
!EOP
!-------------------------------------------------------------------------
include 'stdio.h'
write(stderr,*) 'It have NOT been written yet.'
return
end