Ignore:
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

Location:
XIOS/trunk/src/interface
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • XIOS/trunk/src/interface/c/icfield.cpp

    r591 r943  
    2525   typedef xios::CFieldGroup * XFieldGroupPtr; 
    2626    
     27   typedef xios::CDomain     * XDomainPtr; 
     28   typedef xios::CAxis       * XAxisPtr; 
     29   typedef xios::CScalar     * XScalarPtr; 
     30 
    2731// --------------------------------------------------------------------------    
    2832// ------------------------ Création des handle ----------------------------- 
    … …  
    6973 
    7074// -----------------------------------------------------------------------------------------------------    
    71 // ------------------------- other function----------  --------------------------- --------------------- 
     75// ------------------------- Other functions 
    7276// -----------------------------------------------------------------------------------------------------    
    7377 
    … …  
    7882    CTimer::get("XIOS").suspend() ; 
    7983  } 
    80     
     84 
     85// ----------------------------------------------------------------------------------------------------- 
     86// ------------------------- Retrieving information of grid associated to field 
     87// ----------------------------------------------------------------------------------------------------- 
     88  void cxios_field_get_domain_handle(XDomainPtr * domain_hdl_ret, XFieldPtr field_hdl, int domainIndex) 
     89  { 
     90     CTimer::get("XIOS").resume() ; 
     91     *domain_hdl_ret = field_hdl->grid->getDomain(domainIndex); 
     92     CTimer::get("XIOS").suspend(); 
     93  } 
     94 
     95  void cxios_field_get_axis_handle(XAxisPtr * axis_hdl_ret, XFieldPtr field_hdl, int axisIndex) 
     96  { 
     97     CTimer::get("XIOS").resume() ; 
     98     *axis_hdl_ret = field_hdl->grid->getAxis(axisIndex); 
     99     CTimer::get("XIOS").suspend(); 
     100  } 
     101 
     102  void cxios_field_get_scalar_handle(XScalarPtr * scalar_hdl_ret, XFieldPtr field_hdl, int scalarIndex) 
     103  { 
     104     CTimer::get("XIOS").resume() ; 
     105     *scalar_hdl_ret = field_hdl->grid->getScalar(scalarIndex); 
     106     CTimer::get("XIOS").suspend(); 
     107  } 
    81108} // extern "C" 

  • XIOS/trunk/src/interface/fortran/field_interface.f90

    r501 r943  
    2424      END SUBROUTINE cxios_field_is_active 
    2525 
     26      SUBROUTINE cxios_field_get_domain_handle(ret, field_hdl, idx) BIND(C) 
     27         USE ISO_C_BINDING 
     28         INTEGER  (kind = C_INTPTR_T)               :: ret 
     29         INTEGER  (kind = C_INTPTR_T), VALUE        :: field_hdl 
     30         INTEGER  (kind = C_INT)     , VALUE        :: idx 
     31      END SUBROUTINE cxios_field_get_domain_handle 
     32 
     33      SUBROUTINE cxios_field_get_axis_handle(ret, field_hdl, idx) BIND(C) 
     34         USE ISO_C_BINDING 
     35         INTEGER  (kind = C_INTPTR_T)               :: ret 
     36         INTEGER  (kind = C_INTPTR_T), VALUE        :: field_hdl 
     37         INTEGER  (kind = C_INT)     , VALUE        :: idx 
     38      END SUBROUTINE cxios_field_get_axis_handle 
     39 
     40      SUBROUTINE cxios_field_get_scalar_handle(ret, field_hdl, idx) BIND(C) 
     41         USE ISO_C_BINDING 
     42         INTEGER  (kind = C_INTPTR_T)               :: ret 
     43         INTEGER  (kind = C_INTPTR_T), VALUE        :: field_hdl 
     44         INTEGER  (kind = C_INT)     , VALUE        :: idx 
     45      END SUBROUTINE cxios_field_get_scalar_handle 
     46 
    2647   END INTERFACE 
    2748    

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

    r545 r943  
    88!   USE IFIELDGROUP_ATTR 
    99   USE IDURATION 
     10   USE IDOMAIN 
     11   USE IAXIS 
     12   USE ISCALAR 
    1013    
    1114   TYPE txios(field) 
    … …  
    1619      INTEGER(kind = C_INTPTR_T) :: daddr 
    1720   END TYPE txios(fieldgroup) 
    18     
     21 
    1922   CONTAINS ! Fonctions disponibles pour les utilisateurs. 
    2023 
    … …  
    3538   END SUBROUTINE xios(get_fieldgroup_handle) 
    3639    
     40   SUBROUTINE xios(field_id_get_domain_handle)(field_id, ret, idx) 
     41      IMPLICIT NONE 
     42      CHARACTER(len  = *)     , INTENT(IN)  :: field_id 
     43      INTEGER, OPTIONAL       , INTENT(IN)  :: idx 
     44      TYPE(txios(domain)), INTENT(OUT) :: ret 
     45      TYPE(txios(field))      :: field_hdl 
     46      INTEGER                 :: index 
     47      index = 0 
     48      IF (PRESENT(idx)) THEN 
     49        index = idx 
     50      ENDIF 
     51 
     52      CALL xios(get_field_handle)(field_id,field_hdl) 
     53      CALL xios(field_get_domain_handle(field_hdl, ret, index)) 
     54   END SUBROUTINE xios(field_id_get_domain_handle) 
     55 
     56   SUBROUTINE xios(field_get_domain_handle)(field_hdl, ret, idx) 
     57      IMPLICIT NONE 
     58      TYPE(txios(field))      , INTENT(IN)  :: field_hdl 
     59      INTEGER, OPTIONAL       , INTENT(IN)  :: idx 
     60      TYPE(txios(domain)), INTENT(OUT) :: ret 
     61      INTEGER :: index 
     62      index = 0 
     63      IF (PRESENT(idx)) THEN 
     64        index = idx 
     65      ENDIF 
     66      CALL cxios_field_get_domain_handle(ret%daddr, field_hdl%daddr, index) 
     67   END SUBROUTINE xios(field_get_domain_handle) 
     68 
     69   SUBROUTINE xios(field_id_get_axis_handle)(field_id, ret, idx) 
     70      IMPLICIT NONE 
     71      CHARACTER(len  = *)   , INTENT(IN)  :: field_id 
     72      INTEGER, OPTIONAL     , INTENT(IN)  :: idx 
     73      TYPE(txios(axis)), INTENT(OUT) :: ret 
     74      TYPE(txios(field))     :: field_hdl 
     75      INTEGER                :: index 
     76      index = 0 
     77      IF (PRESENT(idx)) THEN 
     78        index = idx 
     79      ENDIF 
     80      CALL xios(get_field_handle)(field_id,field_hdl) 
     81      CALL xios(field_get_axis_handle(field_hdl, ret, index)) 
     82   END SUBROUTINE xios(field_id_get_axis_handle) 
     83 
     84   SUBROUTINE xios(field_get_axis_handle)(field_hdl, ret, idx) 
     85      IMPLICIT NONE 
     86      TYPE(txios(field))    , INTENT(IN)  :: field_hdl 
     87      INTEGER, OPTIONAL     , INTENT(IN)  :: idx 
     88      TYPE(txios(axis)), INTENT(OUT) :: ret 
     89      INTEGER :: index 
     90      index = 0 
     91      IF (PRESENT(idx)) THEN 
     92        index = idx 
     93      ENDIF 
     94      CALL cxios_field_get_axis_handle(ret%daddr, field_hdl%daddr, index) 
     95   END SUBROUTINE xios(field_get_axis_handle) 
     96 
     97   SUBROUTINE xios(field_id_get_scalar_handle)(field_id, ret, idx) 
     98      IMPLICIT NONE 
     99      CHARACTER(len  = *)     , INTENT(IN)  :: field_id 
     100      INTEGER, OPTIONAL       , INTENT(IN)  :: idx 
     101      TYPE(txios(scalar)), INTENT(OUT) :: ret 
     102      TYPE(txios(field))     :: field_hdl 
     103      INTEGER                :: index 
     104      index = 0 
     105      IF (PRESENT(idx)) THEN 
     106        index = idx 
     107      ENDIF 
     108      CALL xios(get_field_handle)(field_id,field_hdl) 
     109      CALL xios(field_get_scalar_handle(field_hdl, ret, index)) 
     110   END SUBROUTINE xios(field_id_get_scalar_handle) 
     111 
     112   SUBROUTINE xios(field_get_scalar_handle)(field_hdl, ret, idx) 
     113      IMPLICIT NONE 
     114      TYPE(txios(field))      , INTENT(IN)  :: field_hdl 
     115      INTEGER, OPTIONAL       , INTENT(IN)  :: idx 
     116      TYPE(txios(scalar)), INTENT(OUT) :: ret 
     117      INTEGER :: index 
     118      index = 0 
     119      IF (PRESENT(idx)) THEN 
     120        index = idx 
     121      ENDIF 
     122      CALL cxios_field_get_axis_handle(ret%daddr, field_hdl%daddr, index) 
     123   END SUBROUTINE xios(field_get_scalar_handle) 
    37124 
    38125   LOGICAL FUNCTION xios(is_valid_field)(idt) 

  • XIOS/trunk/src/interface/fortran/ixios.F90

    r934 r943  
    113113 
    114114USE XIOS_INTERFACES, ONLY : xios(set_attr), xios(get_attr), xios(is_defined_attr), xios(get_handle), xios(add_child), & 
    115                             xios(send_field), xios(recv_field), xios(field_is_active), xios(getVar), xios(setVar) 
     115                            xios(send_field), xios(recv_field), xios(field_is_active), xios(getVar), xios(setVar), & 
     116                            xios(field_get_domain),xios(field_get_axis),xios(field_get_scalar) 
    116117 
    117118END MODULE XIOS 

  • XIOS/trunk/src/interface/fortran/ixios_interfaces.F90

    r934 r943  
    2424USE idomaingroup_attr, ONLY : xios(set_domaingroup_attr_hdl), xios(get_domaingroup_attr_hdl), xios(is_defined_domaingroup_attr_hdl) 
    2525 
    26 USE ifield, ONLY : xios(get_field_handle), xios(get_fieldgroup_handle), xios(field_is_active_id), xios(field_is_active_hdl) 
     26USE ifield, ONLY : xios(get_field_handle), xios(get_fieldgroup_handle), xios(field_is_active_id), xios(field_is_active_hdl), & 
     27                   xios(field_get_domain_handle), xios(field_get_axis_handle), xios(field_get_scalar_handle), & 
     28                   xios(field_id_get_domain_handle), xios(field_id_get_axis_handle), xios(field_id_get_scalar_handle) 
    2729 
    2830USE ifield_attr, ONLY : xios(set_field_attr_hdl), xios(get_field_attr_hdl), xios(is_defined_field_attr_hdl) 
    … …  
    195197END INTERFACE xios(field_is_active) 
    196198 
     199INTERFACE xios(field_get_domain) 
     200  MODULE PROCEDURE xios(field_get_domain_handle), xios(field_id_get_domain_handle) 
     201END INTERFACE xios(field_get_domain) 
     202 
     203INTERFACE xios(field_get_axis) 
     204  MODULE PROCEDURE xios(field_get_axis_handle),xios(field_id_get_axis_handle) 
     205END INTERFACE xios(field_get_axis) 
     206 
     207INTERFACE xios(field_get_scalar) 
     208  MODULE PROCEDURE xios(field_get_scalar_handle),xios(field_id_get_scalar_handle) 
     209END INTERFACE xios(field_get_scalar) 
     210 
    197211INTERFACE xios(getVar) 
    198212  MODULE PROCEDURE xios(getVar_k8), xios(getVar_k4), xios(getVar_int), xios(getVar_logic), xios(getVar_char) 
    … …  
    204218 
    205219PUBLIC :: xios(set_attr), xios(get_attr), xios(is_defined_attr), xios(get_handle), xios(add_child), & 
    206           xios(send_field), xios(recv_field), xios(field_is_active), xios(getVar), xios(setVar) 
     220          xios(send_field), xios(recv_field), xios(field_is_active), xios(getVar), xios(setVar), & 
     221          xios(field_get_domain),xios(field_get_axis),xios(field_get_scalar) 
    207222 
    208223END MODULE XIOS_INTERFACES 
Note: See TracChangeset for help on using the changeset viewer.