module Map projection More...

Data Types
interface	lonlat2xy_s

interface	mapfactor_s

interface	rotcoef_s

Functions/Subroutines
subroutine, public	mapprojection_setup (DOMAIN_CENTER_X, DOMAIN_CENTER_Y)
	Setup. More...

subroutine	mapprojection_xy2lonlat_0d (x, y, lon, lat)
	(x,y) -> (lon,lat) More...

subroutine, public	mapprojection_mapfactor (IA, IS, IE, JA, JS, JE, lat, m1, m2)
	(x,y) -> (lon,lat) More...

subroutine, public	mapprojection_rotcoef (IA, IS, IE, JA, JS, JE, lon, lat, rotc)
	u(lat,lon) = cos u(x,y) - sin v(x,y) v(lat,lon) = sin u(x,y) + cos v(x,y) More...

subroutine, public	mapprojection_get_attributes (mapping, false_easting, false_northing, longitude_of_central_meridian, longitude_of_projection_origin, latitude_of_projection_origin, straight_vertical_longitude_from_pole, standard_parallel)
	Get mapping attributes. More...

Variables
real(rp), public	mapprojection_basepoint_lon = 135.221_RP

real(rp), public	mapprojection_basepoint_lat = 34.653_RP

procedure(xy2lonlat_s), pointer	xy2lonlat => NULL()

Detailed Description

module Map projection

Description: Map projection module

Author: Team SCALE

NAMELIST

PARAM_MAPPROJECTION

name	type	default value	comment
MAPPROJECTION_BASEPOINT_LON	real(RP)	135.221_RP	position of base point (domain center) in real world [deg]
MAPPROJECTION_BASEPOINT_LAT	real(RP)	34.653_RP	position of base point (domain center) in real world [deg]
MAPPROJECTION_BASEPOINT_X	real(DP)		position of base point in the model [m]
MAPPROJECTION_BASEPOINT_Y	real(DP)		position of base point in the model [m]
MAPPROJECTION_TYPE	character(len=H_SHORT)	'NONE'	map projection type
MAPPROJECTION_ROTATION	real(DP)	0.0_DP	rotation factor
MAPPROJECTION_LC_LAT1	real(DP)	30.0_DP	standard latitude1 for L.C. projection [deg]
MAPPROJECTION_LC_LAT2	real(DP)	60.0_DP	standard latitude2 for L.C. projection [deg]
MAPPROJECTION_PS_LAT	real(DP)		standard latitude1 for P.S. projection [deg]
MAPPROJECTION_M_LAT	real(DP)	0.0_DP	standard latitude1 for Mer. projection [deg]
MAPPROJECTION_EC_LAT	real(DP)	0.0_DP	standard latitude1 for E.C. projection [deg]

History Output: No history output

Function/Subroutine Documentation

◆ mapprojection_setup()

subroutine, public scale_mapprojection::mapprojection_setup	(	real(rp), intent(in)	DOMAIN_CENTER_X,
		real(rp), intent(in)	DOMAIN_CENTER_Y
	)

Setup.

Parameters

[in]	domain_center_x	center position of global domain [m]: x
[in]	domain_center_y	center position of global domain [m]: y

Definition at line 167 of file scale_mapprojection.F90.

References scale_const::const_d2r, scale_const::const_pi, scale_const::const_radius, scale_const::const_undef, scale_io::io_fid_conf, mapprojection_basepoint_lat, mapprojection_basepoint_lon, scale_prc::prc_abort(), and xy2lonlat.

Referenced by scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_setup().

     use scale_prc, only: &
        prc_abort
     use scale_const, only: &
        undef     => const_undef,  &
        pi_rp     => const_pi,     &
        d2r_rp    => const_d2r,    &
        radius_rp => const_radius
     implicit none
 
     real(RP), intent(in) :: domain_center_x
     real(RP), intent(in) :: domain_center_y
 
     namelist / param_mapprojection / &
        mapprojection_basepoint_lon, &
        mapprojection_basepoint_lat, &
        mapprojection_basepoint_x,   &
        mapprojection_basepoint_y,   &
        mapprojection_type,          &
        mapprojection_rotation,      &
        mapprojection_lc_lat1,       &
        mapprojection_lc_lat2,       &
        mapprojection_ps_lat,        &
        mapprojection_m_lat,         &
        mapprojection_ec_lat
 
     integer :: ierr
     !---------------------------------------------------------------------------
 
     log_newline
     log_info("MAPPROJECTION_setup",*) 'Setup'
 
     pi     = real(pi_rp,     kind=dp)
     d2r    = real(d2r_rp,    kind=dp)
     radius = real(radius_rp, kind=dp)
 
     mapprojection_basepoint_x = undef
     mapprojection_basepoint_y = undef
     mapprojection_ps_lat      = undef
 
     mapprojection_basepoint_x = domain_center_x
     mapprojection_basepoint_y = domain_center_y
 
     !--- read namelist
     rewind(io_fid_conf)
     read(io_fid_conf,nml=param_mapprojection,iostat=ierr)
 
     if( ierr < 0 ) then !--- missing
        log_info("MAPPROJECTION_setup",*) 'Not found namelist. Default used.'
     elseif( ierr > 0 ) then !--- fatal error
        log_error("MAPPROJECTION_setup",*) 'Not appropriate names in namelist PARAM_MAPPROJECTION. Check!'
        call prc_abort
     endif
     log_nml(param_mapprojection)
 
     log_newline
     log_info("MAPPROJECTION_setup",*) 'Map projection information'
     log_info_cont('(1x,A,F15.3)') 'Basepoint(x)    [m] : ', mapprojection_basepoint_x
     log_info_cont('(1x,A,F15.3)') 'Basepoint(y)    [m] : ', mapprojection_basepoint_y
     log_info_cont(*)              'Map projection type : ', trim(mapprojection_type)
 
     mapprojection_mapping = ""
     mapprojection_false_easting = undef
     mapprojection_false_northing = undef
     mapprojection_longitude_of_central_meridian = undef
     mapprojection_longitude_of_projection_origin = undef
     mapprojection_latitude_of_projection_origin = undef
     mapprojection_straight_vertical_longitude_from_pole = undef
     mapprojection_standard_parallel(:) = undef
 
     if( mapprojection_ps_lat == undef ) mapprojection_ps_lat = mapprojection_basepoint_lat
 
     select case(trim(mapprojection_type))
     case('NONE')
        log_info_cont(*) '=> NO map projection'
        call mapprojection_none_setup
        xy2lonlat => mapprojection_none_xy2lonlat
        lonlat2xy => mapprojection_none_lonlat2xy
        mapfactor => mapprojection_none_mapfactor
        rotcoef   => mapprojection_none_rotcoef
     case('LC')
        log_info_cont(*) '=> Lambert Conformal projection'
        call mapprojection_lambertconformal_setup
        xy2lonlat => mapprojection_lambertconformal_xy2lonlat
        lonlat2xy => mapprojection_lambertconformal_lonlat2xy
        mapfactor => mapprojection_lambertconformal_mapfactor
        rotcoef   => mapprojection_lambertconformal_rotcoef
     case('PS')
        log_info_cont(*) '=> Polar Stereographic projection'
        call mapprojection_polarstereographic_setup
        xy2lonlat => mapprojection_polarstereographic_xy2lonlat
        lonlat2xy => mapprojection_polarstereographic_lonlat2xy
        mapfactor => mapprojection_polarstereographic_mapfactor
        rotcoef   => mapprojection_polarstereographic_rotcoef
     case('MER')
        log_info_cont(*) '=> Mercator projection'
        call mapprojection_mercator_setup
        xy2lonlat => mapprojection_mercator_xy2lonlat
        lonlat2xy => mapprojection_mercator_lonlat2xy
        mapfactor => mapprojection_mercator_mapfactor
        rotcoef   => mapprojection_mercator_rotcoef
     case('EC')
        log_info_cont(*) '=> Equidistant Cylindrical projection'
        call mapprojection_equidistantcylindrical_setup
        xy2lonlat => mapprojection_equidistantcylindrical_xy2lonlat
        lonlat2xy => mapprojection_equidistantcylindrical_lonlat2xy
        mapfactor => mapprojection_equidistantcylindrical_mapfactor
        rotcoef   => mapprojection_equidistantcylindrical_rotcoef
     case default
        log_error("MAPPROJECTION_setup",*) 'Unsupported MAPPROJECTION_type. STOP'
        call prc_abort
     endselect
 
     mapprojection_rotation = mapprojection_rotation * d2r
     mapprojection_rot_fact_sin = sin(mapprojection_rotation)
     mapprojection_rot_fact_cos = cos(mapprojection_rotation)
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ mapprojection_xy2lonlat_0d()

subroutine scale_mapprojection::mapprojection_xy2lonlat_0d	(	real(rp), intent(in)	x,
		real(rp), intent(in)	y,
		real(rp), intent(out)	lon,
		real(rp), intent(out)	lat
	)

(x,y) -> (lon,lat)

Definition at line 291 of file scale_mapprojection.F90.

References xy2lonlat.

     implicit none
     real(RP), intent(in)  :: x
     real(RP), intent(in)  :: y
     real(RP), intent(out) :: lon ! [rad]
     real(RP), intent(out) :: lat ! [rad]
 
     real(DP) :: xx, yy
     !---------------------------------------------------------------------------
 
     xx = mapprojection_basepoint_x &
        + ( x - mapprojection_basepoint_x ) * mapprojection_rot_fact_cos &
        - ( y - mapprojection_basepoint_y ) * mapprojection_rot_fact_sin
     yy = mapprojection_basepoint_y &
        + ( y - mapprojection_basepoint_y ) * mapprojection_rot_fact_cos &
        + ( x - mapprojection_basepoint_x ) * mapprojection_rot_fact_sin
 
     call xy2lonlat( xx, yy, lon, lat )
 
     return

◆ mapprojection_mapfactor()

subroutine, public scale_mapprojection::mapprojection_mapfactor	(	integer, intent(in)	IA,
		integer, intent(in)	IS,
		integer, intent(in)	IE,
		integer, intent(in)	JA,
		integer, intent(in)	JS,
		integer, intent(in)	JE,
		real(rp), dimension(ia,ja), intent(in)	lat,
		real(rp), dimension (ia,ja), intent(out)	m1,
		real(rp), dimension (ia,ja), intent(out)	m2
	)

(x,y) -> (lon,lat)

Definition at line 397 of file scale_mapprojection.F90.

Referenced by scale_atmos_grid_cartesc_metric::atmos_grid_cartesc_metric_setup().

     implicit none
     integer, intent(in) :: ia, is, ie
     integer, intent(in) :: ja, js, je
 
     real(RP), intent(in)  :: lat(ia,ja) ! [rad]
     real(RP), intent(out) :: m1 (ia,ja)
     real(RP), intent(out) :: m2 (ia,ja)
 
     !real(RP) :: mm1, mm2
     integer :: i, j
     !---------------------------------------------------------------------------
 
     !$omp parallel do
     do j = js, je
     do i = is, ie
        call mapfactor( lat(i,j), m1(i,j), m2(i,j) )
 !!$       call mapfactor( lat(i,j), mm1, mm2 )
 !!$       m1(i,j) = sqrt( (MAPPROJECTION_rot_fact_cos * mm1)**2 + (MAPPROJECTION_rot_fact_sin * mm2)**2 )
 !!$       m2(i,j) = sqrt( (MAPPROJECTION_rot_fact_cos * mm2)**2 + (MAPPROJECTION_rot_fact_sin * mm1)**2 )
     end do
     end do
 
     return

Here is the caller graph for this function:

◆ mapprojection_rotcoef()

subroutine, public scale_mapprojection::mapprojection_rotcoef	(	integer, intent(in)	IA,
		integer, intent(in)	IS,
		integer, intent(in)	IE,
		integer, intent(in)	JA,
		integer, intent(in)	JS,
		integer, intent(in)	JE,
		real(rp), dimension(ia,ja), intent(in)	lon,
		real(rp), dimension(ia,ja), intent(in)	lat,
		real(rp), dimension(ia,ja,2), intent(out)	rotc
	)

u(lat,lon) = cos u(x,y) - sin v(x,y) v(lat,lon) = sin u(x,y) + cos v(x,y)

Parameters

[out] rotc rotc(:,:,1)->cos, rotc(:,:,2)->sin

Definition at line 429 of file scale_mapprojection.F90.

Referenced by scale_atmos_grid_cartesc_metric::atmos_grid_cartesc_metric_rotcoef().

     implicit none
     integer, intent(in) :: ia, is, ie
     integer, intent(in) :: ja, js, je
 
     real(RP), intent(in)  :: lon(ia,ja)    ! [rad]
     real(RP), intent(in)  :: lat(ia,ja)    ! [rad]
     real(RP), intent(out) :: rotc(ia,ja,2)
 
     integer :: i, j
     !---------------------------------------------------------------------------
 
     !$omp parallel do
     do j = js, je
     do i = is, ie
        call rotcoef( lon(i,j), lat(i,j), rotc(i,j,1), rotc(i,j,2) )
     end do
     end do
 
     return

Here is the caller graph for this function:

◆ mapprojection_get_attributes()

subroutine, public scale_mapprojection::mapprojection_get_attributes	(	character(len=*), intent(out)	mapping,
		real(dp), intent(out), optional	false_easting,
		real(dp), intent(out), optional	false_northing,
		real(dp), intent(out), optional	longitude_of_central_meridian,
		real(dp), intent(out), optional	longitude_of_projection_origin,
		real(dp), intent(out), optional	latitude_of_projection_origin,
		real(dp), intent(out), optional	straight_vertical_longitude_from_pole,
		real(dp), dimension(2), intent(out), optional	standard_parallel
	)

Get mapping attributes.

Definition at line 461 of file scale_mapprojection.F90.

References mapprojection_basepoint_lat, mapprojection_basepoint_lon, and scale_prc::prc_abort().

Referenced by scale_file_cartesc::file_cartesc_def_var(), scale_file_history_cartesc::file_history_cartesc_set_pres(), scale_file_history_cartesc::file_history_cartesc_truncate_3d(), and scale_file_cartesc::set_dimension_informations().

     implicit none
 
     character(len=*), intent(out) :: mapping
 
     real(DP), intent(out), optional :: false_easting
     real(DP), intent(out), optional :: false_northing
     real(DP), intent(out), optional :: longitude_of_central_meridian
     real(DP), intent(out), optional :: longitude_of_projection_origin
     real(DP), intent(out), optional :: latitude_of_projection_origin
     real(DP), intent(out), optional :: straight_vertical_longitude_from_pole
     real(DP), intent(out), optional :: standard_parallel(2)
     !---------------------------------------------------------------------------
 
     mapping = mapprojection_mapping
 
     if ( present(false_easting)                         ) &
                  false_easting                          = mapprojection_false_easting
     if ( present(false_northing)                        ) &
                  false_northing                         = mapprojection_false_northing
     if ( present(longitude_of_central_meridian)         ) &
                  longitude_of_central_meridian          = mapprojection_longitude_of_central_meridian
     if ( present(longitude_of_projection_origin)        ) &
                  longitude_of_projection_origin         = mapprojection_longitude_of_projection_origin
     if ( present(latitude_of_projection_origin)         ) &
                  latitude_of_projection_origin          = mapprojection_latitude_of_projection_origin
     if ( present(straight_vertical_longitude_from_pole) ) &
                  straight_vertical_longitude_from_pole  = mapprojection_straight_vertical_longitude_from_pole
     if ( present(standard_parallel)                     ) &
                  standard_parallel(:)                   = mapprojection_standard_parallel(:)
 
     return