module rof_comp_mct ! !USES: use shr_sys_mod use shr_kind_mod , only: IN=>SHR_KIND_IN, R8=>SHR_KIND_R8, CS=>SHR_KIND_CS use shr_file_mod , only: shr_file_getunit, shr_file_getlogunit, shr_file_getloglevel, & shr_file_setlogunit, shr_file_setloglevel, shr_file_setio, & shr_file_freeunit use shr_mpi_mod , only: shr_mpi_bcast use mct_mod use esmf use dead_data_mod use dead_mct_mod use dead_mod use seq_flds_mod use seq_cdata_mod use seq_infodata_mod use seq_timemgr_mod use seq_comm_mct , only: seq_comm_inst, seq_comm_name, seq_comm_suffix use seq_flds_mod , only: flds_d2x => seq_flds_r2x_fields, & flds_x2d => seq_flds_x2r_fields ! ! !PUBLIC TYPES: implicit none save private ! except !-------------------------------------------------------------------------- ! Public interfaces !-------------------------------------------------------------------------- public :: rof_init_mct public :: rof_run_mct public :: rof_final_mct !-------------------------------------------------------------------------- ! Private data interfaces !-------------------------------------------------------------------------- !--- stdin input stuff --- character(CS) :: str ! cpp defined model name !--- other --- integer(IN) :: dbug = 0 ! debug level (higher is more) real(R8), pointer :: gbuf(:,:) ! grid info buffer integer(IN) :: nproc_x ! num of i pes (type 3) integer(IN) :: seg_len ! length of segs (type 4) integer(IN) :: nxg ! global dim i-direction integer(IN) :: nyg ! global dim j-direction integer(IN) :: decomp_type ! data decomp type: character(CS) :: myModelName = 'rof' ! user defined model name integer(IN) :: ncomp = 6 ! component index integer(IN) :: my_task ! my task in mpi communicator mpicom integer(IN) :: master_task=0 ! task number of master task integer(IN) :: logunit ! logging unit number integer :: inst_index ! number of current instance (ie. 1) character(len=16) :: inst_name ! fullname of current instance (ie. "lnd_0001") character(len=16) :: inst_suffix ! char string associated with instance ! (ie. "_0001" or "") !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ CONTAINS !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ !=============================================================================== !BOP =========================================================================== ! ! !IROUTINE: rof_init_mct ! ! !DESCRIPTION: ! initialize dead rof model ! ! !REVISION HISTORY: ! ! !INTERFACE: ------------------------------------------------------------------ subroutine rof_init_mct( EClock, cdata, x2d, d2x, NLFilename ) ! !INPUT/OUTPUT PARAMETERS: type(ESMF_Clock) , intent(in) :: EClock type(seq_cdata) , intent(inout) :: cdata type(mct_aVect) , intent(inout) :: x2d, d2x character(len=*), optional , intent(in) :: NLFilename ! Namelist filename !EOP !--- local variables --- integer(IN) :: unitn ! Unit for namelist file integer(IN) :: ierr ! error code integer(IN) :: local_comm ! local communicator integer(IN) :: mype ! pe info integer(IN) :: totpe ! total number of pes integer(IN), allocatable :: gindex(:) ! global index integer(IN) :: COMPID integer(IN) :: mpicom integer(IN) :: lsize type(mct_gsMap) , pointer :: gsMap type(mct_gGrid) , pointer :: dom type(seq_infodata_type), pointer :: infodata ! Input init object integer(IN) :: shrlogunit, shrloglev ! original log unit and level !--- formats --- character(*), parameter :: F00 = "('(rof_init_mct) ',8a)" character(*), parameter :: F01 = "('(rof_init_mct) ',a,4i8)" character(*), parameter :: F02 = "('(rof_init_mct) ',a,4es13.6)" character(*), parameter :: F03 = "('(rof_init_mct) ',a,i8,a)" character(*), parameter :: F90 = "('(rof_init_mct) ',73('='))" character(*), parameter :: F91 = "('(rof_init_mct) ',73('-'))" character(*), parameter :: subName = "(rof_init_mct) " !------------------------------------------------------------------------------- ! !------------------------------------------------------------------------------- ! !------------------------------------------------------------------------------- ! Set cdata pointers call seq_cdata_setptrs(cdata, ID=COMPID, mpicom=mpicom, & gsMap=gsMap, dom=dom, infodata=infodata) call mpi_comm_rank(mpicom, my_task, ierr) inst_name = seq_comm_name(COMPID) inst_index = seq_comm_inst(COMPID) inst_suffix = seq_comm_suffix(COMPID) !--- open log file --- if (my_task == master_task) then logUnit = shr_file_getUnit() call shr_file_setIO('rof_modelio.nml'//trim(inst_suffix),logUnit) else logUnit = 6 endif !---------------------------------------------------------------------------- ! Reset shr logging to my log file !---------------------------------------------------------------------------- call shr_file_getLogUnit (shrlogunit) call shr_file_getLogLevel(shrloglev) call shr_file_setLogUnit (logUnit) ! read the namelist input (used to configure model) nxg = -9999 nyg = -9999 nproc_x = -9999 seg_len = -9999 decomp_type = -9999 if (my_task == master_task) then unitn = shr_file_getUnit() open( unitn, file='xrof_in'//trim(inst_suffix), status='old' ) read(unitn,*) nxg read(unitn,*) nyg read(unitn,*) decomp_type read(unitn,*) nproc_x read(unitn,*) seg_len close (unitn) call shr_file_freeunit(unitn) endif call shr_mpi_bcast(nxg ,mpicom,'xrof nxg') call shr_mpi_bcast(nyg ,mpicom,'xrof nyg') call shr_mpi_bcast(decomp_type,mpicom,'xrof decomp_type') call shr_mpi_bcast(nproc_x ,mpicom,'xrof nproc_x') call shr_mpi_bcast(seg_len ,mpicom,'xrof seg_len') if (my_task == master_task) then write(logunit,*) ' Read in Xrof input from file= xrof_in' write(logunit,F00) write(logunit,F00) ' Model : ',trim(myModelName) write(logunit,F01) ' NGX : ',nxg write(logunit,F01) ' NGY : ',nyg write(logunit,F01) ' Decomposition : ',decomp_type write(logunit,F03) ' Num pes in X : ',nproc_x,' (type 3 only)' write(logunit,F03) ' Segment Length : ',seg_len,' (type 11 only)' write(logunit,F01) ' inst_index : ',inst_index write(logunit,F00) ' inst_name : ',trim(inst_name) write(logunit,F00) ' inst_suffix : ',trim(inst_suffix) write(logunit,F00) call shr_sys_flush(logunit) end if ! Set flag to specify dead components call seq_infodata_PutData(infodata, dead_comps=.true., rof_present=.true., & rof_prognostic = .true., flood_present = .false., rof_nx=nxg, rof_ny=nyg) ! Determine communicator groups and sizes local_comm = mpicom call MPI_COMM_RANK(local_comm,mype ,ierr) call MPI_COMM_SIZE(local_comm,totpe,ierr) ! Initialize grid call dead_setNewGrid(decomp_type,nxg,nyg,totpe,mype,lsize,gbuf,seg_len,nproc_x) ! Initialize MCT global seg map allocate(gindex(lsize)) gindex(:) = gbuf(:,dead_grid_index) call mct_gsMap_init( gsMap, gindex, mpicom, COMPID, lsize, nxg*nyg) deallocate(gindex) ! Initialize MCT domain call dead_domain_mct(mpicom, gbuf, gsMap, dom) ! Initialize MCT attribute vectors call mct_aVect_init(d2x, rList=flds_d2x, lsize=lsize) call mct_aVect_zero(d2x) call mct_aVect_init(x2d, rList=flds_x2d, lsize=lsize) call mct_aVect_zero(x2d) ! Get relevant data to send back to driver !! call rof_run_mct(EClock, cdata, x2d, d2x) !---------------------------------------------------------------------------- ! Reset shr logging to original values !---------------------------------------------------------------------------- call shr_file_setLogUnit (shrlogunit) call shr_file_setLogLevel(shrloglev) call shr_sys_flush(logunit) end subroutine rof_init_mct !=============================================================================== !BOP =========================================================================== ! ! !IROUTINE: rof_run_mct ! ! !DESCRIPTION: ! run method for dead rof model ! ! !REVISION HISTORY: ! ! !INTERFACE: ------------------------------------------------------------------ subroutine rof_run_mct( EClock, cdata, x2d, d2x) implicit none ! !INPUT/OUTPUT PARAMETERS: type(ESMF_Clock) ,intent(in) :: EClock type(seq_cdata) ,intent(inout) :: cdata type(mct_aVect) ,intent(inout) :: x2d ! driver -> dead type(mct_aVect) ,intent(inout) :: d2x ! dead -> driver !EOP !--- local --- integer(IN) :: CurrentYMD ! model date integer(IN) :: CurrentTOD ! model sec into model date integer(IN) :: n ! index integer(IN) :: nf ! fields loop index integer(IN) :: ki ! index of ifrac integer(IN) :: lsize ! size of AttrVect real(R8) :: lat ! latitude real(R8) :: lon ! longitude integer(IN) :: shrlogunit, shrloglev ! original log unit and level integer :: nflds_d2x, nflds_x2d character(*), parameter :: F00 = "('(rof_run_mct) ',8a)" character(*), parameter :: F04 = "('(rof_run_mct) ',2a,2i8,'s')" character(*), parameter :: subName = "(rof_run_mct) " !------------------------------------------------------------------------------- !---------------------------------------------------------------------------- ! Reset shr logging to my log file !---------------------------------------------------------------------------- call shr_file_getLogUnit (shrlogunit) call shr_file_getLogLevel(shrloglev) call shr_file_setLogUnit (logUnit) lsize = mct_avect_lsize(x2d) nflds_d2x = mct_avect_nRattr(d2x) nflds_x2d = mct_avect_nRattr(x2d) ! UNPACK do nf=1,nflds_x2d do n=1,lsize ! ?? = x2d%rAttr(nf,n) enddo enddo ! PACK do nf=1,nflds_d2x do n=1,lsize d2x%rAttr(nf,n) = (nf+1) * 1.0_r8 ! lon = gbuf(n,dead_grid_lon) ! lat = gbuf(n,dead_grid_lat) ! d2x%rAttr(nf,n) = (nf*100) & ! * cos (SHR_CONST_PI*lat/180.0_R8) & ! * sin((SHR_CONST_PI*lon/180.0_R8) & ! - (ncomp-1)*(SHR_CONST_PI/3.0_R8) ) & ! + (ncomp*10.0_R8) enddo enddo ! log output for model date if (my_task == master_task) then call seq_timemgr_EClockGetData(EClock,curr_ymd=CurrentYMD, curr_tod=CurrentTOD) write(logunit,F04) trim(myModelName),': model date ', CurrentYMD,CurrentTOD call shr_sys_flush(logunit) end if !---------------------------------------------------------------------------- ! Reset shr logging to original values !---------------------------------------------------------------------------- call shr_file_setLogUnit (shrlogunit) call shr_file_setLogLevel(shrloglev) call shr_sys_flush(logunit) end subroutine rof_run_mct !=============================================================================== !BOP =========================================================================== ! ! !IROUTINE: rof_final_mct ! ! !DESCRIPTION: ! finalize method for dead rof model ! ! !REVISION HISTORY: ! ! !INTERFACE: ------------------------------------------------------------------ ! subroutine rof_final_mct( EClock, cdata, x2d, d2x) implicit none ! !INPUT/OUTPUT PARAMETERS: type(ESMF_Clock) ,intent(inout) :: EClock type(seq_cdata) ,intent(inout) :: cdata type(mct_aVect) ,intent(inout) :: x2d type(mct_aVect) ,intent(inout) :: d2x !EOP !--- formats --- character(*), parameter :: F00 = "('(rof_final_mct) ',8a)" character(*), parameter :: F91 = "('(rof_final_mct) ',73('-'))" character(*), parameter :: subName = "(rof_final_mct) " !------------------------------------------------------------------------------- ! !------------------------------------------------------------------------------- if (my_task == master_task) then write(logunit,F91) write(logunit,F00) trim(myModelName),': end of main integration loop' write(logunit,F91) end if end subroutine rof_final_mct !=============================================================================== end module rof_comp_mct