psmile_trf_lonlat_1d_real.F90

Go to the documentation of this file.

!-----------------------------------------------------------------------
! Copyright 2006-2010, NEC Europe Ltd., London, UK.
! All rights reserved. Use is subject to OASIS4 license terms.
!-----------------------------------------------------------------------
!BOP
!
! !ROUTINE: PSMILe_Trf_lonlat_1d_real
!
! !INTERFACE:

      subroutine psmile_trf_lonlat_1d_real (x_coords, y_coords,             &
                                            coords_shape, grid_valid_shape, &
                                            sin_values_lon, cos_values_lon, &
                                            sin_values_lat, cos_values_lat, &
                                            ierror)
!
! !USES:
!
      use PRISM_constants
!
      use PSMILe, dummy_interface => PSMILe_Trf_lonlat_1d_real

      implicit none
!
! !INPUT PARAMETERS:
!
      Integer, Intent (In)            :: coords_shape (2, ndim_2d)

!     Dimension of coordinates (method) x_coords, ...

      Real, Intent (In)               :: x_coords(coords_shape(1,1): 
                                                  coords_shape(2,1))
      Real, Intent (In)               :: y_coords(coords_shape(1,2): 
                                                  coords_shape(2,2))

!     Coordinates of the method

      Integer, Intent (In)            :: grid_valid_shape (2, ndim_2d)

!     Specifies the valid block shape for the "ndim_3d"-dimensional block

!
! !OUTPUT PARAMETERS:
!
      Real, Intent (Out)              :: sin_values_lon (grid_valid_shape(1,1): 
                                                         grid_valid_shape(2,1))
      Real, Intent (Out)              :: sin_values_lat (grid_valid_shape(1,2): 
                                                         grid_valid_shape(2,2))

!     Arrays returning the sin values for lon and lat

      Real, Intent (Out)              :: cos_values_lon (grid_valid_shape(1,1): 
                                                         grid_valid_shape(2,1))
      Real, Intent (Out)              :: cos_values_lat (grid_valid_shape(1,2): 
                                                         grid_valid_shape(2,2))

!     Arrays returning the cos values for lon and lat
!
      Integer, Intent (Out)           :: ierror

!     Returns the error code of PSMILE_Trf_lonlat_1d_real;
!             ierror = 0 : No error
!             ierror > 0 : Severe error
!
! !LOCAL VARIABLES
!
!     ... for error handling
!
!     Integer, parameter              :: nerrp = 2
!     Integer                         :: ierrp (nerrp)
!
!
! !DESCRIPTION:
!
! Subroutine "PSMILe_Trf_lonlat_1d_real" computes sin and cos values of
! the Longitudes and Latitudes of a 2d-dimensional grid of type
! "PRISM_Reglonlatvrt".
!
! !REVISION HISTORY:
!
!   Date      Programmer   Description
! ----------  ----------   -----------
! 03.07.21    H. Ritzdorf  created
!
!EOP
!----------------------------------------------------------------------
!
!  $Id: psmile_trf_lonlat_1d_real.F90 2325 2010-04-21 15:00:07Z valcke $
!  $Author: valcke $
!
   Character(len=len_cvs_string), save :: mycvs = 
       '$Id: psmile_trf_lonlat_1d_real.F90 2325 2010-04-21 15:00:07Z valcke $'
!----------------------------------------------------------------------
!
!  Initialization
!
#ifdef VERBOSE
      print 9990, trim(ch_id)

      call psmile_flushstd
#endif /* VERBOSE */

      ierror = 0
!
#ifdef PRISM_ASSERTION
#endif
!
!===> Transform from degrees into radients
!
      sin_values_lon (:) = &
        x_coords (grid_valid_shape(1,1):grid_valid_shape(2,1)) *      &
        real_deg2rad
!
      sin_values_lat (:) = &
        y_coords (grid_valid_shape(1,2):grid_valid_shape(2,2)) *      &
        real_deg2rad
!
!===> Compute cos and sin values
!
      cos_values_lon = cos (sin_values_lon)
      cos_values_lat = cos (sin_values_lat)
!
      sin_values_lon = sin (sin_values_lon)
      sin_values_lat = sin (sin_values_lat)
!
!===> All done
!
#ifdef VERBOSE
      print 9980, trim(ch_id), ierror

      call psmile_flushstd
#endif /* VERBOSE */
!
!  Formats:
!
9990 format (1x, a, ': psmile_trf_lonlat_1d_real')
9980 format (1x, a, ': psmile_trf_lonlat_1d_real: eof, ierror =', i3)

      end subroutine PSMILe_Trf_lonlat_1d_real