PROCESS.f

Go to the documentation of this file.

      SUBROUTINE process(kth,kpv,th,pv,iostatusD3D)
!
!----------------------------------------------------------------------------
!
      use mpi, only: mpi_wtime

      use ctlblk_mod, only: cfld, etafld2_tim, eta2p_tim, mdl2sigma_tim, surfce2_tim,&
                            mdl2agl_tim, mdl2std_tim, mdl2thandpv_tim, calrad_wcloud_tim,&
                            cldrad_tim, miscln_tim, fixed_tim, ntlfld, me
!- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
      implicit none
!
!------ DECLARE VARIABLES.
!
      integer,intent(in) :: kth
      integer,intent(in) :: kpv
      integer,intent(in) :: iostatusd3d
      real,intent(in)    :: th(kth)
      real,intent(in)    :: pv(kpv)
      real(kind=8)       :: btim
      CHARACTER*6           datset,proj
      LOGICAL               north
!
!
!****************************************************************************
!     START SUBROUTINE PROCESS.
!
      cfld=0
      if(me==0) write(0,*) "PROCESS starts"
!
!     COMPUTE/POST FIELDS ON MDL SURFACES.
!
      btim = mpi_wtime()
      CALL mdlfld
      if(me==0) write(0,*) "PROCESS MDLFLD done"
      etafld2_tim = etafld2_tim +(mpi_wtime() - btim)
!
!     COMPUTE/POST FIELDS ON PRESSURE SURFACES.
      btim = mpi_wtime()
      CALL mdl2p(iostatusd3d)
      if(me==0) write(0,*) "PROCESS MDL2P done"
      eta2p_tim = eta2p_tim +(mpi_wtime() - btim)
!
!     COMPUTE/POST FIELDS ON SIGMA SURFACES.
      btim = mpi_wtime()
      CALL mdl2sigma
      if(me==0) write(0,*) "PROCESS MDL2SIGMA done"
      CALL mdl2sigma2
      if(me==0) write(0,*) "PROCESS MDL2SIGMA2 done"
      mdl2sigma_tim = mdl2sigma_tim +(mpi_wtime() - btim)
!
!     COMPUTE/POST FIELDS ON AGL SURFCES.
      btim = mpi_wtime()
      CALL mdl2agl
      if(me==0) write(0,*) "PROCESS MDL2AGL done"
      mdl2agl_tim = mdl2agl_tim +(mpi_wtime() - btim)
!
!     COMPUTE/POST SURFACE RELATED FIELDS.
      btim = mpi_wtime()
      CALL surfce
      if(me==0) write(0,*) "PROCESS SURFCE done"
      surfce2_tim = surfce2_tim +(mpi_wtime() - btim)
!
!     COMPUTE/POST SOUNDING AND CLOUD RELATED FIELDS.
      btim = mpi_wtime()
      CALL cldrad
      if(me==0) write(0,*) "PROCESS CLDRAD done"
      cldrad_tim = cldrad_tim +(mpi_wtime() - btim)
!
!     COMPUTE/POST TROPOPAUSE DATA, FD LEVEL FIELDS,
!     FREEZING LEVEL HEIGHT AND RH, BOUNDARY LAYER FIELDS,
!     AND LFM-NGM LOOK-ALIKE FIELDS.
      btim = mpi_wtime()
      CALL miscln
      if(me==0) write(0,*) "PROCESS MISCLN done"
      miscln_tim = miscln_tim +(mpi_wtime() - btim)

!     COMPUTE/POST TROPOPAUSE DATA, FD LEVEL FIELDS,
!     FREEZING LEVEL HEIGHT AND RH, BOUNDARY LAYER FIELDS,
!     AND LFM-NGM LOOK-ALIKE FIELDS.
      btim = mpi_wtime()
      CALL mdl2std_p
      if(me==0) write(0,*) "PROCESS MDL2STD_P done"
      mdl2std_tim = mdl2std_tim +(mpi_wtime() - btim)
!
!     POST FIXED FIELDS.
      btim = mpi_wtime()
      CALL fixed
      if(me==0) write(0,*) "PROCESS FIXED done"
      fixed_tim =  fixed_tim +(mpi_wtime() - btim)
!
!     COMPUTE/POST FIELDS ON SIGMA SURFACES.
      btim = mpi_wtime()
      CALL mdl2thandpv(kth,kpv,th,pv)
      if(me==0) write(0,*) "PROCESS MDL2THANDPV done"
      mdl2thandpv_tim = mdl2thandpv_tim +(mpi_wtime() - btim)
!
!     POST RADIANCE AND BRIGHTNESS FIELDS.
      btim = mpi_wtime()
      CALL calrad_wcloud
      if(me==0) write(0,*) "PROCESS CALRAD_WCLOUD done"
      calrad_wcloud_tim = calrad_wcloud_tim +(mpi_wtime() - btim)
!
!     END OF ROUTINE.
!
      ntlfld=cfld
      if(me==0)print *,'nTLFLD=',ntlfld
      if(me==0) write(0,*) "PROCESS done"
!
      RETURN
      END