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