trunk/libf/test_inter_barxy.f90

program test_inter_barxy

  use inter_barxy_m, only: inter_barxy
  use comconst, only: pi, initialize, dtvr, daysec
  use comgeom, only: rlonu, rlatv
  use dimens_m, only: iim, jjm
  use conf_gcm_m, only: conf_gcm, day_step
  use comvert, only: pa
  use inigeom_m, only: inigeom

  implicit none

  integer:: iml_dyn = 3, jml_dyn = 6
  real, allocatable:: lon_rad(:), lon_ini(:)
  real, allocatable:: lat_rad(:), lat_ini(:)
  integer i
  real, allocatable:: var_ana3d(:, :)
  REAL var_tmp2d(iim, jjm + 1)

  LOGICAL:: interbar = .true. ! barycentric interpolation
  namelist /main/iml_dyn, jml_dyn, interbar

  !------------------------

  print *, "Enter namelist 'main'."
  read (unit=*, nml=main)
  write(unit=*, nml=main)

  allocate(lon_rad(iml_dyn), lon_ini(iml_dyn))
  allocate(lat_rad(jml_dyn-1), lat_ini(jml_dyn))
  allocate(var_ana3d(iml_dyn, jml_dyn))

  call initialize
  CALL conf_gcm
  dtvr = daysec / real(day_step)
  print *, 'dtvr = ', dtvr
  pa = 5e4
  CALL iniconst
  CALL inigeom

  lon_ini = - pi + 2 * pi / iml_dyn * (/(i, i = 0, iml_dyn - 1)/)

  forall (i = 1: iml_dyn - 1) lon_rad(i) = (lon_ini(i) + lon_ini(i + 1)) / 2
  lon_rad(iml_dyn) = (lon_ini(iml_dyn) + pi) / 2

  lat_ini = pi / 2 - pi / (jml_dyn - 1) * (/(i, i = 0, jml_dyn - 1)/)
  forall (i = 1: jml_dyn - 1) lat_rad(i) = (lat_ini(i) + lat_ini(i + 1)) / 2

  !!forall (i = 1: jml_dyn) var_ana3d(:, i) = i
  forall (i = 1: jml_dyn) var_ana3d(:, i) = 240. + real(i) / 2

  call inter_barxy(lon_rad, lat_rad, var_ana3d, rlonu(:iim), rlatv, var_tmp2d)

  !!print *, "lon_rad * 180. / pi = ", lon_rad * 180. / pi
  print *, "lat_rad * 180. / pi = ", lat_rad * 180. / pi
  !!print *, "var_tmp2d = ", var_tmp2d

  print *, "minval(var_tmp2d, dim=1) = ", minval(var_tmp2d, dim=1)
  print *, "maxval(var_tmp2d, dim=1) = ", maxval(var_tmp2d, dim=1)

end program test_inter_barxy
1	guez	3	program test_inter_barxy
2
3			use inter_barxy_m, only: inter_barxy
4			use comconst, only: pi, initialize, dtvr, daysec
5			use comgeom, only: rlonu, rlatv
6			use dimens_m, only: iim, jjm
7			use conf_gcm_m, only: conf_gcm, day_step
8			use comvert, only: pa
9	guez	25	use inigeom_m, only: inigeom
10	guez	3
11			implicit none
12
13			integer:: iml_dyn = 3, jml_dyn = 6
14			real, allocatable:: lon_rad(:), lon_ini(:)
15			real, allocatable:: lat_rad(:), lat_ini(:)
16			integer i
17			real, allocatable:: var_ana3d(:, :)
18			REAL var_tmp2d(iim, jjm + 1)
19
20			LOGICAL:: interbar = .true. ! barycentric interpolation
21			namelist /main/iml_dyn, jml_dyn, interbar
22
23			!------------------------
24
25			print *, "Enter namelist 'main'."
26			read (unit=*, nml=main)
27			write(unit=*, nml=main)
28
29			allocate(lon_rad(iml_dyn), lon_ini(iml_dyn))
30			allocate(lat_rad(jml_dyn-1), lat_ini(jml_dyn))
31			allocate(var_ana3d(iml_dyn, jml_dyn))
32
33			call initialize
34			CALL conf_gcm
35			dtvr = daysec / real(day_step)
36			print *, 'dtvr = ', dtvr
37			pa = 5e4
38			CALL iniconst
39			CALL inigeom
40
41			lon_ini = - pi + 2 * pi / iml_dyn * (/(i, i = 0, iml_dyn - 1)/)
42
43			forall (i = 1: iml_dyn - 1) lon_rad(i) = (lon_ini(i) + lon_ini(i + 1)) / 2
44			lon_rad(iml_dyn) = (lon_ini(iml_dyn) + pi) / 2
45
46			lat_ini = pi / 2 - pi / (jml_dyn - 1) * (/(i, i = 0, jml_dyn - 1)/)
47			forall (i = 1: jml_dyn - 1) lat_rad(i) = (lat_ini(i) + lat_ini(i + 1)) / 2
48
49			!!forall (i = 1: jml_dyn) var_ana3d(:, i) = i
50			forall (i = 1: jml_dyn) var_ana3d(:, i) = 240. + real(i) / 2
51
52			call inter_barxy(lon_rad, lat_rad, var_ana3d, rlonu(:iim), rlatv, var_tmp2d)
53
54			!!print , "lon_rad 180. / pi = ", lon_rad * 180. / pi
55			print , "lat_rad 180. / pi = ", lat_rad * 180. / pi
56			!!print *, "var_tmp2d = ", var_tmp2d
57
58			print *, "minval(var_tmp2d, dim=1) = ", minval(var_tmp2d, dim=1)
59			print *, "maxval(var_tmp2d, dim=1) = ", maxval(var_tmp2d, dim=1)
60
61			end program test_inter_barxy