psmile_info_coords_irreg2_dble.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_Info_coords_irreg2_dble
!
! !INTERFACE:

      subroutine psmile_info_coords_irreg2_dble (              &
                               x_coords, y_coords, z_coords,   &
                               coords_shape, grid_valid_shape, &
                               sinvec, cosvec, ierror)
!
! !USES:
!
      use PRISM_constants
!
      use PSMILe, dummy_interface => PSMILe_Info_coords_irreg2_dble

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

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

      Double Precision, Intent (In)   :: x_coords(coords_shape(1,1): 
                                                  coords_shape(2,1), 
                                                  coords_shape(1,2): 
                                                  coords_shape(2,2))
      Double Precision, Intent (In)   :: y_coords(coords_shape(1,1): 
                                                  coords_shape(2,1), 
                                                  coords_shape(1,2): 
                                                  coords_shape(2,2))
      Double Precision, Intent (In)   :: z_coords(coords_shape(1,3): 
                                                  coords_shape(2,3))

!     Coordinates of the method


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

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

!
! !INPUT/OUTPUT PARAMETERS:
!
      Type (dble_vector)              :: sinvec
!
!     Vector containing the sin values
!
      Type (dble_vector)              :: cosvec
!
!     Vector containing the cos values
!
! !OUTPUT PARAMETERS:
!
      Integer, Intent (Out)           :: ierror

!     Returns the error code of PSMILE_Info__coords_irreg2_dble;
!             ierror = 0 : No error
!             ierror > 0 : Severe error
!
! !DEFINED PARAMETERS:
!
!  lon   = Index of Longitudes in arrays "sin_values" and "cos_values"
!  lat   = Index of Latitudes  in arrays "sin_values" and "cos_values"
!
      Integer, Parameter              :: lon = 1
      Integer, Parameter              :: lat = 2
!
! !LOCAL VARIABLES
!
!     ... For sin and cos values
!
      Integer                         :: len
!
!     ... for error handling
!
      Integer, parameter              :: nerrp = 2
      Integer                         :: ierrp (nerrp)
!
!
! !DESCRIPTION:
!
! Subroutine "PSMILe_Info_coords_irreg2_dble" computes sin and cos values of
! the Longitudes and Latitudes of a 3d-dimensional grid of type
! "PRISM_Irrlonlat_Regvrt".
!
! !REVISION HISTORY:
!
!   Date      Programmer   Description
! ----------  ----------   -----------
! 03.07.21    H. Ritzdorf  created
!
!EOP
!----------------------------------------------------------------------
!
!  $Id: psmile_info_coords_irreg2_dble.F90 2325 2010-04-21 15:00:07Z valcke $
!  $Author: valcke $
!
   Character(len=len_cvs_string), save :: mycvs = 
       '$Id: psmile_info_coords_irreg2_dble.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
!
!-----------------------------------------------------------------------
!     Allocate and set temporary array which is used for the 
!     computation of distances.
!     (*) Longitudes and Latitudes are transformed
!         from degrees into radients.
!     (*) The z-values are currently not transformed
!         ??? Whats about PRISM_Irrlonlat_sigmavrt, ...
!
!     ??? Lohnt es sich den Shape weiter einzuschraenken ?
!-----------------------------------------------------------------------
!
      len = lat                                                 &
          * (grid_valid_shape(2,1) - grid_valid_shape(1,1) + 1) &
          * (grid_valid_shape(2,2) - grid_valid_shape(1,2) + 1)
!
      Allocate (sinvec%vector(len),                &
                cosvec%vector(len), STAT = ierror)
      if ( ierror > 0 ) then
         ierrp (1) = ierror
         ierrp (2) = len * 2

         ierror = PRISM_Error_Alloc
         call psmile_error ( ierror, 'sinvec%vector, cosvec%vector', &
                             ierrp, 2, __FILE__, __LINE__ )
         return
      endif
!
!===> Transform lon and lat values
!
      call psmile_trf_lonlat_2d_dble (x_coords, y_coords,             &
                                      coords_shape    (:, 1:ndim_2d), &
                                      grid_valid_shape(:, 1:ndim_2d), &
                                      sinvec%vector, cosvec%vector, ierror)
!
!===> All done
!
#ifdef VERBOSE
      print 9980, trim(ch_id), ierror

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

      end subroutine PSMILe_Info_coords_irreg2_dble