Changeset 492


Ignore:
Timestamp:
10/18/16 11:15:43 (8 years ago)
Author:
ymipsl
Message:

Some bugs fixed in XIOS...
=> Some modifciation in restart and read fonctionnalities.

YM

Location:
codes/icosagcm/trunk/src
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • codes/icosagcm/trunk/src/restart.f90

    r488 r492  
    2828        CALL xios_set_file_attr("restart",name=TRIM(ADJUSTL(restart_file_name)),output_freq=itaumax*xios_timestep) 
    2929        CALL xios_set_fieldgroup_attr("group_restart", freq_op=itaumax*xios_timestep) 
    30         CALL xios_set_axis_attr("lev_read", n_glo=llm) 
     30!        CALL xios_set_axis_attr("lev_read", n_glo=llm) 
    3131      ENDIF 
    3232    ENDIF 
     
    420420  INTEGER    :: status 
    421421  REAL(rstd),ALLOCATABLE :: lon(:),lat(:),bounds_lon(:,:),bounds_lat(:,:) 
    422   REAL(rstd) :: it_temp(1) 
     422  REAL(rstd) :: it_real 
    423423   
    424424    IF (no_io) RETURN 
     
    450450      IF (PRESENT(field18))  THEN ; CALL  xios_read_field(TRIM(field18(1)%name)//'_start',field18)  ; ENDIF 
    451451      IF (PRESENT(field19))  THEN ; CALL  xios_read_field(TRIM(field19(1)%name)//'_start',field19)  ; ENDIF 
    452 ! doesn't work for now, to be decomment when xios is ok 
    453 !      CALL xios_recv_field("it_start",it_temp) 
    454 !      it=it_temp(1) 
    455       it=0 
     452 
     453      CALL xios_recv_field("it_start",it_real) 
     454      it=it_real 
    456455    ELSE 
    457456 
  • codes/icosagcm/trunk/src/xios_mod.F90

    r483 r492  
    5252  REAL(rstd) :: lev_value(llm) 
    5353  REAL(rstd) :: lev_valuep1(llm+1) 
     54  REAL(rstd) :: nq_value(nqtot) 
    5455  INTEGER :: ncell, ncell_tot, ncell_glo(0:mpi_size-1), displ 
    5556  INTEGER :: ind, i,j,k,l,ij 
     
    6566   lev_value(:) = (/ (l,l=1,llm) /)      
    6667   lev_valuep1(:) = (/ (l,l=1,llm+1) /)      
     68   nq_value(:) = (/ (l,l=1,nqtot) /)      
    6769   CALL xios_set_axis_attr("lev",n_glo=llm ,value=lev_value) ; 
    6870   CALL xios_set_axis_attr("levp1",n_glo=llm+1 ,value=lev_valuep1) ; 
    69    CALL xios_set_axis_attr("nq",n_glo=nqtot) ; 
     71   CALL xios_set_axis_attr("nq",n_glo=nqtot, value=nq_value) ; 
    7072    
    7173   ncell=0 
Note: See TracChangeset for help on using the changeset viewer.