Changeset 1076 for XIOS/trunk/src/test

Timestamp:

03/16/17 15:04:53 (7 years ago)

Author:

mhnguyen

Message:

Correcting behavior of detecting_missing_value:

• Missing value detection is activated only when detecting_missing_value = true

and a default_value is defined.

• By default, undefined value by the computation of vertical (horizontal) interpolation will be NaN (not a number).

They are only converted to default_value if missing value detection is activated

Test

• On Curie
• Work

File:

• XIOS/trunk/src/test/test_remap.f90 (modified) (2 diffs)

XIOS/trunk/src/test/test_remap.f90

@@ -90 +90 @@
   DO i=1,src_ni
     src_field_3D(i,:) = src_field_2D(i)
-    IF ((23.5 < src_lat(i)) .AND. (src_lat(i) < 65.5) .AND. (0 < src_lon(i)) .AND. (src_lon(i) < 30)) THEN
+    IF ((0 < src_lat(i)) .AND. (src_lat(i) < 25.5) .AND. (30 < src_lon(i)) .AND. (src_lon(i) < 60)) THEN
       src_mask_2D(i)=.FALSE.
-      src_field_2D(i) = missing_value
+      ! src_field_2D(i) = missing_value
     ELSE
       src_mask_2D(i)=.TRUE.      
@@ -192 +192 @@
   CALL xios_get_domain_attr("src_domain_unstructured_read", ni=src_tmp_ni, nj=src_tmp_nj)
   ALLOCATE(tmp_field_2(src_tmp_ni*src_tmp_nj))
-
-  DO ts=1,1
-    CALL xios_recv_field("src_field_regular", tmp_field_0)
-    CALL xios_recv_field("src_field_curvilinear", tmp_field_1)
-    CALL xios_recv_field("src_field_unstructured", tmp_field_2)
+  
+  CALL xios_recv_field("src_field_regular", tmp_field_0)
+  CALL xios_recv_field("src_field_curvilinear", tmp_field_1)
+  CALL xios_recv_field("src_field_unstructured", tmp_field_2)
+
+  DO ts=1,10
     CALL xios_update_calendar(ts)
     CALL xios_send_field("src_field_2D",src_field_2D)