Changeset 1980 for XIOS/trunk/src/test/generic_testcase.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/generic_testcase.f90 (modified) (8 diffs)

XIOS/trunk/src/test/generic_testcase.f90

@@ -1 @@
-PROGRAM generic_testcase
+PROGRAM generic_testcase 
   USE xios
   USE mod_wait
@@ -161 +161 @@
     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(:,:)
@@ -188 +188 @@
     INTEGER :: ierr      
 
-    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
@@ -290 +290 @@
     w=size(ensemble_value)
 
+
     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))
@@ -364 +366 @@
             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
@@ -382 +385 @@
     field_YZW(:,:,0) = field_YZ(:,:)*(1+0.1*ensemble_value(0))
     
-    
     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") ;
@@ -628 +631 @@
       CALL xios_update_calendar(ts)
 
+
       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)
@@ -2384 +2389 @@
   END SUBROUTINE get_decomposition
 
-END PROGRAM generic_testcase
-
-
+END PROGRAM generic_testcase 
+
+