gen_be_read_regcoeffs.f90
References to this file elsewhere.
1 program gen_be_read_regcoeffs
2
3 !-------------------------------------------------------------------------------------------
4 ! Purpose: Read regression coefficients from BE file and write out in format
5 ! required for input to gen_be_stage2 (used for ep calculation).
6 !
7 ! Owner: Dale Barker
8 !
9 !-------------------------------------------------------------------------------------------
10
11 usa da_control, only : filename_len
12
13 implicit none
14
15 integer, parameter :: unit = 10 ! I/O unit.
16
17 character(len=filename_len) :: filename ! Input filename.
18 character*3 :: be_method ! Be method ('NMC', or 'ENS')
19 integer :: ni, nj, nk ! Dimensions read in.
20 integer :: bin_type ! Type of bin to average over. !!!DALE ADD.
21 integer :: num_bins ! Number of 3D bins.
22 integer :: num_bins2d ! Number of 2D bins.
23 real :: lat_min, lat_max ! Used if bin_type = 2 (degrees).
24 real :: binwidth_lat ! Used if bin_type = 2 (degrees). !!!DALE ADD..
25 real :: binwidth_hgt ! Used if bin_type = 2 (m). !!!DALE ADD..
26 real :: hgt_min, hgt_max ! Used if bin_type = 2 (m).
27
28 integer, allocatable:: bin(:,:,:) ! Bin assigned to each 3D point.
29 integer, allocatable:: bin2d(:,:) ! Bin assigned to each 2D point.
30 real, allocatable :: regcoeff1(:) ! psi/chi regression cooefficient.
31 real, allocatable :: regcoeff2(:,:) ! psi/ps regression cooefficient.
32 real, allocatable :: regcoeff3(:,:,:) ! psi/T regression cooefficient.
33
34 stderr = 0
35 stdout = 6
36
37 be_method = "ENS" ! Hardwired for now!
38
39 !----------------------------------------------------------------------------
40 ! [1] Read regression coefficients.
41 !----------------------------------------------------------------------------
42
43 filename = 'gen_be.'//trim(be_method)//'.dat'
44 open (unit, file = filename, form='unformatted')
45
46 ! Read the dimensions:
47 read(unit)ni, nj, nk
48
49 allocate( bin(1:ni,1:nj,1:nk) )
50 allocate( bin2d(1:ni,1:nj) )
51
52 ! Read bin info:
53 read(unit)bin_type
54 read(unit)lat_min, lat_max, binwidth_lat
55 read(unit)hgt_min, hgt_max, binwidth_hgt
56 read(unit)num_bins, num_bins2d
57 read(unit)bin(1:ni,1:nj,1:nk)
58 read(unit)bin2d(1:ni,1:nj)
59
60 ! Read the regression coefficients:
61 allocate( regcoeff1(1:num_bins) )
62 allocate( regcoeff2(1:nk,1:num_bins2d) )
63 allocate( regcoeff3(1:nk,1:nk,1:num_bins2d) )
64
65 read(unit)regcoeff1
66 read(unit)regcoeff2
67 read(unit)regcoeff3
68
69 close(unit)
70
71 !----------------------------------------------------------------------------
72 ! [2] Write regression coefficients.
73 !----------------------------------------------------------------------------
74
75 filename = 'gen_be_stage2.'//trim(be_method)//'.dat'
76 open (unit, file = filename, form='unformatted')
77 write(unit)ni, nj, nk
78 write(unit)num_bins, num_bins2d
79 write(unit)regcoeff1
80 write(unit)regcoeff2
81 write(unit)regcoeff3
82 close(unit)
83
84 deallocate( bin )
85 deallocate( bin2d )
86 deallocate( regcoeff1 )
87 deallocate( regcoeff2 )
88 deallocate( regcoeff3 )
89
90 end program gen_be_read_regcoeffs
91