gsl_oro_data.f90

Go to the documentation of this file.

 
 
program gsl_oro_data
 
use omp_lib
 
use gsl_oro_data_sm_scale, only: calc_gsl_oro_data_sm_scale
use gsl_oro_data_lg_scale, only: calc_gsl_oro_data_lg_scale
 
implicit none
 
character(len=2) :: tile_num   ! tile number entered by user
character(len=7) :: res_indx   ! grid-resolution index, e.g., 96, 192, 384, 768,
                               ! etc. entered by user
character(len=4) :: halo       ! halo value entered by user (for input grid data)
 
logical :: duplicate_oro_data_file   ! flag for whether oro_data_ls file is a duplicate
                   ! of oro_data_ss due to minimum grid size being less than 7.5km
 
integer :: tid, nthreads
 
! Read in FV3GFS grid info
print *
print *, "Enter tile number:"
read (5,*) tile_num
print *
print *, "Enter grid-resolution index:"
read (5,*) res_indx
print *
print *, "Enter halo number (-999 for no halo):"
read (5,*) halo
print *
print *, "Creating tile oro_data for tile number: ", tile_num
print *, "Grid resolution = ", res_indx
print *, "Halo = ", halo
print *
 
!$OMP PARALLEL PRIVATE(TID)
  tid = omp_get_thread_num()
  if (tid==0) then
    nthreads = omp_get_num_threads()
    print*,'Number of threads = ', nthreads
  endif
!$OMP END PARALLEL
 
call calc_gsl_oro_data_sm_scale(tile_num,res_indx,halo,duplicate_oro_data_file)
 
print *, "duplicate_oro_data_file =", duplicate_oro_data_file
print *
 
if ( .not.duplicate_oro_data_file ) then
   call calc_gsl_oro_data_lg_scale(tile_num,res_indx,halo)
end if
 
 
print *
print *, "End program gsl_oro_data"
print *
 
 
end program gsl_oro_data