readfig.f

水文模型的原始代码

                if (inum3s(idum) == 0) then
                  write (40+inum1s(idum),5400) title
                  write (40+inum1s(idum),5500)
                else
                  iida = 0
                  iida = idaf
                  call xmon
                  write (40+inum1s(idum),5501) iyr, i_mo,               &
     &                                              (iida - ndays(i_mo))
                end if
              end if

            case (10) !! icode = 10  RECDAY command: read in daily values
                      !! with water in cms and rest in tons
              day_in = ""
              read (102,5100) day_in
              call caps(day_in)
              open (555+inum1s(idum),file=day_in,recl=350)
              do ii = 1, 6
                read (555+inum1s(idum),5200) titldum
              end do

            case (11) !! icode = 11  RECCNST command: read in average 
                      !! annual values with water in m^3, sed in t, and 
                      !! the nutrients in kg
            annual_in = ""
            read (102,5100) annual_in
            reccnstps(ihouts(idum))=annual_in(1:index(annual_in,'.')-1)
            call caps(annual_in)
            i = 0
            i = inum1s(idum)
            open (109,file=annual_in,recl=350)
            call readcnst

            case (14) !! icode = 14 SAVECONC command: saves hourly or 
                      !! daily reach constituent concentrations to an 
                      !! output file from a hydrograph node on the 
                      !! channel network
              day_in = ""
              read (102,5100) day_in
              call caps(day_in)
              if (inum1s(idum) <= 50 .and. inum1s(idum) > 0) then
                open (50+inum1s(idum),file=day_in,recl=350)
                write (50+inum1s(idum),5400) title
                write (50+inum1s(idum),5600)
              end if

            case (16) !! icode = 16 AUTOCAL command: identifies node
                      !! on channel network to compare to measured
                      !! data for autocalibration
              auto_in = ""
              read (102,5100) auto_in
              call caps(auto_in)
              close(9000+inum1s(idum))
              close(8000+inum1s(idum))
              if (inum1s(idum) <= 70 .and. inum1s(idum) > 0) then
                open (9000+inum1s(idum),file=auto_in,recl=350)
                select case (inum1s(idum))
                  case(1)
                   open (8000+inum1s(idum), file='autocal1.out')
                  case(2)
                   open (8000+inum1s(idum), file='autocal2.out')
                  case(3)
                   open (8000+inum1s(idum), file='autocal3.out')
                  case(4)
                   open (8000+inum1s(idum), file='autocal4.out')
                  case(5)
                   open (8000+inum1s(idum), file='autocal5.out')
                  case(6)
                   open (8000+inum1s(idum), file='autocal6.out')
                  case(7)
                   open (8000+inum 1s(idum), file='autocal7.out')
                  case(8)
                   open (8000+inum1s(idum), file='autocal8.out')
                  case(9)
                   open (8000+inum1s(idum), file='autocal9.out')
                  case(10)
                   open (8000+inum1s(idum), file='autocal10.out')
                  case(11)
                   open (8000+inum1s(idum), file='autocal11.out')
                  case(12)
                   open (8000+inum1s(idum), file='autocal12.out')
                  case(13)
                   open (8000+inum1s(idum), file='autocal13.out')
                  case(14)
                   open (8000+inum1s(idum), file='autocal14.out')
                  case(15)
                   open (8000+inum1s(idum), file='autocal15.out')
                  case(16)
                   open (8000+inum1s(idum), file='autocal16.out')
                  case(17)
                   open (8000+inum1s(idum), file='autocal17.out')
                  case(18)
                   open (8000+inum1s(idum), file='autocal18.out')
                  case(19)
                   open (8000+inum1s(idum), file='autocal19.out')
                  case(20)
                   open (8000+inum1s(idum), file='autocal20.out')
                  case(21)
                   open (8000+inum1s(idum), file='autocal21.out')
                  case(22)
                   open (8000+inum1s(idum), file='autocal22.out')
                  case(23)
                   open (8000+inum1s(idum), file='autocal23.out')
                  case(24)
                   open (8000+inum1s(idum), file='autocal24.out')
                  case(25)
                   open (8000+inum1s(idum), file='autocal25.out')
                  case(26)
                   open (8000+inum1s(idum), file='autocal26.out')
                  case(27)
                   open (8000+inum1s(idum), file='autocal27.out')
                  case(28)
                   open (8000+inum1s(idum), file='autocal28.out')
                  case(29)
                   open (8000+inum1s(idum), file='autocal29.out')
                  case(30)
                   open (8000+inum1s(idum), file='autocal30.out')
                  case(31)
                   open (8000+inum1s(idum), file='autocal31.out')
                  case(32)
                   open (8000+inum1s(idum), file='autocal32.out')
                  case(33)
                   open (8000+inum1s(idum), file='autocal33.out')
                  case(34)
                   open (8000+inum1s(idum), file='autocal34.out')
                  case(35)
                   open (8000+inum1s(idum), file='autocal35.out')
                  case(36)
                   open (8000+inum1s(idum), file='autocal36.out')
                  case(37)
                   open (8000+inum1s(idum), file='autocal37.out')
                  case(38)
                   open (8000+inum1s(idum), file='autocal38.out')
                  case(39)
                   open (8000+inum1s(idum), file='autocal39.out')
                  case(40)
                   open (8000+inum1s(idum), file='autocal40.out')
                  case(41)
                   open (8000+inum1s(idum), file='autocal41.out')
                  case(42)
                   open (8000+inum1s(idum), file='autocal42.out')
                  case(43)
                   open (8000+inum1s(idum), file='autocal43.out')
                  case(44)
                   open (8000+inum1s(idum), file='autocal44.out')
                  case(45)
                   open (8000+inum1s(idum), file='autocal45.out')
                  case(46)
                   open (8000+inum1s(idum), file='autocal46.out')
                  case(47)
                   open (8000+inum1s(idum), file='autocal47.out')
                  case(48)
                   open (8000+inum1s(idum), file='autocal48.out')
                  case(49)
                   open (8000+inum1s(idum), file='autocal49.out')
                  case(50)
                   open (8000+inum1s(idum), file='autocal50.out')
                  case(51)
                   open (8000+inum1s(idum), file='autocal51.out')
                  case(52)
                   open (8000+inum1s(idum), file='autocal52.out')
                  case(53)
                   open (8000+inum1s(idum), file='autocal53.out')
                  case(54)
                   open (8000+inum1s(idum), file='autocal54.out')
                  case(55)
                   open (8000+inum1s(idum), file='autocal55.out')
                  case(56)
                   open (8000+inum1s(idum), file='autocal56.out')
                  case(57)
                   open (8000+inum1s(idum), file='autocal57.out')
                  case(58)
                   open (8000+inum1s(idum), file='autocal58.out')
                  case(59)
                   open (8000+inum1s(idum), file='autocal59.out')
                  case(60)
                   open (8000+inum1s(idum), file='autocal60.out')
                  case(61)
                   open (8000+inum1s(idum), file='autocal61.out')
                  case(62)
                   open (8000+inum1s(idum), file='autocal62.out')
                  case(63)
                   open (8000+inum1s(idum), file='autocal63.out')
                  case(64)
                   open (8000+inum1s(idum), file='autocal64.out')
                  case(65)
                   open (8000+inum1s(idum), file='autocal65.out')
                  case(66)
                   open (8000+inum1s(idum), file='autocal66.out')
                  case(67)
                   open (8000+inum1s(idum), file='autocal67.out')
                  case(68)
                   open (8000+inum1s(idum), file='autocal68.out')
                  case(69)
                   open (8000+inum1s(idum), file='autocal69.out')
                  case(70)
                   open (8000+inum1s(idum), file='autocal70.out')
            end select
            end if
          end select

        end if
      end do

!! close .fig file
      close (102)


      return
!! isproj = 0
 5000 format (a1,9x,5i6,f6.1,i9)
 5001 format (7x,i3,4x,6f12.3)
 5002 format(a)
!! isproj = 2 (CEAP)
 5003 format (a1,9x,4i6,i5,f8.0,i8)
 5100 format (10x,2a13)
 5200 format (a80)
 5300 format (2i6)
 5400 format (20a4)
 5500 format (//," DAY YEAR HR ","  WATER m^3 ","   SED tons ",         &
     &         "    ORGN kg ","    ORGP kg ","   NO3-N kg ",            &
     &         "   NH3-N kg ","   NO2-N kg ","    MINP kg ",            &
     &         "    CBOD kg ","   DISOX kg ","    CHLA kg ",            &
     &         "  SOLPST mg ","  SORPST mg ","    BACTP # ",            &
     &         "   BACTLP # ","CMETAL#1 kg ","CMETAL#2 kg ",            &
     &         "CMETAL#3 kg ","  TEMP degC ")
 5501 format (i4,2i2,", Point Discharge File created with Save Comand")
 5600 format (//," Year  Day  Hr FLOWm^3/s    SEDmg/L   ORGNmg/L",      &
     &         "   ORGPmg/L    NO3mg/L    NH3mg/L    NO2mg/L",          &
     &         "   MINPmg/L   CBODmg/L  DISOXmg/L   CHLAug/L",          &
     &         " SOLPSTmg/L SORPSTmg/L  BACTPct/L BACTLPct/L",          &
     &         "CMETAL1mg/LCMETAL2mg/LCMETAL3mg/L   TEMPdegC")
      end