File: C:\NOAA\NEMS_11731\src\atmos\gfs\phys\info.f

1     c
2     c***********************************************************************
3     c
4           SUBROUTINE OUT_PARA(dt)
5     c
6     c***********************************************************************
7     c
8           use resol_def
9           use layout1
10           implicit none
11     c
12           real dt
13     c
14     c@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
15     c
16           IF (me.eq.0) THEN
17             write(6,*)
18             write(6,*)'This is the reduced-grid code with new phys'
19             write(6,9100) JCAP,LEVS
20             write(6,9110) NODES
21     cjfe    write(6,9115) NUM_PARTHDS()
22             write(6,9120) LEVS,LATG,LONF
23             write(6,9140) dt
24     c        write(6,9130)
25           ENDIF
26      
27      9100 format( /5x,'   T',i3,'-',i2,' FORECAST MODEL',
28          .        /5x,' ===========================')
29      9110 format(  5x,' Number of MPI Tasks:   ',i4)
30      9115 format(  5x,' Number of Threads:     ',i4)
31      9120 format( /5x,' Number of levels:      ',i4,
32          .        /5x,' Number of latitudes:   ',i4,
33          .        /5x,' Number of longitudes:  ',i4)
34      9130 format(//'                               PERFORMANCE'
35          .        /'                             Times in seconds'
36          .        /'                        Rates in Mbyte/s or Mflop/s'
37          .       //'         TRANSPOSE    GDSUM',
38          .         '         FFT       FLNS2FG      FL2FLN    FIDI/PHYS'
39          .        /' Step   Time  Rate  Time  Rate ',
40          .         ' Time  Rate  Time  Rate  Time  Rate  Time  Rate '
41          .        /' --------------------------------------',
42          .         '---------------------------------------')
43      9140 format( /5x,' Timestep:   ',F6.1,' seconds')
44      
45           RETURN
46           END
47     c
48     c***********************************************************************
49     c
50           SUBROUTINE BAR3(FE,FO,name,levl)
51     c
52     c***********************************************************************
53     c
54           use resol_def
55           use layout1
56           implicit none
57           real FE(len_trie_ls,2,levs),ffbar
58           real FO(len_trio_ls,2,levs)
59           integer n,l,joff,i,jlev,levl
60           character*3 name
61     c
62           do jlev=1,levl
63             ffbar=0.
64             DO i=1,len_trie_ls
65               ffbar=ffbar+FE(i,1,jlev)*FE(i,1,jlev)
66             ENDDO
67             DO i=1,len_trie_ls
68               ffbar=ffbar+FE(i,2,jlev)*FE(i,2,jlev)
69             ENDDO
70             ffbar=ffbar/2.
71             DO i=1,len_trio_ls
72               ffbar=ffbar+FO(i,1,jlev)*FO(i,1,jlev)
73             ENDDO
74             DO i=1,len_trio_ls
75               ffbar=ffbar+FO(i,2,jlev)*FO(i,2,jlev)
76             ENDDO
77             ffbar=sqrt(ffbar)
78             write(*,101)'rms ',name,' lev ',jlev,' = ',ffbar
79           enddo
80     c
81      101  format(3A,i3,A,ES24.17)
82     c
83           RETURN
84           END
85