Changeset 2064 for XIOS/xios_training/hands-on-3/test_tp3.f90
- Timestamp:
- 03/09/21 19:25:28 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
XIOS/xios_training/hands-on-3/test_tp3.f90
r2045 r2064 29 29 30 30 INTEGER :: ts 31 32 INTEGER :: ni, nj, ibegin, jbegin33 31 34 32 CALL MPI_INIT(ierr) … … 67 65 if(rank.EQ.0) print*, "domain size = ", ni_glo, "*", nj_glo 68 66 69 ni=ni_glo/size70 nj=nj_glo67 ALLOCATE(lonvalue(ni_glo)) 68 ALLOCATE(latvalue(nj_glo)) 71 69 72 ibegin=ni*rank 73 jbegin=0 74 75 CALL xios_set_domain_attr("domain_A", ni=ni, nj=nj, ibegin=ibegin,jbegin=jbegin) 76 print*, "Proc ",rank, "occupies subdomain",ibegin,"->",ibegin+ni,"X",jbegin,"->",jbegin+nj 77 78 79 ALLOCATE(lonvalue(ni)) 80 ALLOCATE(latvalue(nj)) 81 82 DO i=1,ni 83 lonvalue(i) = -180 + (rank*ni+i) * 360/ni_glo 70 DO i=1,ni_glo 71 lonvalue(i) = -180 + i * 360/ni_glo 84 72 ENDDO 85 73 86 DO j=1,nj 74 DO j=1,nj_glo 87 75 latvalue(j) = -90 + j * 180/nj_glo 88 76 ENDDO 89 77 90 78 CALL xios_set_domain_attr("domain_A", lonvalue_1d=lonvalue,latvalue_1d=latvalue) 91 92 !--- Case 2 ---! 93 ! CALL xios_set_domain_attr("domain_A", data_dim=2, data_ibegin=-1, data_ni=ni+2, data_jbegin=-1, data_nj=nj+2) 94 95 !--- Case 3 ---! 96 CALL xios_set_domain_attr("domain_A", data_dim=2, data_ibegin=1, data_ni=ni-2, data_jbegin=0, data_nj=nj/2) 97 79 98 80 CALL xios_close_context_definition() 99 81 100 !--- Case 1 ---! 101 ! ALLOCATE(field_A(ni, nj, axis_size)) 102 ! field_A(:,:,:)=rank 103 104 !--- Case 2 ---! 105 ! ALLOCATE(field_A(ni+2, nj+2, axis_size)) 106 ! field_A(:,:,:)=rank 107 108 !--- Case 3 ---! 109 ALLOCATE(field_A(ni-2, nj/2, axis_size)) 82 ALLOCATE(field_A(ni_glo, nj_glo, axis_size)) 110 83 field_A(:,:,:)=rank 111 84
Note: See TracChangeset
for help on using the changeset viewer.