mod_cnv2d Module Reference

module Convert 2D data More...

Functions/Subroutines
subroutine, public	cnv2d_setup
	Setup. More...
subroutine, public	cnv2d_tile_init (dtype, dlat, dlon, dir, catalogue, interp_type, interp_level, nmax)
subroutine, public	cnv2d_grads_init (FILE_NAME, VAR_NAME, LAT_NAME, LON_NAME, interp_type, interp_level, search_limit, POSTFIX)
subroutine, public	cnv2d_exec (var, step, min_value, yrevers)

Detailed Description

module Convert 2D data

Description

subroutines for preparing 2D data (convert from external file)

Author

Team SCALE

Function/Subroutine Documentation

cnv2d_setup()

subroutine, public mod_cnv2d::cnv2d_setup

Setup.

Definition at line 92 of file mod_cnv2d.F90.

    use scale_prc, only: &
       prc_abort
    implicit none
 
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("CNV2D_setup",*) 'Setup'
 
    !--- read namelist
!!$    rewind(IO_FID_CONF)
!!$    read(IO_FID_CONF,nml=PARAM_CNV2D,iostat=ierr)
!!$    if( ierr < 0 ) then !--- missing
!!$       LOG_INFO("CNV2D_setup",*) 'Not found namelist. Default used.'
!!$    elseif( ierr > 0 ) then !--- fatal error
!!$       LOG_ERROR("CNV2D_setup",*) 'Not appropriate names in namelist PARAM_CNV2D. Check!'
!!$       call PRC_abort
!!$    endif
!!$    LOG_NML(PARAM_CNV2D)
 
 
    return

References scale_prc::prc_abort().

Here is the call graph for this function:

cnv2d_tile_init()

subroutine, public mod_cnv2d::cnv2d_tile_init	(	character(len=*), intent(in)	dtype,
		real(rp), intent(in)	dlat,
		real(rp), intent(in)	dlon,
		character(len=*), intent(in)	dir,
		character(len=*), intent(in)	catalogue,
		character(len=*), intent(in)	interp_type,
		integer, intent(in), optional	interp_level,
		integer, intent(in), optional	nmax
	)

Definition at line 124 of file mod_cnv2d.F90.

    use scale_prc, only: &
       prc_abort
    use scale_const, only: &
       undef  => const_undef, &
       d2r    => const_d2r
    use scale_file_tiledata, only: &
       file_tiledata_get_info,   &
       file_tiledata_get_latlon, &
       file_tiledata_get_data
    use scale_atmos_grid_cartesc_real, only: &
       latxv => atmos_grid_cartesc_real_latxv, &
       lonuy => atmos_grid_cartesc_real_lonuy
    implicit none
    character(len=*), intent(in) :: dtype
    real(RP),         intent(in) :: dlat, dlon
    character(len=*), intent(in) :: dir
    character(len=*), intent(in) :: catalogue
    character(len=*), intent(in) :: interp_type
 
    integer, intent(in), optional :: interp_level
    integer, intent(in), optional :: nmax
 
    real(RP) :: DOMAIN_LATS, DOMAIN_LATE
    real(RP) :: DOMAIN_LONS, DOMAIN_LONE
 
    character(len=H_LONG) :: fname
 
    integer :: ierr
    integer :: i, j
    !---------------------------------------------------------------------------
 
    if ( present(nmax) ) then
       tile_nlim = nmax
    else
       tile_nlim = 1000
    end if
 
    domain_lats = minval( latxv(:,:) )
    domain_late = maxval( latxv(:,:) )
    domain_lons = minval( lonuy(:,:) )
    domain_lone = maxval( lonuy(:,:) )
 
    log_info("CNV2D_setup",*) 'Domain Information'
    log_info_cont(*) 'Domain (LAT)    :', domain_lats/d2r, domain_late/d2r
    log_info_cont(*) '       (LON)    :', domain_lons/d2r, domain_lone/d2r
 
    tile_dlat = dlat * d2r
    tile_dlon = dlon * d2r
 
    if (allocated(tile_fname)) deallocate(tile_fname)
    if (allocated(tile_hit)) deallocate(tile_hit)
    if (allocated(tile_js)) deallocate(tile_js)
    if (allocated(tile_je)) deallocate(tile_je)
    if (allocated(tile_is)) deallocate(tile_is)
    if (allocated(tile_ie)) deallocate(tile_ie)
    if (allocated(lat_1d)) deallocate(lat_1d)
    if (allocated(lon_1d)) deallocate(lon_1d)
 
    allocate( tile_fname(tile_nlim), tile_hit(tile_nlim) )
    allocate( tile_js(tile_nlim), tile_je(tile_nlim), tile_is(tile_nlim), tile_ie(tile_nlim) )
 
    ! catalogue file
    fname = trim(dir)//'/'//trim(catalogue)
 
    call file_tiledata_get_info( tile_nlim,                                          & ! [IN]
                                 tile_dlat, tile_dlon,                               & ! [IN]
                                 domain_lats, domain_late, domain_lons, domain_lone, & ! [IN]
                                 fname,                                              & ! [IN]
                                 global_ia,                                          & ! [OUT]
                                 tile_nmax,                                          & ! [OUT]
                                 tile_fname(:), tile_hit(:),                         & ! [OUT]
                                 tile_js(:), tile_je(:), tile_is(:), tile_ie(:),     & ! [OUT]
                                 nlat, nlon, dom_js, dom_je, dom_is, dom_ie,         & ! [OUT]
                                 zonal, pole                                         ) ! [OUT]
 
    allocate( lat_1d(nlat) )
    allocate( lon_1d(nlon) )
 
    call file_tiledata_get_latlon( nlat, nlon,           & ! [IN]
                                   dom_js, dom_is,       & ! [IN]
                                   tile_dlat, tile_dlon, & ! [IN]
                                   lat_1d(:), lon_1d(:)  ) ! [OUT]
 
    cnv2d_ftype = i_tile
    cnv2d_tile_dtype = dtype
    cnv2d_tile_dir   = dir
 
    call cnv2d_init( interp_type,                 &
                     interp_level = interp_level, &
                     ll_struct = .true.           )
 
    return

References scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_latxv, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_lonuy, scale_const::const_d2r, scale_const::const_undef, scale_file_tiledata::file_tiledata_get_info(), scale_file_tiledata::file_tiledata_get_latlon(), and scale_prc::prc_abort().

Referenced by mod_cnvtopo::cnvtopo(), and mod_cnvuser::cnvuser().

Here is the call graph for this function:

Here is the caller graph for this function:

cnv2d_grads_init()

subroutine, public mod_cnv2d::cnv2d_grads_init	(	character(len=*), intent(in)	FILE_NAME,
		character(len=*), intent(in)	VAR_NAME,
		character(len=*), intent(in)	LAT_NAME,
		character(len=*), intent(in)	LON_NAME,
		character(len=*), intent(in)	interp_type,
		integer, intent(in)	interp_level,
		real(rp), intent(in), optional	search_limit,
		character(len=*), intent(in), optional	POSTFIX
	)

Definition at line 228 of file mod_cnv2d.F90.

    use scale_file_grads, only: &
       file_grads_open,      &
       file_grads_get_shape, &
       file_grads_varid,     &
       file_grads_isoned,    &
       file_grads_read
    use scale_const, only: &
       d2r => const_d2r
    implicit none
    character(len=*), intent(in) :: FILE_NAME
    character(len=*), intent(in) :: VAR_NAME
    character(len=*), intent(in) :: LAT_NAME
    character(len=*), intent(in) :: LON_NAME
    character(len=*), intent(in) :: INTERP_TYPE
    integer,          intent(in) :: interp_level
    real(RP),         intent(in), optional :: search_limit
    character(len=*), intent(in), optional :: POSTFIX
 
    integer :: file_id, var_id
    integer :: shape(2)
 
    integer :: i, j
 
    call file_grads_open( file_name, & ! [IN]
                          file_id    ) ! [OUT]
 
    call file_grads_get_shape( file_id, var_name, &
                               shape(:)           )
    nlon = shape(1)
    nlat = shape(2)
 
    if (allocated(lat_org)) deallocate(lat_org)
    if (allocated(lon_org)) deallocate(lon_org)
    if (allocated(lat_1d)) deallocate(lat_1d)
    if (allocated(lon_1d)) deallocate(lon_1d)
    allocate( lat_org(nlon,nlat), lon_org(nlon,nlat) )
    allocate( lat_1d(nlat), lon_1d(nlon) )
 
    ! lat
    call file_grads_varid( file_id, lat_name, & ! (in)
                           var_id             ) ! (out)
    if ( file_grads_isoned( file_id, var_id ) ) then
       call file_grads_read( file_id, var_id, & ! (in)
                             lat_1d(:)        ) ! (out)
       !$omp parallel do collapse(2)
       do j = 1, nlat
       do i = 1, nlon
          lat_org(i,j) = lat_1d(j) * d2r
       end do
       end do
    else
       call file_grads_read( file_id, var_id,  & ! (in)
                             lat_org(:,:),     & ! (out)
                             postfix = postfix ) ! (in)
       !$omp parallel do collapse(2)
       do j = 1, nlat
       do i = 1, nlon
          lat_org(i,j) = lat_org(i,j) * d2r
       end do
       end do
    end if
 
    ! lon
    call file_grads_varid( file_id, lon_name, & ! (in)
                           var_id             ) ! (out)
    if ( file_grads_isoned( file_id, var_id ) ) then
       call file_grads_read( file_id, var_id, & ! (in)
                             lon_1d(:)         ) ! (out)
       !$omp parallel do collapse(2)
       do j = 1, nlat
       do i = 1, nlon
          lon_org(i,j) = lon_1d(i) * d2r
       end do
       end do
    else
       call file_grads_read( file_id, var_id,  & ! (in)
                             lon_org(:,:),     & ! (out)
                             postfix = postfix ) ! (in)
       !$omp parallel do collapse(2)
       do j = 1, nlat
       do i = 1, nlon
          lon_org(i,j) = lon_org(i,j) * d2r
       end do
       end do
    end if
 
    call file_grads_varid( file_id, var_name, & ! [IN]
                           var_id             ) ! [OUT]
 
    cnv2d_ftype = i_grads
    cnv2d_grads_fid = file_id
    cnv2d_grads_vid = var_id
 
    call cnv2d_init( interp_type,                 &
                     interp_level = interp_level, &
                     search_limit = search_limit, &
                     ll_struct    = .false.       )
 
    return

References scale_const::const_d2r, scale_file_grads::file_grads_isoned(), scale_file_grads::file_grads_open(), and scale_file_grads::file_grads_varid().

Referenced by mod_cnvtopo::cnvtopo(), and mod_cnvuser::cnvuser().

Here is the call graph for this function:

Here is the caller graph for this function:

cnv2d_exec()

subroutine, public mod_cnv2d::cnv2d_exec	(	real(rp), dimension(ia,ja), intent(out)	var,
		integer, intent(in), optional	step,
		real(rp), intent(in), optional	min_value,
		logical, intent(in), optional	yrevers
	)

Definition at line 334 of file mod_cnv2d.F90.

    use scale_file_tiledata, only: &
       file_tiledata_get_data
    use scale_file_grads, only: &
       file_grads_read
    use scale_interp, only: &
       interp_interp2d
    implicit none
    real(RP), intent(out) :: var(IA,JA)
    integer,  intent(in), optional :: step
    real(RP), intent(in), optional :: min_value
    logical,  intent(in), optional :: yrevers
 
    select case ( cnv2d_ftype )
    case ( i_tile )
       call file_tiledata_get_data( nlat, nlon,                                     & ! [IN]
                                    cnv2d_tile_dir,                                 & ! [IN]
                                    global_ia,                                      & ! [IN]
                                    tile_nmax,                                      & ! [IN]
                                    tile_dlat, tile_dlon,                           & ! [IN]
                                    tile_fname(:), tile_hit(:),                     & ! [IN]
                                    tile_js(:), tile_je(:), tile_is(:), tile_ie(:), & ! [IN]
                                    dom_js, dom_je, dom_is, dom_ie,                 & ! [IN]
                                    cnv2d_tile_dtype,                               & ! [IN]
                                    data_org(:,:),                                  & ! [OUT]
                                    step = step,                                    & ! [IN]
                                    min_value = min_value, yrevers = yrevers        ) ! [IN]
    case ( i_grads )
       call file_grads_read( cnv2d_grads_fid, cnv2d_grads_vid, & ! [IN]
                             data_org(:,:),                    & ! [OUT]
                             step = step                       ) ! [IN]
    end select
 
 
    call interp_interp2d( itp_lev,                    & ! [IN]
                          nlon,nlat,                  & ! [IN]
                          ia, ja,                     & ! [IN]
                          idx_i(:,:,:), idx_j(:,:,:), & ! [IN]
                          hfact(:,:,:),               & ! [IN]
                          data_org(:,:),              & ! [IN]
                          var(:,:)                    ) ! [OUT]
 
    return

References scale_atmos_grid_cartesc::atmos_grid_cartesc_cx, scale_atmos_grid_cartesc::atmos_grid_cartesc_cy, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_lat, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_lon, scale_const::const_pi, scale_atmos_grid_cartesc_index::ia, scale_interp::interp_factor2d_linear_latlon(), scale_interp::interp_factor2d_linear_xy(), scale_interp::interp_factor2d_weight(), scale_interp::interp_interp2d(), scale_atmos_grid_cartesc_index::ja, and scale_prc::prc_abort().

Referenced by mod_cnvtopo::cnvtopo(), and mod_cnvuser::cnvuser().

Here is the call graph for this function:

Here is the caller graph for this function: