enkftv_allocate.f90

CLM集合卡曼滤波数据同化算法

  SUBROUTINE enkftv_allocate (lpt, kpt, maxsnl, msl,Ne)

!  ------------------------------------------------------------------------
! |                                                                        |
! |                                                                        |
! |                                                                        |
! |   Allocates memory space to ensemble time-varying variables            |
! |   For physical meanings of each variable, read clmtv_ensemble_module   |
! |                                                                        |
! |   source file: enkftv_allocate.f90                                     |
! |   first version: Qinjun,  Dai Yongjiu 5/26/2004                        |
!  -----------------------------------------------------------------------

  USE enkftv_module
  

  IMPLICIT NONE

  integer, INTENT(in) :: &
           lpt,       &!number of land points on clmlon x clmlat grid
           kpt,       &!total number of clm land points, including subgrid points
  !         msub,      &!number of submesh
           msl,       &!number of soil layers
           maxsnl,    &!number of snow layers
		   Ne          ! ensemble size used to perform data assimilation



!  --------------------------------
!  ** CLM time-varying variables **
!  --------------------------------

      allocate (en_snl(1:kpt, 1:Ne))

      allocate (en_dz     (maxsnl+1:msl, 1:kpt, 1:Ne))
      allocate (en_z      (maxsnl+1:msl, 1:kpt, 1:Ne))
      allocate (en_zi     (maxsnl+0:msl, 1:kpt, 1:Ne))

      allocate (en_tss    (maxsnl+1:msl, 1:kpt, 1:Ne))
      allocate (en_wliq   (maxsnl+1:msl, 1:kpt, 1:Ne))
      allocate (en_wice   (maxsnl+1:msl, 1:kpt, 1:Ne))
      allocate (en_rootr  (    1:msl, 1:kpt, 1:Ne))

      allocate (en_fveg   (1:kpt, 1:Ne))
      allocate (en_green  (1:kpt, 1:Ne))
      allocate (en_fsno   (1:kpt, 1:Ne))
      allocate (en_sigf   (1:kpt, 1:Ne))
      allocate (en_tlsun  (1:kpt, 1:Ne))
      allocate (en_tlsha  (1:kpt, 1:Ne))
      allocate (en_ldew   (1:kpt, 1:Ne))
      allocate (en_sag    (1:kpt, 1:Ne))
      allocate (en_scv    (1:kpt, 1:Ne))
      allocate (en_snowdp (1:kpt, 1:Ne))
      allocate (en_etrc   (1:kpt, 1:Ne))
      allocate (en_tg     (1:kpt, 1:Ne))

      allocate (en_lai    (1:kpt, 1:Ne))
      allocate (en_sai    (1:kpt, 1:Ne))

! --------------------------------------------
! **  error  for ensemble kalman filtering  **
! --------------------------------------------

      allocate (en_xerr   (1:kpt, 1:Ne))
      allocate (en_zerr   (1:kpt, 1:Ne))

! -------------------------------------------
! **  time  for ensemble kalman filtering  **
! -------------------------------------------

      allocate (en_iyear  ( 1:Ne))
      allocate (en_jday   ( 1:Ne))
      allocate (en_mcsec  ( 1:Ne))


! -----------------------------------------------------------------
! Albedo calculation for next time step
! -----------------------------------------------------------------

      allocate (en_albg   (1:2, 1:2, 1:kpt, 1:Ne))
      allocate (en_albv   (1:2, 1:2, 1:kpt, 1:Ne))
      allocate (en_alb    (1:2, 1:2, 1:kpt, 1:Ne))
      allocate (en_ssun   (1:2, 1:2, 1:kpt, 1:Ne))
      allocate (en_ssha   (1:2, 1:2, 1:kpt, 1:Ne))
      allocate (en_tranc  (1:2, 1:2, 1:kpt, 1:Ne))
      allocate (en_thermk           (1:kpt, 1:Ne))
      allocate (en_extkb            (1:kpt, 1:Ne))
      allocate (en_extkd            (1:kpt, 1:Ne))
      allocate (en_cosz             (1:kpt, 1:Ne))



  END SUBROUTINE enkftv_allocate