Page 1 Source Listing FLIGHTS_ALLNORMS 2022-05-02 15:46 flights_allnorms.f 1 c ********************************************************************** 2 c ********************************************************************** 3 c ******* ETKF code to calculate total signal variance in the VR ******* 4 c ********** for any flight tracks 'track.01', 'track.02' etc ********** 5 c ********************************************************************** 6 c ********************************************************************** 7 c ***** (c) S.J.Majumdar, C.H.Bishop, B.J.Etherton, December 1999. ***** 8 c *************** Code developed at Penn State University ************** 9 c ********************************************************************** 10 c ********* Thanks to Bob Kohler of HRD for eigenvalue solvers ********* 11 c ********************************************************************** 12 c ********************************************************************** 13 14 15 c 3456789012345678901234567890123456789012345678901234567890123456789012 16 17 program flights_allnorms 18 19 allocatable sistr(:) 20 allocatable sigtop(:) 21 double precision,allocatable::xvvest(:,:) 22 double precision,allocatable::cmat(:,:) 23 double precision,allocatable::ceval(:) 24 double precision,allocatable::yvvest(:,:) 25 double precision,allocatable::ytvec(:,:) 26 double precision,allocatable::hy(:,:) 27 double precision,allocatable::hyt(:,:) 28 double precision,allocatable::obs(:,:,:) 29 double precision,allocatable::bsig(:,:) 30 double precision,allocatable::betasig(:) 31 double precision,allocatable::bsig2(:,:) 32 double precision,allocatable::svpvr(:,:) 33 double precision,allocatable::work3(:,:) 34 double precision sig(100),fev_routine 35 double precision sigold(100) 36 integer eid(100),ira(100) 37 character*2 tfl 38 real uvsd,uv850sd,prcpsd,mslpsd,z500sd 39 40 read(5,*)idim,jdim,jstr,mem,vrlon,vrlat,radvr,nd,nd9, 41 &nd12,nflights,georad,ne1,n2,nvr,nvr12,nv,ne9,ltcode, 42 &casecode,ivnorm,iphase 43 44 allocate(sistr(jstr)) 45 allocate(sigtop(jstr)) 46 allocate(xvvest(ne1,mem)) 47 allocate(cmat(mem,mem)) 48 allocate(ceval(mem)) 49 allocate(yvvest(nvr12,mem)) 50 allocate(ytvec(ne9,mem)) 51 allocate(hy(nd9,mem)) 52 allocate(hyt(mem,nd9)) 53 allocate(obs(idim,jdim,12)) 54 allocate(bsig(mem,mem)) 55 allocate(betasig(mem)) 56 allocate(bsig2(mem,mem)) 57 allocate(svpvr(nvr12,mem)) Page 2 Source Listing FLIGHTS_ALLNORMS 2022-05-02 15:46 flights_allnorms.f 58 allocate(work3(mem,mem)) 59 60 ivr=nint(vrlon/2.5)+1 61 jvr=nint((90.-vrlat)/2.5)+1 62 63 knum1=ltcode+8200 64 knum2=ltcode+8300 65 knum3=ltcode+8400 66 knum5=casecode+8600 67 68 read(knum1) ytvec 69 read(knum2) ceval 70 read(knum3) cmat 71 read(knum5) xvvest 72 73 pi=acos(-1.0) 74 dpsieq=1./36. 75 do j=1,jstr 76 jeq=37-j 77 sistr(j)=asin(jeq*dpsieq) 78 sigtop(j)=(pi/2.)-sistr(j) 79 end do 80 81 knum0=1234 82 open(knum0,form='formatted') 83 read(knum0,*) uvsd 84 read(knum0,*) uv850sd 85 read(knum0,*) prcpsd 86 read(knum0,*) mslpsd 87 read(knum0,*) z500sd 88 close(knum0) 89 90 call calcvrci(ivr,jvr,sistr,sigtop,cmat,ivm,xvvest,yvvest, 91 &jstr,ne1,mem,nvr12,idim,jdim,radvr,georad) 92 call calc_fev_routine(ivm,yvvest,fev_routine, 93 & uvsd,uv850sd,prcpsd,mslpsd,z500sd, 94 & nvr12,mem,ivnorm) 95 96 do 125 i=1,idim 97 do 125 j=1,jdim 98 obs(i,j,1)=2.4**2 99 obs(i,j,2)=2.8**2 100 obs(i,j,3)=2.95**2 101 obs(i,j,4)=2.4**2 102 obs(i,j,5)=2.8**2 103 obs(i,j,6)=2.95**2 104 obs(i,j,7)=0.8**2 105 obs(i,j,8)=0.8**2 106 obs(i,j,9)=1.2**2 107 125 continue 108 109 sigmax=0.0 110 sigmin=1.0e5 111 112 * -------------------------------------------------------- 113 * START INTRODUCING DIFFERENT SETS OF OBSERVATIONS 114 * -------------------------------------------------------- Page 3 Source Listing FLIGHTS_ALLNORMS 2022-05-02 15:46 flights_allnorms.f 115 116 do 10101 iflight=1,nflights 117 118 hy=0.0 119 hyt=0.0 120 betasig=0.0 121 bsig=0.0 122 123 call flight2(iflight,obs,hy,ytvec,nd,ne9,nd9,mem,idim, 124 &jdim,nv) 125 * ------------------------------------------------------------- 126 * Time to calculate the MEM eigenvalues [BETA_sigma] and 127 * eigenvectors [B_sigma] of Hstar*P*Hstar(T) 128 * ------------------------------------------------------------- 129 CALL DTRANS(nd9,MEM,hy,hyt) 130 CALL DMRRRR(MEM,nd9,hyt,MEM,nd9,MEM,hy,nd9,MEM,MEM,WORK3,MEM) 131 CALL DEVCSF(MEM,WORK3,MEM,BETASIG,BSIG,MEM) 132 133 do 1231 nm1=1,mem 134 if (BETASIG(nm1).le.1.0e-12)BETASIG(nm1)=1.0e-15 135 1231 continue 136 137 do 123 nm1=1,mem 138 do 123 nm2=1,mem 139 BSIG2(nm1,nm2)= 140 & BSIG(nm1,nm2)*sqrt(BETASIG(nm2)/(1.+BETASIG(nm2))) 141 123 continue 142 143 * ---------------------------------------------------------------- 144 * Calculate signal variance variance as function of target region 145 * ---------------------------------------------------------------- 146 CALL DMRRRR(nvr12,MEM,yvvest,nvr12,MEM,MEM,BSIG2,MEM, 147 & nvr12,MEM,svpvr,nvr12) 148 sig(iflight)=0.0 149 150 151 * ---------------------------------------------------------------- 152 * HERE IS WHERE WE DEFINE THE NORM !! 153 * ---------------------------------------------------------------- 154 sig(iflight)=0.0 155 do 61 ivcount=1,ivm 156 do 61 k=1,mem 157 iu850=(ivcount-1)*12 + 1 158 iv850=(ivcount-1)*12 + 2 159 iu500=(ivcount-1)*12 + 3 160 iv500=(ivcount-1)*12 + 4 161 iu200=(ivcount-1)*12 + 5 162 iv200=(ivcount-1)*12 + 6 163 it850=(ivcount-1)*12 + 7 164 it500=(ivcount-1)*12 + 8 165 it200=(ivcount-1)*12 + 9 166 iprcp=(ivcount-1)*12 + 10 167 imslp=(ivcount-1)*12 + 11 168 iz500=(ivcount-1)*12 + 12 169 170 * ****** (u,v) verifying norm ****** 171 if (ivnorm.eq.1) then Page 4 Source Listing FLIGHTS_ALLNORMS 2022-05-02 15:46 flights_allnorms.f 172 sig(iflight)=sig(iflight) + 173 & svpvr(iu850,k)**2 + svpvr(iv850,k)**2 + 174 & svpvr(iu500,k)**2 + svpvr(iv500,k)**2 + 175 & svpvr(iu200,k)**2 + svpvr(iv200,k)**2 176 endif 177 178 * ****** Total energy verifying norm ****** 179 if (ivnorm.eq.2) then 180 sig(iflight)=sig(iflight) + 181 & 0.5*( svpvr(iu850,k)**2 + svpvr(iv850,k)**2 + 182 & svpvr(iu500,k)**2 + svpvr(iv500,k)**2 + 183 & svpvr(iu200,k)**2 + svpvr(iv200,k)**2 ) + 184 & 3.346*( svpvr(it850,k)**2 + svpvr(it500,k)**2 + 185 & svpvr(it200,k)**2 ) 186 endif 187 188 * ****** (u850,v850,precip,mslp) verifying norm ****** 189 if (ivnorm.eq.3) then 190 sig(iflight)=sig(iflight) + 191 & (svpvr(iu850,k)**2 + svpvr(iv850,k)**2)/uv850sd**2 + 192 & (svpvr(iprcp,k)/prcpsd)**2 + (svpvr(imslp,k)/mslpsd)**2 193 endif 194 195 61 continue 196 197 if (ivnorm.eq.1) then 198 sig(iflight)=sig(iflight)/(mem*ivm*3.0) 199 endif 200 201 if (ivnorm.eq.2) then 202 sig(iflight)=sig(iflight)/(mem*ivm*3.0) 203 endif 204 205 if (ivnorm.eq.3) then 206 sig(iflight)=sig(iflight)/(mem*ivm) 207 endif 208 209 write(6,99) 'Signal Variance(',iflight,')=',sig(iflight), 210 & ' fcst.err.var of routine network=',fev_routine 211 212 if (sig(iflight).gt.sigmax) then 213 sigmax=sig(iflight) 214 endif 215 216 if (sig(iflight).lt.sigmin.and.sig(iflight).gt.0.1) then 217 sigmin=sig(iflight) 218 endif 219 220 10101 continue 221 222 open(1010,form='formatted') 223 do 602 iflight=1,nflights 224 write(1010,666)iflight,sig(iflight) 225 602 continue 226 close(1010) 227 228 call bubble_sort(sig,sigold,eid,ira,nflights) Page 5 Source Listing FLIGHTS_ALLNORMS 2022-05-02 15:46 flights_allnorms.f 229 C Y.C.Song changed for P-3 flight that we will select the track 230 C from each set of tracks (1-30,31-56,57-70) 231 C 1-30 from Alaska, 31-56 from Haiwaii, 61 - 66 for Gulf 232 C 67 - 90 from Japan 233 234 flag1=0 235 flag2=0 236 flag3=0 237 do 49 i=1,nflights 238 if ( eid(i).le.30.and.flag1.eq.0) then 239 eid(1)=eid(i) 240 flag1=1 241 endif 242 243 if (iphase.eq.2) then 244 if (eid(i).gt.30.and.eid(i).lt.57.and.flag2.eq.0) then 245 eid(2)=eid(i) 246 flag2=1 247 endif 248 else 249 if (eid(i).gt.67.and.flag2.eq.0) then 250 eid(2)=eid(i) 251 flag2=1 252 endif 253 endif 254 255 if (eid(i).gt.60.and.eid(i).lt.67.and.flag3.eq.0) then 256 eid(3)=eid(i) 257 flag3=1 258 endif 259 260 49 continue 261 262 263 499 write(1040,*) 5* int(sigmax/5.0) + 5 264 write(1041,*) 5* int(sigmin/5.0) 265 266 do 55 i=1,3 267 tfl='00' 268 if (eid(i).ge.10) then 269 write(tfl(1:2),'(i2)')eid(i) 270 else 271 write(tfl(2:2),'(i1)')eid(i) 272 endif 273 knum=1016+i 274 write(knum,*) tfl 275 55 continue 276 277 99 format (A16,I2,A2,F7.2,A35,F7.2) 278 666 format(i5,f8.3) 279 280 STOP 281 END Page 6 Source Listing FLIGHTS_ALLNORMS 2022-05-02 15:46 Entry Points flights_allnorms.f ENTRY POINTS Name MAIN__ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 10101 Label 220 116 123 Label 141 137,138 1231 Label 135 133 125 Label 107 96,97 49 Label 260 237 499 Label 263 55 Label 275 266 602 Label 225 223 61 Label 195 155,156 666 Label 278 224 99 Label 277 209 ACOS Func 73 scalar 73 ASIN Func 77 scalar 77 BETASIG Local 30 R(8) 8 1 1 ALC 55,120,131,134,140 BSIG Local 29 R(8) 8 2 1 ALC 54,121,131,140 BSIG2 Local 31 R(8) 8 2 1 ALC 56,139,146 BUBBLE_SORT Subr 228 228 CALCVRCI Subr 90 90 CALC_FEV_ROUTINE Subr 92 92 CASECODE Local 42 R(4) 4 scalar 42,66 CEVAL Local 23 R(8) 8 1 1 ALC 48,69 CMAT Local 22 R(8) 8 2 1 ALC 47,70,90 DEVCSF Subr 131 131 DMRRRR Subr 130 130,146 DPSIEQ Local 74 R(4) 4 scalar 74,77 DTRANS Subr 129 129 EID Local 36 I(4) 4 1 100 228,238,239,244,245,249,250,255,25 6,268,269,271 FEV_ROUTINE Local 34 R(8) 8 scalar 92,210 FLAG1 Local 234 R(4) 4 scalar 234,238,240 FLAG2 Local 235 R(4) 4 scalar 235,244,246,249,251 FLAG3 Local 236 R(4) 4 scalar 236,255,257 FLIGHT2 Subr 123 123 FLIGHTS_ALLNORMS Prog 17 GEORAD Local 41 R(4) 4 scalar 41,91 HY Local 26 R(8) 8 2 1 ALC 51,118,123,129,130 HYT Local 27 R(8) 8 2 1 ALC 52,119,129,130 I Local 96 I(4) 4 scalar 96,98,99,100,101,102,103,104,105,1 06,237,238,239,244,245,249,250,255 ,256,266,268,269,271,273 IDIM Local 40 I(4) 4 scalar 40,53,91,96,123 IFLIGHT Local 116 I(4) 4 scalar 116,123,148,154,172,180,190,198,20 2,206,209,212,213,216,217,223,224 IMSLP Local 167 I(4) 4 scalar 167,192 Page 7 Source Listing FLIGHTS_ALLNORMS 2022-05-02 15:46 Symbol Table flights_allnorms.f Name Object Declared Type Bytes Dimen Elements Attributes References INT Func 263 scalar 263,264 IPHASE Local 42 I(4) 4 scalar 42,243 IPRCP Local 166 I(4) 4 scalar 166,192 IRA Local 36 I(4) 4 1 100 228 IT200 Local 165 I(4) 4 scalar 165,185 IT500 Local 164 I(4) 4 scalar 164,184 IT850 Local 163 I(4) 4 scalar 163,184 IU200 Local 161 I(4) 4 scalar 161,175,183 IU500 Local 159 I(4) 4 scalar 159,174,182 IU850 Local 157 I(4) 4 scalar 157,173,181,191 IV200 Local 162 I(4) 4 scalar 162,175,183 IV500 Local 160 I(4) 4 scalar 160,174,182 IV850 Local 158 I(4) 4 scalar 158,173,181,191 IVCOUNT Local 155 I(4) 4 scalar 155,157,158,159,160,161,162,163,16 4,165,166,167,168 IVM Local 90 I(4) 4 scalar 90,92,155,198,202,206 IVNORM Local 42 I(4) 4 scalar 42,94,171,179,189,197,201,205 IVR Local 60 I(4) 4 scalar 60,90 IZ500 Local 168 I(4) 4 scalar 168 J Local 75 I(4) 4 scalar 75,76,77,78,97,98,99,100,101,102,1 03,104,105,106 JDIM Local 40 I(4) 4 scalar 40,53,91,97,124 JEQ Local 76 I(4) 4 scalar 76,77 JSTR Local 40 I(4) 4 scalar 40,44,45,75,91 JVR Local 61 I(4) 4 scalar 61,90 K Local 156 I(4) 4 scalar 156,173,174,175,181,182,183,184,18 5,191,192 KNUM Local 273 I(4) 4 scalar 273,274 KNUM0 Local 81 I(4) 4 scalar 81,82,83,84,85,86,87,88 KNUM1 Local 63 I(4) 4 scalar 63,68 KNUM2 Local 64 I(4) 4 scalar 64,69 KNUM3 Local 65 I(4) 4 scalar 65,70 KNUM5 Local 66 I(4) 4 scalar 66,71 LTCODE Local 41 I(4) 4 scalar 41,63,64,65 MEM Local 40 I(4) 4 scalar 40,46,47,48,49,50,51,52,54,55,56,5 7,58,91,94,123,129,130,131,133,137 ,138,146,147,156,198,202,206 MSLPSD Local 38 R(4) 4 scalar 86,93,192 N2 Local 41 I(4) 4 scalar 41 ND Local 40 I(4) 4 scalar 40,123 ND12 Local 41 I(4) 4 scalar 41 ND9 Local 40 I(4) 4 scalar 40,51,52,123,129,130 NE1 Local 41 I(4) 4 scalar 41,46,91 NE9 Local 41 I(4) 4 scalar 41,50,123 NFLIGHTS Local 41 I(4) 4 scalar 41,116,223,228,237 NINT Func 60 scalar 60,61 NM1 Local 133 I(4) 4 scalar 133,134,137,139,140 NM2 Local 138 I(4) 4 scalar 138,139,140 NV Local 41 I(4) 4 scalar 41,124 NVR Local 41 I(4) 4 scalar 41 NVR12 Local 41 I(4) 4 scalar 41,49,57,91,94,146,147 OBS Local 28 R(8) 8 3 1 ALC 53,98,99,100,101,102,103,104,105,1 06,123 PI Local 73 R(4) 4 scalar 73,78 PRCPSD Local 38 R(4) 4 scalar 85,93,192 Page 8 Source Listing FLIGHTS_ALLNORMS 2022-05-02 15:46 Symbol Table flights_allnorms.f Name Object Declared Type Bytes Dimen Elements Attributes References RADVR Local 40 R(4) 4 scalar 40,91 SIG Local 34 R(8) 8 1 100 148,154,172,180,190,198,202,206,20 9,212,213,216,217,224,228 SIGMAX Local 109 R(4) 4 scalar 109,212,213,263 SIGMIN Local 110 R(4) 4 scalar 110,216,217,264 SIGOLD Local 35 R(8) 8 1 100 228 SIGTOP Local 20 R(4) 4 1 1 ALC 45,78,90 SISTR Local 19 R(4) 4 1 1 ALC 44,77,78,90 SQRT Func 140 scalar 140 SVPVR Local 32 R(8) 8 2 1 ALC 57,147,173,174,175,181,182,183,184 ,185,191,192 TFL Local 37 CHAR 2 scalar 267,269,271,274 UV850SD Local 38 R(4) 4 scalar 84,93,191 UVSD Local 38 R(4) 4 scalar 83,93 VRLAT Local 40 R(4) 4 scalar 40,61 VRLON Local 40 R(4) 4 scalar 40,60 WORK3 Local 33 R(8) 8 2 1 ALC 58,130,131 XVVEST Local 21 R(8) 8 2 1 ALC 46,71,90 YTVEC Local 25 R(8) 8 2 1 ALC 50,68,123 YVVEST Local 24 R(8) 8 2 1 ALC 49,90,92,146 Z500SD Local 38 R(4) 4 scalar 87,93 Page 9 Source Listing FLIGHTS_ALLNORMS 2022-05-02 15:46 flights_allnorms.f 282 283 ************************************************************************ 284 ************************************************************************ 285 286 287 ************************************************************************ 288 * ----------------------------------------------------------------- 289 * Noting that keeping a bunch of extra variables around, the time 290 * has come to make matrix transposes a calculates, rather than a 291 * stored quantity. 292 * ----------------------------------------------------------------- 293 ************************************************************************ 294 SUBROUTINE DTRANS(M,N,A,B) 295 ************************************************************************ 296 297 double precision A(m,n) 298 double precision B(n,m) 299 do 1111 i=1,n 300 do 1111 j=1,m 301 1111 B(i,j)=A(j,i) 302 303 return 304 end ENTRY POINTS Name dtrans_ Page 10 Source Listing DTRANS 2022-05-02 15:46 Symbol Table flights_allnorms.f SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 1111 Label 301 299,300 A Dummy 294 R(8) 8 2 0 ARG,INOUT 301 B Dummy 294 R(8) 8 2 0 ARG,INOUT 301 DTRANS Subr 294 I Local 299 I(4) 4 scalar 299,301 J Local 300 I(4) 4 scalar 300,301 M Dummy 294 I(4) 4 scalar ARG,INOUT 297,298,300 N Dummy 294 I(4) 4 scalar ARG,INOUT 297,298,299 Page 11 Source Listing DTRANS 2022-05-02 15:46 flights_allnorms.f 305 306 307 308 309 ************************************************************************ 310 * ----------------------------------------------------------------- 311 * Prescribed flight tracks -- set matrix HY 312 * ----------------------------------------------------------------- 313 subroutine flight2(iflight,obs,hy,ytvec, 314 &nd,ne9,nd9,mem,idim,jdim,nv) 315 * ----------------------------------------------------------------- 316 ************************************************************************ 317 318 real*4 obslon(nd),obslat(nd) 319 integer lat1(nd),lat2(nd),lon1(nd),lon2(nd) 320 double precision ytvec(ne9,mem),hmat(nd,4),hy(nd9,mem) 321 double precision obs(idim,jdim,nv) 322 character*12 droptrack 323 324 * ---------------------------------------------------------- 325 * Create the translation matrix from gridpoint space to 326 * observation space known as H. 327 * ---------------------------------------------------------- 328 droptrack='wsr_track.00' 329 * print*, 'iflight=',iflight 330 if (iflight .ge. 10) then 331 write(droptrack(11:12),'(i2)')iflight 332 else 333 * write(droptrack(11:11),'(i1)')0 334 write(droptrack(12:12),'(i1)')iflight 335 endif 336 * print*, 'droptrack=',droptrack 337 open(7,file=droptrack,form='formatted',status='unknown') 338 339 pi=acos(-1.0) 340 read(7,*) ndd 341 342 do 20 iobs=1,ndd 343 read(7,*) obslon(iobs), obslat(iobs) 344 C Y Song added if statement to allow positive longitude 345 if(obslon(iobs).le.0) then 346 obslon(iobs)=360.+obslon(iobs) 347 endif 348 xlo1=int(obslon(iobs)*144./360.)*2.5 349 xlo2=xlo1+2.5 350 yla1=int(obslat(iobs)*72./180.)*2.5 351 yla2=yla1+2.5 352 zquot=(xlo2-xlo1)*(yla2-yla1) 353 c print*, 'obslon(iobs), obslat(iobs), xlo1,xlo2,yla1,yla2' 354 c print*, obslon(iobs), obslat(iobs), xlo1,xlo2,yla1,yla2 355 hmat(iobs,1)= (obslon(iobs)-xlo2)*(obslat(iobs)-yla2)/zquot 356 hmat(iobs,2)=-(obslon(iobs)-xlo2)*(obslat(iobs)-yla1)/zquot 357 hmat(iobs,3)=-(obslon(iobs)-xlo1)*(obslat(iobs)-yla2)/zquot 358 hmat(iobs,4)= (obslon(iobs)-xlo1)*(obslat(iobs)-yla1)/zquot 359 c print*, 'h1=',hmat(iobs,1),' h2=',hmat(iobs,2) 360 c print*, 'h3=',hmat(iobs,3),' h4=',hmat(iobs,4) 361 lon1(iobs)=nint(xlo1/2.5)+1 Page 12 Source Listing FLIGHT2 2022-05-02 15:46 flights_allnorms.f 362 lon2(iobs)=nint(xlo2/2.5)+1 363 lat1(iobs)=37-nint(yla1/2.5) 364 lat2(iobs)=37-nint(yla2/2.5) 365 c print*, 'lon1=',lon1(iobs),' lon2=',lon2(iobs) 366 c print*, 'lat1=',lat1(iobs),' lat2=',lat2(iobs) 367 20 continue 368 close(7) 369 370 * -------------------------------------------------------- 371 * Generate the nice small matrix H*Y, using our brains! 372 * -------------------------------------------------------- 373 do 33 nm=1,mem 374 do 33 iobs=1,ndd 375 do 33 iv=1,nv 376 mm1=9*(lat1(iobs)-1)*idim + 9*(lon1(iobs)-1) + iv 377 mm2=9*(lat2(iobs)-1)*idim + 9*(lon1(iobs)-1) + iv 378 mm3=9*(lat1(iobs)-1)*idim + 9*(lon2(iobs)-1) + iv 379 mm4=9*(lat2(iobs)-1)*idim + 9*(lon2(iobs)-1) + iv 380 io=9*(iobs-1) + iv 381 hy(io,nm)=hmat(iobs,1)*ytvec(mm1,nm)+hmat(iobs,2)*ytvec(mm2,nm)+ 382 & hmat(iobs,3)*ytvec(mm3,nm)+hmat(iobs,4)*ytvec(mm4,nm) 383 hy(io,nm)=hy(io,nm)/sqrt(obs(1,1,iv)) 384 33 continue 385 386 return 387 end ENTRY POINTS Name flight2_ Page 13 Source Listing FLIGHT2 2022-05-02 15:46 Symbol Table flights_allnorms.f SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 20 Label 367 342 33 Label 384 373,374,375 ACOS Func 339 scalar 339 DROPTRACK Local 322 CHAR 12 scalar 328,331,334,337 FLIGHT2 Subr 313 HMAT Local 320 R(8) 8 2 0 355,356,357,358,381,382 HY Dummy 313 R(8) 8 2 0 ARG,INOUT 381,383 IDIM Dummy 314 I(4) 4 scalar ARG,INOUT 321,376,377,378,379 IFLIGHT Dummy 313 I(4) 4 scalar ARG,INOUT 330,331,334 INT Func 348 scalar 348,350 IO Local 380 I(4) 4 scalar 380,381,383 IOBS Local 342 I(4) 4 scalar 342,343,345,346,348,350,355,356,35 7,358,361,362,363,364,374,376,377, 378,379,380,381,382 IV Local 375 I(4) 4 scalar 375,376,377,378,379,380,383 JDIM Dummy 314 I(4) 4 scalar ARG,INOUT 321 LAT1 Local 319 I(4) 4 1 0 363,376,378 LAT2 Local 319 I(4) 4 1 0 364,377,379 LON1 Local 319 I(4) 4 1 0 361,376,377 LON2 Local 319 I(4) 4 1 0 362,378,379 MEM Dummy 314 I(4) 4 scalar ARG,INOUT 320,373 MM1 Local 376 I(4) 4 scalar 376,381 MM2 Local 377 I(4) 4 scalar 377,381 MM3 Local 378 I(4) 4 scalar 378,382 MM4 Local 379 I(4) 4 scalar 379,382 ND Dummy 314 I(4) 4 scalar ARG,INOUT 318,319,320 ND9 Dummy 314 I(4) 4 scalar ARG,INOUT 320 NDD Local 340 I(4) 4 scalar 340,342,374 NE9 Dummy 314 I(4) 4 scalar ARG,INOUT 320 NINT Func 361 scalar 361,362,363,364 NM Local 373 I(4) 4 scalar 373,381,382,383 NV Dummy 314 I(4) 4 scalar ARG,INOUT 321,375 OBS Dummy 313 R(8) 8 3 0 ARG,INOUT 383 OBSLAT Local 318 R(4) 4 1 0 343,350,355,356,357,358 OBSLON Local 318 R(4) 4 1 0 343,345,346,348,355,356,357,358 PI Local 339 R(4) 4 scalar 339 SQRT Func 383 scalar 383 XLO1 Local 348 R(4) 4 scalar 348,349,352,357,358,361 XLO2 Local 349 R(4) 4 scalar 349,352,355,356,362 YLA1 Local 350 R(4) 4 scalar 350,351,352,356,358,363 YLA2 Local 351 R(4) 4 scalar 351,352,355,357,364 YTVEC Dummy 313 R(8) 8 2 0 ARG,INOUT 381,382 ZQUOT Local 352 R(4) 4 scalar 352,355,356,357,358 Page 14 Source Listing FLIGHT2 2022-05-02 15:46 flights_allnorms.f 388 389 390 391 392 393 394 ************************************************************************ 395 * ----------------------------------------------------------------- 396 * Circular verification region, calculate Y matrix in VR 397 * ----------------------------------------------------------------- 398 subroutine calcvrci(ivr,jvr,sistr,sigtop,cmat,ivm,xvvest,yvvest, 399 &jstr,ne1,mem,nvr12,idim,jdim,radvr,georad) 400 * ----------------------------------------------------------------- 401 ************************************************************************ 402 403 dimension sistr(jstr),sigtop(jstr) 404 integer ivra(500),jvra(500) 405 double precision xvvest(ne1,mem),cmat(mem,mem),yvvest(nvr12,mem) 406 407 pi=acos(-1.0) 408 409 c compute coordinates of the verification area 410 c ============================================ 411 radjvr=float(jvr-1)*pi/72. 412 radivr=float(ivr-1)*2.0*pi/144. 413 zvr=cos(radjvr) 414 xvr=sin(radjvr)*cos(radivr) 415 yvr=sin(radjvr)*sin(radivr) 416 417 c ================================================================ 418 c define distance of grid point from center of verification region 419 c ================================================================ 420 ie2=0 421 c print*, 'j, sistr(j)*180./pi, sigtop(j)*180./pi' 422 do 70 j=1,jstr 423 c print*, j, sistr(j)*180./pi, sigtop(j)*180./pi 424 radj=sigtop(j) 425 z=cos(radj) 426 hr=sin(radj) 427 do 60 i=1,idim 428 radi=float(i-1)*2.0*pi/144. 429 x=hr*cos(radi) 430 y=hr*sin(radi) 431 costh=(xvr*x+yvr*y+zvr*z) 432 th=acos(costh) 433 dist=georad*th 434 if (dist.le.radvr) then 435 ie2=ie2+1 436 ivra(ie2)=i 437 jvra(ie2)=j 438 c print*, ie2,float(i-1)*2.5,sistr(j)*180./pi 439 end if 440 60 continue 441 70 continue 442 ivm=ie2 443 c print*, 'ivmemo=',ivm 444 Page 15 Source Listing CALCVRCI 2022-05-02 15:46 flights_allnorms.f 445 do 61 ivcount=1,ivm 446 mmm=12*((jvra(ivcount)-1)*idim) + 12*(ivra(ivcount)-1) 447 do 61 k=1,mem 448 do 61 ivar=1,12 449 num=mmm+ivar 450 iv=(ivcount-1)*12+ivar 451 yvvest(iv,k)=0.0 452 do 62 nm2=1,mem 453 yvvest(iv,k)=yvvest(iv,k)+xvvest(num,nm2)*cmat(nm2,k) 454 62 continue 455 61 continue 456 457 return 458 end ENTRY POINTS Name calcvrci_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 60 Label 440 427 61 Label 455 445,447,448 62 Label 454 452 70 Label 441 422 ACOS Func 407 scalar 407,432 CALCVRCI Subr 398 CMAT Dummy 398 R(8) 8 2 0 ARG,INOUT 453 COS Func 413 scalar 413,414,425,429 COSTH Local 431 R(4) 4 scalar 431,432 DIST Local 433 R(4) 4 scalar 433,434 FLOAT Func 411 scalar 411,412,428 GEORAD Dummy 399 R(4) 4 scalar ARG,INOUT 433 HR Local 426 R(4) 4 scalar 426,429,430 I Local 427 I(4) 4 scalar 427,428,436 IDIM Dummy 399 I(4) 4 scalar ARG,INOUT 427,446 IE2 Local 420 I(4) 4 scalar 420,435,436,437,442 IV Local 450 I(4) 4 scalar 450,451,453 IVAR Local 448 I(4) 4 scalar 448,449,450 IVCOUNT Local 445 I(4) 4 scalar 445,446,450 IVM Dummy 398 I(4) 4 scalar ARG,INOUT 442,445 IVR Dummy 398 I(4) 4 scalar ARG,INOUT 412 IVRA Local 404 I(4) 4 1 500 436,446 J Local 422 I(4) 4 scalar 422,424,437 JDIM Dummy 399 I(4) 4 scalar ARG,INOUT JSTR Dummy 399 I(4) 4 scalar ARG,INOUT 403,422 JVR Dummy 398 I(4) 4 scalar ARG,INOUT 411 JVRA Local 404 I(4) 4 1 500 437,446 K Local 447 I(4) 4 scalar 447,451,453 MEM Dummy 399 I(4) 4 scalar ARG,INOUT 405,447,452 MMM Local 446 I(4) 4 scalar 446,449 Page 16 Source Listing CALCVRCI 2022-05-02 15:46 Symbol Table flights_allnorms.f Name Object Declared Type Bytes Dimen Elements Attributes References NE1 Dummy 399 I(4) 4 scalar ARG,INOUT 405 NM2 Local 452 I(4) 4 scalar 452,453 NUM Local 449 I(4) 4 scalar 449,453 NVR12 Dummy 399 I(4) 4 scalar ARG,INOUT 405 PI Local 407 R(4) 4 scalar 407,411,412,428 RADI Local 428 R(4) 4 scalar 428,429,430 RADIVR Local 412 R(4) 4 scalar 412,414,415 RADJ Local 424 R(4) 4 scalar 424,425,426 RADJVR Local 411 R(4) 4 scalar 411,413,414,415 RADVR Dummy 399 R(4) 4 scalar ARG,INOUT 434 SIGTOP Dummy 398 R(4) 4 1 0 ARG,INOUT 424 SIN Func 414 scalar 414,415,426,430 SISTR Dummy 398 R(4) 4 1 0 ARG,INOUT TH Local 432 R(4) 4 scalar 432,433 X Local 429 R(4) 4 scalar 429,431 XVR Local 414 R(4) 4 scalar 414,431 XVVEST Dummy 398 R(8) 8 2 0 ARG,INOUT 453 Y Local 430 R(4) 4 scalar 430,431 YVR Local 415 R(4) 4 scalar 415,431 YVVEST Dummy 398 R(8) 8 2 0 ARG,INOUT 451,453 Z Local 425 R(4) 4 scalar 425,431 ZVR Local 413 R(4) 4 scalar 413,431 Page 17 Source Listing CALCVRCI 2022-05-02 15:46 flights_allnorms.f 459 460 461 ************************************************************************ 462 * ----------------------------------------------------------------- 463 * Calculate forecast error variance for routine obs. network 464 * ----------------------------------------------------------------- 465 subroutine calc_fev_routine(ivm,yvvest,fev_routine, 466 & uvsd,uv850sd,prcpsd,mslpsd,z500sd, 467 & nvr12,mem,ivnorm) 468 * ----------------------------------------------------------------- 469 ************************************************************************ 470 471 integer ivm 472 double precision yvvest(nvr12,mem),fev_routine 473 474 real uvsd,uv850sd,prcpsd,mslpsd,z500sd 475 476 fev_routine=0.0 477 do 161 ivcount=1,ivm 478 do 161 k=1,mem 479 iu850=(ivcount-1)*12 + 1 480 iv850=(ivcount-1)*12 + 2 481 iu500=(ivcount-1)*12 + 3 482 iv500=(ivcount-1)*12 + 4 483 iu200=(ivcount-1)*12 + 5 484 iv200=(ivcount-1)*12 + 6 485 it850=(ivcount-1)*12 + 7 486 it500=(ivcount-1)*12 + 8 487 it200=(ivcount-1)*12 + 9 488 iprcp=(ivcount-1)*12 + 10 489 imslp=(ivcount-1)*12 + 11 490 iz500=(ivcount-1)*12 + 12 491 492 * ****** (u,v) verifying norm ****** 493 if (ivnorm.eq.1) then 494 fev_routine = fev_routine + 495 & yvvest(iu850,k)**2 + yvvest(iv850,k)**2 + 496 & yvvest(iu500,k)**2 + yvvest(iv500,k)**2 + 497 & yvvest(iu200,k)**2 + yvvest(iv200,k)**2 498 endif 499 500 * ****** Total energy verifying norm ****** 501 if (ivnorm.eq.2) then 502 fev_routine = fev_routine + 503 & 0.5*( yvvest(iu850,k)**2 + yvvest(iv850,k)**2 + 504 & yvvest(iu500,k)**2 + yvvest(iv500,k)**2 + 505 & yvvest(iu200,k)**2 + yvvest(iv200,k)**2 ) + 506 & 3.346*( yvvest(it850,k)**2 + yvvest(it500,k)**2 + 507 & yvvest(it200,k)**2 ) 508 endif 509 510 * ****** (u850,v850,precip,mslp) verifying norm ****** 511 if (ivnorm.eq.3) then 512 fev_routine = fev_routine + 513 & (yvvest(iu850,k)**2 + yvvest(iv850,k)**2)/uv850sd**2 + 514 & (yvvest(iprcp,k)/prcpsd)**2 + (yvvest(imslp,k)/mslpsd)**2 515 endif Page 18 Source Listing CALC_FEV_ROUTINE 2022-05-02 15:46 flights_allnorms.f 516 517 161 continue 518 519 if (ivnorm.eq.1) then 520 fev_routine = fev_routine/(mem*ivm*3.0) 521 endif 522 523 if (ivnorm.eq.2) then 524 fev_routine = fev_routine/(mem*ivm*3.0) 525 endif 526 527 if (ivnorm.eq.3) then 528 fev_routine = fev_routine/(mem*ivm) 529 endif 530 531 return 532 end ENTRY POINTS Name calc_fev_routine_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 161 Label 517 477,478 CALC_FEV_ROUTINE Subr 465 FEV_ROUTINE Dummy 465 R(8) 8 scalar ARG,INOUT 476,494,502,512,520,524,528 IMSLP Local 489 I(4) 4 scalar 489,514 IPRCP Local 488 I(4) 4 scalar 488,514 IT200 Local 487 I(4) 4 scalar 487,507 IT500 Local 486 I(4) 4 scalar 486,506 IT850 Local 485 I(4) 4 scalar 485,506 IU200 Local 483 I(4) 4 scalar 483,497,505 IU500 Local 481 I(4) 4 scalar 481,496,504 IU850 Local 479 I(4) 4 scalar 479,495,503,513 IV200 Local 484 I(4) 4 scalar 484,497,505 IV500 Local 482 I(4) 4 scalar 482,496,504 IV850 Local 480 I(4) 4 scalar 480,495,503,513 IVCOUNT Local 477 I(4) 4 scalar 477,479,480,481,482,483,484,485,48 6,487,488,489,490 IVM Dummy 465 I(4) 4 scalar ARG,INOUT 477,520,524,528 IVNORM Dummy 467 I(4) 4 scalar ARG,INOUT 493,501,511,519,523,527 IZ500 Local 490 I(4) 4 scalar 490 K Local 478 I(4) 4 scalar 478,495,496,497,503,504,505,506,50 7,513,514 MEM Dummy 467 I(4) 4 scalar ARG,INOUT 472,478,520,524,528 MSLPSD Dummy 466 R(4) 4 scalar ARG,INOUT 514 NVR12 Dummy 467 I(4) 4 scalar ARG,INOUT 472 PRCPSD Dummy 466 R(4) 4 scalar ARG,INOUT 514 UV850SD Dummy 466 R(4) 4 scalar ARG,INOUT 513 UVSD Dummy 466 R(4) 4 scalar ARG,INOUT Page 19 Source Listing CALC_FEV_ROUTINE 2022-05-02 15:46 Symbol Table flights_allnorms.f Name Object Declared Type Bytes Dimen Elements Attributes References YVVEST Dummy 465 R(8) 8 2 0 ARG,INOUT 495,496,497,503,504,505,506,507,51 3,514 Z500SD Dummy 466 R(4) 4 scalar ARG,INOUT Page 20 Source Listing CALC_FEV_ROUTINE 2022-05-02 15:46 flights_allnorms.f 533 534 535 c ********************************************************************** 536 subroutine bubble_sort(signal,sigold,iid,ira,nflights) 537 c ********************************************************************** 538 539 integer iid(500), iid2(500), ira(500), nflights 540 double precision signal(500), sigold(500), sigve2(500), maxsig 541 542 print *,nflights 543 do 10 k=1,nflights 544 iid(k)=k 545 iid2(k)=iid(k) 546 sigve2(k)=signal(k) 547 sigold(k)=signal(k) 548 10 continue 549 550 do 18 k=1,nflights 551 maxsig=0.0 552 do 19 kk=k,nflights 553 if (signal(kk).gt.maxsig) then 554 maxsig=signal(kk) 555 kkk=kk 556 endif 557 19 continue 558 print *,kkk 559 do 20 kk=k+1,kkk 560 signal(kk)=sigve2(kk-1) 561 iid(kk)=iid2(kk-1) 562 20 continue 563 signal(k)=maxsig 564 iid(k)=iid2(kkk) 565 do 21 kk=1,nflights 566 sigve2(kk)=signal(kk) 567 iid2(kk)=iid(kk) 568 21 continue 569 18 continue 570 close(17) 571 572 do 30 ir=1,nd 573 ira(iid(ir))=ir 574 30 continue 575 return 576 end Page 21 Source Listing BUBBLE_SORT 2022-05-02 15:46 Entry Points flights_allnorms.f ENTRY POINTS Name bubble_sort_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 10 Label 548 543 18 Label 569 550 19 Label 557 552 20 Label 562 559 21 Label 568 565 30 Label 574 572 BUBBLE_SORT Subr 536 IID Dummy 536 I(4) 4 1 500 ARG,INOUT 544,545,561,564,567,573 IID2 Local 539 I(4) 4 1 500 545,561,564,567 IR Local 572 I(4) 4 scalar 572,573 IRA Dummy 536 I(4) 4 1 500 ARG,INOUT 573 K Local 543 I(4) 4 scalar 543,544,545,546,547,550,552,559,56 3,564 KK Local 552 I(4) 4 scalar 552,553,554,555,559,560,561,565,56 6,567 KKK Local 555 I(4) 4 scalar 555,558,559,564 MAXSIG Local 540 R(8) 8 scalar 551,553,554,563 ND Local 572 I(4) 4 scalar 572 NFLIGHTS Dummy 536 I(4) 4 scalar ARG,INOUT 542,543,550,552,565 SIGNAL Dummy 536 R(8) 8 1 500 ARG,INOUT 546,547,553,554,560,563,566 SIGOLD Dummy 536 R(8) 8 1 500 ARG,INOUT 547 SIGVE2 Local 540 R(8) 8 1 500 546,560,566 Page 22 Source Listing BUBBLE_SORT 2022-05-02 15:46 Subprograms/Common Blocks flights_allnorms.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References BUBBLE_SORT Subr 536 CALCVRCI Subr 398 CALC_FEV_ROUTINE Subr 465 DTRANS Subr 294 FLIGHT2 Subr 313 FLIGHTS_ALLNORMS Prog 17 COMPILER OPTIONS BEING USED -align noall -align nonone -align nocommons -align nodcommons -align noqcommons -align nozcommons -align records -align sequence -align norec1byte -align norec2byte -align norec4byte -align norec8byte -align norec16byte -align norec32byte -align norec64byte -align noarray8byte -align noarray16byte -align noarray32byte -align noarray64byte -align noarray128byte -align noarray256byte -altparam -assume accuracy_sensitive -assume nobscc -assume nobuffered_io -assume nobuffered_stdout -assume byterecl -assume nocontiguous_assumed_shape -assume nocontiguous_pointer -assume nocc_omp -assume nocstring -assume nodummy_aliases -assume nofpe_summary -assume noieee_fpe_flags -assume nominus0 -assume noold_boz -assume old_complex_align -assume old_unit_star -assume old_inquire_recl -assume old_ldout_format -assume old_ldout_zero -assume noold_logical_assign -assume noold_logical_ldio -assume old_maxminloc -assume old_xor -assume noprotect_allocates -assume protect_constants -assume noprotect_parens -assume split_common -assume source_include -assume nostd_intent_in -assume std_minus0_rounding -assume nostd_mod_proc_name -assume std_value -assume realloc_lhs -assume underscore -assume no2underscores -assume norecursion no -auto -auto_scalar no -bintext -ccdefault default -check noarg_temp_created -check noassume -check nobounds -check nocontiguous -check noformat -check nooutput_conversion -check nooverflow -check nopointers -check noshape -check nostack -check nouninitialized -check noudio_iostat -coarray-num-procs 0 no -coarray-config-file -convert big_endian -cross_reference -D __INTEL_COMPILER=1910 -D __INTEL_COMPILER_UPDATE=3 -D __unix__ -D __unix -D __linux__ -D __linux Page 23 Source Listing BUBBLE_SORT 2022-05-02 15:46 flights_allnorms.f -D __gnu_linux__ -D unix -D linux -D __ELF__ -D __x86_64 -D __x86_64__ -D __amd64 -D __amd64__ -D __INTEL_COMPILER_BUILD_DATE=20200925 -D __INTEL_OFFLOAD -D __MMX__ -D __SSE__ -D __SSE_MATH__ -D __SSE2__ -D __SSE2_MATH__ -D __SSE3__ -D __SSSE3__ -D __SSE4_1__ -D __SSE4_2__ -D __POPCNT__ -D __PCLMUL__ -D __AES__ -D __AVX__ -D __F16C__ -D __AVX_I__ -D __RDRND__ -D __FMA__ -D __FP_FAST_FMA -D __FP_FAST_FMAF -D __BMI__ -D __LZCNT__ -D __AVX2__ -D __haswell -D __haswell__ -D __tune_haswell__ -D __core_avx2 -D __core_avx2__ -D __tune_core_avx2__ -double_size 64 no -d_lines no -Qdyncom -error_limit 30 no -f66 no -f77rtl no -fast -fpscomp nofilesfromcmd -fpscomp nogeneral -fpscomp noioformat -fpscomp noldio_spacing -fpscomp nologicals -fixed no -fpconstant -fpe3 -fprm nearest -ftz -fp_model precise -fp_model nofast -fp_model nostrict -fp_model nosource -fp_model nodouble -fp_model noextended -fp_model novery_fast -fp_model noexcept -fp_model nono_except -fp_modbits nofp_contract -fp_modbits nono_fp_contract -fp_modbits nofenv_access -fp_modbits nono_fenv_access -fp_modbits nocx_limited_range -fp_modbits nono_cx_limited_range -fp_modbits noprec_div -fp_modbits no_prec_div -fp_modbits noprec_sqrt -fp_modbits no_prec_sqrt -fp_modbits ftz -fp_modbits nono_ftz -fp_modbits nointrin_limited_range -fp_modbits nono_intrin_limited_range -fp_modbits notrunc_compares -fp_modbits nono_trunc_compares -fp_modbits noieee_nan_compares -fp_modbits nono_ieee_nan_compares -fp_modbits nohonor_f32_conversion -fp_modbits nono_honor_f32_conversion -fp_modbits nohonor_f64_conversion -fp_modbits nono_honor_f64_conversion -fp_modbits nono_x87_copy -fp_modbits nono_no_x87_copy -fp_modbits noexception_semantics -fp_modbits nono_exception_semantics -fp_modbits noprecise_libm_functions -fp_modbits no_precise_libm_functions -heap_arrays 0 no -threadprivate_compat -g0 -iface nomixed_str_len_arg -iface nono_mixed_str_len_arg -init noarrays -init nohuge -init noinfinity -init nominus_huge -init nominus_infinity -init nominus_tiny -init nonan -init nosnan -init notiny -init nozero no -intconstant -integer_size 32 no -mixed_str_len_arg no -module -names lowercase no -noinclude no -o Page 24 Source Listing BUBBLE_SORT 2022-05-02 15:46 flights_allnorms.f -offload-build=host -openmp-simd -O3 no -pad_source -real_size 32 no -recursive -reentrancy threaded -vec=simd -show nofullpath -show noinclude -show map -show options no -syntax_only no -threadcom no -U no -vms -w noall -w nonone -w alignments -w nodeclarations -w noexternals -w general -w noignore_bounds -w noignore_loc -w nointerfaces -w noshape -w notruncated_source -w uncalled -w uninitialized -w nounused -w usage no -wrap-margins -includepath : /pe/intel/compilers_and_libraries_2020.4.304/linux/ipp/include/,/pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/intel64/, /usr/include/,.FOR,./.FOR,/opt/cray/pe/mpich/8.1.9/ofi/intel/19.0/include/.FOR,/pe/intel/compilers_and_libraries_2020.4.304/linux/mkl/include/intel64/lp64/.FOR, /pe/intel/compilers_and_libraries_2020.4.304/linux/mkl/include/.FOR,/pe/intel/compilers_and_libraries_2020.4.304/linux/ipp/include/.FOR, /pe/intel/compilers_and_libraries_2020.4.304/linux/mkl/include/.FOR,/pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/include/.FOR, /pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/stdlib/.FOR,/pe/intel/compilers_and_libraries_2020.4.304/linux/tbb/include/.FOR, /pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/intel64/.FOR,/pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/icc/.FOR, /pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/.FOR,/usr/lib64/gcc/x86_64-suse-linux/7/include/.FOR, /usr/lib64/gcc/x86_64-suse-linux/7/include-fixed/.FOR,/usr/include/.FOR,/usr/include/.FOR,/usr/include/.FOR -list filename : flights_allnorms.lst no -o COMPILER: Intel(R) Fortran 19.1-1655