Changeset 943 for XIOS/trunk/src/interface/fortran/ifield.F90

Timestamp:

09/22/16 10:59:27 (8 years ago)

Author:

mhnguyen

Message:

Adding Fortran interfaces for retrieving domains, axis, scalars of a field

+) Add some functions in Grid to get its domain, axis, scalar

Test
+) On Curie
+) Work

File:

• XIOS/trunk/src/interface/fortran/ifield.F90 (modified) (3 diffs)

XIOS/trunk/src/interface/fortran/ifield.F90

@@ -8 +8 @@
 !   USE IFIELDGROUP_ATTR
    USE IDURATION
+   USE IDOMAIN
+   USE IAXIS
+   USE ISCALAR
    
    TYPE txios(field)
@@ -16 +19 @@
       INTEGER(kind = C_INTPTR_T) :: daddr
    END TYPE txios(fieldgroup)
-   
+
    CONTAINS ! Fonctions disponibles pour les utilisateurs.
 
@@ -35 +38 @@
    END SUBROUTINE xios(get_fieldgroup_handle)
    
+   SUBROUTINE xios(field_id_get_domain_handle)(field_id, ret, idx)
+      IMPLICIT NONE
+      CHARACTER(len  = *)     , INTENT(IN)  :: field_id
+      INTEGER, OPTIONAL       , INTENT(IN)  :: idx
+      TYPE(txios(domain)), INTENT(OUT) :: ret
+      TYPE(txios(field))      :: field_hdl
+      INTEGER                 :: index
+      index = 0
+      IF (PRESENT(idx)) THEN
+        index = idx
+      ENDIF
+
+      CALL xios(get_field_handle)(field_id,field_hdl)
+      CALL xios(field_get_domain_handle(field_hdl, ret, index))
+   END SUBROUTINE xios(field_id_get_domain_handle)
+
+   SUBROUTINE xios(field_get_domain_handle)(field_hdl, ret, idx)
+      IMPLICIT NONE
+      TYPE(txios(field))      , INTENT(IN)  :: field_hdl
+      INTEGER, OPTIONAL       , INTENT(IN)  :: idx
+      TYPE(txios(domain)), INTENT(OUT) :: ret
+      INTEGER :: index
+      index = 0
+      IF (PRESENT(idx)) THEN
+        index = idx
+      ENDIF
+      CALL cxios_field_get_domain_handle(ret%daddr, field_hdl%daddr, index)
+   END SUBROUTINE xios(field_get_domain_handle)
+
+   SUBROUTINE xios(field_id_get_axis_handle)(field_id, ret, idx)
+      IMPLICIT NONE
+      CHARACTER(len  = *)   , INTENT(IN)  :: field_id
+      INTEGER, OPTIONAL     , INTENT(IN)  :: idx
+      TYPE(txios(axis)), INTENT(OUT) :: ret
+      TYPE(txios(field))     :: field_hdl
+      INTEGER                :: index
+      index = 0
+      IF (PRESENT(idx)) THEN
+        index = idx
+      ENDIF
+      CALL xios(get_field_handle)(field_id,field_hdl)
+      CALL xios(field_get_axis_handle(field_hdl, ret, index))
+   END SUBROUTINE xios(field_id_get_axis_handle)
+
+   SUBROUTINE xios(field_get_axis_handle)(field_hdl, ret, idx)
+      IMPLICIT NONE
+      TYPE(txios(field))    , INTENT(IN)  :: field_hdl
+      INTEGER, OPTIONAL     , INTENT(IN)  :: idx
+      TYPE(txios(axis)), INTENT(OUT) :: ret
+      INTEGER :: index
+      index = 0
+      IF (PRESENT(idx)) THEN
+        index = idx
+      ENDIF
+      CALL cxios_field_get_axis_handle(ret%daddr, field_hdl%daddr, index)
+   END SUBROUTINE xios(field_get_axis_handle)
+
+   SUBROUTINE xios(field_id_get_scalar_handle)(field_id, ret, idx)
+      IMPLICIT NONE
+      CHARACTER(len  = *)     , INTENT(IN)  :: field_id
+      INTEGER, OPTIONAL       , INTENT(IN)  :: idx
+      TYPE(txios(scalar)), INTENT(OUT) :: ret
+      TYPE(txios(field))     :: field_hdl
+      INTEGER                :: index
+      index = 0
+      IF (PRESENT(idx)) THEN
+        index = idx
+      ENDIF
+      CALL xios(get_field_handle)(field_id,field_hdl)
+      CALL xios(field_get_scalar_handle(field_hdl, ret, index))
+   END SUBROUTINE xios(field_id_get_scalar_handle)
+
+   SUBROUTINE xios(field_get_scalar_handle)(field_hdl, ret, idx)
+      IMPLICIT NONE
+      TYPE(txios(field))      , INTENT(IN)  :: field_hdl
+      INTEGER, OPTIONAL       , INTENT(IN)  :: idx
+      TYPE(txios(scalar)), INTENT(OUT) :: ret
+      INTEGER :: index
+      index = 0
+      IF (PRESENT(idx)) THEN
+        index = idx
+      ENDIF
+      CALL cxios_field_get_axis_handle(ret%daddr, field_hdl%daddr, index)
+   END SUBROUTINE xios(field_get_scalar_handle)
 
    LOGICAL FUNCTION xios(is_valid_field)(idt)