Changeset 1025 for XIOS/dev/dev_olga/src/test/test_client.f90
- Timestamp:
- 01/11/17 15:14:22 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
XIOS/dev/dev_olga/src/test/test_client.f90
r987 r1025 17 17 CHARACTER(len=15) :: calendar_type 18 18 TYPE(xios_context) :: ctx_hdl 19 INTEGER,PARAMETER :: ni_glo= 220 INTEGER,PARAMETER :: nj_glo= 221 INTEGER,PARAMETER :: llm= 119 INTEGER,PARAMETER :: ni_glo=4 20 INTEGER,PARAMETER :: nj_glo=4 21 INTEGER,PARAMETER :: llm=5 22 22 DOUBLE PRECISION :: lval(llm)=1 23 23 TYPE(xios_field) :: field_hdl … … 28 28 DOUBLE PRECISION,DIMENSION(ni_glo,nj_glo) :: lon_glo,lat_glo 29 29 DOUBLE PRECISION :: field_A_glo(ni_glo,nj_glo,llm) 30 DOUBLE PRECISION,ALLOCATABLE :: lon(:,:),lat(:,:),field_A(:,:,:), lonvalue(:,:) ;30 DOUBLE PRECISION,ALLOCATABLE :: lon(:,:),lat(:,:),field_A(:,:,:), lonvalue(:,:), axisValue(:), field_domain(:,:) ; 31 31 INTEGER :: ni,ibegin,iend,nj,jbegin,jend 32 32 INTEGER :: i,j,l,ts,n … … 66 66 iend=ibegin+ni-1 ; jend=jbegin+nj-1 67 67 68 ALLOCATE(lon(ni,nj),lat(ni,nj),field_A(0:ni+1,-1:nj+2,llm),lonvalue(ni,nj) )68 ALLOCATE(lon(ni,nj),lat(ni,nj),field_A(0:ni+1,-1:nj+2,llm),lonvalue(ni,nj), axisValue(nj), field_domain(0:ni+1,-1:nj+2)) 69 69 lon(:,:)=lon_glo(ibegin+1:iend+1,jbegin+1:jend+1) 70 70 lat(:,:)=lat_glo(ibegin+1:iend+1,jbegin+1:jend+1) 71 71 field_A(1:ni,1:nj,:)=field_A_glo(ibegin+1:iend+1,jbegin+1:jend+1,:) 72 field_domain(1:ni,1:nj) = field_A_glo(ibegin+1:iend+1,jbegin+1:jend+1,1) 73 axisValue(1:nj)=field_A(1,1:nj,1); 72 74 73 75 CALL xios_context_initialize("test",comm) … … 78 80 PRINT *, "calendar_type = ", calendar_type 79 81 82 ! CALL xios_set_axis_attr("axis_A",n_glo=nj_glo ,value=axisValue, n=nj, begin=jbegin) ; 80 83 CALL xios_set_axis_attr("axis_A",n_glo=llm ,value=lval) ; 81 84 CALL xios_set_domain_attr("domain_A",ni_glo=ni_glo, nj_glo=nj_glo, ibegin=ibegin, ni=ni,jbegin=jbegin,nj=nj,type='curvilinear') … … 128 131 129 132 PRINT*,"field field_A is active ? ",xios_field_is_active("field_A") 130 DO ts=1, 24133 DO ts=1,40 131 134 CALL xios_update_calendar(ts) 132 CALL xios_send_field("field_A",field_A) 135 ! CALL xios_send_field("field_A",field_A) 136 ! CALL xios_send_field("field_Axis",axisValue) 137 ! CALL xios_send_field("field_Axis",lval) 138 CALL xios_send_field("field_Domain",field_domain) 133 139 CALL wait_us(5000) ; 134 140 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.