topo.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	topo

Function/Subroutine Documentation

subroutine topo

(

)

Definition at line 3 of file topo.f.

Referenced by wocss().

C*********************************************************************
C
C   READ  AND  COMPUTE TOPOGRAPHY AT GRID POINTS.  GET
C   FIRST GUESS FLOW SURFACES.  RESIG WILL BE USED LATER TO DEFINE THEM
C   IN TERMS OF THE CRITICAL STREAMLINE PARAMETERS. 
C
    include 'NGRIDS.PAR'       
C
    include 'ANCHOR.CMM'
    include 'FLOWER.CMM'
    include 'LIMITS.CMM'
C
C  MAKE SURE THAT  THINGS READ ON 1ST PASS ARE THERE ON LATER CALLS
C
        SAVE
C
C  READ TERRAIN HEIGHT VALUES AT GRID POINTS IN METERS, ALL GRIDS
C
    WRITE (*,*) 'READ TOPORAPHY'
    sfclow=1.e6
    sfchi=-1.e6
C
C  IN DATA FILE, NORTHERN ROW IS FIRST.
C  READ HEIGHTS AT GRID POINTS & SET INITIAL BOUNDARY LAYER HT.
C  USE TERRAIN FOLLOWING 1ST; WILL RESET IN RESIG LATER.
C
    DO 28 jinvy=1,nrow
       jy=1+nrow-jinvy

       READ (11,6004) (sfcht(ix,jy),ix=1,ncol)
       IF (mod(jinvy,10).EQ.1) 
     $          WRITE (*,6004) (sfcht(ix,jy),ix=1,ncol,10)
C
C  GET ELEVATION AT ORIGIN
C
      zorig=sfcht(1,1)
          DO 18 ix=1,ncol
            IF(sfcht(ix,jy).LT.sfclow) sfclow=sfcht(ix,jy)
            IF(sfcht(ix,jy).GT. sfchi) sfchi =sfcht(ix,jy)
18    CONTINUE
28  CONTINUE
C
        zrise=sfchi-sfclow
        relht=avthk
C
C GET HEIGHT OF EACH 1ST GUESS SURFACE RELATIVE TO TERRAIN FOR EACH 
C GRID PT.  CHANGE IN HT RELATIVE TO THE LOWEST HT (MSL) IS ASSUMED 
C PROPORTIONAL TO SIGMA FOR THAT SFC.
C
    DO 67 jy=1,nrow
           DO 67 ix=1,ncol
              DO 65 kz=1,nlvl
C
         rhs(ix,jy,kz) =sigma(kz)*avthk- 
     $                   (sfcht(ix,jy)-sfclow)*(1.0-slfac)
65         CONTINUE
67  CONTINUE
C
6001    FORMAT(f3.0,107f4.0)
6002    FORMAT (1x, 45f5.0)
6003    FORMAT (1x, 161f5.0)
6004    FORMAT (1x, 180f5.0)
C
        RETURN
C

Here is the caller graph for this function: