mod_realinput_scale Module Reference

module REAL input SCALE More...

Functions/Subroutines
subroutine, public	parentatmossetupscale (dims)
	Atmos Setup. More...
subroutine, public	parentatmosopenscale (lon_org, lat_org, cz_org, basename_org, dims)
subroutine, public	parentatmosinputscale (velz_org, velx_org, vely_org, pres_org, dens_org, pott_org, qv_org, qtrc_org, cz_org, basename_org, sfc_diagnoses, same_mptype, dims, it)
subroutine, public	parentlandsetupscale (ldims)
	Land Setup. More...
subroutine, public	parentlandinputscale (tg_org, strg_org, lst_org, ust_org, albg_org, topo_org, lmask_org, llon_org, llat_org, lz_org, basename_land, ldims, use_file_landwater, it)
subroutine, public	parentoceansetupscale (odims)
	Ocean Setup. More...
subroutine, public	parentoceanopenscale (olon_org, olat_org, omask_org, basename_ocean)
subroutine, public	parentoceaninputscale (tw_org, sst_org, albw_org, z0w_org, omask_org, basename_ocean, it)

Detailed Description

module REAL input SCALE

Description

read data from SCALE file for real atmospheric simulations

Author

Team SCALE

Function/Subroutine Documentation

parentatmossetupscale()

subroutine, public mod_realinput_scale::parentatmossetupscale

(

integer, dimension(6), intent(out)

dims

)

Atmos Setup.

Definition at line 64 of file mod_realinput_scale.F90.

    use scale_comm_cartesc_nest, only: &
       parent_kmax,                                     &
       parent_imax,                                     &
       parent_jmax,                                     &
       nest_tile_num_x => comm_cartesc_nest_tile_num_x, &
       nest_tile_num_y => comm_cartesc_nest_tile_num_y
    implicit none
 
    integer, intent(out) :: dims(6)
    !---------------------------------------------------------------------------
 
    log_info("ParentAtmosSetupSCALE",*) 'Real Case/Atmos Input File Type: SCALE-RM'
 
    ! full level
    dims(1) = parent_kmax(handle)
    dims(2) = parent_imax(handle) * nest_tile_num_x
    dims(3) = parent_jmax(handle) * nest_tile_num_y
    ! half level
    dims(4) = dims(1)
    dims(5) = dims(2)
    dims(6) = dims(3)
 
    allocate( read2d(        parent_imax(handle),parent_jmax(handle)) )
    allocate( read3d(dims(1),parent_imax(handle),parent_jmax(handle)) )
 
    allocate( rotc_cos(dims(2),dims(3)), rotc_sin(dims(2),dims(3)) )
 
    return

References scale_comm_cartesc_nest::comm_cartesc_nest_tile_num_x, scale_comm_cartesc_nest::comm_cartesc_nest_tile_num_y, scale_comm_cartesc_nest::parent_imax, scale_comm_cartesc_nest::parent_jmax, and scale_comm_cartesc_nest::parent_kmax.

Referenced by mod_realinput::realinput_surface().

Here is the caller graph for this function:

parentatmosopenscale()

subroutine, public mod_realinput_scale::parentatmosopenscale	(	real(rp), dimension(:,:), intent(out)	lon_org,
		real(rp), dimension(:,:), intent(out)	lat_org,
		real(rp), dimension (:,:,:), intent(out)	cz_org,
		character(len=*), intent(in)	basename_org,
		integer, dimension(6), intent(in)	dims
	)

Definition at line 101 of file mod_realinput_scale.F90.

    use scale_const, only: &
       d2r => const_d2r
    use scale_prc, only: &
       prc_abort
    use scale_comm_cartesc_nest, only: &
       parent_imax,                                     &
       parent_jmax,                                     &
       nest_tile_num_x => comm_cartesc_nest_tile_num_x, &
       nest_tile_id    => comm_cartesc_nest_tile_id
    use scale_file, only: &
       file_open, &
       file_get_attribute
    use scale_file_cartesc, only: &
       file_cartesc_read
    use scale_mapprojection, only: &
       mappinginfo, &
       mappingparam, &
       mapprojection_get_param, &
       mapprojection_rotcoef
    implicit none
 
    real(RP),         intent(out) :: lon_org(:,:)
    real(RP),         intent(out) :: lat_org(:,:)
    real(RP),         intent(out) :: cz_org (:,:,:)
    character(len=*), intent(in)  :: basename_org
    integer,          intent(in)  :: dims(6)
 
    ! mapping information
    character(len=H_SHORT) :: vname
    type(mappinginfo)  :: mapping_info
    type(mappingparam) :: mapping_param
 
    integer :: rank
    integer :: xloc, yloc
    integer :: xs, xe, ys, ye
 
    logical :: existed
 
    integer :: fid
    integer :: i, j
    !---------------------------------------------------------------------------
 
    do i = 1, size( nest_tile_id(:) )
       ! read data from split files
       rank = nest_tile_id(i)
 
       xloc = mod( i-1, nest_tile_num_x ) + 1
       yloc = int( real(i-1) / real(nest_tile_num_x) ) + 1
 
       xs = parent_imax(handle) * (xloc-1) + 1
       xe = parent_imax(handle) * xloc
       ys = parent_jmax(handle) * (yloc-1) + 1
       ye = parent_jmax(handle) * yloc
 
       call file_open( basename_org,                  & ! (in)
                       fid,                           & ! (out)
                       aggregate=.false., rankid=rank ) ! (in)
 
       call file_cartesc_read( fid, "lon", read2d(:,:) )
       lon_org(xs:xe,ys:ye) = read2d(:,:) * d2r
 
       call file_cartesc_read( fid, "lat", read2d(:,:) )
       lat_org(xs:xe,ys:ye) = read2d(:,:) * d2r
 
       call file_cartesc_read( fid, "height", read3d(:,:,:) )
       cz_org(3:dims(1)+2,xs:xe,ys:ye) = read3d(:,:,:)
 
       call file_cartesc_read( fid, "topo", read2d(:,:), existed=existed )
       cz_org(2,xs:xe,ys:ye) = read2d(:,:)
 
       if ( i == 1 ) then
          call file_get_attribute( fid, "DENS", "grid_mapping", vname, existed=existed )
          if ( existed ) then
             call file_get_attribute( fid, vname, "grid_mapping_name", mapping_info%mapping_name )
             call file_get_attribute( fid, vname, "false_easting", mapping_info%false_easting, existed=existed )
             call file_get_attribute( fid, vname, "false_northing", mapping_info%false_northing, existed=existed )
             call file_get_attribute( fid, vname, "longitude_of_central_meridian", mapping_info%longitude_of_central_meridian, existed=existed )
             call file_get_attribute( fid, vname, "longitude_of_projection_origin", mapping_info%longitude_of_projection_origin, existed=existed )
             call file_get_attribute( fid, vname, "latitude_of_projection_origin", mapping_info%latitude_of_projection_origin, existed=existed )
             call file_get_attribute( fid, vname, "straight_vertical_longitude_from_pole", mapping_info%straight_vertical_longitude_from_pole, existed=existed )
             call file_get_attribute( fid, vname, "standard_parallel", mapping_info%standard_parallel(:), existed=existed )
             if ( .not. existed ) then
                call file_get_attribute( fid, vname, "standard_parallel", mapping_info%standard_parallel(1), existed=existed )
             end if
             call file_get_attribute( fid, vname, "rotation", mapping_info%rotation, existed=existed )
             if ( .not. existed ) mapping_info%rotation = 0.0_dp
          else
             mapping_info%mapping_name = ""
          end if
       end if
 
    end do
 
    cz_org(1,:,:) = 0.0_rp
 
    if ( mapping_info%mapping_name .ne. "" ) then
       call mapprojection_get_param( mapping_info, & ! (in)
                                     mapping_param ) ! (out)
       call mapprojection_rotcoef( dims(2), 1, dims(2), dims(3), 1, dims(3), &
                                   lon_org(:,:), lat_org(:,:),  & ! (in)
                                   mapping_info%mapping_name,   & ! (in)
                                   mapping_param,               & ! (in)
                                   rotc_cos(:,:), rotc_sin(:,:) ) ! (out)
    else
!OCL XFILL
       !$omp parallel do
       do j = 1, dims(3)
       do i = 1, dims(2)
          rotc_cos(i,j) = 1.0_rp
          rotc_sin(i,j) = 0.0_rp
       end do
       end do
    end if
 
    return

References scale_comm_cartesc_nest::comm_cartesc_nest_tile_id, scale_comm_cartesc_nest::comm_cartesc_nest_tile_num_x, scale_const::const_d2r, scale_file::file_open(), scale_mapprojection::mapprojection_get_param(), scale_comm_cartesc_nest::parent_imax, scale_comm_cartesc_nest::parent_jmax, and scale_prc::prc_abort().

Referenced by mod_realinput::realinput_surface().

Here is the call graph for this function:

Here is the caller graph for this function:

parentatmosinputscale()

subroutine, public mod_realinput_scale::parentatmosinputscale	(	real(rp), dimension(:,:,:), intent(out)	velz_org,
		real(rp), dimension(:,:,:), intent(out)	velx_org,
		real(rp), dimension(:,:,:), intent(out)	vely_org,
		real(rp), dimension(:,:,:), intent(out)	pres_org,
		real(rp), dimension(:,:,:), intent(out)	dens_org,
		real(rp), dimension(:,:,:), intent(out)	pott_org,
		real(rp), dimension (:,:,:), intent(out)	qv_org,
		real(rp), dimension(:,:,:,:), intent(out)	qtrc_org,
		real(rp), dimension (:,:,:), intent(in)	cz_org,
		character(len=*), intent(in)	basename_org,
		logical, intent(in)	sfc_diagnoses,
		logical, intent(in)	same_mptype,
		integer, dimension(6), intent(in)	dims,
		integer, intent(in)	it
	)

Definition at line 234 of file mod_realinput_scale.F90.

    use scale_const, only: &
       undef => const_undef, &
       p00   => const_pre00, &
       cpdry => const_cpdry, &
       rdry  => const_rdry,  &
       grav  => const_grav,  &
       laps  => const_laps
    use scale_comm_cartesc, only: &
       comm_vars8, &
       comm_wait
    use scale_comm_cartesc_nest, only: &
       parent_imax,                                     &
       parent_jmax,                                     &
       nest_tile_num_x => comm_cartesc_nest_tile_num_x, &
       nest_tile_id    => comm_cartesc_nest_tile_id
    use scale_file, only: &
       file_open
    use scale_file_cartesc, only: &
       file_cartesc_read
    use scale_atmos_grid_cartesc_metric, only: &
       rotc => atmos_grid_cartesc_metric_rotc
    use scale_atmos_thermodyn, only: &
       thermodyn_specific_heat  => atmos_thermodyn_specific_heat, &
       thermodyn_rhot2temp_pres => atmos_thermodyn_rhot2temp_pres
    use scale_atmos_hydrometeor, only: &
       n_hyd,    &
       hyd_name, &
       num_name
    use mod_atmos_phy_mp_driver, only: &
       atmos_phy_mp_driver_qhyd2qtrc
    use mod_atmos_phy_mp_vars, only: &
       qa_mp, &
       qs_mp, &
       qe_mp
    implicit none
 
    real(RP),         intent(out) :: velz_org(:,:,:)
    real(RP),         intent(out) :: velx_org(:,:,:)
    real(RP),         intent(out) :: vely_org(:,:,:)
    real(RP),         intent(out) :: pres_org(:,:,:)
    real(RP),         intent(out) :: dens_org(:,:,:)
    real(RP),         intent(out) :: pott_org(:,:,:)
    real(RP),         intent(out) :: qv_org  (:,:,:)
    real(RP),         intent(out) :: qtrc_org(:,:,:,:)
    real(RP),         intent(in)  :: cz_org  (:,:,:)
    character(len=*), intent(in)  :: basename_org
    logical,          intent(in)  :: sfc_diagnoses
    logical,          intent(in)  :: same_mptype
    integer,          intent(in)  :: dims(6)
    integer,          intent(in)  :: it
 
    real(RP) :: momz_org(dims(1)+2,dims(2),dims(3))
    real(RP) :: momx_org(dims(1)+2,dims(2),dims(3))
    real(RP) :: momy_org(dims(1)+2,dims(2),dims(3))
    real(RP) :: rhot_org(dims(1)+2,dims(2),dims(3))
    real(RP) :: tsfc_org(          dims(2),dims(3))
    real(RP) :: qhyd_org(dims(1)+2,dims(2),dims(3),N_HYD)
    real(RP) :: qnum_org(dims(1)+2,dims(2),dims(3),N_HYD)
    real(RP) :: qdry    (dims(1)+2)
    real(RP) :: Rtot    (dims(1)+2)
    real(RP) :: CVtot   (dims(1)+2)
    real(RP) :: CPtot   (dims(1)+2)
    real(RP) :: temp_org(dims(1)+2)
    real(RP) :: dz
 
    real(RP) :: u, v
 
    integer  :: xs, xe
    integer  :: ys, ye
    integer  :: nx, ny
    integer  :: xloc, yloc
    integer  :: rank
 
    logical :: qnum_flag
 
    integer  :: fid
    logical  :: existed, existed_t2, existed_mslp
    integer  :: k, i, j, iq
    !---------------------------------------------------------------------------
 
    qnum_flag = .false.
 
    do i = 1, size( nest_tile_id(:) )
       ! read data from split files
       rank = nest_tile_id(i)
 
       xloc = mod( i-1, nest_tile_num_x ) + 1
       yloc = int( real(i-1) / real(nest_tile_num_x) ) + 1
 
       xs = parent_imax(handle) * (xloc-1) + 1
       xe = parent_imax(handle) * xloc
       ys = parent_jmax(handle) * (yloc-1) + 1
       ye = parent_jmax(handle) * yloc
 
       nx = xe - xs + 1
       ny = ye - ys + 1
 
       call file_open( basename_org,                  & ! (in)
                       fid,                           & ! (out)
                       aggregate=.false., rankid=rank ) ! (in)
 
       if ( sfc_diagnoses ) then
          call file_cartesc_read( fid, "T2", read2d(:,:), step=it, existed=existed_t2 )
          if ( existed_t2 ) then
!OCL XFILL
             tsfc_org(xs:xe,ys:ye) = read2d(:,:)
          end if
 
          call file_cartesc_read( fid, "MSLP", read2d(:,:), step=it, existed=existed_mslp )
          if ( existed_mslp ) then
!OCL XFILL
             pres_org(1,xs:xe,ys:ye) = read2d(:,:)
          end if
       end if
 
       call file_cartesc_read( fid, "DENS", read3d(:,:,:), step=it )
!OCL XFILL
       dens_org(3:dims(1)+2,xs:xe,ys:ye) = read3d(:,:,:)
 
       call file_cartesc_read( fid, "MOMZ", read3d(:,:,:), step=it )
!OCL XFILL
       momz_org(3:dims(1)+2,xs:xe,ys:ye) = read3d(:,:,:)
 
       call file_cartesc_read( fid, "MOMX", read3d(:,:,:), step=it )
!OCL XFILL
       momx_org(3:dims(1)+2,xs:xe,ys:ye) = read3d(:,:,:)
 
       call file_cartesc_read( fid, "MOMY", read3d(:,:,:), step=it )
!OCL XFILL
       momy_org(3:dims(1)+2,xs:xe,ys:ye) = read3d(:,:,:)
 
 
       call file_cartesc_read( fid, "RHOT", read3d(:,:,:), step=it )
!OCL XFILL
       rhot_org(3:dims(1)+2,xs:xe,ys:ye) = read3d(:,:,:)
 
 
       if ( same_mptype ) then
 
          do iq = qs_mp, qe_mp
             call file_cartesc_read( fid, tracer_name(iq), read3d(:,:,:), step=it )
!OCL XFILL
             qtrc_org(3:dims(1)+2,xs:xe,ys:ye,iq) = read3d(:,:,:)
!OCL XFILL
             qtrc_org(2,xs:xe,ys:ye,iq) = qtrc_org(3,xs:xe,ys:ye,iq)
!OCL XFILL
             qtrc_org(1,xs:xe,ys:ye,iq) = qtrc_org(3,xs:xe,ys:ye,iq)
          enddo
 
       else
 
!OCL XFILL
          do iq = 1, n_hyd
             qhyd_org(:,xs:xe,ys:ye,iq) = 0.0_rp
             qnum_org(:,xs:xe,ys:ye,iq) = 0.0_rp
          end do
 
          call file_cartesc_read( fid, "QV", read3d(:,:,:), step=it, existed=existed )
          if ( existed ) then
!OCL XFILL
             qv_org(3:dims(1)+2,xs:xe,ys:ye) = read3d(:,:,:)
!OCL XFILL
             qv_org(2,xs:xe,ys:ye) = qv_org(3,xs:xe,ys:ye)
!OCL XFILL
             qv_org(1,xs:xe,ys:ye) = qv_org(3,xs:xe,ys:ye)
          else
!OCL XFILL
             qv_org(:,:,:) = 0.0_rp
          end if
 
          do iq = 1, n_hyd
             call file_cartesc_read( fid, hyd_name(iq), read3d(:,:,:), step=it, existed=existed )
             if ( existed ) then
!OCL XFILL
                qhyd_org(3:dims(1)+2,xs:xe,ys:ye,iq) = read3d(:,:,:)
!OCL XFILL
                qhyd_org(2,xs:xe,ys:ye,iq) = qhyd_org(3,xs:xe,ys:ye,iq)
!OCL XFILL
                qhyd_org(1,xs:xe,ys:ye,iq) = qhyd_org(3,xs:xe,ys:ye,iq)
             else
!OCL XFILL
                qhyd_org(:,:,:,iq) = 0.0_rp
             end if
 
             ! number density
             call file_cartesc_read( fid, num_name(iq), read3d(:,:,:), step=it, existed=existed )
             if ( existed ) then
!OCL XFILL
                qnum_org(3:dims(1)+2,xs:xe,ys:ye,iq) = read3d(:,:,:)
!OCL XFILL
                qnum_org(2,xs:xe,ys:ye,iq) = qnum_org(3,xs:xe,ys:ye,iq)
!OCL XFILL
                qnum_org(1,xs:xe,ys:ye,iq) = qnum_org(3,xs:xe,ys:ye,iq)
                qnum_flag = .true.
             else
!OCL XFILL
                qnum_org(:,:,:,iq) = undef
             end if
          end do
 
       endif
 
!       call FILE_CARTESC_read( fid, "Q2", read2D(:,:), step=it )
!       qv_org(2,xs:xe,ys:ye) = read2D(:,:)
 
    end do
 
    ! convert from momentum to velocity
!OCL XFILL
    do j = 1, dims(3)
    do i = 1, dims(2)
    do k = 4, dims(1)+2
       velz_org(k,i,j) = ( momz_org(k-1,i,j) + momz_org(k,i,j) ) / dens_org(k,i,j) * 0.5_rp
    end do
    end do
    end do
!OCL XFILL
    do j = 1, dims(3)
    do i = 1, dims(2)
       velz_org(1:3    ,i,j) = 0.0_rp
       velz_org(dims(1)+2,i,j) = 0.0_rp
    end do
    end do
 
    ! convert from momentum to velocity
!OCL XFILL
    do j = 1, dims(3)
    do i = 2, dims(2)
    do k = 3, dims(1)+2
       velx_org(k,i,j) = ( momx_org(k,i-1,j) + momx_org(k,i,j) ) / dens_org(k,i,j) * 0.5_rp
    end do
    end do
    end do
!OCL XFILL
    do j = 1, dims(3)
    do k = 3, dims(1)+2
       velx_org(k,1,j) = momx_org(k,1,j) / dens_org(k,1,j)
    end do
    end do
    if ( sfc_diagnoses ) then
!OCL XFILL
       velx_org(1:2,:,:) = 0.0_rp
    else
!OCL XFILL
       velx_org(1:2,:,:) = undef
    end if
 
    ! convert from momentum to velocity
!OCL XFILL
    do j = 2, dims(3)
    do i = 1, dims(2)
    do k = 3, dims(1)+2
       vely_org(k,i,j) = ( momy_org(k,i,j-1) + momy_org(k,i,j) ) / dens_org(k,i,j) * 0.5_rp
    end do
    end do
    end do
!OCL XFILL
    do i = 1, dims(2)
    do k = 3, dims(1)+2
       vely_org(k,i,1) = momy_org(k,i,1) / dens_org(k,i,1)
    end do
    end do
    if ( sfc_diagnoses ) then
!OCL XFILL
       vely_org(1:2,:,:) = 0.0_rp
    else
!OCL XFILL
       vely_org(1:2,:,:) = undef
    end if
 
 
    ! rotation
    !$omp parallel do &
    !$omp private(u,v)
    do j = 1, dims(3)
    do i = 1, dims(2)
    do k = 1, dims(1)+2
       u = velx_org(k,i,j)
       v = vely_org(k,i,j)
       if ( u .ne. undef .and. v .ne.undef ) then
          velx_org(k,i,j) = u * rotc_cos(i,j) - v * rotc_sin(i,j)
          vely_org(k,i,j) = u * rotc_sin(i,j) + v * rotc_cos(i,j)
       else
          velx_org(k,i,j) = undef
          vely_org(k,i,j) = undef
       end if
    end do
    end do
    end do
 
 
    do iq = 1, size(qtrc_org,4)
       if ( iq >= qs_mp .and. iq <= qe_mp ) cycle
!OCL XFILL
       qtrc_org(:,:,:,iq) = 0.0_rp
    end do
 
    if ( qa_mp > 0 .AND. .NOT. same_mptype ) then
       if ( qnum_flag ) then
          call atmos_phy_mp_driver_qhyd2qtrc( dims(1)+2, 1, dims(1)+2, dims(2), 1, dims(2), dims(3), 1, dims(3), &
                                              qv_org(:,:,:), qhyd_org(:,:,:,:), & ! [IN]
                                              qtrc_org(:,:,:,qs_mp:qe_mp),      & ! [OUT]
                                              qnum=qnum_org(:,:,:,:)            ) ! [IN]
       else
          call atmos_phy_mp_driver_qhyd2qtrc( dims(1)+2, 1, dims(1)+2, dims(2), 1, dims(2), dims(3), 1, dims(3), &
                                              qv_org(:,:,:), qhyd_org(:,:,:,:), & ! [IN]
                                              qtrc_org(:,:,:,qs_mp:qe_mp)       ) ! [OUT]
       end if
    end if
 
    !$omp parallel do default(none) &
    !$omp private(dz,temp_org,qdry,rtot,cvtot,cptot) &
    !$omp shared(dims,QA,UNDEF,LAPS,P00,Rdry,CPdry,GRAV,TRACER_MASS,TRACER_R,TRACER_CV,TRACER_CP, &
    !$omp        dens_org,rhot_org,qtrc_org,pres_org,pott_org,tsfc_org,cz_org, &
    !$omp        existed_t2,existed_mslp,sfc_diagnoses)
    do j = 1, dims(3)
    do i = 1, dims(2)
 
       ! diagnose temp and pres
       call thermodyn_specific_heat( dims(1)+2, 3, dims(1)+2, qa,                             & ! [IN]
                                     qtrc_org(:,i,j,:),                                       & ! [IN]
                                     tracer_mass(:), tracer_r(:), tracer_cv(:), tracer_cp(:), & ! [IN]
                                     qdry(:), rtot(:), cvtot(:), cptot(:)                     ) ! [OUT]
 
       call thermodyn_rhot2temp_pres( dims(1), 1, dims(1),                                       & ! [IN]
                                      dens_org(3:dims(1)+2,i,j), rhot_org(3:dims(1)+2,i,j),      & ! [IN]
                                      rtot(3:dims(1)+2), cvtot(3:dims(1)+2), cptot(3:dims(1)+2), & ! [IN]
                                      temp_org(3:dims(1)+2), pres_org(3:dims(1)+2,i,j)           ) ! [OUT]
 
!OCL XFILL
       do k = 3, dims(1)+2
          pott_org(k,i,j) = rhot_org(k,i,j) / dens_org(k,i,j)
       end do
 
       if ( sfc_diagnoses ) then
          ! at the surface
          dz = cz_org(3,i,j) - cz_org(2,i,j)
 
          if ( .not. existed_t2 ) then
             tsfc_org(i,j) = temp_org(3) + laps * dz
          end if
          dens_org(2,i,j) = ( pres_org(3,i,j) + grav * dens_org(3,i,j) * dz * 0.5_rp ) &
                          / ( rdry * tsfc_org(i,j) - grav * dz * 0.5_rp )
          pres_org(2,i,j) = dens_org(2,i,j) * rdry * tsfc_org(i,j)
          pott_org(2,i,j) = tsfc_org(i,j) * ( p00 / pres_org(2,i,j) )**(rdry/cpdry)
 
          ! at the sea-level
          temp_org(1) = tsfc_org(i,j) + laps * cz_org(2,i,j)
          if ( existed_mslp ) then
             pott_org(1,i,j) = temp_org(1) * ( p00 / pres_org(1,i,j) )**(rdry/cpdry)
             dens_org(1,i,j) = pres_org(1,i,j) / ( rdry * temp_org(1) )
          else
             dens_org(1,i,j) = ( pres_org(2,i,j) + grav * dens_org(2,i,j) * cz_org(2,i,j) * 0.5_rp ) &
                             / ( rdry * temp_org(1) - grav * cz_org(2,i,j) * 0.5_rp )
             pres_org(1,i,j) = dens_org(1,i,j) * rdry * temp_org(1)
             pott_org(1,i,j) = temp_org(1) * ( p00 / pres_org(1,i,j) )**(rdry/cpdry)
          end if
       else
          dens_org(1:2,i,j) = undef
          pres_org(1:2,i,j) = undef
          pott_org(1:2,i,j) = undef
       end if
 
    end do
    end do
 
    return

References scale_atmos_grid_cartesc_metric::atmos_grid_cartesc_metric_rotc, mod_atmos_phy_mp_driver::atmos_phy_mp_driver_qhyd2qtrc(), scale_comm_cartesc_nest::comm_cartesc_nest_tile_id, scale_comm_cartesc_nest::comm_cartesc_nest_tile_num_x, scale_const::const_cpdry, scale_const::const_grav, scale_const::const_laps, scale_const::const_pre00, scale_const::const_rdry, scale_const::const_undef, scale_file::file_open(), scale_atmos_hydrometeor::hyd_name, scale_tracer::k, scale_atmos_hydrometeor::n_hyd, scale_atmos_hydrometeor::num_name, scale_comm_cartesc_nest::parent_imax, scale_comm_cartesc_nest::parent_jmax, scale_tracer::qa, mod_atmos_phy_mp_vars::qa_mp, mod_atmos_phy_mp_vars::qe_mp, mod_atmos_phy_mp_vars::qs_mp, scale_tracer::tracer_cp, scale_tracer::tracer_cv, scale_tracer::tracer_mass, scale_tracer::tracer_name, and scale_tracer::tracer_r.

Referenced by mod_realinput::realinput_surface().

Here is the call graph for this function:

Here is the caller graph for this function:

parentlandsetupscale()

subroutine, public mod_realinput_scale::parentlandsetupscale

(

integer, dimension(3), intent(out)

ldims

)

Land Setup.

Definition at line 607 of file mod_realinput_scale.F90.

    use scale_comm_cartesc_nest, only: &
       parent_imax,                                     &
       parent_jmax,                                     &
       parent_lkmax,                                    &
       nest_tile_num_x => comm_cartesc_nest_tile_num_x, &
       nest_tile_num_y => comm_cartesc_nest_tile_num_y
    implicit none
 
    integer, intent(out) :: ldims(3)
    !---------------------------------------------------------------------------
 
    log_info("ParentLandSetupSCALE",*) 'Real Case/Land Input File Type: SCALE-RM'
 
    ldims(1) = parent_lkmax(handle)
    ldims(2) = parent_imax(handle) * nest_tile_num_x
    ldims(3) = parent_jmax(handle) * nest_tile_num_y
 
    if ( .NOT. allocated(read2d) ) then
       allocate( read2d(parent_imax(handle),parent_jmax(handle)) )
    endif
 
    allocate( read3dl(ldims(1),parent_imax(handle),parent_jmax(handle)) )
 
    return

References scale_comm_cartesc_nest::comm_cartesc_nest_tile_num_x, scale_comm_cartesc_nest::comm_cartesc_nest_tile_num_y, scale_comm_cartesc_nest::parent_imax, scale_comm_cartesc_nest::parent_jmax, and scale_comm_cartesc_nest::parent_lkmax.

Referenced by mod_realinput::realinput_surface().

Here is the caller graph for this function:

parentlandinputscale()

subroutine, public mod_realinput_scale::parentlandinputscale	(	real(rp), dimension (:,:,:), intent(out)	tg_org,
		real(rp), dimension (:,:,:), intent(out)	strg_org,
		real(rp), dimension (:,:), intent(out)	lst_org,
		real(rp), dimension (:,:), intent(out)	ust_org,
		real(rp), dimension (:,:,:,:), intent(out)	albg_org,
		real(rp), dimension (:,:), intent(out)	topo_org,
		real(rp), dimension(:,:), intent(out)	lmask_org,
		real(rp), dimension (:,:), intent(out)	llon_org,
		real(rp), dimension (:,:), intent(out)	llat_org,
		real(rp), dimension (:), intent(out)	lz_org,
		character(len=*), intent(in)	basename_land,
		integer, dimension(3), intent(in)	ldims,
		logical, intent(in)	use_file_landwater,
		integer, intent(in)	it
	)

Definition at line 649 of file mod_realinput_scale.F90.

    use scale_const, only: &
       undef => const_undef, &
       d2r   => const_d2r
    use scale_comm_cartesc_nest, only: &
       parent_imax,                                     &
       parent_jmax,                                     &
       nest_tile_num_x => comm_cartesc_nest_tile_num_x, &
       nest_tile_id    => comm_cartesc_nest_tile_id
    use scale_file, only: &
       file_open, &
       file_read, &
       file_get_datainfo
    use scale_file_cartesc, only: &
       file_cartesc_read
    implicit none
 
    real(RP),         intent(out) :: tg_org   (:,:,:)
    real(RP),         intent(out) :: strg_org (:,:,:)
    real(RP),         intent(out) :: lst_org  (:,:)
    real(RP),         intent(out) :: ust_org  (:,:)
    real(RP),         intent(out) :: albg_org (:,:,:,:)
    real(RP),         intent(out) :: topo_org (:,:)
    real(RP),         intent(out) :: lmask_org(:,:)
    real(RP),         intent(out) :: llon_org (:,:)
    real(RP),         intent(out) :: llat_org (:,:)
    real(RP),         intent(out) :: lz_org   (:)
    character(len=*), intent(in)  :: basename_land
    integer,          intent(in)  :: ldims(3)
    logical,          intent(in)  :: use_file_landwater ! use land water data from files
    integer,          intent(in)  :: it
 
    integer :: rank
    integer :: fid
    integer :: xloc, yloc
    integer :: xs, xe
    integer :: ys, ye
    logical :: existed
 
    integer :: i
    !---------------------------------------------------------------------------
 
    do i = 1, size( nest_tile_id(:) )
       ! read data from split files
       rank = nest_tile_id(i)
 
       xloc = mod( i-1, nest_tile_num_x ) + 1
       yloc = int( real(i-1) / real(nest_tile_num_x) ) + 1
 
       xs = parent_imax(handle) * (xloc-1) + 1
       xe = parent_imax(handle) * xloc
       ys = parent_jmax(handle) * (yloc-1) + 1
       ye = parent_jmax(handle) * yloc
 
       call file_open( basename_land,                 & ! (in)
                       fid,                           & ! (out)
                       aggregate=.false., rankid=rank ) ! (in)
 
       call file_cartesc_read( fid, "LAND_TEMP",  read3dl(:,:,:), step=it )
       tg_org(1:ldims(1),xs:xe,ys:ye) = read3dl(:,:,:)
 
       if( use_file_landwater )then
          call file_cartesc_read( fid, "LAND_WATER", read3dl(:,:,:), step=it )
          strg_org(1:ldims(1),xs:xe,ys:ye) = read3dl(:,:,:)
       endif
 
       call file_cartesc_read( fid, "lon", read2d(:,:) )
       llon_org(xs:xe,ys:ye)  = read2d(:,:) * d2r
 
       call file_cartesc_read( fid, "lat", read2d(:,:) )
       llat_org(xs:xe,ys:ye)  = read2d(:,:) * d2r
 
       call file_cartesc_read( fid, "LAND_SFC_TEMP", read2d(:,:), step=it )
       lst_org(xs:xe,ys:ye) = read2d(:,:)
 
       call file_get_datainfo( fid, "URBAN_SFC_TEMP", istep=it, existed=existed )
       if ( existed ) then
          call file_cartesc_read( fid, "URBAN_SFC_TEMP", read2d(:,:), step=it )
          ust_org(xs:xe,ys:ye) = read2d(:,:)
       else
          ust_org(xs:xe,ys:ye) = undef
       end if
 
       call file_cartesc_read( fid, "LAND_SFC_ALB_IR_dir", read2d(:,:), step=it )
       albg_org(xs:xe,ys:ye,i_r_direct ,i_r_ir ) = read2d(:,:)
 
       call file_cartesc_read( fid, "LAND_SFC_ALB_IR_dif", read2d(:,:), step=it )
       albg_org(xs:xe,ys:ye,i_r_diffuse,i_r_ir ) = read2d(:,:)
 
       call file_cartesc_read( fid, "LAND_SFC_ALB_NIR_dir", read2d(:,:), step=it )
       albg_org(xs:xe,ys:ye,i_r_direct ,i_r_nir) = read2d(:,:)
 
       call file_cartesc_read( fid, "LAND_SFC_ALB_NIR_dif", read2d(:,:), step=it )
       albg_org(xs:xe,ys:ye,i_r_diffuse,i_r_nir) = read2d(:,:)
 
       call file_cartesc_read( fid, "LAND_SFC_ALB_VIS_dir", read2d(:,:), step=it )
       albg_org(xs:xe,ys:ye,i_r_direct ,i_r_vis) = read2d(:,:)
 
       call file_cartesc_read( fid, "LAND_SFC_ALB_VIS_dif", read2d(:,:), step=it )
       albg_org(xs:xe,ys:ye,i_r_diffuse,i_r_vis) = read2d(:,:)
 
       call file_cartesc_read( fid, "topo", read2d(:,:) )
       topo_org(xs:xe,ys:ye) = read2d(:,:)
 
       call file_cartesc_read( fid, "lsmask", read2d(:,:) )
       lmask_org(xs:xe,ys:ye) = read2d(:,:)
 
    end do
 
    call file_read( fid, "lz", lz_org(:) )
 
    return

References scale_comm_cartesc_nest::comm_cartesc_nest_tile_id, scale_comm_cartesc_nest::comm_cartesc_nest_tile_num_x, scale_const::const_d2r, scale_const::const_undef, scale_file::file_open(), scale_cpl_sfc_index::i_r_diffuse, scale_cpl_sfc_index::i_r_direct, scale_cpl_sfc_index::i_r_ir, scale_cpl_sfc_index::i_r_nir, scale_cpl_sfc_index::i_r_vis, scale_comm_cartesc_nest::parent_imax, and scale_comm_cartesc_nest::parent_jmax.

Referenced by mod_realinput::realinput_surface().

Here is the call graph for this function:

Here is the caller graph for this function:

parentoceansetupscale()

subroutine, public mod_realinput_scale::parentoceansetupscale

(

integer, dimension(2), intent(out)

odims

)

Ocean Setup.

Definition at line 766 of file mod_realinput_scale.F90.

    use scale_comm_cartesc_nest, only: &
       parent_imax,                                     &
       parent_jmax,                                     &
       nest_tile_num_x => comm_cartesc_nest_tile_num_x, &
       nest_tile_num_y => comm_cartesc_nest_tile_num_y
    implicit none
 
    integer, intent(out) :: odims(2)
    !---------------------------------------------------------------------------
 
    log_info("ParentOceanSetupSCALE",*) 'Real Case/Ocean Input File Type: SCALE-RM'
 
    odims(1) = parent_imax(handle) * nest_tile_num_x
    odims(2) = parent_jmax(handle) * nest_tile_num_y
 
    if ( .NOT. allocated(read2d) ) then
       allocate( read2d(parent_imax(handle),parent_jmax(handle)) )
    end if
 
    allocate( read3do(1,parent_imax(handle),parent_jmax(handle)) )
 
    return

References scale_comm_cartesc_nest::comm_cartesc_nest_tile_num_x, scale_comm_cartesc_nest::comm_cartesc_nest_tile_num_y, scale_comm_cartesc_nest::parent_imax, and scale_comm_cartesc_nest::parent_jmax.

Referenced by mod_realinput::realinput_surface().

Here is the caller graph for this function:

parentoceanopenscale()

subroutine, public mod_realinput_scale::parentoceanopenscale	(	real(rp), dimension (:,:), intent(out)	olon_org,
		real(rp), dimension (:,:), intent(out)	olat_org,
		real(rp), dimension(:,:), intent(out)	omask_org,
		character(len=*), intent(in)	basename_ocean
	)

Definition at line 796 of file mod_realinput_scale.F90.

    use scale_const, only: &
       d2r => const_d2r
    use scale_comm_cartesc_nest, only: &
       parent_imax,                                     &
       parent_jmax,                                     &
       nest_tile_num_x => comm_cartesc_nest_tile_num_x, &
       nest_tile_id    => comm_cartesc_nest_tile_id
    use scale_file, only: &
       file_open
    use scale_file_cartesc, only: &
       file_cartesc_read
    implicit none
 
    real(RP),         intent(out) :: olon_org (:,:)
    real(RP),         intent(out) :: olat_org (:,:)
    real(RP),         intent(out) :: omask_org(:,:)
    character(len=*), intent(in)  :: basename_ocean
 
    integer :: rank
    integer :: xloc, yloc
    integer :: xs, xe, ys, ye
 
    integer :: fid
    integer :: i
    !---------------------------------------------------------------------------
 
    do i = 1, size( nest_tile_id(:) )
       ! read data from split files
       rank = nest_tile_id(i)
 
       xloc = mod( i-1, nest_tile_num_x ) + 1
       yloc = int( real(i-1) / real(nest_tile_num_x) ) + 1
 
       xs = parent_imax(handle) * (xloc-1) + 1
       xe = parent_imax(handle) * xloc
       ys = parent_jmax(handle) * (yloc-1) + 1
       ye = parent_jmax(handle) * yloc
 
       call file_open( basename_ocean,                & ! (in)
                       fid,                           & ! (out)
                       aggregate=.false., rankid=rank ) ! (in)
 
       call file_cartesc_read( fid, "lon", read2d(:,:) )
       olon_org(xs:xe,ys:ye) = read2d(:,:) * d2r
 
       call file_cartesc_read( fid, "lat", read2d(:,:) )
       olat_org(xs:xe,ys:ye) = read2d(:,:) * d2r
 
       call file_cartesc_read( fid, "lsmask", read2d(:,:) )
       omask_org(xs:xe,ys:ye) = read2d(:,:)
 
    end do
 
    return

References scale_comm_cartesc_nest::comm_cartesc_nest_tile_id, scale_comm_cartesc_nest::comm_cartesc_nest_tile_num_x, scale_const::const_d2r, scale_file::file_open(), scale_comm_cartesc_nest::parent_imax, and scale_comm_cartesc_nest::parent_jmax.

Referenced by mod_realinput::realinput_surface().

Here is the call graph for this function:

Here is the caller graph for this function:

parentoceaninputscale()

subroutine, public mod_realinput_scale::parentoceaninputscale	(	real(rp), dimension (:,:), intent(out)	tw_org,
		real(rp), dimension (:,:), intent(out)	sst_org,
		real(rp), dimension (:,:,:,:), intent(out)	albw_org,
		real(rp), dimension (:,:), intent(out)	z0w_org,
		real(rp), dimension(:,:), intent(out)	omask_org,
		character(len=*), intent(in)	basename_ocean,
		integer, intent(in)	it
	)

Definition at line 861 of file mod_realinput_scale.F90.

    use scale_comm_cartesc_nest, only: &
       parent_imax,                                     &
       parent_jmax,                                     &
       nest_tile_num_x => comm_cartesc_nest_tile_num_x, &
       nest_tile_id    => comm_cartesc_nest_tile_id
    use scale_file, only: &
       file_open, &
       file_get_datainfo
    use scale_file_cartesc, only: &
       file_cartesc_read
    implicit none
 
    real(RP),         intent(out) :: tw_org   (:,:)
    real(RP),         intent(out) :: sst_org  (:,:)
    real(RP),         intent(out) :: albw_org (:,:,:,:)
    real(RP),         intent(out) :: z0w_org  (:,:)
    real(RP),         intent(out) :: omask_org(:,:)
    character(len=*), intent(in)  :: basename_ocean
    integer,          intent(in)  :: it
 
    integer :: rank
 
    integer :: fid
    integer :: ndim
    integer :: i
    integer :: xloc, yloc
    integer :: xs, xe
    integer :: ys, ye
    !---------------------------------------------------------------------------
 
    do i = 1, size( nest_tile_id(:) )
       ! read data from split files
       rank = nest_tile_id(i)
 
       xloc = mod( i-1, nest_tile_num_x ) + 1
       yloc = int( real(i-1) / real(nest_tile_num_x) ) + 1
 
       xs = parent_imax(handle) * (xloc-1) + 1
       xe = parent_imax(handle) * xloc
       ys = parent_jmax(handle) * (yloc-1) + 1
       ye = parent_jmax(handle) * yloc
 
       call file_open( basename_ocean,                & ! (in)
                       fid,                           & ! (out)
                       aggregate=.false., rankid=rank ) ! (in)
 
       call file_get_datainfo( fid, "OCEAN_TEMP", dim_rank=ndim )
       select case (ndim)
       case ( 2 ) ! for old file
          call file_cartesc_read( fid, "OCEAN_TEMP", read2d(:,:), step=it )
          tw_org(xs:xe,ys:ye) = read2d(:,:)
       case ( 3 )
          call file_cartesc_read( fid, "OCEAN_TEMP", read3do(:,:,:), step=it )
          tw_org(xs:xe,ys:ye) = read3do(1,:,:)
       end select
 
       call file_cartesc_read( fid, "OCEAN_SFC_TEMP", read2d(:,:), step=it )
       sst_org(xs:xe,ys:ye) = read2d(:,:)
 
       call file_cartesc_read( fid, "OCEAN_SFC_ALB_IR_dir", read2d(:,:), step=it )
       albw_org(xs:xe,ys:ye,i_r_direct ,i_r_ir ) = read2d(:,:)
 
       call file_cartesc_read( fid, "OCEAN_SFC_ALB_IR_dif", read2d(:,:), step=it )
       albw_org(xs:xe,ys:ye,i_r_diffuse,i_r_ir ) = read2d(:,:)
 
       call file_cartesc_read( fid, "OCEAN_SFC_ALB_NIR_dir", read2d(:,:), step=it )
       albw_org(xs:xe,ys:ye,i_r_direct ,i_r_nir) = read2d(:,:)
 
       call file_cartesc_read( fid, "OCEAN_SFC_ALB_NIR_dif", read2d(:,:), step=it )
       albw_org(xs:xe,ys:ye,i_r_diffuse,i_r_nir) = read2d(:,:)
 
       call file_cartesc_read( fid, "OCEAN_SFC_ALB_VIS_dir", read2d(:,:), step=it )
       albw_org(xs:xe,ys:ye,i_r_direct ,i_r_vis) = read2d(:,:)
 
       call file_cartesc_read( fid, "OCEAN_SFC_ALB_VIS_dif", read2d(:,:), step=it )
       albw_org(xs:xe,ys:ye,i_r_diffuse,i_r_vis) = read2d(:,:)
 
       call file_cartesc_read( fid, "OCEAN_SFC_Z0M", read2d(:,:), step=it )
       z0w_org(xs:xe,ys:ye) = read2d(:,:)
 
       call file_cartesc_read( fid, "lsmask", read2d(:,:) )
       omask_org(xs:xe,ys:ye) = read2d(:,:)
 
    end do
 
    return

References scale_comm_cartesc_nest::comm_cartesc_nest_tile_id, scale_comm_cartesc_nest::comm_cartesc_nest_tile_num_x, scale_file::file_open(), scale_cpl_sfc_index::i_r_diffuse, scale_cpl_sfc_index::i_r_direct, scale_cpl_sfc_index::i_r_ir, scale_cpl_sfc_index::i_r_nir, scale_cpl_sfc_index::i_r_vis, scale_comm_cartesc_nest::parent_imax, and scale_comm_cartesc_nest::parent_jmax.

Referenced by mod_realinput::realinput_surface().

Here is the call graph for this function:

Here is the caller graph for this function: