Ignore:
Timestamp:
08/24/15 14:53:17 (9 years ago)
Author:
mhnguyen
Message:

Changing some domain attributes

+) Longitude and latitude from now on can be 2 dimension array

Test
+) On Curie
+) test_client, test_complete pass and are correct

File:
1 edited

Legend:

Unmodified
Added
Removed
  • XIOS/trunk/src/test/test_complete.f90

    r655 r664  
    2525  DOUBLE PRECISION,DIMENSION(ni_glo,nj_glo) :: lon_glo,lat_glo 
    2626  DOUBLE PRECISION :: field_A_glo(ni_glo,nj_glo,llm) 
    27   DOUBLE PRECISION,ALLOCATABLE :: lon(:,:),lat(:,:),field_A_atm(:,:,:), field_A_srf(:,:), lonvalue(:) 
     27  DOUBLE PRECISION,ALLOCATABLE :: lon(:,:),lat(:,:),field_A_atm(:,:,:), field_A_srf(:,:), lonvalue(:,:) 
    2828  INTEGER, ALLOCATABLE :: kindex(:) 
    2929  INTEGER :: ni,ibegin,iend,nj,jbegin,jend 
     
    7171  iend=ibegin+ni-1 ; jend=jbegin+nj-1 
    7272 
    73   ALLOCATE(lon(ni,nj),lat(ni,nj),field_A_atm(0:ni+1,-1:nj+2,llm),lonvalue(ni*nj)) 
     73  ALLOCATE(lon(ni,nj),lat(ni,nj),field_A_atm(0:ni+1,-1:nj+2,llm),lonvalue(ni,nj)) 
    7474  lon(:,:)=lon_glo(ibegin+1:iend+1,jbegin+1:jend+1) 
    7575  lat(:,:)=lat_glo(ibegin+1:iend+1,jbegin+1:jend+1) 
     
    8989  CALL xios_set_axis_attr("axis_atm",size=llm ,value=lval) ; 
    9090 
    91   CALL xios_set_domain_attr("domain_atm",ni_glo=ni_glo, nj_glo=nj_glo, ibegin=ibegin, ni=ni,jbegin=jbegin,nj=nj) 
     91  CALL xios_set_domain_attr("domain_atm",ni_glo=ni_glo, nj_glo=nj_glo, ibegin=ibegin, ni=ni,jbegin=jbegin,nj=nj, type='curvilinear') 
    9292  CALL xios_set_domain_attr("domain_atm",data_dim=2, data_ibegin=-1, data_ni=ni+2, data_jbegin=-2, data_nj=nj+4) 
    93   CALL xios_set_domain_attr("domain_atm",lonvalue=RESHAPE(lon,(/ni*nj/)),latvalue=RESHAPE(lat,(/ni*nj/))) 
     93  CALL xios_set_domain_attr("domain_atm",lonvalue_2D=lon,latvalue_2D=lat) 
    9494 
    9595  CALL xios_set_domain_attr("domain_atm_zoom",ni_glo=ni_glo, nj_glo=nj_glo, ibegin=ibegin, ni=ni,jbegin=jbegin,nj=nj) 
    9696  CALL xios_set_domain_attr("domain_atm_zoom",data_dim=2, data_ibegin=-1, data_ni=ni+2, data_jbegin=-2, data_nj=nj+4) 
    97   CALL xios_set_domain_attr("domain_atm_zoom",lonvalue=RESHAPE(lon,(/ni*nj/)),latvalue=RESHAPE(lat,(/ni*nj/))) 
     97  CALL xios_set_domain_attr("domain_atm_zoom",lonvalue_2D=lon,latvalue_2D=lat) 
    9898!  CALL xios_set_domain_attr("domain_atm_zoom",zoom_ibegin=39, zoom_ni=20, zoom_jbegin=39, zoom_nj=5) 
    9999 
     
    124124!!! Recupration des valeurs des longitudes et de taille des domaines locaux (pour test de fonctionnalité) 
    125125 
    126   ni=0 ; lonvalue(:)=0 
    127   CALL xios_get_domain_attr("domain_atm",ni=ni,lonvalue=lonvalue) 
     126  ni=0 ; lonvalue(:,:)=0 
     127  CALL xios_get_domain_attr("domain_atm",ni=ni,lonvalue_2D=lonvalue) 
    128128 
    129129  PRINT *,"ni",ni 
    130   PRINT *,"lonvalue",lonvalue ; 
     130  PRINT *,"lonvalue",lonvalue; 
    131131 
    132132!!! Fin de la definition du contexte 
     
    177177 
    178178  CALL xios_set_axis_attr("axis_srf",size=llm ,value=lval) ; 
    179   CALL xios_set_domain_attr("domain_srf",ni_glo=ni_glo, nj_glo=nj_glo, ibegin=ibegin, ni=ni,jbegin=jbegin,nj=nj) 
     179  CALL xios_set_domain_attr("domain_srf",ni_glo=ni_glo, nj_glo=nj_glo, ibegin=ibegin, ni=ni,jbegin=jbegin,nj=nj, type='curvilinear') 
    180180  CALL xios_set_domain_attr("domain_srf",data_dim=1, data_ibegin=0, data_ni=nb_pt) 
    181181  CALL xios_set_domain_attr("domain_srf",data_n_index=nb_pt, data_i_index=kindex) 
    182   CALL xios_set_domain_attr("domain_srf",lonvalue=RESHAPE(lon,(/ni*nj/)),latvalue=RESHAPE(lat,(/ni*nj/))) 
     182  CALL xios_set_domain_attr("domain_srf",lonvalue_2D=lon,latvalue_2D=lat) 
    183183 
    184184!!! Création d un nouveau champ 
     
    204204!!! Recupration des valeurs des longitudes et de taille des domaines locaux (pour test de fonctionnalité) 
    205205 
    206   ni=0 ; lonvalue(:)=0 
    207   CALL xios_get_domain_attr("domain_srf",ni=ni,lonvalue=lonvalue) 
     206  ni=0 ; lonvalue(:,:)=0 
     207  CALL xios_get_domain_attr("domain_srf",ni=ni,lonvalue_2D=lonvalue) 
    208208 
    209209  PRINT *,"ni",ni 
Note: See TracChangeset for help on using the changeset viewer.