Changeset 1980 for XIOS/trunk/src/test/test_interpolate.f90

Timestamp:

11/24/20 11:33:59 (4 years ago)

Author:

yushan

Message:

trunk : axis interpolate can have coordinate source (coordinate_src) and coordinate destination (coordinate_dst), while previous attribute coordinate compatible to source

File:

• XIOS/trunk/src/test/test_interpolate.f90 (modified) (6 diffs)

XIOS/trunk/src/test/test_interpolate.f90

@@ -162 +162 @@
     DOUBLE PRECISION, POINTER :: field3D_sub(:,:), other_field3D_sub(:,:)
     DOUBLE PRECISION, POINTER :: field3D_recv(:,:), other_field3D_recv(:,:)
-    DOUBLE PRECISION, POINTER :: pressure(:,:), other_pressure(:,:)
+    DOUBLE PRECISION, POINTER :: pressure(:,:), other_pressure(:,:), pressure_shifted(:,:)
 
     DOUBLE PRECISION, POINTER :: field2D_W(:,:), other_field2D_W(:,:)
@@ -190 +190 @@
 
     LOGICAL :: ok_field_in
-    LOGICAL :: ok_field2D, ok_field3D, ok_pressure, ok_field2D_sub, ok_field3D_sub,ok_field3D_recv, ok_field3D_send
+    LOGICAL :: ok_field2D, ok_field3D, ok_pressure_shifted, ok_pressure, ok_field2D_sub, ok_field3D_sub,ok_field3D_recv, ok_field3D_send
     LOGICAL :: ok_field_X, ok_field_Y, ok_field_XY, ok_field_Z, ok_field_XYZ, ok_field_XZ, ok_field_YZ
     LOGICAL :: ok_field2D_W, ok_field3D_W, ok_pressure_W, ok_field2D_sub_W, ok_field3D_sub_W,ok_field3D_recv_W, ok_field3D_send_W
@@ -296 +296 @@
     ALLOCATE(field2D(0:xy-1))
     ALLOCATE(field3D(0:xy-1,0:z-1))
+    ALLOCATE(pressure_shifted(0:xy-1,0:z-1))
     ALLOCATE(pressure(0:xy-1,0:z-1))
     ALLOCATE(field3D_recv(0:xy-1,0:z-1))
@@ -368 +369 @@
             dist=sqrt((lat(k)/90.)**2+(lon(k)/180.)**2) ;
             pressure(i,j)=pressure(i,j)*(1+params%pressure_factor*exp(-(4*dist)**2))
+            pressure_shifted(i,j)=pressure(i,j)+5000
           ENDIF
         ENDIF
@@ -386 +388 @@
     field_YZW(:,:,0) = field_YZ(:,:)*(1+0.1*ensemble_value(0))
     
-    ok_field_in=xios_is_valid_field("field_in") ;
+    ok_field_in=xios_is_valid_field("field_input") ;
+    ok_field_in=.FALSE.
     ok_field2D=xios_is_valid_field("field2D") ;
     ok_field3D=xios_is_valid_field("field3D") ;
+    ok_pressure_shifted=xios_is_valid_field("pressure_shifted") ;    
     ok_pressure=xios_is_valid_field("pressure") ;
     ok_field2D_sub=xios_is_valid_field("field2D_sub") ;
@@ -632 +636 @@
       CALL xios_update_calendar(ts)
 
-      if (ok_field_in) CALL xios_recv_field("field_in", field_in)
+      if (ok_field_in) CALL xios_recv_field("field_input", field_in)
       if (ok_field_in) CALL xios_send_field("field_out", field_in)
 
       IF (ok_field2D) CALL xios_send_field("field2D",field2D)
       IF (ok_field3D) CALL xios_send_field("field3D",field3D)
+      IF (ok_pressure_shifted) CALL xios_send_field("pressure_shifted",pressure_shifted)
       IF (ok_pressure) CALL xios_send_field("pressure",pressure)
       IF (ok_field_X) CALL xios_send_field("field_X",field_X)