c subroutine getinp include "tgcmsat.h" namelist/inp/ satfile,histvols,mtimes,iden,fields,inttime, + height,lineplt,iwrascii c c Initialize inputs not set in blkdat: c mtimes(:,:) = ispval do i=1,mxvols call clearstr(histvols(i)) enddo c c Read inputs: c read(5,inp) c c Number of model times: c ntms = 0 do i=1,mxtms if (mtimes(1,i).ne.ispval) then ntms = ntms+1 if (mtimes(2,i).eq.ispval.or.mtimes(3,i).eq.ispval) then write(6,"('>>> mtimes must be given in triplets')") stop 'mtimes' endif else goto 100 endif enddo 100 continue c c Number of model volumes: c nhvols = 0 do i=1,mxvols if (lenstr(histvols(i)).gt.0) nhvols = nhvols+1 enddo c c Number of fields to get from history: c (always get z) c nfget = 0 do i=1,nfhist if (fields(i).gt.0.or.i.eq.ixz) then nfget = nfget+1 ifget(i) = nfget endif enddo c c Height at which to interpolate: c if (height.eq.0.) then write(6,"('>>> Please specify height -- height (km)', + ' at which to interpolate tgcm data')") stop 'height' endif return end