Vgrid/exemples/vgd new/new read

From Wiki
Jump to: navigation, search
  1. program use_new_read
  2.    use vGrid_Descriptors, only: vgrid_descriptor,vgd_new,vgd_print,VGD_OK
  3.  
  4.    implicit none
  5.    
  6.    integer :: status,fnom,fstouv,iun=10
  7.    type(vgrid_descriptor) :: vgd
  8.  
  9.    status = fnom(iun,'../tests/data_Linux/dm_5002_from_model_run','RND+R/O',0)
  10.    status = fstouv(iun,'RND')
  11.  
  12.    status = vgd_new(vgd,iun,'fst')
  13.    
  14.    if (status /= VGD_OK)then
  15.       print*,'ERROR'
  16.       call exit(1)
  17.    endif
  18.  
  19.    status = vgd_print(vgd,6)
  20.  
  21. end program use_new_read