write_res0.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Subroutine: WRITE_RES0                                              C
!  Purpose: write out the initial restart records (namelist data)      C
!                                                                      C
!  Author: P. Nicoletti                               Date: 13-DEC-91  C
!  Reviewer: P. Nicoletti, W. Rogers, M. Syamlal      Date: 24-JAN-92  C
!                                                                      C
!  Revision Number:                                                    C
!  Purpose:                                                            C
!  Author:                                            Date: dd-mmm-yy  C
!  Reviewer:                                          Date: dd-mmm-yy  C
!                                                                      C
! TODO:                                                                C
!    this file may need work for GHD which has internally incremented  C
!    user mmax to mmax+1 to represent a mixture solids phase and the   C
!    number of real soldis phases is now smax... consequently several  C
!    physical properties will not be set at 'mmax' the mixture phase   C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C

      SUBROUTINE write_res0

!-----------------------------------------------
! Modules
!-----------------------------------------------
      USE bc
      USE cdist
      USE compar
      USE constant
      USE funits
      USE geometry
      USE ic
      USE in_binary_512i
      USE is
      USE leqsol
      USE machine
      USE mpi_utility      ! for gather
      USE output
      USE param
      USE param1
      USE physprop
      USE run
      USE rxns
      USE scalars
      USE scales
      USE sendrecv         ! for filling the boundary information
      USE stiff_chem
      USE toleranc
      USE ur_facs

      IMPLICIT NONE
!-----------------------------------------------
! Local variables
!-----------------------------------------------
      integer, allocatable :: arr1(:)
      integer, allocatable :: arr2(:)
      integer :: work_around(100)
! loop counters
      INTEGER :: LC, L, NN, IDX
! Pointer to the next record
      INTEGER :: NEXT_RECA
! file version id

! Place holder for deprecated variables:
      LOGICAL, PARAMETER :: lCALL_ISAT = .false.

      CHARACTER(LEN=512) :: VERSION
!-----------------------------------------------

      next_reca = 5

! note: the value below for version must be the same as the value
!       of version in mfix.f
!       If you change it in this subroutine, you must change it in
!       mfix.f also.

!       The value should also be consistent with the check in
!       read_res0

      version = 'RES = 01.8'

! Add new data entries at the end of the file and identify version no.

      if (mype.eq.pe_io) then
         allocate (arr1(ijkmax3))
         allocate (arr2(ijkmax2))
      else
         allocate (arr1(1))
         allocate (arr2(1))
         goto 1100
      end if

!------------------------------------------------------------------------
      WRITE (unit_res, rec=1) version
      WRITE (unit_res, rec=2) run_name, id_month, id_day, id_year, id_hour, &
         id_minute, id_second
      WRITE (unit_res, rec=3) next_reca
      WRITE (unit_res, rec=4) imin1, jmin1, kmin1, imax, jmax, kmax, imax1, &
         jmax1, kmax1, imax2, jmax2, kmax2, ijmax2, ijkmax2, mmax, dimension_ic&
         , dimension_bc, dimension_c, dimension_is, dt, xmin, xlength, ylength&
         , zlength, c_e, c_f, phi, phi_w
      CALL out_bin_512 (unit_res, c, dimension_c, next_reca)
      next_reca = 1 + next_reca                  ! work around for -O3 compiler bug
      next_reca = next_reca - 1
      DO lc = 1, dimension_c
         WRITE (unit_res, rec=next_reca) c_name(lc)
         next_reca = next_reca + 1
      END DO
      write (*,*) 'next_reca = ' , next_reca
      do l = 0,mmax
         write (*,*) l,nmax(l)
         work_around(l+1) = nmax(l)
      end do

!      next_reca = next_reca + 1
!      CALL OUT_BIN_512I (UNIT_RES, work_around, MMAX+1, NEXT_RECA)
!      WRITE (UNIT_RES, REC=NEXT_RECA) (work_around(L),L=1,MMAX+1)

      WRITE (unit_res, rec=next_reca) (nmax(l),l=0,mmax)

!      do l =0,mmax
!         write (unit_res,rec=next_reca) nmax(l)
!         next_reca = next_reca + 1
!      end do
!      write (unit_res,rec=next_reca) nmax(0) , nmax(1)

      next_reca = next_reca + 1
      write (*,*) ' next_reca = ' , next_reca

      CALL out_bin_512 (unit_res, dx(1), imax2, next_reca)
      CALL out_bin_512 (unit_res, dy(1), jmax2, next_reca)
      CALL out_bin_512 (unit_res, dz(1), kmax2, next_reca)
      WRITE (unit_res, rec=next_reca) run_name, description, units, run_type, &
         coordinates
      next_reca = next_reca + 1
      WRITE (unit_res, rec=next_reca) (d_p0(l),l=1,mmax), (ro_s0(l),l=1,mmax), &
         ep_star, ro_g0, mu_g0, mw_avg
      next_reca = next_reca + 1
      CALL out_bin_512 (unit_res, mw_g, nmax(0), next_reca)
      DO lc = 1, mmax
         WRITE (unit_res, rec=next_reca) (mw_s(lc,nn),nn=1,nmax(lc))
         next_reca = next_reca + 1
      END DO
      CALL out_bin_512 (unit_res, ic_x_w, dimension_ic, next_reca)
      CALL out_bin_512 (unit_res, ic_x_e, dimension_ic, next_reca)
      CALL out_bin_512 (unit_res, ic_y_s, dimension_ic, next_reca)
      CALL out_bin_512 (unit_res, ic_y_n, dimension_ic, next_reca)
      CALL out_bin_512 (unit_res, ic_z_b, dimension_ic, next_reca)
      CALL out_bin_512 (unit_res, ic_z_t, dimension_ic, next_reca)
      CALL out_bin_512i (unit_res, ic_i_w, dimension_ic, next_reca)
      CALL out_bin_512i (unit_res, ic_i_e, dimension_ic, next_reca)
      CALL out_bin_512i (unit_res, ic_j_s, dimension_ic, next_reca)
      CALL out_bin_512i (unit_res, ic_j_n, dimension_ic, next_reca)
      CALL out_bin_512i (unit_res, ic_k_b, dimension_ic, next_reca)
      CALL out_bin_512i (unit_res, ic_k_t, dimension_ic, next_reca)
      CALL out_bin_512 (unit_res, ic_ep_g, dimension_ic, next_reca)
      CALL out_bin_512 (unit_res, ic_p_g, dimension_ic, next_reca)
      CALL out_bin_512 (unit_res, ic_t_g, dimension_ic, next_reca)
      DO nn = 1, nmax(0)
         CALL out_bin_512 (unit_res, ic_x_g(1,nn), dimension_ic, next_reca)
      END DO
      CALL out_bin_512 (unit_res, ic_u_g, dimension_ic, next_reca)
      CALL out_bin_512 (unit_res, ic_v_g, dimension_ic, next_reca)
      CALL out_bin_512 (unit_res, ic_w_g, dimension_ic, next_reca)

      DO lc = 1, mmax
         CALL out_bin_512 (unit_res, ic_rop_s(1,lc), dimension_ic, next_reca)
         CALL out_bin_512 (unit_res, ic_u_s(1,lc), dimension_ic, next_reca)
         CALL out_bin_512 (unit_res, ic_v_s(1,lc), dimension_ic, next_reca)
         CALL out_bin_512 (unit_res, ic_w_s(1,lc), dimension_ic, next_reca)
         CALL out_bin_512 (unit_res, ic_t_s(1,lc), dimension_ic, next_reca)
         DO nn = 1, nmax(lc)
            CALL out_bin_512(unit_res,ic_x_s(1,lc,nn),dimension_ic,next_reca)
         END DO
      END DO
      CALL out_bin_512 (unit_res, bc_x_w, dimension_bc, next_reca)
      CALL out_bin_512 (unit_res, bc_x_e, dimension_bc, next_reca)
      CALL out_bin_512 (unit_res, bc_y_s, dimension_bc, next_reca)
      CALL out_bin_512 (unit_res, bc_y_n, dimension_bc, next_reca)
      CALL out_bin_512 (unit_res, bc_z_b, dimension_bc, next_reca)
      CALL out_bin_512 (unit_res, bc_z_t, dimension_bc, next_reca)

      CALL out_bin_512i (unit_res, bc_i_w, dimension_bc, next_reca)
      CALL out_bin_512i (unit_res, bc_i_e, dimension_bc, next_reca)
      CALL out_bin_512i (unit_res, bc_j_s, dimension_bc, next_reca)
      CALL out_bin_512i (unit_res, bc_j_n, dimension_bc, next_reca)
      CALL out_bin_512i (unit_res, bc_k_b, dimension_bc, next_reca)
      CALL out_bin_512i (unit_res, bc_k_t, dimension_bc, next_reca)
      CALL out_bin_512 (unit_res, bc_ep_g, dimension_bc, next_reca)
      CALL out_bin_512 (unit_res, bc_p_g, dimension_bc, next_reca)
      CALL out_bin_512 (unit_res, bc_t_g, dimension_bc, next_reca)
      DO nn = 1, nmax(0)
         CALL out_bin_512 (unit_res, bc_x_g(1,nn), dimension_bc, next_reca)
      END DO
      CALL out_bin_512 (unit_res, bc_u_g, dimension_bc, next_reca)
      CALL out_bin_512 (unit_res, bc_v_g, dimension_bc, next_reca)
      CALL out_bin_512 (unit_res, bc_w_g, dimension_bc, next_reca)
      CALL out_bin_512 (unit_res, bc_ro_g, dimension_bc, next_reca)
      CALL out_bin_512 (unit_res, bc_rop_g, dimension_bc, next_reca)
      CALL out_bin_512 (unit_res, bc_volflow_g, dimension_bc, next_reca)
      CALL out_bin_512 (unit_res, bc_massflow_g, dimension_bc, next_reca)
      DO lc = 1, mmax
         CALL out_bin_512 (unit_res, bc_rop_s(1,lc), dimension_bc, next_reca)
         CALL out_bin_512 (unit_res, bc_u_s(1,lc), dimension_bc, next_reca)
         CALL out_bin_512 (unit_res, bc_v_s(1,lc), dimension_bc, next_reca)
         CALL out_bin_512 (unit_res, bc_w_s(1,lc), dimension_bc, next_reca)
         CALL out_bin_512 (unit_res, bc_t_s(1,lc), dimension_bc, next_reca)
         DO nn = 1, nmax(lc)
            CALL out_bin_512(unit_res,bc_x_s(1,lc,nn),dimension_bc,next_reca)
         END DO
         CALL out_bin_512 (unit_res, bc_volflow_s(1,lc), dimension_bc, &
            next_reca)
         CALL out_bin_512 (unit_res, bc_massflow_s(1,lc), dimension_bc, &
            next_reca)
      END DO
      DO lc = 1, dimension_bc
         WRITE (unit_res, rec=next_reca) bc_type(lc)
         next_reca = next_reca + 1
      END DO


 1100 continue

!      call MPI_Barrier(MPI_COMM_WORLD,mpierr)  !//PAR_I/O enforce barrier here
!      call send_recv (flag,2)
      call gather (flag,arr1,root)
! To take care of filling the processor ghost layers with the correct values
      call scatter (flag,arr1,root)
!      call MPI_Barrier(MPI_COMM_WORLD,mpierr)  !//PAR_I/O enforce barrier here
      if (mype .ne. pe_io) goto 1200
      call convert_to_io_i(arr1,arr2,ijkmax2)
      CALL out_bin_512i (unit_res, arr2, ijkmax2, next_reca)

      CALL out_bin_512 (unit_res, is_x_w, dimension_is, next_reca)
      CALL out_bin_512 (unit_res, is_x_e, dimension_is, next_reca)
      CALL out_bin_512 (unit_res, is_y_s, dimension_is, next_reca)
      CALL out_bin_512 (unit_res, is_y_n, dimension_is, next_reca)
      CALL out_bin_512 (unit_res, is_z_b, dimension_is, next_reca)
      CALL out_bin_512 (unit_res, is_z_t, dimension_is, next_reca)
      CALL out_bin_512i (unit_res, is_i_w, dimension_is, next_reca)
      CALL out_bin_512i (unit_res, is_i_e, dimension_is, next_reca)
      CALL out_bin_512i (unit_res, is_j_s, dimension_is, next_reca)
      CALL out_bin_512i (unit_res, is_j_n, dimension_is, next_reca)
      CALL out_bin_512i (unit_res, is_k_b, dimension_is, next_reca)
      CALL out_bin_512i (unit_res, is_k_t, dimension_is, next_reca)
      CALL out_bin_512 (unit_res, is_pc(1,1), dimension_is, next_reca)
      CALL out_bin_512 (unit_res, is_pc(1,2), dimension_is, next_reca)
      DO lc = 1, mmax
         CALL out_bin_512 (unit_res, is_vel_s(1,lc), dimension_is, next_reca)
      END DO
      DO lc = 1, dimension_is
         WRITE (unit_res, rec=next_reca) is_type(lc)
         next_reca = next_reca + 1
      END DO
      WRITE (unit_res, rec=next_reca) cyclic_x, cyclic_y, cyclic_z, cyclic_x_pd&
         , cyclic_y_pd, cyclic_z_pd, delp_x, delp_y, delp_z, u_g0, u_s0, v_g0, &
         v_s0, w_g0, w_s0
      next_reca = next_reca + 1

! Version 01.09
! ------------------------------------------------------------------------
      WRITE (unit_res, rec=next_reca) time, tstop, energy_eq, res_dt, out_dt, &
         nlog, l_scale0, no_i, no_j, no_k, call_usr
      next_reca = next_reca + 1
      write (unit_res,rec=next_reca) n_spx
      next_reca = next_reca + 1
      DO lc = 1, n_spx
         WRITE (unit_res, rec=next_reca) spx_dt(lc)
         next_reca = next_reca + 1
      END DO
      DO lc = 0, mmax
         WRITE (unit_res, rec=next_reca) species_eq(lc)
         next_reca = next_reca + 1
      END DO
      CALL out_bin_512 (unit_res, usr_dt, dimension_usr, next_reca)
      CALL out_bin_512 (unit_res, usr_x_w, dimension_usr, next_reca)
      CALL out_bin_512 (unit_res, usr_x_e, dimension_usr, next_reca)
      CALL out_bin_512 (unit_res, usr_y_s, dimension_usr, next_reca)
      CALL out_bin_512 (unit_res, usr_y_n, dimension_usr, next_reca)
      CALL out_bin_512 (unit_res, usr_z_b, dimension_usr, next_reca)
      CALL out_bin_512 (unit_res, usr_z_t, dimension_usr, next_reca)
      DO lc = 1, dimension_usr
         WRITE (unit_res, rec=next_reca) usr_format(lc), usr_ext(lc), usr_type(&
            lc), usr_var(lc)
         next_reca = next_reca + 1
      END DO
      CALL out_bin_512 (unit_res, ic_p_star, dimension_ic, next_reca)
      CALL out_bin_512 (unit_res, ic_l_scale, dimension_ic, next_reca)
      DO lc = 1, dimension_ic
         WRITE (unit_res, rec=next_reca) ic_type(lc)
         next_reca = next_reca + 1
      END DO
      CALL out_bin_512 (unit_res, bc_dt_0, dimension_bc, next_reca)
      CALL out_bin_512 (unit_res, bc_jet_g0, dimension_bc, next_reca)
      CALL out_bin_512 (unit_res, bc_dt_h, dimension_bc, next_reca)
      CALL out_bin_512 (unit_res, bc_jet_gh, dimension_bc, next_reca)
      CALL out_bin_512 (unit_res, bc_dt_l, dimension_bc, next_reca)
      CALL out_bin_512 (unit_res, bc_jet_gl, dimension_bc, next_reca)


! Version 01.10
! ------------------------------------------------------------------------
      WRITE (unit_res, rec=next_reca) mu_gmax
      next_reca = next_reca + 1

! Version 01.11
! ------------------------------------------------------------------------
      WRITE (unit_res, rec=next_reca) v_ex, model_b
      next_reca = next_reca + 1

! Version 01.12
! ------------------------------------------------------------------------
      WRITE (unit_res, rec=next_reca) p_ref, p_scale, ur_fac, tol_resid, dt_max&
         , dt_min, dt_fac, close_packed, gravity, mu_s0(1)
      next_reca = next_reca + 1
      WRITE (unit_res, rec=next_reca) leq_it, leq_method
      next_reca = next_reca + 1
      CALL out_bin_512 (unit_res, bc_hw_g, dimension_bc, next_reca)
      CALL out_bin_512 (unit_res, bc_uw_g, dimension_bc, next_reca)
      CALL out_bin_512 (unit_res, bc_vw_g, dimension_bc, next_reca)
      CALL out_bin_512 (unit_res, bc_ww_g, dimension_bc, next_reca)
      DO lc = 1, mmax
         CALL out_bin_512 (unit_res, bc_hw_s(1,lc), dimension_bc, next_reca)
         CALL out_bin_512 (unit_res, bc_uw_s(1,lc), dimension_bc, next_reca)
         CALL out_bin_512 (unit_res, bc_vw_s(1,lc), dimension_bc, next_reca)
         CALL out_bin_512 (unit_res, bc_ww_s(1,lc), dimension_bc, next_reca)
      END DO
      WRITE (unit_res, rec=next_reca) momentum_x_eq, momentum_y_eq, &
         momentum_z_eq, tol_diverge, discretize, full_log
      next_reca = next_reca + 1

! Version 01.14
! ------------------------------------------------------------------------
      WRITE (unit_res, rec=next_reca) detect_stall
      next_reca = next_reca + 1

! Version 01.15
! ------------------------------------------------------------------------
      WRITE (unit_res, rec=next_reca) &
         k_g0, k_s0(1), c_pg0, c_ps0(1), tol_resid_t, tol_resid_x
      next_reca = next_reca + 1
      CALL out_bin_512 (unit_res, ic_gama_rg, dimension_ic, next_reca)
      CALL out_bin_512 (unit_res, ic_t_rg, dimension_ic, next_reca)
      DO lc = 1, mmax
         CALL out_bin_512 (unit_res, ic_gama_rs(1,lc), dimension_ic, next_reca)
         CALL out_bin_512 (unit_res, ic_t_rs(1,lc), dimension_ic, next_reca)
      END DO

! Version 01.2
! ------------------------------------------------------------------------
      WRITE (unit_res, rec=next_reca) norm_g, norm_s
      next_reca = next_reca + 1

! Version 01.3
! ------------------------------------------------------------------------
      WRITE (unit_res, rec=next_reca) nscalar, tol_resid_scalar, dim_scalar
      next_reca = next_reca + 1
      CALL out_bin_512i (unit_res, phase4scalar, dim_scalar, next_reca)

! Version 1.4 -- write radiation variables in write_res1
! ------------------------------------------------------------------------

! Version 1.5 -- write nRR
! ------------------------------------------------------------------------
      WRITE (unit_res, rec=next_reca) nrr
      next_reca = next_reca + 1

! Version 1.6 -- write k and epsilon in write_res1 and spx1
! ------------------------------------------------------------------------
      WRITE (unit_res, rec=next_reca) k_epsilon
      next_reca = next_reca + 1

! Version 1.7 -- write STIFF_CHEMISTRY and lCALL_ISAT
! ------------------------------------------------------------------------
      WRITE (unit_res, rec=next_reca) stiff_chemistry, lcall_isat
      next_reca = next_reca + 1

! Version 1.8 -- write densities of each solids species
! ------------------------------------------------------------------------
      WRITE (unit_res, rec=next_reca) (solve_ros(lc),lc=1,mmax)
      next_reca = next_reca + 1
      DO lc = 1, mmax
         idx = inert_species(lc)
         WRITE (unit_res, rec=next_reca) idx, ro_s0(lc),&
            (x_s0(lc,nn),nn=1,nmax(lc))
         next_reca = next_reca + 1
      ENDDO



!  Add new write statements above this line.  Remember to update NEXT_RECA.
!  Remember to change the version number near beginning of this subroutine.
!  Also modify READ_RES0.  The routines such as OUT_BIN_512 etc. writes
!  arrays dimensioned ARRAY(DIM).  So arrays dimensioned ARRAY(DIM1:DIM2)
!  should be passed as ARRAY(DIM1) and array length as DIM2-DIM1+1.
!---------------------------------------------------------------------------
      WRITE (unit_res, rec=3) next_reca
      FLUSH (unit_res)

 1200 continue

      if (bdist_io .and. mype.ne.pe_io) then
         WRITE (unit_res, rec=1) version
         WRITE (unit_res, rec=2) run_name, id_month, id_day, id_year, id_hour, &
             id_minute, id_second
         WRITE (unit_res, rec=3) 4
         FLUSH (unit_res)
      end if

      deallocate (arr1)
      deallocate (arr2)

      RETURN
      END SUBROUTINE write_res0