!-------------------------------------- LICENCE BEGIN ------------------------------------
!Environment Canada - Atmospheric Science and Technology License/Disclaimer,
!                     version 3; Last Modified: May 7, 2008.
!This is free but copyrighted software; you can use/redistribute/modify it under the terms
!of the Environment Canada - Atmospheric Science and Technology License/Disclaimer
!version 3 or (at your option) any later version that should be found at:
!http://collaboration.cmc.ec.gc.ca/science/rpn.comm/license.html
!
!This software is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
!without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
!See the above mentioned License/Disclaimer for more details.
!You should have received a copy of the License/Disclaimer along with this software;
!if not, you can write to: EC-RPN COMM Group, 2121 TransCanada, suite 500, Dorval (Quebec),
!CANADA, H9P 1J3; or send e-mail to service.rpn@ec.gc.ca
!-------------------------------------- LICENCE END --------------------------------------
!

      subroutine stdevout_la 1,8
#if defined (DOC)
*
***s/r stdevout_la  - Write LAM background-error st-dev used in the current experiment on the 
!                     analysis increment file (analysis resolution).
*
*Author  : L. Fillion - ARMA/EC - 21 Feb 2007.
*Revision: 
*

*Arguments
*
#endif
C
      IMPLICIT NONE
*implicits
#include "taglam4d.cdk"
#include "pardim.cdk"
#include "comdim.cdk"
#include "comlun.cdk"
#include "comcst.cdk"
#include "comgem.cdk"
#include "cvcord.cdk"
#include "rpnstd.cdk"
#include "comgdpar.cdk"
#include "comct0.cdk"
#include "comgrd_param.cdk"
#include "comgrd.cdk"
#include "comgemla.cdk"
#include "comrgsigla.cdk"
*
      integer initer
      character*8 clletiket
!
      logical lltictac
      character*2 clvar
      integer VFSTLUK
      integer jlev,imode,ilev,ji,jj,jk,jpt
      integer imaxlev
      parameter(imaxlev=200)
      integer itrlgid,iip1s(imaxlev),iip2,iip3
      integer itrlnlev
      integer iig2,iig3,iig4,ikey
      integer ikind,ibrpstamp
      integer ezgprm
      integer idum1,idum2,idum3,idum4
!
      real*8 zmin,zmax
      real*8 zwork
      real*8 zfldin(mni_in,mnj_in,nflev)
      real*8 pgdpsi(ni,nj,nflev)
      real*8 pgdchiu(ni,nj,nflev)
      real*8 pgdtu(ni,nj,nflev)
      real*8 pgdpsu(ni,nj)
!
!
      character*1 cltypinc,clstring

      integer jvar,idim,iig1,iunit,FNOM,FSTOUV,IERR
      integer fstopc
      real*8 zcon

!
!!
!      iunit=700
!      IERR = FNOM(iunit,'aincr.fst','RND',0)
!      IERR = FSTOUV(iunit,'RND')

      lltictac = .false.
!
      ibrpstamp=nbrpstamp
      IPAK = -32
      IDATYP = 5
      cltypvar = 'E'
!
!     Writing positional parameters
!     -----------------------------
!
      if(lltictac) then
!
! Write >>:
        cltypvar = 'X'
        clnomvar = '>>'
        clgrtyp = 'E'
        ierr = vfstecr(grd_x_8,zwork,-inbits,nulinclr,ibrpstamp,
     &          ideet,inpas,mni_in,1,1,mig1in,mig2in,mig3in,cltypvar
     &          ,clnomvar,cetikinc,clgrtyp,mig1tic,mig2tic,mig3tic,mig4tic
     &          ,ndtinc,.true.)
!
! Write ^^
        cltypvar = 'X'
        clnomvar = '^^'
        clgrtyp = 'E'
        ierr = vfstecr(grd_y_8,zwork,-inbits,nulinclr,ibrpstamp,
     &          ideet,inpas,1,mnj_in,1,mig1in,mig2in,mig3in,cltypvar
     &          ,clnomvar,cetikinc,clgrtyp,mig1tic,mig2tic,mig3tic,mig4tic
     &          ,ndtinc,.true.)
      endif
!
!     Writing Fields
!     --------------
!
! Psi
      do jlev=1,nflev
        do jj=1,mnj_in
          do ji=1,mni_in
            zfldin(ji,jj,jlev) = rgsiglapp(ji,jlev,jj)
          enddo
        enddo
      enddo
      clletiket='PP_STDEV'
      call ecr_la(zfldin,'PP',niter,clletiket,ipak,ibrpstamp,ideet,inpas,
     &            idatyp,cltypvar,.false.,nulinclr,mni_in,mnj_in,nflev)
!
! Chiu
      do jlev=1,nflev
        do jj=1,mnj_in
          do ji=1,mni_in
            zfldin(ji,jj,jlev) = rgsiglacu(ji,jlev,jj)
          enddo
        enddo
      enddo
      clletiket='CU_STDEV'
      call ecr_la(zfldin,'CC',niter,clletiket,ipak,ibrpstamp,ideet,inpas,
     &            idatyp,cltypvar,.false.,nulinclr,mni_in,mnj_in,nflev)
!
! Tu
      do jlev=1,nflev
        do jj=1,mnj_in
          do ji=1,mni_in
            zfldin(ji,jj,jlev) = rgsiglatu(ji,jlev,jj)
          enddo
        enddo
      enddo
      clletiket='TU_STDEV'
      call ecr_la(zfldin,'TT',niter,clletiket,ipak,ibrpstamp,ideet,inpas,
     &            idatyp,cltypvar,.false.,nulinclr,mni_in,mnj_in,nflev)
!
! Lnq
      do jlev=1,nflev
        do jj=1,mnj_in
          do ji=1,mni_in
            zfldin(ji,jj,jlev) = rgsiglalq(ji,jlev,jj)
          enddo
        enddo
      enddo
      clletiket='LQ_STDEV'
      call ecr_la(zfldin,'HU',niter,clletiket,ipak,ibrpstamp,ideet,inpas,
     &            idatyp,cltypvar,.false.,nulinclr,mni_in,mnj_in,nflev)
!
! Psu
      do jj=1,mnj_in
        do ji=1,mni_in
          zfldin(ji,jj,1) = rgsiglapu(ji,1,jj)
        enddo
      enddo
      clletiket='PS_STDEV'
      call ecr_la(zfldin,'P0',niter,clletiket,ipak,ibrpstamp,ideet,inpas,
     &            idatyp,cltypvar,.false.,nulinclr,mni_in,mnj_in,1)
!
! Coriol
      do jj=1,mnj_in
        do ji=1,mni_in
          zfldin(ji,jj,1) = coriol(ji,jj)
        enddo
      enddo
      clletiket='CORIOLIS'
      call ecr_la(zfldin,'QQ',niter,clletiket,ipak,ibrpstamp,ideet,inpas,
     &            idatyp,cltypvar,.false.,nulinclr,mni_in,mnj_in,1)
!
      return
      end