enkf_clm.f90

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

    
      PROGRAM enkf_clm 
!
!     Common Land Model, initial version managered by Yongjiu Dai
! ======================================================================
!
      use enkftv_module
      use clmtv_module
      IMPLICIT NONE

      integer &
           lon,             &!longitude points
           lat,             &!latitude points
           msub,            &!number of submesh
           mstep,           &!model step for simulation [-]

           mbyear,          &!base year of run
           mbday,           &!base julian day of run
           mbsec,           &!base seconds of starting date (e.g., 0)
           nsrest,          &!0 = initial run, > 0 = re-start run

           lun_ini,         &!logical unit number of initial data
           lun_rst,         &!logical unit number of re-start file
           lun_srf,         &!logical unit number of surface data
           lun_lai,         &!logical unit number of LAI data
           lun_met,         &!logical unit number of meteorological data
           lun_out,         &!logical unit number of output file

           istep             !loop number

      real dt_step           ! time step of model (second)

      logical &
           rdlsf,           &!true if read in initial data set
           rdlai,           &!true if read in LAI data set
           greenwich         !true if the input starting time is greenwich time
                             !false -> local time
      integer::loc_yea,loc_day,loc_hou    


! ======================================================================
! (1) Define the run and open files (for off-line use)
! ======================================================================

      call  data (lon     ,lat     ,msub    ,&
                  nsrest  ,rdlsf   ,rdlai   ,dt_step   ,mstep ,&
                  mbyear  ,mbday   ,mbsec   ,greenwich ,& 
                  lun_ini ,lun_rst ,lun_srf ,lun_lai   ,lun_met, lun_out)

! ======================================================================
! (2) Set the time-invariant and time-varying variables
! ======================================================================

      call clm_ini (lon     ,lat     ,msub ,&
                    nsrest  ,rdlsf   ,rdlai ,&
                    dt_step ,mbyear  ,mbday   ,mbsec ,greenwich ,&
                    lun_ini ,lun_rst ,lun_srf ,lun_lai)
      print*,"clm_ini"
!  -------------------------------------------------------------------------------
!  **  initilizing ensemble state variables and derived time-varying parameters **
!  -------------------------------------------------------------------------------
     
     call enkf_ini( )
     print*," enkf_ensemble_ini "
!===========================================================================

      open(8,file='output_txt',form='formatted')!open output file 
!      open(777,file='Tb_txt',form='formatted')
!      open(unit=666,file='date_list_30.txt',form='formatted')!open local time file
!==========================================================================

! ======================================================================
! (3) Time loop 
! ======================================================================

      do istep = 1, mstep

       ! Read in the atmospheric forcing
         call getmet (lon ,lat ,lun_met)

       ! Main driver for CLM
         call enkf_drv (lon ,lat ,istep ,lun_rst, rdlai)


       !  call clm_drv(lon ,lat ,istep ,lun_rst, rdlai)

       ! Write-out
       !  call clm_out (lun_out ,lon ,lat)

      end do
!
! ======================================================================

      write(6,*) 'Execution Completed'
      print*,'game over'

      END PROGRAM enkf_clm