CTLBLK.f

  module ctlblk_mod
!-----------------------------------------------------------------------
! module: CTLBLK
!
! ABSTRACT: 
! this module is replacing the CTLBLK.comm, all the comm block is 
!    removed.
! Revision Log:
!  2011-02    Jun Wang  - ADD variables for grib2
!  2011-12-14 SARAH LU  - ADD AER FILENAME
!  2011-12-23 SARAH LU  - ADD NBIN FOR DU, SS, OC, BC, SU
!  2021-09-30 JESSE MENG- 2D DECOMPOSITION
!-----------------------------------------------------------------------
!
  implicit none
!
  type field_info
    integer ifld
    integer lvl
    integer lvl1,lvl2
    integer ntrange
    integer tinvstat
  end type
  integer, parameter :: komax=70
  integer, parameter :: lsmdef=46             ! default number of p levels
  integer,PARAMETER  :: nfd=18,nbnd=6
  REAL,  PARAMETER   :: qmin = 1.e-15
!
  integer :: novegtype ! max number of veg type
!
  character(len=256) :: filename,filenameflux,filenamed3d,filenameaer, &
                        filenameflat
  character(len=19)  :: datestr
  character(len=4)   :: modelname, submodelname
  character(len=8)   :: fullmodelname
  character(len=20)  :: ioform
  character(len=4)   :: vtimeunits
!      
  character(5) :: grib
  type(field_info),allocatable :: fld_info(:)
  integer :: cfld,ntlfld,npset
  real*8 :: gdsdegr
  real,allocatable :: datapd(:,:,:)
!
  logical :: gocart_on, d3d_on, hyb_sigp, rdaod, aqfcmaq_on
  logical :: sigma,run,first,restrt
  logical :: global
  logical :: smflag
  integer :: idat(5),ihrst, nfcst,nbc,list,iout,ifhr,ntstm,            &
             nddamp,nprec,idtad,nboco,nshde,ncp,imdlty,nphs,           &
             nrads,nradl,imin,ifmin,datahandle,imp_physics,            &
             icu_physics,isf_surface_physics,isec,icount_calmict,      &
             ivegsrc
  real :: dt,sdat(3),avrain,avcnvc,dtq2,pt,pdtop,                      &
          spl(komax),alsl(komax),prec_acc_dt,pt_tbl,prec_acc_dt1,spval
! real :: SPVAL=9.9e10                                     ! Moorthi
!
  integer :: num_procs,me,jsta,jend,ista,iend,                         &
               jsta_m,jend_m, jsta_m2,jend_m2,                         &
               ista_m,iend_m,ista_m2,iend_m2,                          &
             iup,idn,icnt(0:1023),idsp(0:1023),  icnt2(0:1023),idsp2(0:1023),                      & 
             jsta_2l, jend_2u,jvend_2u,                                &
             ista_2l, iend_2u,ivend_2u,                                &
             num_servers, mpi_comm_inter,    &
             mpi_comm_comp, im,jm,lm,nsoil,lp1,lm1,im_jm,              &
             ileft,iright,                                             &
                ileftb,irightb ,                                       & 
             ibsize,ibsum,                                             &
             lsm,lsmp1                                    !comm mpi
        integer, allocatable :: icoords(:,:),ibcoords(:,:)
        real   , allocatable :: rcoords(:,:),rbcoords(:,:)
        real, allocatable :: bufs(:),buff(:)                           
        integer , allocatable :: isxa(:),iexa(:),jsxa(:),jexa(:)       
        integer numx
        integer, allocatable :: ibufs(:)
        real, allocatable :: rbufs(:)
!
  real :: ardsw, ardlw, asrfc, tsrfc,trdlw,trdsw,tclod,theat,          &
          tprec,tmaxmin,td3d                              !comm rad
!
  real pthresh ! moved from params because it is defined differently for NAM
!  
  real(kind=8) :: etafld2_tim=0.,eta2p_tim=0.,surfce2_tim=0.,          &
                  cldrad_tim=0.,miscln_tim=0.,fixed_tim=0.,            &
                  mdl2sigma_tim=0.,readxml_tim=0.,mdl2agl_tim=0.,      &
                  mdl2std_tim=0.,mdl2thandpv_tim=0.,calrad_wcloud_tim=0.!comm tim_info
!
  real(kind=8) :: time_output=0., time_e2out=0.           !comm jjt
!
  real :: spldef(lsmdef) =                                             &
      (/200.,500.,700.,1000.,2000.,3000.                               &
      ,5000.,7000.,7500.,10000.,12500.,15000.,17500.,20000.,22500.     &
      ,25000.,27500.,30000.,32500.,35000.,37500.,40000.,42500.,45000.  &
      ,47500.,50000.,52500.,55000.,57500.,60000.,62500.,65000.         &
      ,67500.,70000.,72500.,75000.,77500.,80000.,82500.,85000.         &
      ,87500.,90000.,92500.,95000.,97500.,100000./)
!
  REAL htfd(nfd),petabnd(nbnd),sigbnd(nbnd)

! Add GOCART aerosol specification
  integer, parameter :: nbin_du = 5         ! dust
  integer, parameter :: nbin_ss = 5         ! sea salt
  integer, parameter :: nbin_oc = 2         ! organic carbon
  integer, parameter :: nbin_bc = 2         ! black carbon
  integer, parameter :: nbin_su = 1         ! sulfate
  integer, parameter :: nbin_sm = 1             ! smoke
!
!     SET FD LEVEL HEIGHTS IN GEOPOTENTAL METERS.
      DATA htfd  / 20.e0,30.e0,40.e0,50.e0,80.e0,100.e0,305.e0,457.e0,610.e0,   &
           914.e0,1524.e0,1829.e0,2134.e0,2743.e0,3658.e0,4572.e0, &
       6000.e0,7010.e0/
!
!     SET MIDPOINT "SIGMA" VALUES FOR ETA BOUNDARY LAYERS.
      DATA sigbnd / 0.985,0.955,0.925,0.895,0.865,0.835 /
      DATA petabnd / 15.,45.,75.,105.,135.,165./
!
!-----------------------------------------------------------------------
  end module ctlblk_mod