Page 1 Source Listing NGMFLD2 2025-03-12 18:23 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 2025-03-12 18:23 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 2025-03-12 18:23 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 2025-03-12 18:23 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 2025-03-12 18:23 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 2025-03-12 18:23 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 2025-03-12 18:23 Symbol Table NGMFLD2.f SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References A2 Param 143 R(4) 4 scalar 331 A3 Param 143 R(4) 4 scalar 331 A4 Param 143 R(4) 4 scalar 331 ABS Func 418 scalar 418 AD05 Param 135 R(4) 4 scalar AI Local 324 R(4) 4 scalar 324,327,332 ALPM Local 299 R(4) 4 scalar 299,301 AMAX1 Func 304 scalar 304 AMIN1 Func 323 scalar 323,332 BETA Param 139 R(4) 4 scalar 139 BI Local 325 R(4) 4 scalar 325,328,332 BTG Param 139 R(4) 4 scalar CAPA Param 130 R(4) 4 scalar 137 CFT0 Param 135 R(4) 4 scalar CLDWTR Common 167 17610408 SAVE CLIMIT Local 253 R(4) 4 scalar 253,307 CM1 Param 125 R(4) 4 scalar CM2 Param 125 R(4) 4 scalar CM3 Param 125 R(4) 4 scalar CP Param 136 R(4) 4 scalar 138,142 CPDR Param 138 R(4) 4 scalar CTLBLK Common 240 96 SAVE D00 Param 127 R(4) 4 scalar 259,260,261,262,263,264,265,266,26 7,268,269,362 D001 Param 127 R(4) 4 scalar D0065 Param 131 R(4) 4 scalar D01 Param 127 R(4) 4 scalar 356,357 D05 Param 131 R(4) 4 scalar 135 D1 Param 127 R(4) 4 scalar D115 Param 131 R(4) 4 scalar D125 Param 129 R(4) 4 scalar D25 Param 129 R(4) 4 scalar 365,366,367,368 D50 Param 127 R(4) 4 scalar 299 D608 Param 131 R(4) 4 scalar DP10M Param 137 R(4) 4 scalar DTR Param 136 R(4) 4 scalar 136 DYNAM Common 218 3675964 SAVE DZ Local 300 R(4) 4 scalar 300,378,379,384,385,390,391,396,39 7,402,403 ELDR Param 138 R(4) 4 scalar ELWV Param 138 R(4) 4 scalar 138 EPS Param 141 R(4) 4 scalar 141 EPSQ2 Param 142 R(4) 4 scalar ERAD Param 136 R(4) 4 scalar EXCH2 Subr 282 282 EXP Func 301 scalar 301,331 EXTRA Common 182 71506268 SAVE FMD Param 141 R(4) 4 scalar 141 FMW Param 141 R(4) 4 scalar 141 G Param 134 R(4) 4 scalar 134 GI Param 134 R(4) 4 scalar 138 Page 8 Source Listing NGMFLD2 2025-03-12 18:23 Symbol Table NGMFLD2.f Name Object Declared Type Bytes Dimen Elements Attributes References GOT1847 Local 150 L(4) 4 scalar GOT4710 Local 150 L(4) 4 scalar GOT4796 Local 150 L(4) 4 scalar GOT8498 Local 150 L(4) 4 scalar GOT8510 Local 150 L(4) 4 scalar H1 Param 128 R(4) 4 scalar 135,352,353 H100 Param 129 R(4) 4 scalar H1000 Param 129 R(4) 4 scalar H10000 Param 130 R(4) 4 scalar H10E5 Param 130 R(4) 4 scalar H1E2 Param 132 R(4) 4 scalar H1E4 Param 132 R(4) 4 scalar H1M12 Param 128 R(4) 4 scalar 304,418 H1M2 Param 132 R(4) 4 scalar H1M5 Param 129 R(4) 4 scalar H2 Param 128 R(4) 4 scalar H4 Param 128 R(4) 4 scalar H999 Param 127 R(4) 4 scalar H99999 Param 128 R(4) 4 scalar I Local 258 I(4) 4 scalar 258,259,260,261,262,263,264,265,26 6,267,268,269,284,287,288,289,290, 291,292,293,299,300,302,303,306,30 7,309,311,313,316,339,363,364,365, 366,367,368,369,372,374,378,379,38 4,385,390,391,396,397,402,403,411, 413,414,416,418,420,421,423,426,42 7,429,432,433,435,438,439,441 IDIM1 Param 97 I(4) 4 scalar IDIM2 Param 97 I(4) 4 scalar IE Local 363 I(4) 4 scalar 363,365,367,368,371 IGSTL Param 92 I(4) 4 scalar 97 IGSTR Param 92 I(4) 4 scalar 97 IM Param 86 I(4) 4 scalar 94,97,106,108,109,111,152,153,154, 155,161,163,168,174,176,181,183,18 4,185,186,187,188,200,224,225,258, 284,411 IM1 Param 108 I(4) 4 scalar IMJM Param 181 I(4) 4 scalar 184 IMT Param 106 I(4) 4 scalar INDX Common 230 6192 SAVE INPES Param 89 I(4) 4 scalar 94,97 ITAIL Param 94 I(4) 4 scalar IW Local 364 I(4) 4 scalar 364,366,367,368,371 IWL Local 153 R(4) 4 2 91719 269,306,309,311,313,316,339 IWM1 Local 153 R(4) 4 scalar 306,313 J Local 257 I(4) 4 scalar 257,259,260,261,262,263,264,265,26 6,267,268,269,283,287,288,289,290, 291,292,293,299,300,302,303,306,30 7,309,311,313,316,339,363,364,365, 366,367,368,369,371,372,374,378,37 9,384,385,390,391,396,397,402,403, 410,413,414,416,418,420,421,423,42 6,427,429,432,433,435,438,439,441 JAM Param 108 I(4) 4 scalar 198,222 JDIM1 Param 98 I(4) 4 scalar Page 9 Source Listing NGMFLD2 2025-03-12 18:23 Symbol Table NGMFLD2.f Name Object Declared Type Bytes Dimen Elements Attributes References JDIM2 Param 98 I(4) 4 scalar JGSTL Param 93 I(4) 4 scalar 98 JGSTR Param 93 I(4) 4 scalar 98 JM Param 86 I(4) 4 scalar 95,98,106,108,111,152,153,154,155, 161,163,168,174,176,181,183,184,18 5,186,187,188,200,224,225,231 JM2 Param 108 I(4) 4 scalar JMT Param 106 I(4) 4 scalar JNPES Param 89 I(4) 4 scalar 95,98 JTAIL Param 95 I(4) 4 scalar KSAV Param 116 I(4) 4 scalar KSFI Param 115 I(4) 4 scalar KSLI Param 116 I(4) 4 scalar KSLP Param 116 I(4) 4 scalar KSLPD Param 115 I(4) 4 scalar KSMUD Param 115 I(4) 4 scalar KSRH Param 116 I(4) 4 scalar KSSLP Param 115 I(4) 4 scalar KST Param 115 I(4) 4 scalar KSUV Param 116 I(4) 4 scalar L Local 281 I(4) 4 scalar 281,282,299,300,302,303,307,365,36 6,367,368,371,372 LB Param 111 I(4) 4 scalar LHEAT Param 142 I(4) 4 scalar 142 LHTOCP Param 142 I(4) 4 scalar LM Param 86 I(4) 4 scalar 110,121,163,168,176,184,211,220,28 1 LM1 Param 110 I(4) 4 scalar LOOPS Common 197 748952 SAVE LP1 Param 110 I(4) 4 scalar 183,211,220 LSM Param 86 I(4) 4 scalar 210 MAPOT Common 205 1448 SAVE MASKS Common 173 34853220 SAVE NGMFLD2 Subr 1 NINC Param 109 I(4) 4 scalar NSTAT Param 121 I(4) 4 scalar NSUMD Param 117 I(4) 4 scalar NWORD Param 121 I(4) 4 scalar 121 NWORD3 Param 121 I(4) 4 scalar NWORD6 Param 121 I(4) 4 scalar ONEPS Param 141 R(4) 4 scalar OPTIONS Common 193 24 SAVE OVERRC Param 135 R(4) 4 scalar 135 P100 Local 287 R(4) 4 scalar 287,377,383 P1000 Param 137 R(4) 4 scalar P18 Local 293 R(4) 4 scalar 293,401 P47 Local 292 R(4) 4 scalar 292,383,395,401 P84 Local 291 R(4) 4 scalar 291,389 P85 Local 290 R(4) 4 scalar 290,377 P96 Local 289 R(4) 4 scalar 289,395 P98 Local 288 R(4) 4 scalar 288,389 PARALLEL Common 247 8240 PM Local 301 R(4) 4 scalar 301,330,377,383,389,395,401 PQ0 Param 143 R(4) 4 scalar 330 QDIV Local 373 R(4) 4 scalar 373,374 Page 10 Source Listing NGMFLD2 2025-03-12 18:23 Symbol Table NGMFLD2.f Name Object Declared Type Bytes Dimen Elements Attributes References QI Local 332 R(4) 4 scalar 332,335,340 QINT Local 333 R(4) 4 scalar 333,337,342 QK Local 365 R(4) 4 scalar 365,371 QKHN Local 367 R(4) 4 scalar 367,372 QKHS Local 368 R(4) 4 scalar 368,372 QKM1 Local 366 R(4) 4 scalar 366,371 QM Local 303 R(4) 4 scalar 303,304,351,354,358 QM8510 Dummy 1 R(4) 4 2 91719 ARG,INOUT 259,379,414,416,418 QMCVG Local 362 R(4) 4 scalar 362,374,379 QS Local 335 R(4) 4 scalar 335,337,340,342,347,351,354,358 QUDX Local 371 R(4) 4 scalar 371,373 QVDY Local 372 R(4) 4 scalar 372,373 QW Local 330 R(4) 4 scalar 330,332,333,347 R2DX Local 369 R(4) 4 scalar 369,371 R2DY Local 370 R(4) 4 scalar 370,372 RCAPA Param 137 R(4) 4 scalar RD Param 134 R(4) 4 scalar 134,138,142 RDOCP Param 142 R(4) 4 scalar RG Param 138 R(4) 4 scalar RH Local 351 R(4) 4 scalar 351,352,353,354,356,357,358,385,39 1,397,403 RH1847 Dummy 1 R(4) 4 2 91719 ARG,INOUT 263,403,439,441 RH4710 Dummy 1 R(4) 4 2 91719 ARG,INOUT 260,385,421,423 RH4796 Dummy 1 R(4) 4 2 91719 ARG,INOUT 262,397,433,435 RH8498 Dummy 1 R(4) 4 2 91719 ARG,INOUT 261,391,427,429 ROG Param 134 R(4) 4 scalar RTD Param 136 R(4) 4 scalar SIG100 Param 144 R(4) 4 scalar 287 SIG18 Param 146 R(4) 4 scalar 293 SIG47 Param 146 R(4) 4 scalar 292 SIG78 Param 146 R(4) 4 scalar SIG84 Param 145 R(4) 4 scalar 291 SIG85 Param 145 R(4) 4 scalar 290 SIG89 Param 145 R(4) 4 scalar SIG96 Param 144 R(4) 4 scalar 289 SIG98 Param 144 R(4) 4 scalar 288 SMALL Param 147 R(4) 4 scalar 418 SMTHA Param 117 R(4) 4 scalar SMTHB Param 117 R(4) 4 scalar TFRZ Param 141 R(4) 4 scalar TM Local 302 R(4) 4 scalar 302,308,310,322,331 TMT0 Local 322 R(4) 4 scalar 322,323,326,332,334,336 TMT15 Local 323 R(4) 4 scalar 323,333 VRBLS Common 160 67138308 Z1847 Local 155 R(4) 4 2 91719 268,402,438,439 Z4710 Local 154 R(4) 4 2 91719 266,384,420,421 Z4796 Local 155 R(4) 4 2 91719 267,396,432,433 Z8498 Local 154 R(4) 4 2 91719 265,390,426,427 Z8510 Local 154 R(4) 4 2 91719 264,378,413,414 Page 11 Source Listing NGMFLD2 2025-03-12 18:23 Symbol Table NGMFLD2.f TYPE COMPONENTS/COMMON VARIABLES Name Type Bytes Offset Dimen Elements Attributes References AETA R(4) 4 808 1 45 COM AETAD R(4) 4 396 1 45 COM ALPINT R(4) 4 16876296 3 4219074 COM 299 ALSL R(4) 4 472 1 39 COM CMLD R(4) 4 292 scalar COM CPGFU R(4) 4 740956 2 91719 COM CPGFV R(4) 4 4 scalar COM CPHI0 R(4) 4 1420 scalar COM CURV R(4) 4 1107832 2 91719 COM CWM R(4) 4 0 3 4127355 COM 307 DDMPU R(4) 4 2942212 2 91719 COM DDMPV R(4) 4 3309088 2 91719 COM DETA R(4) 4 628 1 45 COM DETAD R(4) 4 36 1 45 COM DFL R(4) 4 940 1 46 COM DISLP R(4) 4 308 scalar COM DLMD R(4) 4 284 scalar COM DP30 R(4) 4 296 scalar COM DPHD R(4) 4 288 scalar COM DT R(4) 4 48 scalar COM DX R(4) 4 7204 2 91719 COM 369 DY R(4) 4 0 scalar COM 370 EF4T R(4) 4 32 scalar COM EM R(4) 4 1124 1 760 COM EMT R(4) 4 4164 1 760 COM EN R(4) 4 8 scalar COM ENT R(4) 4 12 scalar COM ERLAM0 R(4) 4 1416 scalar COM ETA R(4) 4 988 1 46 COM ETAD R(4) 4 756 1 46 COM F R(4) 4 2575336 2 91719 COM F4D R(4) 4 24 scalar COM F4Q R(4) 4 28 scalar COM F4Q2 R(4) 4 576 1 45 COM FAD R(4) 4 2208460 2 91719 COM FCP R(4) 4 1474708 2 91719 COM FDIV R(4) 4 1841584 2 91719 COM FIRST R(4) 4 24 scalar COM FIS R(4) 4 733752 2 91719 COM FSL R(4) 4 70405632 2 91719 COM HBM2 R(4) 4 0 2 91719 COM 374 HTM R(4) 4 1834380 3 4127355 COM IBESSL I(4) 4 4 scalar COM ICNT I(4) 4 40 1 1024 COM IDAT I(4) 4 8 1 3 COM IDN I(4) 4 36 scalar COM IDSP I(4) 4 4136 1 1024 COM IDTAD I(4) 4 68 scalar COM IFLAG I(4) 4 16 scalar COM IHE I(4) 4 0 1 387 COM 363 IHHA I(4) 4 3040 1 760 COM IHLA I(4) 4 0 1 760 COM IHRST I(4) 4 20 scalar COM Page 12 Source Listing NGMFLD2 2025-03-12 18:23 Symbol Table NGMFLD2.f Name Type Bytes Offset Dimen Elements Attributes References IHW I(4) 4 1548 1 387 COM 364 IOFFS I(4) 4 12 scalar COM IOUT I(4) 4 44 scalar COM ISHDE I(4) 4 12 1 61 COM IUNIT I(4) 4 1440 scalar COM IUP I(4) 4 32 scalar COM IVE I(4) 4 3096 1 387 COM IVHA I(4) 4 9120 1 760 COM IVLA I(4) 4 6080 1 760 COM IVW I(4) 4 4644 1 387 COM IXM I(4) 4 4 scalar COM IYM I(4) 4 8 scalar COM JEND I(4) 4 12 scalar COM 257 JEND_2U I(4) 4 8236 scalar COM JEND_M I(4) 4 20 scalar COM JEND_M2 I(4) 4 28 scalar COM 283,410 JRA I(4) 4 12160 1 760 COM JSTA I(4) 4 8 scalar COM 257 JSTA_2L I(4) 4 8232 scalar COM JSTA_M I(4) 4 16 scalar COM JSTA_M2 I(4) 4 24 scalar COM 283,410 KMNT I(4) 4 50629068 2 4118670 COM KMNTM I(4) 4 50628888 1 45 COM KSB I(4) 4 8 scalar COM KUNIT I(4) 4 1428 scalar COM LC I(4) 4 16876656 2 91719 COM LIST I(4) 4 40 scalar COM LMH I(4) 4 15200 2 91719 COM LMV I(4) 4 382076 2 91719 COM LSL I(4) 4 0 scalar COM LUNIT I(4) 4 1432 scalar COM ME I(4) 4 4 scalar COM MUNIT I(4) 4 1436 scalar COM NBC I(4) 4 36 scalar COM NBOCO I(4) 4 72 scalar COM NCP I(4) 4 80 scalar COM NDDAMP I(4) 4 60 scalar COM NFCST I(4) 4 32 scalar COM NPHS I(4) 4 84 scalar COM NPREC I(4) 4 64 scalar COM NRADL I(4) 4 92 scalar COM NRADS I(4) 4 88 scalar COM NSHDE I(4) 4 76 scalar COM NTSD I(4) 4 52 scalar COM NTSTM I(4) 4 56 scalar COM NUM_PROCS I(4) 4 0 scalar COM NUNIT I(4) 4 1444 scalar COM OLDRD L(4) 4 71506260 scalar COM PD R(4) 4 0 2 91719 COM 287,288,289,290,291,292,293 PDSL R(4) 4 67837500 2 91719 COM PDVP1 R(4) 4 68204376 2 91719 COM PFSL R(4) 4 69305004 2 91719 COM PINT R(4) 4 0 3 4219074 COM PSLP R(4) 4 68571252 2 91719 COM PT R(4) 4 260 scalar COM 287,288,289,290,291,292,293 Page 13 Source Listing NGMFLD2 2025-03-12 18:23 Symbol Table NGMFLD2.f Name Type Bytes Offset Dimen Elements Attributes References PTDYN R(4) 4 20 scalar COM PTSL R(4) 4 68938128 2 91719 COM Q R(4) 4 50628888 3 4127355 COM 282,303,365,366,367,368 QSL R(4) 4 70772508 2 91719 COM R R(4) 4 256 scalar COM RDETA R(4) 4 216 1 45 COM RDYN R(4) 4 16 scalar COM RES R(4) 4 366876 2 91719 COM RESTRT R(4) 4 28 scalar COM RUN R(4) 4 4 scalar COM SATDEL R(4) 4 20 scalar COM SBD R(4) 4 272 scalar COM SICE R(4) 4 1467504 2 91719 COM SIGMA R(4) 4 0 scalar COM SLP R(4) 4 67103748 2 91719 COM SM R(4) 4 1100628 2 91719 COM SPHI0 R(4) 4 1424 scalar COM SPL R(4) 4 316 1 39 COM SPVAL R(4) 4 0 scalar COM 416,418,423,429,435,441 SR R(4) 4 17243532 2 91719 COM STDRD L(4) 4 71506264 scalar COM T R(4) 4 34119468 3 4127355 COM 302 T500 R(4) 4 67470624 2 91719 COM TLM0D R(4) 4 276 scalar COM TPH0D R(4) 4 280 scalar COM TSHDE R(4) 4 1172 1 61 COM TSL R(4) 4 69671880 2 91719 COM TSPH R(4) 4 264 scalar COM TTV R(4) 4 70038756 2 91719 COM U R(4) 4 1100628 3 4127355 COM 371 U00 R(4) 4 16509420 2 91719 COM UL R(4) 4 16876296 1 90 COM V R(4) 4 17610048 3 4127355 COM 372 VBM2 R(4) 4 366876 2 91719 COM VBM3 R(4) 4 733752 2 91719 COM VTM R(4) 4 18343800 3 4127355 COM WBD R(4) 4 268 scalar COM WPDAR R(4) 4 374080 2 91719 COM X1P R(4) 4 300 scalar COM Y1P R(4) 4 304 scalar COM Z0SLP R(4) 4 312 scalar COM Z1000 R(4) 4 71139384 2 91719 COM ZINT R(4) 4 33752592 3 4219074 COM 300 Page 14 Source Listing NGMFLD2 2025-03-12 18:23 NGMFLD2.f 451 Page 15 Source Listing NGMFLD2 2025-03-12 18:23 Subprograms/Common Blocks NGMFLD2.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References CLDWTR Common 167 17610408 SAVE CTLBLK Common 240 96 SAVE DYNAM Common 218 3675964 SAVE EXTRA Common 182 71506268 SAVE INDX Common 230 6192 SAVE LOOPS Common 197 748952 SAVE MAPOT Common 205 1448 SAVE MASKS Common 173 34853220 SAVE NGMFLD2 Subr 1 OPTIONS Common 193 24 SAVE PARALLEL Common 247 8240 VRBLS Common 160 67138308 COMPILER OPTIONS BEING USED -align noall -align nonone -align nocommons -align nodcommons -align noqcommons -align nozcommons -align records -align nosequence -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 Page 16 Source Listing NGMFLD2 2025-03-12 18:23 NGMFLD2.f -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 -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__ -D __CRAY_X86_ROME -D __CRAYXT_COMPUTE_LINUX_TARGET -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 -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 -offload-build=host -openmp-simd -O2 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 Page 17 Source Listing NGMFLD2 2025-03-12 18:23 NGMFLD2.f 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/pstl/include/,/pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/, .f90,./.f90,/opt/cray/pe/mpich/8.1.12/ofi/intel/19.0/include/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/ipp/include/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/mkl/include/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/include/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/stdlib/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/tbb/include/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/intel64/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/icc/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/.f90,/usr/lib64/gcc/x86_64-suse-linux/7/include/.f90, /usr/lib64/gcc/x86_64-suse-linux/7/include-fixed/.f90,/usr/include/.f90,/usr/include/.f90,/usr/include/.f90 -list filename : NGMFLD2.lst no -o COMPILER: Intel(R) Fortran 19.1-1655