input.f

Go to the documentation of this file.

C*********************************************************************
        SUBROUTINE runrd(NITER,KFYL1,KFYL2)
C*********************************************************************
C
C  ALLEN BECKER ADDED THIS TO READ VALUES OF MOST PARAMETERS FROM FILE 
C  'RUNDAT.DAT'.  THIS REPLACES THE OLD BLOCK DATA FILE.  
C   
C     FLUDWIG REVISED 5/97
C
        include 'NGRIDS.PAR'
C
        include 'ANCHOR.CMM'
        include 'FLOWER.CMM'
        include 'LIMITS.CMM'
        include 'STALOC.CMM'
C
C  OPEN RUN DESCRIPTORS.
C
C
        OPEN (22,file=chpath//'RUNSTF'//chmesh//'.DAT', 
     $          status='UNKNOWN',form='FORMATTED')
        READ (22,*)
        READ (22,6002) chpout
    chpvex=chpout
    WRITE (*,*) chpout
C
    READ (22,*) dbug,dbug2
    READ (22,*) dobin,doflat
C
C  READ NUMBER OF VERTICAL FLOW SURFACES 
C
        READ (22,*) nlvl
    IF (nlvl.GT.nzgrd .OR. nlvl .LT. 1) THEN
       WRITE (*,*) 'BAD NUMBER OF LEVELS', nlvl, nzgrd
       pause
       stop
    END IF
C
C  READ FLOW LEVELS (FRACTION OF TOP LEVEL)
C
        READ(22,*) (sigma(k),k=1,nlvl)
    WRITE (*,*) (sigma(k),k=1,nlvl)
C
        READ (22,*) nflat
    IF (nflat.GT.nhoriz .OR. nflat .LT. 1) THEN
       WRITE (*,*) 'BAD NUMBER OF LEVELS', nhoriz, nflat
       pause
       stop
    END IF
C
C  READ FLOW LEVELS IN M ABOVE LOW POINT
C
        READ(22,*) (zchooz(k),k=1,nflat)
C
C  READ FACTOR FOR CONVERTING TO M/S
C
    READ (22,*) spdcnv
C
C  READ NO. COLS,ROWS FOR GRID
C
        READ (22,*) nrow,ncol
    IF (nrow .GT. nygrd .OR. ncol .GT. nxgrd) THEN
       WRITE (*,*) 'BAD X,Y DIMENSIONS ARE ',ncol,nrow
       WRITE (*,*) 'THEY SHOULD NOT BE > ', nxgrd,nygrd
       pause
       stop
    END IF
C
C  READ GRID PT INTERVAL FOR ASCII OUTPUT -- IS ALL
C  GRID PTS, 2, EVRY OTHER ONE ETC.
C
        READ (22,*) nskip
        READ (22,*) nrytx, nryty
        READ (22,*) nascx0, nascy0
C
C  CHECK FOR OUT-OF-BOUNDS ARRAY STUFF.
C
    IF (nascx0-1+nrytx .GT. ncol) THEN
       WRITE (*,*) 'WRITE GRID RUNS OVER CALC. GRID -- X DIRECT.'
       nrytx=ncol-nascx0+1
       WRITE (*,*) 'NO OF OUTPUT COLUMNS SET TO ',nrytx
    END IF
    IF (nascy0-1+nryty .GT. ncol) THEN
       WRITE (*,*) 'WRITE GRID RUNS OVER CALC. GRID -- Y DIRECT.'
       nryty=nrow-nascy0+1
       WRITE (*,*) 'NO OF OUTPUT ROWS SET TO ',nryty
    END IF
C
C  READ GRID INTERVALS IN KM
C
        READ (22,*) dscrs
C
C  READ UTM COORDINATE OF REFERENCE POINT (1,1)--(X,Y) 
C
        READ (22,*) utmapx,utmapy
    WRITE (*,*) utmapx,utmapy
C
C  READ ANCHOR POINT GRID INDICES (IX,JY) 
C
        READ (22,*) kgridx,kgridy
    xorig=utmapx-dscrs*float(kgridx)
    yorig=utmapy-dscrs*float(kgridy)
C
C  CONSTANTS USED TO DEFINE SURFACE SHAPES SUBROUTINE RESIG
C  REDEFINES SURFACE SHAPES USING CRITICAL STREAMLINE METHODOLOGY.
C
        READ (22,*) avthk
        READ (22,*) slfac
        READ (22,*) cmpres
        READ (22,*) dpotmin
    WRITE (*,*) avthk,slfac,cmpres,dpotmin
C
C   CONSTANTS USED FOR INTERPOLATION -- D2MIN IS THE MINIMUM 
C       DISTANCE ALLOWED IN THE INVERSE WEIGHTING DENOMINATOR.  HT2DIS 
C       DEFINES THE RELATIVE IMPORT OF CHANGES IN ELEVATION VERSUS 
C       HORIZONTAL DISTANCE IN THE INVERSE DISTANCE WEIGHTING.  WHEN 
C       SET TO ZERO, VERTICAL TERRAIN EFFECT IS NOT INCLUDED.
C
C       READ ROUGHNESS LENGTH (ZZERO) AND ANEMOMETER HEIGHT (Z10)
C
        READ (22,*) zzero,z10
    z0=zzero
    zchooz(1)=z10
        READ (22,*) d2min
        READ (22,*) dtwt
C
C       NITER IS THE UPPER LIMIT ON ITERATIONS TOWARD NONDIVERGENCE IN 
C       BAL5 ADJMAX IS THE FRACTION OF THE USUAL ITERATIVE ADJUSTMENT 
C       TOWARD NONDIVERGENCE THAT IS MADE AT GRID POINTS NEAR 
C       OBSERVATIONS IN SUBROUTINE BAL5.
C
        READ (22,*) niter
        READ (22,*) adjmax
    WRITE (*,*) d2min,dtwt,niter,adjmax
C
C  READ INDICES OF LOW POINTS TO BE USED IN CRITICAL STREAMLINE 
C  CALCULATIONS FOR EACH TEMPERATURE SOUNDING.
C
    DO 97 lows=1,5 
       READ (22,*) lowix(lows),lowiy(lows)
97  CONTINUE
C
        READ (22,*) nend
    WRITE (*,*) 'number of files to do ', nend
C
C  EXTRA VARIABLES TO OUTPUT -- VERTICAL MOTION (W), THETA
C
    READ(22,*) dowcmp
    READ(22,*) dothet
    READ(22,*) dopres
    READ(22,*) dobri
    READ(22,*) dobvpd
C
C  DISABLED PROFILES & TIME SERIES INSTEAD OF ARRAYS
C****************************************************
C
C       READ (22,*) DOWVSZ
C       READ (22,*) NPFYLS
    dowvsz=.false.
    npfyls=0
    IF (npfyls.GT.nwndpt .AND. dowvsz) THEN
       WRITE(*,*) npfyls,' IS TOO MANY WIND PROFILE OUTPUTS!'
       stop
    END IF
C****************************************************
    IF (.NOT. dowvsz) THEN
       npfyls=0
    ELSE
       IF (nend.GT.111) THEN
          WRITE (*,*) 'TOO MANY CASES FOR PROFILE PLOTS.  ',
     $             'MUST BE 111 OR FEWER.'
              stop
       END IF
       DO 107 it=1,npfyls
          READ (22,*)  ptutmx,ptutmy
          jprylx(it)=1+nint((ptutmx-xorig)/dscrs)
          jpryly(it)=1+nint((ptutmy-yorig)/dscrs)
107    CONTINUE
    END IF
C
6002    FORMAT (a7)
6006    FORMAT (1x, 20f10.2)
6011    FORMAT (1x, 'SIGMAS: ',20f8.4)
C
    RETURN
C
        END