23 SUBROUTINE set_outflds(kth,th,kpv,pv)
31 use grib2_module, only: num_pset,pset,nrecout,first_grbtbl,grib_info_init
33 use ctlblk_mod
, only: npset, me, fld_info
34 use rqstfld_mod
, only: mxfld, iget, ritehd, lvlsxml, datset, ident, &
41 integer,
intent(in) :: kth,kpv
42 real,
intent(in) :: th(kth),pv(kpv)
44 integer l,ifld,mfld,iavbl,irec,i,j
45 CHARACTER*50 avblgrb_name
71 DO j=1,
size(lvlsxml,2)
72 DO i=1,
size(lvlsxml,1)
77 pset = paramset(npset)
79 if (me==0)print *,
'in SET_OUTFLDS, num_pset=',num_pset,
'datset=',trim(pset%datset),
'npset=',npset
87 mfld =
size(pset%param)
92 post_avblflds%param =>paramset(npset)%param
94 write(0,*)
'Size of pset is: ',mfld
95 write(0,*)
'datset is: ',datset
96 write(0,*)
'MXFLD is: ',mxfld
97 write(0,*)
'size of lvlsxml: ',
size(lvlsxml)
98 write(0,*)
'size of post_avblflds param',
size(post_avblflds%param)
100 if(
size(post_avblflds%param) <= 0)
then
101 write(0,*)
'WRONG: post available fields not ready!!!'
125 iavbl = post_avblflds%param(i)%post_avblfldidx
130 call set_lvlsxml(pset%param(i),ifld,irec,kpv,pv,kth,th)
147 allocate(fld_info(nrecout+100))
153 fld_info(i)%ntrange = 0
154 fld_info(i)%tinvstat = 0
156 if(me==0)
write(0,*)
'in readxml. nfld=',nfld,
'nrecout=',nrecout
182 if(me==0)print *,
'end of read_postcntrl_xml'