libsim Versione 7.2.4
|
◆ geo_proj_get_val()
Definizione alla linea 597 del file geo_proj_class.F90. 603
604IF (c_e(this%xoff) .AND. this%xoff /= 0.0d0) THEN
605 print*,"False easting",this%xoff
606ENDIF
607IF (c_e(this%yoff) .AND. this%yoff /= 0.0d0) THEN
608 print*,"False northing",this%yoff
609ENDIF
610
611IF (c_e(this%lov)) THEN
612 print*,"Central Meridian",this%lov
613ENDIF
614
615IF (this%proj_type == 'rotated_ll' .OR. this%proj_type == 'stretched_rotated_ll') THEN
616 print*,"Rotated projection:"
617 print*,"lon of south pole",this%rotated%longitude_south_pole
618 print*,"lat of south pole",this%rotated%latitude_south_pole
619 print*,"angle of rotation",this%rotated%angle_rotation
620ENDIF
621
622IF (this%proj_type == 'stretched_ll' .OR. this%proj_type == 'stretched_rotated_ll') THEN
623 print*,"Stretched projection:"
624 print*,"lon of stretch pole",this%stretched%longitude_stretch_pole
625 print*,"lat of stretch pole",this%stretched%latitude_stretch_pole
626 print*,"stretching factor",this%stretched%stretch_factor
627ENDIF
628
629IF (this%proj_type == 'lambert' .OR. this%proj_type == 'polar_stereographic') THEN
630 print*,"Polar projection:"
631 IF (c_e(this%polar%latin1) .OR. c_e(this%polar%latin2)) THEN
632 print*,"lat of intersections",this%polar%latin1,this%polar%latin2
633 ENDIF
634 IF (c_e(this%polar%lad)) THEN
635 print*,"isometric latitude",this%polar%lad
636 ENDIF
637 IF (iand(this%polar%projection_center_flag, 128) == 0) THEN
638 print*,"North Pole"
639 ELSE
640 print*,"South Pole"
641 ENDIF
642ENDIF
643
644IF (this%proj_type == 'mercator') THEN
645 IF (c_e(this%polar%lad)) THEN
646 print*,"isometric latitude",this%polar%lad
647 ENDIF
648ENDIF
649
650IF (this%ellips%f == 0.0d0) THEN
651 print*,"Spherical Earth:"
652 print*,"Radius (m)",this%ellips%a
653ELSE
654 print*,"Ellipsoid:"
655 print*,"Flattening",this%ellips%f
656 print*,"Reverse of flattening",1.0d0/this%ellips%f
657 print*,"Semi-major axis (m)",this%ellips%a
658ENDIF
659
660
661END SUBROUTINE geo_proj_display
662
663
666ELEMENTAL SUBROUTINE geo_proj_proj(this, lon, lat, x, y)
667TYPE(geo_proj),INTENT(in) :: this
669DOUBLE PRECISION, INTENT(in) :: lon, lat
671DOUBLE PRECISION, INTENT(out) :: x, y
672
673SELECT CASE(this%proj_type)
674
|