File: C:\NOAA\NEMS_11731\src\chem\gocart\src\GMAO_Shared\MAPL_Base\tests\testhorz.F90

1     use MAPL_HorzTransformMod
2     
3     implicit none
4     
5     integer, parameter :: imIN=144, imOUT=72
6     integer, parameter :: jmIN=91,  jmOUT=46
7     real, parameter :: pi=3.14159
8     integer :: rc
9     
10     real a(imin,jmin), b(imout,jmout)
11     integer i,j
12     
13     type(MAPL_HorzTransform) :: Trans
14     
15     
16     do j=1,jmin
17     do i=1,imin
18     a(i,j) = sin(4.*pi*(i-1)/imin)*sin(pi*(j-1)/(jmin-1))
19     end do
20     end do
21     
22     a(20:80,20:30) = -999.
23     
24     
25     open(10,file='inputarr',form='unformatted',status='unknown')
26     write(10) a
27     close(10)
28     
29     
30     call MAPL_HorzTransformCreate  (Trans,           &
31                                            imin,  jmin,   &
32                                            imout, jmout,  &
33     !                                       .false.,0,      &
34                                                         rc=rc)
35     
36     print *, rc
37     
38     call MAPL_HorzTransformRun(Trans,a,b,undef=-999.,rc=rc)
39     
40     open(10,file='outputarr',form='unformatted',status='unknown')
41     write(10) b
42     close(10)
43     
44     call MAPL_HorzTransformDestroy  (Trans)
45     
46     stop
47     end
48     
49