Changeset 2541


Ignore:
Timestamp:
08/29/23 16:48:44 (9 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:
1 edited

Legend:

Unmodified
Added
Removed
  • XIOS3/trunk/src/test/generic_testcase.f90

    r2520 r2541  
    5959  INTEGER :: nb_procs(4) 
    6060  LOGICAL :: who_i_am(4) 
     61  INTEGER :: provided 
    6162    
    6263  OPEN(unit=unit, file='param.def',status='old',iostat=ierr) 
     
    7273!!! MPI Initialization 
    7374 
    74   CALL MPI_INIT(ierr) 
     75  CALL MPI_INIT_THREAD(MPI_THREAD_SERIALIZED, provided, ierr) 
     76 
    7577  CALL init_wait 
    7678  CALL MPI_COMM_RANK(MPI_COMM_WORLD,rank,ierr) 
     
    20572059    ENDDO 
    20582060     
    2059     bounds_value_glo(2,0)=value_glo(0)-(value_glo(1)-value_glo(0))/2 
     2061    bounds_value_glo(1,0)=value_glo(0)-(value_glo(1)-value_glo(0))/2 
    20602062    DO i=1,nlev_glo-1 
    2061      bounds_value_glo(2,i)=(value_glo(i-1)+value_glo(i)) /2 
     2063     bounds_value_glo(1,i)=(value_glo(i-1)+value_glo(i)) /2 
    20622064    ENDDO 
    20632065 
    20642066    DO i=0,nlev_glo-2 
    2065      bounds_value_glo(1,i)=(value_glo(i)+value_glo(i+1)) /2 
    2066     ENDDO 
    2067     bounds_value_glo(1,nlev_glo-1)=value_glo(nlev_glo-1)-(value_glo(nlev_glo-2)-value_glo(nlev_glo-1))/2 
     2067     bounds_value_glo(2,i)=(value_glo(i)+value_glo(i+1)) /2 
     2068    ENDDO 
     2069    bounds_value_glo(2,nlev_glo-1)=value_glo(nlev_glo-1)-(value_glo(nlev_glo-2)-value_glo(nlev_glo-1))/2 
    20682070 
    20692071    nlev=nlev_glo/axis_proc_size 
     
    20862088    return_mask=.TRUE. 
    20872089    CALL set_axis_mask(params,value,return_mask)     
    2088     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')     
     2090    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')     
    20892091    
    20902092 
Note: See TracChangeset for help on using the changeset viewer.