Changeset 2541

Timestamp:

08/29/23 16:48:44 (17 months ago)

Author:

ymipsl

Message:

for compatibility with ferret :

• change vertical axis units from Pa to mb
• invert lower an upper bounds for vertical axis (decreasing values)

YM

File:

• XIOS3/trunk/src/test/generic_testcase.f90 (modified) (4 diffs)

XIOS3/trunk/src/test/generic_testcase.f90

@@ -59 +59 @@
   INTEGER :: nb_procs(4)
   LOGICAL :: who_i_am(4)
+  INTEGER :: provided
    
   OPEN(unit=unit, file='param.def',status='old',iostat=ierr)
@@ -72 +73 @@
 !!! MPI Initialization
 
-  CALL MPI_INIT(ierr)
+  CALL MPI_INIT_THREAD(MPI_THREAD_SERIALIZED, provided, ierr)
+
   CALL init_wait
   CALL MPI_COMM_RANK(MPI_COMM_WORLD,rank,ierr)
@@ -2057 +2059 @@
     ENDDO
     
-    bounds_value_glo(2,0)=value_glo(0)-(value_glo(1)-value_glo(0))/2
+    bounds_value_glo(1,0)=value_glo(0)-(value_glo(1)-value_glo(0))/2
     DO i=1,nlev_glo-1
-     bounds_value_glo(2,i)=(value_glo(i-1)+value_glo(i)) /2
+     bounds_value_glo(1,i)=(value_glo(i-1)+value_glo(i)) /2
     ENDDO
 
     DO i=0,nlev_glo-2
-     bounds_value_glo(1,i)=(value_glo(i)+value_glo(i+1)) /2
-    ENDDO
-    bounds_value_glo(1,nlev_glo-1)=value_glo(nlev_glo-1)-(value_glo(nlev_glo-2)-value_glo(nlev_glo-1))/2
+     bounds_value_glo(2,i)=(value_glo(i)+value_glo(i+1)) /2
+    ENDDO
+    bounds_value_glo(2,nlev_glo-1)=value_glo(nlev_glo-1)-(value_glo(nlev_glo-2)-value_glo(nlev_glo-1))/2
 
     nlev=nlev_glo/axis_proc_size
@@ -2086 +2088 @@
     return_mask=.TRUE.
     CALL set_axis_mask(params,value,return_mask)    
-    CALL xios_set_axis_attr(axis_id, n_glo=nlev_glo, begin=begin, n=nlev, value=value*100000, mask=return_mask, bounds=bounds_value*100000, unit='Pa', positive='up')    
+    CALL xios_set_axis_attr(axis_id, n_glo=nlev_glo, begin=begin, n=nlev, value=value*1000, mask=return_mask, bounds=bounds_value*1000, unit='mb', positive='up')