Page 1 Source Listing NGMFLD2 2014-12-17 20:48 NGMFLD2.f 1 SUBROUTINE NGMFLD2(RH4710,RH4796,RH1847,RH8498,QM8510) 2 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 3 C . . . 4 C SUBPROGRAM: NGMFLD2 COMPUTES LAYER MEAN NGM FIELDS 5 C PRGRMMR: TREADON ORG: W/NP2 DATE: 92-12-22 6 C 7 C ABSTRACT: 8 C THIS ROUTINE COMPUTES A HANDFUL OF NGM LAYER MEAN 9 C FIELDS. THIS IS DONE TO PROVIDE A FULLY COMPLETE 10 C ETA NGM LOOK-ALIKE OUTPUT FILE. THE SIGMA (LAYER) 11 C FIELDS COMPUTED BY THIS ROUTINE ARE TABULATED BELOW. 12 C 13 C SIGMA (LAYER) FIELD(S) 14 C --------------- -------------- 15 C 0.47191-1.00000 RH 16 C 0.47171-0.96470 RH 17 C 0.18019-0.47191 RH 18 C 0.84368-0.98230 RH 19 C 0.85000-1.00000 MCONV 20 C WHERE 21 C RH = RELATIVE HUMIDITY 22 C MCONV = MOISTURE CONVERGENCE 23 C 24 C LAYER MEANS ARE A SUMMATION OVER ETA LAYERS MAPPING INTO 25 C THE PRESSURE RANGE CORRESPONDING TO THE SIGMA RANGE ABOVE. 26 C THE CALCULATION OF THESE BOUNDING PRESSURES IS DONE AT 27 C EACH HORIZONTAL GRID POINT BASED ON THE SURFACE PRESSURE. 28 C EACH TERM IN THE SUMMATION IS WEIGHTED BY THE THICKNESS OF 29 C THE ETA LAYER. THE FINAL LAYER MEAN IS THIS SUM NORMALIZED 30 C BY THE TOTAL DEPTH OF THE LAYER. 31 32 C 33 C 34 C PROGRAM HISTORY LOG: 35 C 92-12-22 RUSS TREADON 36 C 93-07-27 RUSS TREADON - MODIFIED SUMMATION LIMITS FROM 37 C 0.66*PSFC TO 0.75*PSFC AND 0.33*PSFC 38 C TO 0.50*PSFC, WHERE PSFC IS THE 39 C SURFACES PRESSURE. THE REASON FOR 40 C THIS CHANGE WAS RECOGNITION THAT IN 41 C THE LFM 0.33 AND 0.66 WERE MEASURED 42 C FROM THE SURFACE TO THE TROPOPAUSE, 43 C NOT THE TOP OF THE MODEL. 44 C 93-09-13 RUSS TREADON - RH CALCULATIONS WERE MADE INTERNAL 45 C TO THE ROUTINE. 46 C 98-06-16 T BLACK - CONVERSION FROM 1-D TO 2-D 47 C 98-08-18 MIKE BALDWIN - COMPUTE RH OVER ICE 48 C 98-12-22 MIKE BALDWIN - BACK OUT RH OVER ICE 49 C 00-01-04 JIM TUCCILLO - MPI VERSION 50 C 51 C 52 C USAGE: CALL NGMFLD2(RH4710,RH4796,RH1847,RH8498,QM8510) 53 C INPUT ARGUMENT LIST: 54 C NONE 55 C 56 C OUTPUT ARGUMENT LIST: 57 C RH4710 - SIGMA LAYER 0.47-1.00 MEAN RELATIVE HUMIDITY. Page 2 Source Listing NGMFLD2 2014-12-17 20:48 NGMFLD2.f 58 C RH4796 - SIGMA LAYER 0.47-0.96 MEAN RELATIVE HUMIDITY. 59 C RH1847 - SIGMA LAYER 0.18-0.47 MEAN RELATIVE HUMIDITY. 60 C RH8498 - SIGMA LAYER 0.84-0.98 MEAN RELATIVE HUMIDITY. 61 C QM8510 - SIGMA LAYER 0.85-1.00 MEAN MOISTURE CONVERGENCE. 62 C 63 C OUTPUT FILES: 64 C NONE 65 C 66 C LIBRARY: 67 C COMMON - VRBLS 68 C MASKS 69 C EXTRA 70 C OPTIONS 71 C LOOPS 72 C MAPOT 73 C DYNAMD 74 C INDX 75 C 76 C ATTRIBUTES: 77 C LANGUAGE: FORTRAN 78 C MACHINE : CRAY C-90 79 C$$$ 80 C 81 C 82 C INCLUDE PARAMETERS 83 INCLUDE "parmeta" 101 INCLUDE "params" 102 C 144 PARAMETER (SIG100=1.00000, SIG98=0.98230, SIG96=0.96470) 145 PARAMETER (SIG89 =0.89671, SIG85=0.85000, SIG84=0.84368) 146 PARAMETER (SIG78 =0.78483, SIG47=0.47191, SIG18=0.18018) 147 PARAMETER (SMALL = 1.E-6) 148 C 149 C DECLARE VARIABLES. 150 LOGICAL GOT8510,GOT4710,GOT4796,GOT1847,GOT8498 151 LOGICAL OLDRD, STDRD 152 REAL QM8510(IM,JM),RH4710(IM,JM),RH8498(IM,JM) 153 REAL RH4796(IM,JM),RH1847(IM,JM),IWL(IM,JM),IWM1 154 REAL Z8510(IM,JM),Z4710(IM,JM),Z8498(IM,JM) 155 REAL Z4796(IM,JM),Z1847(IM,JM) 156 C 157 C INCLUDE COMMON BLOCKS. 158 INCLUDE "VRBLS.comm" 165 INCLUDE "CLDWTR.comm" 171 INCLUDE "MASKS.comm" 179 INCLUDE "EXTRA.comm" 191 INCLUDE "OPTIONS.comm" 195 INCLUDE "LOOPS.comm" 203 INCLUDE "MAPOT.comm" 216 INCLUDE "DYNAMD.comm" 228 INCLUDE "INDX.comm" 234 INCLUDE "CTLBLK.comm" 235 C 236 C******************************************************************** 237 C START NGMFLD HERE. 238 C 253 CLIMIT =1.0E-20 Page 3 Source Listing NGMFLD2 2014-12-17 20:48 NGMFLD2.f 254 C 255 C INITIALIZE ARRAYS. 256 !$omp parallel do 257 DO J=JSTA,JEND 258 DO I=1,IM 259 QM8510(I,J) = D00 260 RH4710(I,J) = D00 261 RH8498(I,J) = D00 262 RH4796(I,J) = D00 263 RH1847(I,J) = D00 264 Z8510(I,J) = D00 265 Z8498(I,J) = D00 266 Z4710(I,J) = D00 267 Z4796(I,J) = D00 268 Z1847(I,J) = D00 269 IWL (I,J) = D00 270 ENDDO 271 ENDDO 272 C 273 C LOOP OVER HORIZONTAL GRID. 274 C 275 !$omp parallel do 276 !$omp& private(ai,bi,dz,ie,iw,iwl,iwm1,p100,p18,p47,p84,p85, 277 !$omp& p96,p98,pm,qdiv,qi,qint,qk,qkhn,qkhs,qkm1,qm,qm8510, 278 !$omp& qmcvg,qs,qudx,qvdy,qw,r2dx,r2dy,rh,rh1847,rh4710, 279 !$omp& rh4796,rh8498,tm,tmt0,tmt15,z1847,z4710,z4796, 280 !$omp& z8498,z8510) 281 DO L=1,LM 282 CALL EXCH2(Q(1,1,L)) 283 DO J=JSTA_M2,JEND_M2 284 DO I=2,IM-1 285 C 286 C SET TARGET PRESSURES. 287 P100 = SIG100*(PD(I,J)+PT) 288 P98 = SIG98*(PD(I,J)+PT) 289 P96 = SIG96*(PD(I,J)+PT) 290 P85 = SIG85*(PD(I,J)+PT) 291 P84 = SIG84*(PD(I,J)+PT) 292 P47 = SIG47*(PD(I,J)+PT) 293 P18 = SIG18*(PD(I,J)+PT) 294 C 295 C 296 C COMPUTE LAYER MEAN FIELDS AT THE GIVEN K. 297 C 298 C COMPUTE P, Z, T, AND Q AT THE MIDPOINT OF THE CURRENT ETA LAYER. 299 ALPM = D50*(ALPINT(I,J,L)+ALPINT(I,J,L+1)) 300 DZ = ZINT(I,J,L)-ZINT(I,J,L+1) 301 PM = EXP(ALPM) 302 TM = T(I,J,L) 303 QM = Q(I,J,L) 304 QM = AMAX1(QM,H1M12) 305 C 306 IWM1=IWL(I,J) 307 IF(CWM(I,J,L).GT.CLIMIT) THEN 308 IF(TM.LT.258.15)THEN 309 IWL(I,J)=1. 310 ELSEIF(TM.GE.273.15)THEN Page 4 Source Listing NGMFLD2 2014-12-17 20:48 NGMFLD2.f 311 IWL(I,J)=0. 312 ELSE 313 IF(IWM1.EQ.1.0)IWL(I,J)=1. 314 ENDIF 315 ELSE 316 IWL(I,J)=0. 317 ENDIF 318 C 319 C 320 C COMPUTE RELATIVE HUMIDITY. 321 C 322 TMT0=TM-273.16 323 TMT15=AMIN1(TMT0,-15.) 324 AI=0.008855 325 BI=1. 326 IF(TMT0.LT.-20.)THEN 327 AI=0.007225 328 BI=0.9674 329 ENDIF 330 QW=PQ0/PM 331 1 *EXP(A2*(TM-A3)/(TM-A4)) 332 QI=QW*(BI+AI*AMIN1(TMT0,0.)) 333 QINT=QW*(1.-0.00032*TMT15*(TMT15+15.)) 334 IF(TMT0.LT.-15.)THEN 335 QS=QI 336 ELSEIF(TMT0.GE.0.)THEN 337 QS=QINT 338 ELSE 339 IF(IWL(I,J).GT.0.0) THEN 340 QS=QI 341 ELSE 342 QS=QINT 343 ENDIF 344 ENDIF 345 CMEB 12/22/98 SWITCH TO RH VS WATER NO MATTER WHAT 346 C DELETE THIS LINE TO SWITCH BACK TO RH VS ICE 347 QS=QW 348 CMEB 12/22/98 SWITCH TO RH VS WATER NO MATTER WHAT 349 350 C 351 RH = QM/QS 352 IF (RH.GT.H1) THEN 353 RH = H1 354 QM = RH*QS 355 ENDIF 356 IF (RH.LT.D01) THEN 357 RH = D01 358 QM = RH*QS 359 ENDIF 360 C 361 C COMPUTE MOISTURE CONVERGENCE. WE NEED Q AT V POINTS. 362 QMCVG = D00 363 IE=I+IHE(J) 364 IW=I+IHW(J) 365 QK=D25*(Q(I+1,J,L)+Q(I,J,L)+Q(IE,J+1,L)+Q(IE,J-1,L)) 366 QKM1=D25*(Q(I,J,L)+Q(I-1,J,L)+Q(IW,J+1,L)+Q(IW,J-1,L)) 367 QKHN=D25*(Q(IW,J+1,L)+Q(IE,J+1,L)+Q(I,J+2,L)+Q(I,J,L)) Page 5 Source Listing NGMFLD2 2014-12-17 20:48 NGMFLD2.f 368 QKHS=D25*(Q(IW,J-1,L)+Q(IE,J-1,L)+Q(I,J,L)+Q(I,J-2,L)) 369 R2DX =1./(2.*DX(I,J)) 370 R2DY =1./(2.*DY) 371 QUDX =(QK*U(IE,J,L)-QKM1*U(IW,J,L))*R2DX 372 QVDY =(QKHN*V(I,J+1,L)-QKHS*V(I,J-1,L))*R2DY 373 QDIV =QUDX+QVDY 374 QMCVG=-1.*QDIV*HBM2(I,J) 375 C 376 C SIGMA 0.85-1.00 MOISTURE CONVERGENCE. 377 IF ((PM.LE.P100).AND.(PM.GE.P85)) THEN 378 Z8510(I,J) = Z8510(I,J) + DZ 379 QM8510(I,J) = QM8510(I,J) + QMCVG*DZ 380 ENDIF 381 C 382 C SIGMA 0.47-1.00 RELATIVE HUMIDITY. 383 IF ((PM.LE.P100).AND.(PM.GE.P47)) THEN 384 Z4710(I,J) = Z4710(I,J) + DZ 385 RH4710(I,J) = RH4710(I,J) + RH*DZ 386 ENDIF 387 C 388 C SIGMA 0.84-0.98 RELATIVE HUMIDITY. 389 IF ((PM.LE.P98).AND.(PM.GE.P84)) THEN 390 Z8498(I,J) = Z8498(I,J) + DZ 391 RH8498(I,J) = RH8498(I,J) + RH*DZ 392 ENDIF 393 C 394 C SIGMA 0.47-0.96 RELATIVE HUMIDITY. 395 IF ((PM.LE.P96).AND.(PM.GE.P47)) THEN 396 Z4796(I,J) = Z4796(I,J) + DZ 397 RH4796(I,J) = RH4796(I,J) + RH*DZ 398 ENDIF 399 C 400 C SIGMA 0.18-0.47 RELATIVE HUMIDITY. 401 IF ((PM.LE.P47).AND.(PM.GE.P18)) THEN 402 Z1847(I,J) = Z1847(I,J) + DZ 403 RH1847(I,J) = RH1847(I,J) + RH*DZ 404 ENDIF 405 C 406 ENDDO 407 ENDDO 408 ENDDO 409 C 410 DO J=JSTA_M2,JEND_M2 411 DO I=2,IM-1 412 C NORMALIZE TO GET LAYER MEAN VALUES. 413 IF (Z8510(I,J).GT.0) THEN 414 QM8510(I,J) = QM8510(I,J)/Z8510(I,J) 415 ELSE 416 QM8510(I,J) = SPVAL 417 ENDIF 418 IF (ABS(QM8510(I,J)-SPVAL).LT.SMALL)QM8510(I,J)=H1M12 419 C 420 IF (Z4710(I,J).GT.0) THEN 421 RH4710(I,J) = RH4710(I,J)/Z4710(I,J) 422 ELSE 423 RH4710(I,J) = SPVAL 424 ENDIF Page 6 Source Listing NGMFLD2 2014-12-17 20:48 NGMFLD2.f 425 C 426 IF (Z8498(I,J).GT.0) THEN 427 RH8498(I,J) = RH8498(I,J)/Z8498(I,J) 428 ELSE 429 RH8498(I,J) = SPVAL 430 ENDIF 431 C 432 IF (Z4796(I,J).GT.0) THEN 433 RH4796(I,J) = RH4796(I,J)/Z4796(I,J) 434 ELSE 435 RH4796(I,J) = SPVAL 436 ENDIF 437 C 438 IF (Z1847(I,J).GT.0) THEN 439 RH1847(I,J) = RH1847(I,J)/Z1847(I,J) 440 ELSE 441 RH1847(I,J) = SPVAL 442 ENDIF 443 ENDDO 444 ENDDO 445 C 446 C 447 C END OF ROUTINE. 448 C 449 RETURN 450 END ENTRY POINTS Name ngmfld2_ Page 7 Source Listing NGMFLD2 2014-12-17 20:48 Symbol Table NGMFLD2.f SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References A2 Param 41 R(4) 4 scalar 192 A3 Param 41 R(4) 4 scalar 192 A4 Param 41 R(4) 4 scalar 192 ABS Func 279 scalar 279 AD05 Param 33 R(4) 4 scalar AETA Scalar 8 R(4) 4 1 45 COM AETAD Scalar 4 R(4) 4 1 45 COM AI Local 185 R(4) 4 scalar 185,188,193 ALPINT Scalar 4 R(4) 4 3 4219074 COM 160 ALPM Local 160 R(4) 4 scalar 160,162 ALSL Scalar 7 R(4) 4 1 39 COM AMAX1 Func 165 scalar 165 AMIN1 Func 184 scalar 184,193 BETA Param 37 R(4) 4 scalar 37 BI Local 186 R(4) 4 scalar 186,189,193 BTG Param 37 R(4) 4 scalar CAPA Param 28 R(4) 4 scalar 35 CFT0 Param 33 R(4) 4 scalar CLDWTR Common 2 17610408 CLIMIT Local 114 R(4) 4 scalar 114,168 CM1 Param 23 R(4) 4 scalar CM2 Param 23 R(4) 4 scalar CM3 Param 23 R(4) 4 scalar CMLD Scalar 6 R(4) 4 scalar COM CP Param 34 R(4) 4 scalar 36,40 CPDR Param 36 R(4) 4 scalar CPGFU Scalar 8 R(4) 4 2 91719 COM CPGFV Scalar 3 R(4) 4 scalar COM CPHI0 Scalar 9 R(4) 4 scalar COM CTLBLK Common 2 96 CURV Scalar 8 R(4) 4 2 91719 COM CWM Scalar 3 R(4) 4 3 4127355 COM 168 D00 Param 25 R(4) 4 scalar 120,121,122,123,124,125,126,127,12 8,129,130,223 D001 Param 25 R(4) 4 scalar D0065 Param 29 R(4) 4 scalar D01 Param 25 R(4) 4 scalar 217,218 D05 Param 29 R(4) 4 scalar 33 D1 Param 25 R(4) 4 scalar D115 Param 29 R(4) 4 scalar D125 Param 27 R(4) 4 scalar D25 Param 27 R(4) 4 scalar 226,227,228,229 D50 Param 25 R(4) 4 scalar 160 D608 Param 29 R(4) 4 scalar DDMPU Scalar 9 R(4) 4 2 91719 COM DDMPV Scalar 9 R(4) 4 2 91719 COM DETA Scalar 8 R(4) 4 1 45 COM DETAD Scalar 4 R(4) 4 1 45 COM DFL Scalar 4 R(4) 4 1 46 COM DISLP Scalar 7 R(4) 4 scalar COM DLMD Scalar 6 R(4) 4 scalar COM Page 8 Source Listing NGMFLD2 2014-12-17 20:48 Symbol Table NGMFLD2.f Name Object Declared Type Bytes Dimen Elements Attributes References DP10M Param 35 R(4) 4 scalar DP30 Scalar 6 R(4) 4 scalar COM DPHD Scalar 6 R(4) 4 scalar COM DT Scalar 4 R(4) 4 scalar COM DTR Param 34 R(4) 4 scalar 34 DX Scalar 8 R(4) 4 2 91719 COM 230 DY Scalar 3 R(4) 4 scalar COM 231 DYNAM Common 2 3675964 DZ Local 161 R(4) 4 scalar 161,239,240,245,246,251,252,257,25 8,263,264 EF4T Scalar 3 R(4) 4 scalar COM ELDR Param 36 R(4) 4 scalar ELWV Param 36 R(4) 4 scalar 36 EM Scalar 6 R(4) 4 1 760 COM EMT Scalar 6 R(4) 4 1 760 COM EN Scalar 3 R(4) 4 scalar COM ENT Scalar 3 R(4) 4 scalar COM EPS Param 39 R(4) 4 scalar 39 EPSQ2 Param 40 R(4) 4 scalar ERAD Param 34 R(4) 4 scalar ERLAM0 Scalar 9 R(4) 4 scalar COM ETA Scalar 8 R(4) 4 1 46 COM ETAD Scalar 4 R(4) 4 1 46 COM EXCH2 Subr 143 143 EXP Func 162 scalar 162,192 EXTRA Common 3 71506268 F Scalar 9 R(4) 4 2 91719 COM F4D Scalar 3 R(4) 4 scalar COM F4Q Scalar 3 R(4) 4 scalar COM F4Q2 Scalar 4 R(4) 4 1 45 COM FAD Scalar 9 R(4) 4 2 91719 COM FCP Scalar 8 R(4) 4 2 91719 COM FDIV Scalar 9 R(4) 4 2 91719 COM FIRST Scalar 3 R(4) 4 scalar COM FIS Scalar 3 R(4) 4 2 91719 COM FMD Param 39 R(4) 4 scalar 39 FMW Param 39 R(4) 4 scalar 39 FSL Scalar 8 R(4) 4 2 91719 COM G Param 32 R(4) 4 scalar 32 GI Param 32 R(4) 4 scalar 36 GOT1847 Local 92 L(4) 4 scalar GOT4710 Local 92 L(4) 4 scalar GOT4796 Local 92 L(4) 4 scalar GOT8498 Local 92 L(4) 4 scalar GOT8510 Local 92 L(4) 4 scalar H1 Param 26 R(4) 4 scalar 33,213,214 H100 Param 27 R(4) 4 scalar H1000 Param 27 R(4) 4 scalar H10000 Param 28 R(4) 4 scalar H10E5 Param 28 R(4) 4 scalar H1E2 Param 30 R(4) 4 scalar H1E4 Param 30 R(4) 4 scalar H1M12 Param 26 R(4) 4 scalar 165,279 H1M2 Param 30 R(4) 4 scalar H1M5 Param 27 R(4) 4 scalar Page 9 Source Listing NGMFLD2 2014-12-17 20:48 Symbol Table NGMFLD2.f Name Object Declared Type Bytes Dimen Elements Attributes References H2 Param 26 R(4) 4 scalar H4 Param 26 R(4) 4 scalar H999 Param 25 R(4) 4 scalar H99999 Param 26 R(4) 4 scalar HBM2 Scalar 3 R(4) 4 2 91719 COM 235 HTM Scalar 5 R(4) 4 3 4127355 COM I Local 119 I(4) 4 scalar 119,120,121,122,123,124,125,126,12 7,128,129,130,145,148,149,150,151, 152,153,154,160,161,163,164,167,16 8,170,172,174,177,200,224,225,226, 227,228,229,230,233,235,239,240,24 5,246,251,252,257,258,263,264,272, 274,275,277,279,281,282,284,287,28 8,290,293,294,296,299,300,302 IBESSL Scalar 2 I(4) 4 scalar COM ICNT Scalar 7 I(4) 4 1 1024 COM IDAT Scalar 3 I(4) 4 1 3 COM IDIM1 Param 14 I(4) 4 scalar IDIM2 Param 14 I(4) 4 scalar IDN Scalar 7 I(4) 4 scalar COM IDSP Scalar 7 I(4) 4 1 1024 COM IDTAD Scalar 4 I(4) 4 scalar COM IE Local 224 I(4) 4 scalar 224,226,228,229,232 IFLAG Scalar 2 I(4) 4 scalar COM IGSTL Param 9 I(4) 4 scalar 14 IGSTR Param 9 I(4) 4 scalar 14 IHE Scalar 3 I(4) 4 1 387 COM 224 IHHA Scalar 3 I(4) 4 1 760 COM IHLA Scalar 3 I(4) 4 1 760 COM IHRST Scalar 3 I(4) 4 scalar COM IHW Scalar 3 I(4) 4 1 387 COM 225 IM Param 3 I(4) 4 scalar 2,3,4,5,6,7,8,9,11,14,94,95,96,97, 119,145,272 IM1 Param 6 I(4) 4 scalar IMJM Param 2 I(4) 4 scalar 5 IMT Param 4 I(4) 4 scalar INDX Common 2 6192 INPES Param 6 I(4) 4 scalar 11,14 IOFFS Scalar 2 I(4) 4 scalar COM IOUT Scalar 3 I(4) 4 scalar COM ISHDE Scalar 4 I(4) 4 1 61 COM ITAIL Param 11 I(4) 4 scalar IUNIT Scalar 10 I(4) 4 scalar COM IUP Scalar 7 I(4) 4 scalar COM IVE Scalar 3 I(4) 4 1 387 COM IVHA Scalar 3 I(4) 4 1 760 COM IVLA Scalar 3 I(4) 4 1 760 COM IVW Scalar 3 I(4) 4 1 387 COM IW Local 225 I(4) 4 scalar 225,227,228,229,232 IWL Local 95 R(4) 4 2 91719 130,167,170,172,174,177,200 IWM1 Local 95 R(4) 4 scalar 167,174 IXM Scalar 3 I(4) 4 scalar COM IYM Scalar 3 I(4) 4 scalar COM J Local 118 I(4) 4 scalar 118,120,121,122,123,124,125,126,12 7,128,129,130,144,148,149,150,151, Page 10 Source Listing NGMFLD2 2014-12-17 20:48 Symbol Table NGMFLD2.f Name Object Declared Type Bytes Dimen Elements Attributes References 152,153,154,160,161,163,164,167,16 8,170,172,174,177,200,224,225,226, 227,228,229,230,232,233,235,239,24 0,245,246,251,252,257,258,263,264, 271,274,275,277,279,281,282,284,28 7,288,290,293,294,296,299,300,302 JAM Param 6 I(4) 4 scalar 3,6 JDIM1 Param 15 I(4) 4 scalar JDIM2 Param 15 I(4) 4 scalar JEND Scalar 6 I(4) 4 scalar COM 118 JEND_2U Scalar 8 I(4) 4 scalar COM JEND_M Scalar 7 I(4) 4 scalar COM JEND_M2 Scalar 7 I(4) 4 scalar COM 144,271 JGSTL Param 10 I(4) 4 scalar 15 JGSTR Param 10 I(4) 4 scalar 15 JM Param 3 I(4) 4 scalar 2,3,4,5,6,7,8,9,12,15,94,95,96,97 JM2 Param 6 I(4) 4 scalar JMT Param 4 I(4) 4 scalar JNPES Param 6 I(4) 4 scalar 12,15 JRA Scalar 3 I(4) 4 1 760 COM JSTA Scalar 6 I(4) 4 scalar COM 118 JSTA_2L Scalar 8 I(4) 4 scalar COM JSTA_M Scalar 7 I(4) 4 scalar COM JSTA_M2 Scalar 7 I(4) 4 scalar COM 144,271 JTAIL Param 12 I(4) 4 scalar KMNT Scalar 5 I(4) 4 2 4118670 COM KMNTM Scalar 5 I(4) 4 1 45 COM KSAV Param 14 I(4) 4 scalar KSB Scalar 2 I(4) 4 scalar COM KSFI Param 13 I(4) 4 scalar KSLI Param 14 I(4) 4 scalar KSLP Param 14 I(4) 4 scalar KSLPD Param 13 I(4) 4 scalar KSMUD Param 13 I(4) 4 scalar KSRH Param 14 I(4) 4 scalar KSSLP Param 13 I(4) 4 scalar KST Param 13 I(4) 4 scalar KSUV Param 14 I(4) 4 scalar KUNIT Scalar 10 I(4) 4 scalar COM L Local 142 I(4) 4 scalar 142,143,160,161,163,164,168,226,22 7,228,229,232,233 LB Param 9 I(4) 4 scalar LC Scalar 3 I(4) 4 2 91719 COM LHEAT Param 40 I(4) 4 scalar 40 LHTOCP Param 40 I(4) 4 scalar LIST Scalar 3 I(4) 4 scalar COM LM Param 3 I(4) 4 scalar 3,4,5,8,19,142 LM1 Param 8 I(4) 4 scalar LMH Scalar 5 I(4) 4 2 91719 COM LMV Scalar 5 I(4) 4 2 91719 COM LOOPS Common 2 748952 LP1 Param 8 I(4) 4 scalar 4,8 LSL Scalar 3 I(4) 4 scalar COM LSM Param 3 I(4) 4 scalar 7 LUNIT Scalar 10 I(4) 4 scalar COM Page 11 Source Listing NGMFLD2 2014-12-17 20:48 Symbol Table NGMFLD2.f Name Object Declared Type Bytes Dimen Elements Attributes References MAPOT Common 2 1448 MASKS Common 2 34853220 ME Scalar 6 I(4) 4 scalar COM MUNIT Scalar 10 I(4) 4 scalar COM NBC Scalar 3 I(4) 4 scalar COM NBOCO Scalar 4 I(4) 4 scalar COM NCP Scalar 4 I(4) 4 scalar COM NDDAMP Scalar 4 I(4) 4 scalar COM NFCST Scalar 3 I(4) 4 scalar COM NGMFLD2 Subr 1 NINC Param 7 I(4) 4 scalar NPHS Scalar 5 I(4) 4 scalar COM NPREC Scalar 4 I(4) 4 scalar COM NRADL Scalar 5 I(4) 4 scalar COM NRADS Scalar 5 I(4) 4 scalar COM NSHDE Scalar 4 I(4) 4 scalar COM NSTAT Param 19 I(4) 4 scalar NSUMD Param 15 I(4) 4 scalar NTSD Scalar 4 I(4) 4 scalar COM NTSTM Scalar 4 I(4) 4 scalar COM NUM_PROCS Scalar 6 I(4) 4 scalar COM NUNIT Scalar 10 I(4) 4 scalar COM NWORD Param 19 I(4) 4 scalar 19 NWORD3 Param 19 I(4) 4 scalar NWORD6 Param 19 I(4) 4 scalar OLDRD Scalar 93 L(4) 4 scalar COM ONEPS Param 39 R(4) 4 scalar OPTIONS Common 2 24 OVERRC Param 33 R(4) 4 scalar 33 P100 Local 148 R(4) 4 scalar 148,238,244 P1000 Param 35 R(4) 4 scalar P18 Local 154 R(4) 4 scalar 154,262 P47 Local 153 R(4) 4 scalar 153,244,256,262 P84 Local 152 R(4) 4 scalar 152,250 P85 Local 151 R(4) 4 scalar 151,238 P96 Local 150 R(4) 4 scalar 150,256 P98 Local 149 R(4) 4 scalar 149,250 PARALLEL Common 9 8240 PD Scalar 3 R(4) 4 2 91719 COM 148,149,150,151,152,153,154 PDSL Scalar 6 R(4) 4 2 91719 COM PDVP1 Scalar 6 R(4) 4 2 91719 COM PFSL Scalar 7 R(4) 4 2 91719 COM PINT Scalar 4 R(4) 4 3 4219074 COM PM Local 162 R(4) 4 scalar 162,191,238,244,250,256,262 PQ0 Param 41 R(4) 4 scalar 191 PSLP Scalar 6 R(4) 4 2 91719 COM PT Scalar 5 R(4) 4 scalar COM 148,149,150,151,152,153,154 PTDYN Scalar 3 R(4) 4 scalar COM PTSL Scalar 7 R(4) 4 2 91719 COM Q Scalar 5 R(4) 4 3 4127355 COM 143,164,226,227,228,229 QDIV Local 234 R(4) 4 scalar 234,235 QI Local 193 R(4) 4 scalar 193,196,201 QINT Local 194 R(4) 4 scalar 194,198,203 QK Local 226 R(4) 4 scalar 226,232 QKHN Local 228 R(4) 4 scalar 228,233 Page 12 Source Listing NGMFLD2 2014-12-17 20:48 Symbol Table NGMFLD2.f Name Object Declared Type Bytes Dimen Elements Attributes References QKHS Local 229 R(4) 4 scalar 229,233 QKM1 Local 227 R(4) 4 scalar 227,232 QM Local 164 R(4) 4 scalar 164,165,212,215,219 QM8510 Dummy 1 R(4) 4 2 91719 ARG,INOUT 120,240,275,277,279 QMCVG Local 223 R(4) 4 scalar 223,235,240 QS Local 196 R(4) 4 scalar 196,198,201,203,208,212,215,219 QSL Scalar 9 R(4) 4 2 91719 COM QUDX Local 232 R(4) 4 scalar 232,234 QVDY Local 233 R(4) 4 scalar 233,234 QW Local 191 R(4) 4 scalar 191,193,194,208 R Scalar 5 R(4) 4 scalar COM R2DX Local 230 R(4) 4 scalar 230,232 R2DY Local 231 R(4) 4 scalar 231,233 RCAPA Param 35 R(4) 4 scalar RD Param 32 R(4) 4 scalar 32,36,40 RDETA Scalar 4 R(4) 4 1 45 COM RDOCP Param 40 R(4) 4 scalar RDYN Scalar 3 R(4) 4 scalar COM RES Scalar 3 R(4) 4 2 91719 COM RESTRT Scalar 3 R(4) 4 scalar COM RG Param 36 R(4) 4 scalar RH Local 212 R(4) 4 scalar 212,213,214,215,217,218,219,246,25 2,258,264 RH1847 Dummy 1 R(4) 4 2 91719 ARG,INOUT 124,264,300,302 RH4710 Dummy 1 R(4) 4 2 91719 ARG,INOUT 121,246,282,284 RH4796 Dummy 1 R(4) 4 2 91719 ARG,INOUT 123,258,294,296 RH8498 Dummy 1 R(4) 4 2 91719 ARG,INOUT 122,252,288,290 ROG Param 32 R(4) 4 scalar RTD Param 34 R(4) 4 scalar RUN Scalar 3 R(4) 4 scalar COM SATDEL Scalar 2 R(4) 4 scalar COM SBD Scalar 6 R(4) 4 scalar COM SICE Scalar 3 R(4) 4 2 91719 COM SIG100 Param 86 R(4) 4 scalar 148 SIG18 Param 88 R(4) 4 scalar 154 SIG47 Param 88 R(4) 4 scalar 153 SIG78 Param 88 R(4) 4 scalar SIG84 Param 87 R(4) 4 scalar 152 SIG85 Param 87 R(4) 4 scalar 151 SIG89 Param 87 R(4) 4 scalar SIG96 Param 86 R(4) 4 scalar 150 SIG98 Param 86 R(4) 4 scalar 149 SIGMA Scalar 3 R(4) 4 scalar COM SLP Scalar 5 R(4) 4 2 91719 COM SM Scalar 3 R(4) 4 2 91719 COM SMALL Param 89 R(4) 4 scalar 279 SMTHA Param 15 R(4) 4 scalar SMTHB Param 15 R(4) 4 scalar SPHI0 Scalar 9 R(4) 4 scalar COM SPL Scalar 7 R(4) 4 1 39 COM SPVAL Scalar 2 R(4) 4 scalar COM 277,279,284,290,296,302 SR Scalar 3 R(4) 4 2 91719 COM STDRD Scalar 93 L(4) 4 scalar COM T Scalar 5 R(4) 4 3 4127355 COM 163 T500 Scalar 6 R(4) 4 2 91719 COM Page 13 Source Listing NGMFLD2 2014-12-17 20:48 Symbol Table NGMFLD2.f Name Object Declared Type Bytes Dimen Elements Attributes References TFRZ Param 39 R(4) 4 scalar TLM0D Scalar 6 R(4) 4 scalar COM TM Local 163 R(4) 4 scalar 163,169,171,183,192 TMT0 Local 183 R(4) 4 scalar 183,184,187,193,195,197 TMT15 Local 184 R(4) 4 scalar 184,194 TPH0D Scalar 6 R(4) 4 scalar COM TSHDE Scalar 8 R(4) 4 1 61 COM TSL Scalar 8 R(4) 4 2 91719 COM TSPH Scalar 5 R(4) 4 scalar COM TTV Scalar 8 R(4) 4 2 91719 COM U Scalar 5 R(4) 4 3 4127355 COM 232 U00 Scalar 3 R(4) 4 2 91719 COM UL Scalar 3 R(4) 4 1 90 COM V Scalar 5 R(4) 4 3 4127355 COM 233 VBM2 Scalar 3 R(4) 4 2 91719 COM VBM3 Scalar 3 R(4) 4 2 91719 COM VRBLS Common 2 67138308 VTM Scalar 5 R(4) 4 3 4127355 COM WBD Scalar 6 R(4) 4 scalar COM WPDAR Scalar 8 R(4) 4 2 91719 COM X1P Scalar 6 R(4) 4 scalar COM Y1P Scalar 6 R(4) 4 scalar COM Z0SLP Scalar 7 R(4) 4 scalar COM Z1000 Scalar 9 R(4) 4 2 91719 COM Z1847 Local 97 R(4) 4 2 91719 129,263,299,300 Z4710 Local 96 R(4) 4 2 91719 127,245,281,282 Z4796 Local 97 R(4) 4 2 91719 128,257,293,294 Z8498 Local 96 R(4) 4 2 91719 126,251,287,288 Z8510 Local 96 R(4) 4 2 91719 125,239,274,275 ZINT Scalar 4 R(4) 4 3 4219074 COM 161 Page 14 Source Listing NGMFLD2 2014-12-17 20:48 NGMFLD2.f 451 Page 15 Source Listing NGMFLD2 2014-12-17 20:48 Subprograms/Common Blocks NGMFLD2.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References CLDWTR Common 2 17610408 CTLBLK Common 2 96 DYNAM Common 2 3675964 EXTRA Common 3 71506268 INDX Common 2 6192 LOOPS Common 2 748952 MAPOT Common 2 1448 MASKS Common 2 34853220 NGMFLD2 Subr 1 OPTIONS Common 2 24 PARALLEL Common 9 8240 VRBLS Common 2 67138308 COMPILER OPTIONS BEING USED -align nocommons -align nodcommons -align noqcommons -align records -align nosequence -align norec1byte -align norec2byte -align norec4byte -align norec8byte -align norec16byte -altparam -assume accuracy_sensitive -assume nobscc -assume nobuffered_io -assume byterecl -assume nocc_omp -assume nocstring -assume nodummy_aliases -assume nofpe_summary -assume noieee_fpe_flags -assume nominus0 -assume noold_boz -assume old_unit_star -assume old_ldout_format -assume noold_logical_ldio -assume old_maxminloc -assume old_xor -assume protect_constants -assume noprotect_parens -assume split_common -assume source_include -assume nostd_intent_in -assume nostd_mod_proc_name -assume norealloc_lhs -assume underscore -assume no2underscores no -auto -auto_scalar no -bintext -ccdefault default -check noargs -check noarg_temp_created -check nobounds -check noformat -check nooutput_conversion -check nooverflow -check nopointers -check power -check noshape -check nounderflow -check nouninitialized -coarray-num-procs 0 no -coarray-config-file -convert big_endian -cross_reference -D __INTEL_COMPILER=1210 -D __unix__ -D __unix -D __linux__ -D __linux -D __gnu_linux__ -D unix -D linux -D __ELF__ -D __x86_64 -D __x86_64__ -D _MT -D __INTEL_COMPILER_BUILD_DATE=20120612 -D __i686 -D __i686__ -D __pentiumpro -D __pentiumpro__ Page 16 Source Listing NGMFLD2 2014-12-17 20:48 NGMFLD2.f -D __pentium4 -D __pentium4__ -D __tune_pentium4__ -D __SSE2__ -D __SSE__ -D __MMX__ -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 no -ftz -fp_model noprecise -fp_model fast -fp_model nostrict -fp_model nosource -fp_model nodouble -fp_model noextended -fp_model novery_fast -fp_model noexcept -fp_model nono_except -heap_arrays 0 no -threadprivate_compat -g2 -iface nomixed_str_len_arg -iface nono_mixed_str_len_arg no -intconstant -integer_size 32 no -mixed_str_len_arg no -module -names lowercase no -noinclude -O2 no -pad_source -real_size 32 no -recursive -reentrancy none no -sharable_localsaves -vec=simd -show noinclude -show map -show options no -syntax_only no -threadcom no -U no -vms -w noall -w nonone -w alignments -w noargument_checking -w nodeclarations -w general -w noignore_bounds -w noignore_loc -w nointerfaces -w notruncated_source -w uncalled -w uninitialized -w nounused -w usage -includepath : /gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/,.f,./.f,/usrx/local/intel/composerxe/mkl/include/.f, /usrx/local/intel/composerxe/tbb/include/.f,/gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/intel64/.f, /gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/.f,/usr/local/include/.f,/usr/lib/gcc/x86_64-redhat-linux/4.4.7/include/.f, /usr/include/.f,/usr/include/.f -list filename : NGMFLD2.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100