Ignore:
Timestamp:
07/10/17 15:18:16 (4 years ago)
Author:
yushan
Message:

prep to merge with trunk @1200

File:
1 edited

Legend:

Unmodified
Added
Removed
  • XIOS/dev/branch_yushan_merged/src/test/test_remap_omp.f90

    r1179 r1203  
    183183  ierr=NF90_GET_VAR(ncid,varid, dst_boundslat, start=(/1,dst_ibegin+1/),count=(/dst_nvertex,dst_ni/)) 
    184184 
    185  
     185   !$omp barrier 
     186 
     187  !$omp master  
     188  CALL MPI_barrier(comm, ierr) 
     189  !$omp end master 
     190 
     191  !$omp barrier 
     192 
     193   
    186194  CALL xios_context_initialize("test",comm) 
    187195  CALL xios_get_handle("test",ctx_hdl) 
     
    214222  CALL xios_close_context_definition() 
    215223 
    216   CALL xios_get_domain_attr("src_domain_regular_read", ni=src_tmp_ni, nj=src_tmp_nj) 
    217   ALLOCATE(tmp_field_0(src_tmp_ni*src_tmp_nj)) 
    218  
    219   CALL xios_get_axis_attr("src_axis_curvilinear_read", n=src_tmp_n) 
    220   CALL xios_get_domain_attr("src_domain_curvilinear_read", ni=src_tmp_ni, nj=src_tmp_nj) 
    221   ALLOCATE(tmp_field_1(src_tmp_ni*src_tmp_nj*src_tmp_n)) 
    222  
    223   CALL xios_get_domain_attr("src_domain_unstructured_read", ni=src_tmp_ni, nj=src_tmp_nj) 
    224   ALLOCATE(tmp_field_2(src_tmp_ni*src_tmp_nj)) 
    225    
    226   CALL xios_recv_field("src_field_regular", tmp_field_0) 
    227   CALL xios_recv_field("src_field_curvilinear", tmp_field_1) 
    228   CALL xios_recv_field("src_field_unstructured", tmp_field_2) 
     224!  CALL xios_get_domain_attr("src_domain_regular_read", ni=src_tmp_ni, nj=src_tmp_nj) 
     225!  ALLOCATE(tmp_field_0(src_tmp_ni*src_tmp_nj)) 
     226 
     227!  CALL xios_get_axis_attr("src_axis_curvilinear_read", n=src_tmp_n) 
     228!  CALL xios_get_domain_attr("src_domain_curvilinear_read", ni=src_tmp_ni, nj=src_tmp_nj) 
     229!  ALLOCATE(tmp_field_1(src_tmp_ni*src_tmp_nj*src_tmp_n)) 
     230 
     231!  CALL xios_get_domain_attr("src_domain_unstructured_read", ni=src_tmp_ni, nj=src_tmp_nj) 
     232!  ALLOCATE(tmp_field_2(src_tmp_ni*src_tmp_nj)) 
     233   
     234!  CALL xios_recv_field("src_field_regular", tmp_field_0) 
     235!  CALL xios_recv_field("src_field_curvilinear", tmp_field_1) 
     236!  CALL xios_recv_field("src_field_unstructured", tmp_field_2) 
    229237 
    230238  DO ts=1,10 
     
    232240    CALL xios_send_field("src_field_2D",src_field_2D) 
    233241     
    234     !DO i=1,src_ni 
    235     !  src_field_2D_clone(i) = src_field_2D(i) 
    236     !  IF ((23.5 * ts < src_lat(i)) .AND. (src_lat(i) < 65.5 *ts) .AND. (0 < src_lon(i)) .AND. (src_lon(i) < 30*ts)) THEN       
    237     !    src_field_2D_clone(i) = missing_value     
    238     !  ENDIF 
    239     !ENDDO 
    240  
    241     !CALL xios_send_field("src_field_2D_clone",src_field_2D_clone) 
    242     !CALL xios_send_field("src_field_3D",src_field_3D) 
    243     !CALL xios_send_field("src_field_3D_clone",src_field_3D) 
    244     !CALL xios_send_field("src_field_4D",src_field_4D) 
    245     !CALL xios_send_field("src_field_3D_pression",src_field_pression) 
    246     CALL xios_send_field("tmp_field_0",tmp_field_0) 
    247     CALL xios_send_field("tmp_field_1",tmp_field_1) 
    248     CALL xios_send_field("tmp_field_2",tmp_field_2) 
     242    DO i=1,src_ni 
     243      src_field_2D_clone(i) = src_field_2D(i) 
     244      IF ((23.5 * ts < src_lat(i)) .AND. (src_lat(i) < 65.5 *ts) .AND. (0 < src_lon(i)) .AND. (src_lon(i) < 30*ts)) THEN       
     245        src_field_2D_clone(i) = missing_value     
     246      ENDIF 
     247    ENDDO 
     248 
     249    CALL xios_send_field("src_field_2D_clone",src_field_2D_clone) 
     250    CALL xios_send_field("src_field_3D",src_field_3D) 
     251    CALL xios_send_field("src_field_3D_clone",src_field_3D) 
     252    CALL xios_send_field("src_field_4D",src_field_4D) 
     253    CALL xios_send_field("src_field_3D_pression",src_field_pression) 
     254 !   CALL xios_send_field("tmp_field_0",tmp_field_0) 
     255 !   CALL xios_send_field("tmp_field_1",tmp_field_1) 
     256 !   CALL xios_send_field("tmp_field_2",tmp_field_2) 
    249257    CALL wait_us(5000) ; 
    250258   ENDDO 
     
    254262  DEALLOCATE(src_lon, src_lat, src_boundslon,src_boundslat, src_field_2D) 
    255263  DEALLOCATE(dst_lon, dst_lat, dst_boundslon,dst_boundslat) 
    256   DEALLOCATE(tmp_field_0, tmp_field_1, tmp_field_2) 
     264  !DEALLOCATE(tmp_field_0, tmp_field_1, tmp_field_2) 
    257265   
    258266  CALL xios_finalize() 
Note: See TracChangeset for help on using the changeset viewer.