Page 1 Source Listing NGMSLP2 2014-12-17 20:48 NGMSLP2.f 1 SUBROUTINE NGMSLP2 2 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 3 C . . . 4 C SUBPROGRAM: NGMSLP2 NMC SEA LEVEL PRESSURE REDUCTION 5 C PRGRMMR: TREADON ORG: W/NP2 DATE: 93-02-02 6 C 7 C ABSTRACT: 8 C 9 C THIS ROUTINE COMPUTES SEA LEVEL PRESSURE USING THE 10 C HYDROSTATIC EQUATION WITH THE SHUELL CORRECTION. THE 11 C FOLLOWING IS BASED ON DOCUMENTATION IN SUBROUTINE 12 C OUTHYDRO OF THE NGM: 13 C 14 C THE FUNDAMENTAL HYDROSTATIC EQUATION IS 15 C D(HEIGHT) 16 C --------- = TAU = VIRTUAL TEMPERATURE * (RGAS/GRAVITY) 17 C D (Z) 18 C WHERE 19 C Z = MINUS LOG OF PRESSURE (-LN(P)). 20 C 21 C SEA-LEVEL PRESSURE IS COMPUTED FROM THE FORMULA 22 C PRESS(MSL) = PRESS(GROUND) * EXP( F) 23 C WHERE 24 C F = HEIGHT OF GROUND / MEAN TAU 25 C MEAN TAU = ( TAU(GRND) + TAU(SL) ) / 2 26 C 27 C IN THE NGM TAU(GRND) AND TAU(SL) ARE FIRST SET USING A 28 C 6.5DEG/KM LAPSE RATE FROM SIGMA LEVEL K=1. THIS IS MODIFIED 29 C BY A CORRECTION BASED ON THE CRITICAL TAU OF THE SHUELL 30 C CORRECTION: 31 C TAUCR=(RGASD/GRAVITY) * 290.66 32 C 33 C 1) WHERE ONLY TAU(SL) EXCEEDS TAUCR, CHANGE TAU(SL) TO TAUCR. 34 C 35 C 2) WHERE BOTH TAU(SL) AND TAU(GRND) EXCEED TAUCR, 36 C CHANGE TAU(SL) TO TAUCR-CONST*(TAU(GRND)-TAUCR )**2 37 C WHERE CONST = .005 (GRAVITY/RGASD) 38 C 39 C THE AVERAGE OF TAU(SL) AND TAU(GRND) IS THEN USED TOGETHER 40 C WITH THE GROUND HEIGHT AND PRESSURE TO DERIVE THE PRESSURE 41 C AT SEA LEVEL. 42 C 43 C EXTRAPOLATING FROM TEMPERATURES IN THE FIRST ATMOSPHERIC ETA 44 C LAYER LEAD TO EXCESSIVELY NOISY UNDERGROUND TEMPERATURE FIELDS. 45 C AFTER EXPERIMENTATION WE OPTED TO USE A MEAN TEMPERATURE FROM 46 C THE LOWEST DP=30MB ABOVE THE GROUND. FROM THIS MEAN TEMPERATURE 47 C WE EXTRAPOLATE TO GET THE GROUND AND SEA LEVEL TEMPERATURES. 48 C USING A TRUE MASS WEIGHTED LAYER MEAN INSTEAD OF THE SIMPLE 49 C ARITHMETIC MEAN CODED SHOWED LITTLE QUALITATIVE DIFFERENCE. 50 C 51 C HEIGHT OF THE 1000MB SURFACE IS COMPUTED FROM THE MSL PRESSURE 52 C FIELD USING THE FORMULA: 53 C 54 C P(MSL) - P(1000MB) = MEAN DENSITY * GRAVITY * HGT(1000MBS) 55 C 56 C WHERE P(MSL) IS THE SEA LEVEL PRESSURE FIELD WE HAVE JUST 57 C COMPUTED. Page 2 Source Listing NGMSLP2 2014-12-17 20:48 NGMSLP2.f 58 C 59 C BELOW GROUND VIRTUAL TEMPERATURES ARE COMPUTED ASSUMING A CONSTANT 60 C 6.5DEG/KM LAPSE RATE. TO CONVERT VIRTUAL TEMPERATURE TO A DRY 61 C TEMPERATURE WE USE THE DP=30MB MEAN SPECIFIC HUMIDITY. 62 C 63 C LINE COMMENTED OUT WITH 'CX' ARE LEFT OVER FROM EARLIER 64 C EXPERIMENTATION. THEY REMAIN TO REMIND READERS WHAT 65 C WAS TRIED. 66 C . 67 C 68 C PROGRAM HISTORY LOG: 69 C 93-02-02 RUSS TREADON 70 C 98-06-08 T BLACK - CONVERSION FROM 1-D TO 2-D 71 C 00-01-04 JIM TUCCILLO - MPI VERSION 72 C 73 C USAGE: CALL NGMSLP2 74 C INPUT ARGUMENT LIST: 75 C NONE 76 C 77 C OUTPUT ARGUMENT LIST: 78 C NONE 79 C 80 C OUTPUT FILES: 81 C NONE 82 C 83 C SUBPROGRAMS CALLED: 84 C UTILITIES: 85 C NONE 86 C LIBRARY: 87 C COMMON - VRBLS 88 C EXTRA 89 C MAPOT 90 C LOOPS 91 C DYNAMD 92 C MASKS 93 C PVRBLS 94 C INDX 95 C 96 C ATTRIBUTES: 97 C LANGUAGE: FORTRAN 98 C MACHINE : CRAY C-90 99 C$$$ 100 C 101 C 102 C INCLUDE PARAMETERS. SET LOCAL PARAMETERS. 103 INCLUDE "parmeta" 121 INCLUDE "params" 163 PARAMETER (GAMMA=6.5/1000.,ZSL=0.0) 164 PARAMETER (TAUCR=RD*GI*290.66,CONST=0.005*G/RD) 165 PARAMETER (GORD=G/RD,DP=60.E2) 166 PARAMETER (ISMTHT=4,ISMTHQ=2) 167 C 168 C DECLARE VARIABLES 169 REAL TBAR(IM,JM),QBAR(IM,JM),ALPBAR(IM,JM) 170 REAL TSFC(IM,JM),QSFC(IM,JM),EGRID(IM,JM) 171 C 172 C INCLUDE COMMON BLOCKS Page 3 Source Listing NGMSLP2 2014-12-17 20:48 NGMSLP2.f 173 INCLUDE "VRBLS.comm" 180 INCLUDE "EXTRA.comm" 192 INCLUDE "MAPOT.comm" 205 INCLUDE "LOOPS.comm" 213 INCLUDE "DYNAMD.comm" 225 INCLUDE "MASKS.comm" 233 INCLUDE "PVRBLS.comm" 245 INCLUDE "INDX.comm" 251 INCLUDE "CTLBLK.comm" 252 C 253 C********************************************************************** 254 C START NGMSLP HERE. 255 C 256 C LOOP OVER HORIZONTAL GRID. 257 C 258 !$omp parallel do 259 !$omp& private(alpm,alps,l,llmh,nlev,pm,psfc,ptop,qsum,thsfc,tsum) 274 DO 30 J=JSTA,JEND 275 DO 30 I=1,IM 276 LLMH = LMH(I,J) 277 C 278 C LOCATE TOP OF LAYER OVER WHICH TO COMPUTE MEAN FIELDS. 279 PSFC = PD(I,J)+PT 280 PTOP = PSFC-DP 281 QSFC(I,J) = (1.-SM(I,J))*QS(I,J)+SM(I,J)*QZ0(I,J) 282 QSFC(I,J) = AMAX1(H1M12,QSFC(I,J)) 283 THSFC = (1.-SM(I,J))*THS(I,J)+SM(I,J)*THZ0(I,J) 284 TSFC(I,J) = THSFC*(P1000/PSFC)**(-CAPA) 285 C 286 C COMPUTE MEAN FIELDS. 287 NLEV = 1 288 ALPS = D50*(ALPINT(I,J,LLMH)+ALPINT(I,J,LLMH+1)) 289 TSUM = T(I,J,LLMH) 290 QSUM = Q(I,J,LLMH) 291 TBAR(I,J) = TSUM 292 QBAR(I,J) = QSUM 293 ALPBAR(I,J) = ALPS 294 IF (LLMH.EQ.LM) GOTO 30 295 DO 10 L = LLMH-1,1,-1 296 ALPM = D50*(ALPINT(I,J,L)+ALPINT(I,J,L+1)) 297 PM = EXP(ALPM) 298 IF (PM.LT.PTOP) GOTO 20 299 NLEV = NLEV + 1 300 ALPS = ALPS + ALPM 301 TSUM = TSUM + T(I,J,L) 302 QSUM = QSUM + Q(I,J,L) 303 10 CONTINUE 304 20 CONTINUE 305 RNLEV = 1./NLEV 306 ALPBAR(I,J) = ALPS*RNLEV 307 TBAR(I,J) = TSUM*RNLEV 308 QBAR(I,J) = QSUM*RNLEV 309 30 CONTINUE 310 C 311 C SMOOTH LAYER MEAN TEMPERATURE AND SPECIFIC HUMIDITY. 312 C BOUND SPECIFIC HUMIDITY SO THAT IT IS NON-NEGATIVE. 313 C Page 4 Source Listing NGMSLP2 2014-12-17 20:48 NGMSLP2.f 314 CALL P2FILT(ISMTHT,HBM2,TBAR) 315 CALL P2FILT(ISMTHT,HBM2,TSFC) 316 CALL P2FILT(ISMTHQ,HBM2,QBAR) 317 CALL P2FILT(ISMTHQ,HBM2,QSFC) 318 CALL BOUNDL(QBAR,H1M12,H99999,IM,JM) 319 CALL BOUNDL(QSFC,H1M12,H99999,IM,JM) 320 C 321 C LOOP OVER HORIZONTAL GRID. 322 C 323 !$omp parallel do 324 !$omp& private(alpavg,alpsfc,llmh,pavg,psfc,ptop,qavg,rhoavg,rrhog, 325 !$omp& tau,tauavg,tausfc,tausl,tavg,tvrbar,tvrsfc,tvrsl, 326 !$omp& tvrt,tvrtal,zbar,zl,zsfc) 327 DO 50 J=JSTA,JEND 328 DO 50 I=1,IM 329 LLMH = LMH(I,J) 330 ZSFC = FIS(I,J)*GI 331 PSFC = PD(I,J)+PT 332 ALPSFC = ALOG(PSFC) 333 PTOP = PSFC-DP 334 SLP(I,J) = PSFC 335 CX IF (LLMH.EQ.LM) GOTO 50 336 C 337 C COMPUTE LAYER MEAN TAU (VIRTUAL TEMP*RD/G). 338 TVRT = TBAR(I,J)*(H1+D608*QBAR(I,J)) 339 TAU = TVRT*RD*GI 340 C 341 C COMPUTE HEIGHT OF MEAN ATMOSPHERIC LAYER. 342 CX QSFC = (1.-SM(I,J))*QS(I,J)+SM(I,J)*QZ0(I,J) 343 CX QSFC = AMAX1(H1M12,QSFC) 344 QAVG = D50*(QSFC(I,J)+QBAR(I,J)) 345 CX THSFC = (1.-SM(I,J))*THS(I,J)+SM(I,J)*THZ0(I,J) 346 CX TSFC = THSFC*(P1000/PSFC)**(-CAPA) 347 TAVG = D50*(TSFC(I,J)+TBAR(I,J)) 348 TVRBAR = TAVG*(H1+D608*QAVG) 349 ZBAR = RD*GI*TVRBAR*(ALPSFC-ALPBAR(I,J)) + ZSFC 350 C 351 C COMPUTE TAU AT THE GROUND (Z=ZSFC) AND SEA LEVEL (Z=0) 352 C ASSUMING A CONSTANT LAPSE RATE OF GAMMA=6.5DEG/KM. 353 TVRSFC = TVRT + (ZBAR-ZSFC)*GAMMA 354 TAUSFC = TVRSFC*RD*GI 355 TVRSL = TVRT + (ZBAR- ZSL)*GAMMA 356 TAUSL = TVRSL*RD*GI 357 C 358 C IF NEED BE APPLY SHEULL CORRECTION. 359 IF ((TAUSL.GT.TAUCR).AND.(TAUSFC.LE.TAUCR)) THEN 360 TAUSL=TAUCR 361 ELSEIF ((TAUSL.GT.TAUCR).AND.(TAUSFC.GT.TAUCR)) THEN 362 TAUSL = TAUCR-CONST*(TAUSFC-TAUCR)**2 363 ENDIF 364 C 365 C COMPUTE MEAN TAU. 366 TAUAVG = D50*(TAUSL+TAUSFC) 367 C 368 C COMPUTE SEA LEVEL PRESSURE. 369 IF (LLMH.LT.LM) SLP(I,J) = PSFC*EXP(ZSFC/TAUAVG) 370 C Page 5 Source Listing NGMSLP2 2014-12-17 20:48 NGMSLP2.f 371 C COMPUTE 1000MB HEIGHTS. 372 ALPAVG = D50*(ALOG(PSFC)+ALOG(SLP(I,J))) 373 PAVG = EXP(ALPAVG) 374 RHOAVG = PAVG*GI/TAUAVG 375 RRHOG = H1/(RHOAVG*G) 376 Z1000(I,J) = (SLP(I,J)-P1000)*RRHOG 377 C 378 C COMPUTE TEMPERATURES ON ETA LEVELS BELOW GROUND. 379 IF ((ZSFC-ZSL).EQ.0.) GOTO 50 380 CX TLAPSE = GORD*(TAUSFC-TAUSL)/(ZSFC-ZSL) 381 IF (LLMH.EQ.LM) GOTO 50 382 DO 40 L = LLMH+1,LM 383 ZL = GI*DFL(L) 384 TVRTL = TVRT + (ZBAR-ZL)*GAMMA 385 T(I,J,L) = TVRTL/(H1+D608*QBAR(I,J)) 386 40 CONTINUE 387 C 388 C MOVE TO NEXT HORIZONTAL GRIDPOINT. 389 50 CONTINUE 390 C 391 C WE NOW HAVE THE SHEULL SEA LEVEL PRESSURE FIELD 392 C IN ARRAY SLP. APPLY A FIVE POINT SMOOTHER. 393 C 394 !$omp parallel do 395 DO J=JSTA,JEND 396 DO I=1,IM 397 EGRID(I,J)=SLP(I,J) 398 ENDDO 399 ENDDO 400 C 401 CALL EXCH(SLP) 402 DO 90 KS=1,KSLPD 403 DO 80 J=JSTA_M2,JEND_M2 404 IEND=IM-1-MOD(J+1,2) 405 DO I=2,IEND 406 IF (FIS(I,J).GT.H1000) GO TO 70 407 IF (FIS(I+IHW(J),J-1).GT.H1000)GO TO 70 408 IF (FIS(I+IHE(J),J-1).GT.H1000)GO TO 70 409 IF (FIS(I+IHW(J),J+1).GT.H1000)GO TO 70 410 IF (FIS(I+IHE(J),J+1).GT.H1000)GO TO 70 411 IF (KS.GT.1) GO TO 80 412 70 EGRID(I,J)=D125* 413 1 (SLP(I+IHW(J),J-1)+SLP(I+IHE(J),J-1) 414 1 +SLP(I+IHW(J),J+1)+SLP(I+IHE(J),J+1) 415 2 +H4*SLP(I,J)) 416 ENDDO 417 80 CONTINUE 418 DO J=JSTA,JEND 419 DO I=1,IM 420 SLP(I,J) = EGRID(I,J) 421 ENDDO 422 ENDDO 423 90 CONTINUE 424 C 425 C END OF ROUTINE. 426 C 427 Page 6 Source Listing NGMSLP2 2014-12-17 20:48 NGMSLP2.f 428 RETURN 429 END ENTRY POINTS Name ngmslp2_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 10 Label 161 153 20 Label 162 156 30 Label 167 132,133,152 40 Label 244 240 50 Label 247 185,186,237,239 70 Label 270 264,265,266,267,268 80 Label 275 261,269 90 Label 281 260 A2 Param 41 R(4) 4 scalar A3 Param 41 R(4) 4 scalar A4 Param 41 R(4) 4 scalar ACCLIQ Scalar 7 R(4) 4 2 91719 COM ACPREC Scalar 7 R(4) 4 2 91719 COM 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 AKHS Scalar 3 R(4) 4 2 91719 COM AKMS Scalar 3 R(4) 4 2 91719 COM ALOG Func 190 scalar 190,230 ALPAVG Local 230 R(4) 4 scalar 230,231 ALPBAR Local 111 R(4) 4 2 91719 151,164,207 ALPINT Scalar 4 R(4) 4 3 4219074 COM 146,154 ALPM Local 154 R(4) 4 scalar 154,155,158 ALPS Local 146 R(4) 4 scalar 146,151,158,164 ALPSFC Local 190 R(4) 4 scalar 190,207 ALSL Scalar 7 R(4) 4 1 39 COM AMAX1 Func 140 scalar 140 BETA Param 37 R(4) 4 scalar 37 BOUNDL Subr 176 176,177 BTG Param 37 R(4) 4 scalar CAPA Param 28 R(4) 4 scalar 35,142 CFT0 Param 33 R(4) 4 scalar CLDEFI Scalar 6 R(4) 4 2 91719 COM 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 CONST Param 106 R(4) 4 scalar 220 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 Page 7 Source Listing NGMSLP2 2014-12-17 20:48 Symbol Table NGMSLP2.f Name Object Declared Type Bytes Dimen Elements Attributes References CPHI0 Scalar 9 R(4) 4 scalar COM CTLBLK Common 2 96 CUPREC Scalar 7 R(4) 4 2 91719 COM CURV Scalar 8 R(4) 4 2 91719 COM D00 Param 25 R(4) 4 scalar D001 Param 25 R(4) 4 scalar D0065 Param 29 R(4) 4 scalar D01 Param 25 R(4) 4 scalar 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 270 D25 Param 27 R(4) 4 scalar D50 Param 25 R(4) 4 scalar 146,154,202,205,224,230 D608 Param 29 R(4) 4 scalar 196,206,243 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 241 DISLP Scalar 7 R(4) 4 scalar COM DLMD Scalar 6 R(4) 4 scalar COM DP Param 107 R(4) 4 scalar 138,191 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 DY Scalar 3 R(4) 4 scalar COM DYNAM Common 2 3675964 EF4T Scalar 3 R(4) 4 scalar COM EGRID Local 112 R(4) 4 2 91719 255,270,278 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 EXCH Subr 259 259 EXP Func 155 scalar 155,227,231 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 Page 8 Source Listing NGMSLP2 2014-12-17 20:48 Symbol Table NGMSLP2.f Name Object Declared Type Bytes Dimen Elements Attributes References FIRST Scalar 3 R(4) 4 scalar COM FIS Scalar 3 R(4) 4 2 91719 COM 188,264,265,266,267,268 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,106,107,233 GAMMA Param 105 R(4) 4 scalar 211,213,242 GI Param 32 R(4) 4 scalar 36,106,188,197,207,212,214,232,241 GORD Param 107 R(4) 4 scalar H1 Param 26 R(4) 4 scalar 33,196,206,233,243 H100 Param 27 R(4) 4 scalar H1000 Param 27 R(4) 4 scalar 264,265,266,267,268 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 140,176,177 H1M2 Param 30 R(4) 4 scalar H1M5 Param 27 R(4) 4 scalar H2 Param 26 R(4) 4 scalar H4 Param 26 R(4) 4 scalar 273 H999 Param 25 R(4) 4 scalar H99999 Param 26 R(4) 4 scalar 176,177 HBM2 Scalar 3 R(4) 4 2 91719 COM 172,173,174,175 HTM Scalar 5 R(4) 4 3 4127355 COM I Local 133 I(4) 4 scalar 133,134,137,139,140,141,142,146,14 7,148,149,150,151,154,159,160,164, 165,166,186,187,188,189,192,196,20 2,205,207,227,230,234,243,254,255, 263,264,265,266,267,268,270,271,27 2,273,277,278 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 IEND Local 262 I(4) 4 scalar 262,263 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 266,268,271,272 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 265,267,271,272 IM Param 3 I(4) 4 scalar 2,3,4,5,6,7,8,9,11,14,111,112,133, 176,177,186,254,262,277 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 IOUT Scalar 3 I(4) 4 scalar COM ISHDE Scalar 4 I(4) 4 1 61 COM Page 9 Source Listing NGMSLP2 2014-12-17 20:48 Symbol Table NGMSLP2.f Name Object Declared Type Bytes Dimen Elements Attributes References ISMTHQ Param 108 I(4) 4 scalar 174,175 ISMTHT Param 108 I(4) 4 scalar 172,173 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 IXM Scalar 3 I(4) 4 scalar COM IYM Scalar 3 I(4) 4 scalar COM J Local 132 I(4) 4 scalar 132,134,137,139,140,141,142,146,14 7,148,149,150,151,154,159,160,164, 165,166,185,187,188,189,192,196,20 2,205,207,227,230,234,243,253,255, 261,262,264,265,266,267,268,270,27 1,272,273,276,278 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 132,185,253,276 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 261 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,111,112,176, 177 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 132,185,253,276 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 261 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 KS Local 260 I(4) 4 scalar 260,269 KSAV Param 14 I(4) 4 scalar 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 260 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 153 I(4) 4 scalar 153,154,159,160,240,241,243 LB Param 9 I(4) 4 scalar LHEAT Param 40 I(4) 4 scalar 40 LHTOCP Param 40 I(4) 4 scalar Page 10 Source Listing NGMSLP2 2014-12-17 20:48 Symbol Table NGMSLP2.f Name Object Declared Type Bytes Dimen Elements Attributes References LIST Scalar 3 I(4) 4 scalar COM LLMH Local 134 I(4) 4 scalar 134,146,147,148,152,153,187,227,23 9,240 LM Param 3 I(4) 4 scalar 4,5,8,9,19,152,227,239,240 LM1 Param 8 I(4) 4 scalar LMH Scalar 5 I(4) 4 2 91719 COM 134,187 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 MAPOT Common 2 1448 MASKS Common 2 34853220 ME Scalar 6 I(4) 4 scalar COM MOD Func 262 scalar 262 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 NGMSLP2 Subr 1 NINC Param 7 I(4) 4 scalar NLEV Local 145 I(4) 4 scalar 145,157,163 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 10 R(4) 4 scalar COM ONEPS Param 39 R(4) 4 scalar OVERRC Param 33 R(4) 4 scalar 33 P1000 Param 35 R(4) 4 scalar 142,234 P2FILT Subr 172 172,173,174,175 PARALLEL Common 9 8240 PAVG Local 231 R(4) 4 scalar 231,232 PD Scalar 3 R(4) 4 2 91719 COM 137,189 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 155 R(4) 4 scalar 155,156 PQ0 Param 41 R(4) 4 scalar PREC Scalar 7 R(4) 4 2 91719 COM PSFC Local 137 R(4) 4 scalar 137,138,142,189,190,191,192,227,23 Page 11 Source Listing NGMSLP2 2014-12-17 20:48 Symbol Table NGMSLP2.f Name Object Declared Type Bytes Dimen Elements Attributes References 0 PSLP Scalar 6 R(4) 4 2 91719 COM PT Scalar 5 R(4) 4 scalar COM 137,189 PTDYN Scalar 3 R(4) 4 scalar COM PTOP Local 138 R(4) 4 scalar 138,156,191 PTSL Scalar 7 R(4) 4 2 91719 COM PVRBLS Common 2 23480064 Q Scalar 5 R(4) 4 3 4127355 COM 148,160 Q2 Scalar 9 R(4) 4 3 4127355 COM QAVG Local 202 R(4) 4 scalar 202,206 QBAR Local 111 R(4) 4 2 91719 150,166,174,176,196,202,243 QS Scalar 3 R(4) 4 2 91719 COM 139 QSFC Local 112 R(4) 4 2 91719 139,140,175,177,202 QSL Scalar 9 R(4) 4 2 91719 COM QSUM Local 148 R(4) 4 scalar 148,150,160,166 QWBS Scalar 5 R(4) 4 2 91719 COM QZ0 Scalar 4 R(4) 4 2 91719 COM 139 R Scalar 5 R(4) 4 scalar COM RCAPA Param 35 R(4) 4 scalar RD Param 32 R(4) 4 scalar 32,36,40,106,107,197,207,212,214 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 RF Scalar 5 R(4) 4 2 91719 COM RG Param 36 R(4) 4 scalar RHOAVG Local 232 R(4) 4 scalar 232,233 RNLEV Local 163 R(4) 4 scalar 163,164,165,166 ROG Param 32 R(4) 4 scalar RRHOG Local 233 R(4) 4 scalar 233,234 RTD Param 34 R(4) 4 scalar RUN Scalar 3 R(4) 4 scalar COM SBD Scalar 6 R(4) 4 scalar COM SI Scalar 6 R(4) 4 2 91719 COM SICE Scalar 3 R(4) 4 2 91719 COM SIGMA Scalar 3 R(4) 4 scalar COM SLP Scalar 5 R(4) 4 2 91719 COM 192,227,230,234,255,259,271,272,27 3,278 SM Scalar 3 R(4) 4 2 91719 COM 139,141 SMTHA Param 15 R(4) 4 scalar SMTHB Param 15 R(4) 4 scalar SNO Scalar 6 R(4) 4 2 91719 COM SPHI0 Scalar 9 R(4) 4 scalar COM SPL Scalar 7 R(4) 4 1 39 COM STDRD Scalar 10 R(4) 4 scalar COM T Scalar 5 R(4) 4 3 4127355 COM 147,159,243 T500 Scalar 6 R(4) 4 2 91719 COM TAU Local 197 R(4) 4 scalar 197 TAUAVG Local 224 R(4) 4 scalar 224,227,232 TAUCR Param 106 R(4) 4 scalar 217,218,219,220 TAUSFC Local 212 R(4) 4 scalar 212,217,219,220,224 TAUSL Local 214 R(4) 4 scalar 214,217,218,219,220,224 TAVG Local 205 R(4) 4 scalar 205,206 TBAR Local 111 R(4) 4 2 91719 149,165,172,196,205 Page 12 Source Listing NGMSLP2 2014-12-17 20:48 Symbol Table NGMSLP2.f Name Object Declared Type Bytes Dimen Elements Attributes References TFRZ Param 39 R(4) 4 scalar THS Scalar 3 R(4) 4 2 91719 COM 141 THSFC Local 141 R(4) 4 scalar 141,142 THZ0 Scalar 4 R(4) 4 2 91719 COM 141 TLM0D Scalar 6 R(4) 4 scalar COM TPH0D Scalar 6 R(4) 4 scalar COM TSFC Local 112 R(4) 4 2 91719 142,173,205 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 TSUM Local 147 R(4) 4 scalar 147,149,159,165 TTV Scalar 8 R(4) 4 2 91719 COM TVRBAR Local 206 R(4) 4 scalar 206,207 TVRSFC Local 211 R(4) 4 scalar 211,212 TVRSL Local 213 R(4) 4 scalar 213,214 TVRT Local 196 R(4) 4 scalar 196,197,211,213,242 TVRTL Local 242 R(4) 4 scalar 242,243 TWBS Scalar 5 R(4) 4 2 91719 COM U Scalar 5 R(4) 4 3 4127355 COM UZ0 Scalar 4 R(4) 4 2 91719 COM V Scalar 5 R(4) 4 3 4127355 COM 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 VZ0 Scalar 4 R(4) 4 2 91719 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 Z0 Scalar 3 R(4) 4 2 91719 COM Z0SLP Scalar 7 R(4) 4 scalar COM Z1000 Scalar 9 R(4) 4 2 91719 COM 234 ZBAR Local 207 R(4) 4 scalar 207,211,213,242 ZINT Scalar 4 R(4) 4 3 4219074 COM ZL Local 241 R(4) 4 scalar 241,242 ZSFC Local 188 R(4) 4 scalar 188,207,211,227,237 ZSL Param 105 R(4) 4 scalar 213,237 Page 13 Source Listing NGMSLP2 2014-12-17 20:48 Subprograms/Common Blocks NGMSLP2.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References 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 NGMSLP2 Subr 1 PARALLEL Common 9 8240 PVRBLS Common 2 23480064 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__ -D __pentium4 -D __pentium4__ Page 14 Source Listing NGMSLP2 2014-12-17 20:48 NGMSLP2.f -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 : NGMSLP2.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100