| 125 | == The 700nm problem == |
| 126 | |
| 127 | This is were the data is read from the files |
| 128 | {{{#! |
| 129 | C======================================================================= |
| 130 | |
| 131 | SUBROUTINE TPDINI( OBLECE ) |
| 132 | |
| 133 | C----------------------------------------------------------------------- |
| 134 | [...] |
| 135 | C READ THE TABLE OF ATMOSPHERIC EXTINCTION AND STORE IT IN THE |
| 136 | C ATMABS MATRIX. |
| 137 | [...] |
| 138 | DO I = 1, 105 |
| 139 | READ(MCERABS,21) WLT |
| 140 | 21 FORMAT(I4) |
| 141 | READ(MCERABS,*,ERR=995,END=995) (ATMABS(I,J), J=0,50) |
| 142 | ENDDO |
| 143 | [...] |
| 144 | C CALCULATE THE ATM. EXTINCTION FOR OBS. LEVEL INTERPOLATING THE TABLE |
| 145 | DO IWL = 1, 105 |
| 146 | FX0 = ATMABS(IWL,X0) |
| 147 | FX1 = ATMABS(IWL,X1) |
| 148 | ATEOBS(IWL) = LINEAR(X,X0,X1,FX0,FX1) |
| 149 | ENDDO |
| 150 | ENDIF |
| 151 | [...] |
| 152 | C READ THE QUANTUM EFFECIENCY OF THE PMT |
| 153 | OPEN(UNIT=MCERQEF,FILE=DATDIR(1:INDEX(DATDIR,' ')-1)// |
| 154 | * 'quanteff.dat',STATUS='OLD',FORM='FORMATTED',ERR=996) |
| 155 | READ(MCERQEF,20) TEXTQEF |
| 156 | READ(MCERQEF,23,ERR=997,END=997) (QUAEFF(I),I = 1,105) |
| 157 | 23 FORMAT(8F6.3) |
| 158 | CLOSE( MCERQEF ) |
| 159 | WRITE(MONIOU,22) TEXTQEF |
| 160 | ENDIF |
| 161 | [...] |
| 162 | C READ THE MIRROR REFLECTIVITY |
| 163 | OPEN(UNIT=MCERMIR,FILE=DATDIR(1:INDEX(DATDIR,' ')-1)// |
| 164 | * 'mirreff.dat',STATUS='OLD',FORM='FORMATTED',ERR=998) |
| 165 | READ(MCERMIR,20) TEXTREF |
| 166 | READ(MCERMIR,23,END=999,ERR=999) (MIRREF(I), I=1,105) |
| 167 | CLOSE( MCERMIR ) |
| 168 | WRITE(MONIOU,22) TEXTREF |
| 169 | ENDIF |
| 170 | }}} |
| 171 | |
| 172 | |
| 173 | This is where the index or the y-value and the corresponding x-values are calculated for mirror and pde efficiency. |
| 174 | {{{#!fortran |
| 175 | C===================================================================== |
| 176 | |
| 177 | SUBROUTINE TELEFF( ABSORB ) |
| 178 | |
| 179 | C----------------------------------------------------------------------- |
| 180 | [...] |
| 181 | C CALCULATE THE REFERENCE WL AND INDEX OF WL FOR THE INTERPOLATIONS |
| 182 | RIWL = 1 + INT( (WL-180.D0)/5.D0 ) |
| 183 | WLI0 = RIWL*5 + 175 |
| 184 | WLI1 = RIWL*5 + 175 + 5 |
| 185 | }}} |
| 186 | |
| 187 | This is the atmospheric absorption |
| 188 | {{{#!fortran |
| 189 | *-- Author : V. de Souza Filho, Uni. Campinas 22/06/1999 |
| 190 | C======================================================================= |
| 191 | |
| 192 | SUBROUTINE ATABSO( ABSORB ) |
| 193 | |
| 194 | C----------------------------------------------------------------------- |
| 195 | [...] |
| 196 | C CALCULATE THE REFERENCE WL AND INDEX OF WL FOR THE INTERPOLATIONS |
| 197 | RIWL = 1 + INT( (WL-180.D0)/5.D0 ) |
| 198 | WLI0 = RIWL*5 + 175 |
| 199 | WLI1 = RIWL*5 + 175 + 5 |
| 200 | }}} |
| 201 | |
| 202 | |