Page 1 Source Listing SFCDIF 2014-12-17 20:47 /tmpnwprd/ifortEl06Ya.i 1 # 1 "SFCDIF.F" 2 C&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 3 SUBROUTINE SFCDIF(LMHK,SM,THS,QS,UZ0,VZ0,THZ0,QZ0 4 &, USTAR,WSTAR,Z0,ZEFF,AKMS,AKHS,HPBL,CT 5 &, U10,V10,TH02,TH10,Q02,Q10 6 &, TH30,Q30,U30,V30 7 &, ULM,VLM,T,Q,APE,Z,PD,PT,TLM) 8 C ****************************************************************** 9 C * * 10 C * SURFACE LAYER * 11 C * * 12 C ****************************************************************** 13 C Ammended to use the "Effective roughness" of Mason (1986, see 14 C Georgelin et al., MWR July 1994), by FM, RW, June 1995 15 C----------------------------------------------------------------------- 16 INCLUDE "parmeta" 39 INCLUDE "mpp.h" 40 # 1 "./sp.h" 1 41 # 4 42 43 # 18 "SFCDIF.F" 2 44 C----------------------------------------------------------------------- 103 P A R A M E T E R 104 &(LP1=LM+1) 105 C----------------------------------------------------------------------- 106 P A R A M E T E R 107 &(WWST=1.2,WWST2=WWST*WWST,G=9.8,USTFC=0.018/G 108 &,VKRM=0.40,RIC=0.183,RFC=0.191,FHNEU=0.8 109 &,RRIC=1.0/RIC,RFAC=RIC/(FHNEU*RFC*RFC),EXCM=0.001 110 &,BETA=1./270.,BTG=BETA*G 111 &,ELFC=VKRM*BTG,CNV=0.608*G/BTG 112 &,WOLD=.15,WNEW=1.-WOLD,ITRMX=05 113 &,PIHF=3.14159265/2.,PIFR=3.14159265/4. 114 C----------------------------------------------------------------------- 115 &,EPSU2=1.E-4,EPSUST=0.07,EPSIT=1.E-4,EPSA=1.E-8 116 &,ZTMIN=-5.,ZTMAX=1. 117 C----------------------------------------------------------------------- 118 &,SMALL=0.35, GLKBS=30.0,GLKBR=10.0,GRRS=GLKBR/GLKBS 119 &,CZIV=SMALL*GLKBS 120 &,VISC=1.5E-5, TVISC=2.1E-5, QVISC=2.1E-5 121 &,RVISC=1./VISC,RTVISC=1./TVISC,RQVISC=1./QVISC 122 &,SQPR=0.84,SQSC=0.84,ZQRZT=SQSC/SQPR 123 &,USTR=0.225,USTC=0.7 124 &,FZU1=CZIV*VISC,FZT1=RVISC *TVISC*SQPR, FZQ1=RTVISC*QVISC*ZQRZT 125 &, FZT2=CZIV*GRRS*TVISC*SQPR,FZQ2=RTVISC*QVISC*ZQRZT 126 C----------------------------------------------------------------------- 127 &,ZTFC=1.0 128 c &,CZIL=.1000,SQVISC=258.2,ZILFC=-CZIL*VKRM*SQVISC 129 &,CZIL=.2000,SQVISC=258.2,ZILFC=-CZIL*VKRM*SQVISC 130 &,PQ0=379.90516,A2=17.2693882,A3=273.16,A4=35.86 131 &,CAPA=0.28589641E0,H1M5=1.E-5) 132 C----------------------------------------------------------------------- 133 D I M E N S I O N 134 & T (LM),Q (LM) 135 D I M E N S I O N 136 & APE (LM) 137 &,Z (LP1) Page 2 Source Listing SFCDIF 2014-12-17 20:47 SFCDIF.F 138 C-ZEFF-ZEFF-ZEFF-ZEFF 139 &,ZEFF (4) 140 C-ZEFF-ZEFF-ZEFF-ZEFF 141 C----------------------------------------------------------------------- 142 PSLMU(ZZ)=-0.96*ALOG(1.0-4.5*ZZ) 143 PSLMS(ZZ)=ZZ*RRIC-2.076*(1.-1./(ZZ+1.)) 144 PSLHU(ZZ)=-0.96*ALOG(1.0-4.5*ZZ) 145 PSLHS(ZZ)=ZZ*RFAC-2.076*(1.-1./(ZZ+1.)) 146 C 147 PSPMU(XX)=-2.*ALOG((XX+1.)*0.5)-ALOG((XX*XX+1.)*0.5)+2.*ATAN(XX) 148 & -PIHF 149 PSPMS(YY)=5.*YY 150 PSPHU(XX)=-2.*ALOG((XX*XX+1.)*0.5) 151 PSPHS(YY)=5.*YY 152 C*********************************************************************** 153 LMHP=LMHK+1 154 C 155 THLM=T(LMHK)*APE(LMHK) 156 QLM=Q(LMHK) 157 C----------------------------------------------------------------------- 158 Z0=(1.-SM)*Z0+SM*AMAX1(USTFC*USTAR*USTAR,1.59E-5) 159 C--------------VISCOUS SUBLAYER----------------------------------------- 160 IF(SM.GT.0.5.AND.USTAR.LT.USTC)THEN 161 C----------------------------------------------------------------------- 162 IF(USTAR.LT.USTR)THEN 163 C 164 ZU=FZU1*SQRT(SQRT(Z0*USTAR*RVISC))/USTAR 165 WGHT=AKMS*ZU*RVISC 166 RWGH=WGHT/(WGHT+1.) 167 UZ0=(ULM*RWGH+UZ0)*0.5 168 VZ0=(VLM*RWGH+VZ0)*0.5 169 C 170 ZT=FZT1*ZU 171 WGHT=AKHS*ZT*RTVISC 172 THZ0=((WGHT*THLM+THS)/(WGHT+1.)+THZ0)*0.5 173 C 174 ZQ=FZQ1*ZT 175 WGHT=AKHS*ZQ*RQVISC 176 QZ0 =((WGHT*QLM+QS)/(WGHT+1.)+QZ0)*0.5 177 C 178 ENDIF 179 IF(USTAR.GE.USTR.AND.USTAR.LT.USTC)THEN 180 C 181 ZU=Z0 182 UZ0=0. 183 VZ0=0. 184 C 185 ZT=FZT2*SQRT(SQRT(Z0*USTAR*RVISC))/USTAR 186 WGHT=AKHS*ZT*RTVISC 187 THZ0=((WGHT*THLM+THS)/(WGHT+1.)+THZ0)*0.5 188 C 189 ZQ=FZQ2*ZT 190 WGHT=AKHS*ZQ*RQVISC 191 QZ0 =((WGHT*QLM+QS)/(WGHT+1.)+QZ0)*0.5 192 ENDIF 193 C----------------------------------------------------------------------- 194 ELSE Page 3 Source Listing SFCDIF 2014-12-17 20:47 SFCDIF.F 195 C----------------------------------------------------------------------- 196 ZU=Z0 197 C-ZEFF-ZEFF-ZEFF-ZEFF 198 IF(SM.LE.0.5)THEN 199 IF(ULM.EQ.0.) ULM=EPSU2 200 ALPHA=ABS(ATAN(VLM/ULM)+PIHF-EPSA) 201 X=ALPHA/PIFR 202 ML=1+X 203 C 2012-11 WNE 204 if (ml.lt.1.or.ml.gt.4) then 205 write (*,*)'>>ml problem ml,ulm,vlm=',ml,ulm,vlm 206 ml = 1 207 endif 208 MH=1+MOD(ML,4) 209 210 WLOW=X-ML+1 211 ZU=WLOW*ZEFF(ML)+(1.-WLOW)*ZEFF(MH) 212 ENDIF 213 C-ZEFF-ZEFF-ZEFF-ZEFF 214 UZ0=0. 215 VZ0=0. 216 C 217 ZT=Z0 218 THZ0=THS 219 C 220 ZQ=Z0 221 QZ0=QS 222 C----------------------------------------------------------------------- 223 ENDIF 224 C----------------------------------------------------------------------- 225 ZSL=(Z(LMHK)-Z(LMHP))*0.5 226 C-ZEFF-ZEFF-ZEFF-ZEFF 227 ZU=AMIN1(ZU,0.5*ZSL) 228 C-ZEFF-ZEFF-ZEFF-ZEFF 229 RDZ=1./ZSL 230 CXCH=EXCM*RDZ 231 C----------------------------------------------------------------------- 232 IF(SM.GT.0.5)THEN 233 DTHV=(0.608*QLM+1.)*THLM-(0.608*QZ0+1.)*THZ0 234 ELSE 235 DTHV=(QLM-QZ0)*CNV+THLM-THZ0 236 ZT=Z0*ZTFC 237 ENDIF 238 C 239 DU2=AMAX1((ULM-UZ0)**2+(VLM-VZ0)**2,EPSU2) 240 C----------------------------------------------------------------------- 241 RIB=BTG*DTHV*ZSL/DU2 242 C--------------BELJARS CORRECTION OF USTAR------------------------------ 243 BTGH=BTG*HPBL 244 WSTAR2=WWST2*ABS(BTGH*AKHS*DTHV)**(2./3.) 245 USTAR=AMAX1(SQRT(AKMS*SQRT(DU2+WSTAR2)),EPSUST) 246 C--------------ZILITINKEVITCH FIX FOR ZT-------------------------------- 247 IF(SM.LT.0.5)ZT=EXP(ZILFC*SQRT(USTAR*Z0))*Z0 248 C----------------------------------------------------------------------- 249 IF(SM.GT.0.5.AND.RIB.GE.RIC)THEN 250 C----------------------------------------------------------------------- 251 AKMS=AMAX1( VISC*RDZ,CXCH) Page 4 Source Listing SFCDIF 2014-12-17 20:47 SFCDIF.F 252 AKHS=AMAX1(TVISC*RDZ,CXCH) 253 C----------------------------------------------------------------------- 254 ELSE 255 C----------------------------------------------------------------------- 256 ZSLU=ZSL+ZU 257 ZSLT=ZSL+ZT 258 C 259 RLOGU=ALOG(ZSLU/ZU) 260 RLOGT=ALOG(ZSLT/ZT) 261 C 262 RLMO=ELFC*AKHS*DTHV/USTAR**3 263 C--------------SEA POINTS FIRST ... ------------------------------------ 264 IF(SM.GT.0.5)THEN 265 DO 100 ITR=1,ITRMX 266 C--------------1./MONIN-OBUKKHOV LENGTH-SCALE--------------------------- 267 ZETALT=AMAX1(ZSLT*RLMO,ZTMIN) 268 RLMO=ZETALT/ZSLT 269 ZETALU=ZSLU*RLMO 270 C 271 ZETAU=ZU*RLMO 272 ZETAT=ZT*RLMO 273 C--------------LL FUNCTIONS OVER SEA------------------------------------ 274 IF(RLMO.LT.0.)THEN 275 PSMZ=PSLMU(ZETAU) 276 SIMM= PSLMU(ZETALU)-PSMZ+RLOGU 277 PSHZ=PSLHU(ZETAT) 278 SIMH=FHNEU*(PSLHU(ZETALT)-PSHZ+RLOGT) 279 ELSE 280 PSMZ=PSLMS(ZETAU) 281 SIMM= PSLMS(ZETALU)-PSMZ+RLOGU 282 PSHZ=PSLHS(ZETAT) 283 SIMH=FHNEU*(PSLHS(ZETALT)-PSHZ+RLOGT) 284 ENDIF 285 C--------------BELJAARS CORRECTION FOR USTAR---------------------------- 286 USTAR=AMAX1(SQRT(AKMS*SQRT(DU2+WSTAR2)),EPSUST) 287 C----------------------------------------------------------------------- 288 USTARK=USTAR*VKRM 289 AKMS=AMAX1(USTARK/SIMM,CXCH) 290 AKHS=AMAX1(USTARK/SIMH,CXCH) 291 C----------------------------------------------------------------------- 292 WSTAR2=WWST2*ABS(BTGH*AKHS*DTHV)**(2./3.) 293 RLMN=ELFC*AKHS*DTHV/USTAR**3 294 C----------------------------------------------------------------------- 295 RLMP=RLMO 296 RLMA=RLMO*WOLD+RLMN*WNEW 297 C----------------------------------------------------------------------- 298 C IF(ABS((RLMN-RLMO)/RLMA).LT.EPSIT) GO TO 110 299 C----------------------------------------------------------------------- 300 RLMO=RLMA 301 C----------------------------------------------------------------------- 302 100 CONTINUE 303 C----------------------------------------------------------------------- 304 110 CONTINUE 305 C--------------END OF SEA POINT PROCESSING------------------------------ 306 ELSE 307 C--------------NOW LAND POINTS ...-------------------------------------- 308 DO 200 ITR=1,ITRMX Page 5 Source Listing SFCDIF 2014-12-17 20:47 SFCDIF.F 309 C--------------1./MONIN-OBUKKHOV LENGTH-SCALE--------------------------- 310 ZETALT=AMAX1(ZSLT*RLMO,ZTMIN) 311 RLMO=ZETALT/ZSLT 312 ZETALU=ZSLU*RLMO 313 C 314 ZETAU=ZU*RLMO 315 ZETAT=ZT*RLMO 316 C--------------PAULSON 1970 FUNCTIONS OVER LAND W RAD. SKIN T----------- 317 IF(RLMO.LT.0.)THEN 318 XLU4=1.-16.*ZETALU 319 XLT4=1.-16.*ZETALT 320 XU4 =1.-16.*ZETAU 321 XT4 =1.-16.*ZETAT 322 C 323 XLU=SQRT(SQRT(XLU4)) 324 XLT=SQRT(SQRT(XLT4)) 325 XU =SQRT(SQRT(XU4)) 326 XT =SQRT(SQRT(XT4)) 327 C 328 PSMZ=PSPMU(XU) 329 SIMM=PSPMU(XLU)-PSMZ+RLOGU 330 PSHZ=PSPHU(XT) 331 SIMH=PSPHU(XLT)-PSHZ+RLOGT 332 ELSE 333 ZETAU=AMIN1(ZETAU,ZTMAX) 334 ZETAT=AMIN1(ZETAT,ZTMAX) 335 ZETALU=AMIN1(ZETALU,ZTMAX) 336 ZETALT=AMIN1(ZETALT,ZTMAX) 337 PSMZ=PSPMS(ZETAU) 338 SIMM=PSPMS(ZETALU)-PSMZ+RLOGU 339 PSHZ=PSPHS(ZETAT) 340 SIMH=PSPHS(ZETALT)-PSHZ+RLOGT 341 ENDIF 342 C--------------BELJAARS CORRECTION FOR USTAR---------------------------- 343 USTAR=AMAX1(SQRT(AKMS*SQRT(DU2+WSTAR2)),EPSUST) 344 C--------------ZILITINKEVITCH FIX FOR ZT-------------------------------- 345 ZT=EXP(ZILFC*SQRT(USTAR*Z0))*Z0 346 ZSLT=ZSL+ZT 347 RLOGT=ALOG(ZSLT/ZT) 348 C----------------------------------------------------------------------- 349 USTARK=USTAR*VKRM 350 AKMS=AMAX1(USTARK/SIMM,CXCH) 351 AKHS=AMAX1(USTARK/SIMH,CXCH) 352 C----------------------------------------------------------------------- 353 WSTAR2=WWST2*ABS(BTGH*AKHS*DTHV)**(2./3.) 354 RLMN=ELFC*AKHS*DTHV/USTAR**3 355 C----------------------------------------------------------------------- 356 RLMP=RLMO 357 RLMA=RLMO*WOLD+RLMN*WNEW 358 C----------------------------------------------------------------------- 359 C IF(ABS((RLMN-RLMO)/RLMA).LT.EPSIT) GO TO 210 360 C----------------------------------------------------------------------- 361 RLMO=RLMA 362 C----------------------------------------------------------------------- 363 200 CONTINUE 364 C----------------------------------------------------------------------- 365 210 CONTINUE Page 6 Source Listing SFCDIF 2014-12-17 20:47 SFCDIF.F 366 C--------------END OF LAND POINT PROCESSING AND SEA-LAND BRANCHING------ 367 ENDIF 368 C--------------END OF TURBULENCE-NO TURBULENCE BRANCHING---------------- 369 ENDIF 370 C--------------COUNTERGRADIENT FIX-------------------------------------- 371 C HV=-AKHS*DTHV 372 C IF(HV.GT.0.)THEN 373 C FCT=-10.*(BTG)**(-1./3.) 374 C CT=FCT*(HV/(HPBL*HPBL))**(2./3.) 375 C ELSE 376 CT=0. 377 C ENDIF 378 C--------------DIAGNOSTIC BLOCK----------------------------------------- 379 WSTAR=SQRT(WSTAR2)/WWST 380 C 381 UMFLX=AKMS*(ULM -UZ0 ) 382 VMFLX=AKMS*(VLM -VZ0 ) 383 HSFLX=AKHS*(THLM-THZ0) 384 HLFLX=AKHS*(QLM -QZ0 ) 385 C----------------------------------------------------------------------- 386 IF(SM.GT.0.5.AND.RIB.GE.RIC)THEN 387 C----------------------------------------------------------------------- 388 AKMS10=AMAX1( VISC/10.,CXCH) 389 AKMS30=AMAX1( VISC/30.,CXCH) 390 AKHS02=AMAX1(TVISC/02.,CXCH) 391 AKHS10=AMAX1(TVISC/10.,CXCH) 392 AKHS30=AMAX1(TVISC/30.,CXCH) 393 C----------------------------------------------------------------------- 394 ELSE 395 C----------------------------------------------------------------------- 396 ZU10=ZU+10. 397 ZU30=ZU+30. 398 ZT02=ZT+02. 399 ZT10=ZT+10. 400 ZT30=ZT+30. 401 C 402 RLNU10=ALOG(ZU10/ZU) 403 RLNU30=ALOG(ZU30/ZU) 404 RLNT02=ALOG(ZT02/ZT) 405 RLNT10=ALOG(ZT10/ZT) 406 RLNT30=ALOG(ZT30/ZT) 407 C 408 ZTAU10=ZU10*RLMP 409 ZTAU30=ZU30*RLMP 410 ZTAT02=ZT02*RLMP 411 ZTAT10=ZT10*RLMP 412 ZTAT30=ZT30*RLMP 413 C--------------LL FUNCTIONS OVER SEA------------------------------------ 414 IF(SM.GT.0.5)THEN 415 C----------------------------------------------------------------------- 416 IF(RLMP.LT.0.)THEN 417 SIMM10= PSLMU(ZTAU10)-PSMZ+RLNU10 418 SIMM30= PSLMU(ZTAU30)-PSMZ+RLNU30 419 SIMH02=FHNEU*(PSLHU(ZTAT02)-PSHZ+RLNT02) 420 SIMH10=FHNEU*(PSLHU(ZTAT10)-PSHZ+RLNT10) 421 SIMH30=FHNEU*(PSLHU(ZTAT30)-PSHZ+RLNT30) 422 ELSE Page 7 Source Listing SFCDIF 2014-12-17 20:47 SFCDIF.F 423 SIMM10= PSLMS(ZTAU10)-PSMZ+RLNU10 424 SIMM30= PSLMS(ZTAU30)-PSMZ+RLNU30 425 SIMH02=FHNEU*(PSLHS(ZTAT02)-PSHZ+RLNT02) 426 SIMH10=FHNEU*(PSLHS(ZTAT10)-PSHZ+RLNT10) 427 SIMH30=FHNEU*(PSLHS(ZTAT30)-PSHZ+RLNT30) 428 ENDIF 429 C--------------PAULSON 1970 FUNCTIONS OVER LAND W RAD. SKIN T----------- 430 ELSE 431 C----------------------------------------------------------------------- 432 IF(RLMP.LT.0.)THEN 433 XLU104=1.-16.*ZTAU10 434 XLU304=1.-16.*ZTAU30 435 XLT024=1.-16.*ZTAT02 436 XLT104=1.-16.*ZTAT10 437 XLT304=1.-16.*ZTAT30 438 C 439 XLU10=SQRT(SQRT(XLU104)) 440 XLU30=SQRT(SQRT(XLU304)) 441 XLT02=SQRT(SQRT(XLT024)) 442 XLT10=SQRT(SQRT(XLT104)) 443 XLT30=SQRT(SQRT(XLT304)) 444 C 445 SIMM10=PSPMU(XLU10)-PSMZ+RLNU10 446 SIMM30=PSPMU(XLU30)-PSMZ+RLNU30 447 SIMH02=PSPHU(XLT02)-PSHZ+RLNT02 448 SIMH10=PSPHU(XLT10)-PSHZ+RLNT10 449 SIMH30=PSPHU(XLT30)-PSHZ+RLNT30 450 ELSE 451 ZTAU10=AMIN1(ZTAU10,ZTMAX) 452 ZTAU30=AMIN1(ZTAU30,ZTMAX) 453 ZTAT02=AMIN1(ZTAT02,ZTMAX) 454 ZTAT10=AMIN1(ZTAT10,ZTMAX) 455 ZTAT30=AMIN1(ZTAT30,ZTMAX) 456 C 457 SIMM10=PSPMS(ZTAU10)-PSMZ+RLNU10 458 SIMM30=PSPMS(ZTAU30)-PSMZ+RLNU30 459 SIMH02=PSPHS(ZTAT02)-PSHZ+RLNT02 460 SIMH10=PSPHS(ZTAT10)-PSHZ+RLNT10 461 SIMH30=PSPHS(ZTAT30)-PSHZ+RLNT30 462 ENDIF 463 C----------------------------------------------------------------------- 464 ENDIF 465 C----------------------------------------------------------------------- 466 AKMS10=AMAX1(USTARK/SIMM10,CXCH) 467 AKMS30=AMAX1(USTARK/SIMM30,CXCH) 468 AKHS02=AMAX1(USTARK/SIMH02,CXCH) 469 AKHS10=AMAX1(USTARK/SIMH10,CXCH) 470 AKHS30=AMAX1(USTARK/SIMH30,CXCH) 471 C----------------------------------------------------------------------- 472 ENDIF 473 C----------------------------------------------------------------------- 474 U10 =UMFLX/AKMS10+UZ0 475 V10 =VMFLX/AKMS10+VZ0 476 U30 =UMFLX/AKMS30+UZ0 477 V30 =VMFLX/AKMS30+VZ0 478 TH02=HSFLX/AKHS02+THZ0 479 TH10=HSFLX/AKHS10+THZ0 Page 8 Source Listing SFCDIF 2014-12-17 20:47 SFCDIF.F 480 TH30=HSFLX/AKHS30+THZ0 481 482 C GSM changed this section in response to problem with 2-m 483 C dew point occasionally being greater than 2-m temperature 484 C and similar problem at 10-m. Now, a saturation Q is 485 C calculated at each level, and the Q is constrained to 486 C be no higher than the saturation value. 487 488 PDS=PD+PT 489 TERM1=-0.068283/TLM 490 PSHLTR=PDS*EXP(TERM1) 491 T02=TH02*(PSHLTR*H1M5)**CAPA 492 QSAT2 = PQ0/PSHLTR*EXP(A2*(T02-A3)/(T02-A4)) 493 Q02 =HLFLX/AKHS02+QZ0 494 IF (Q02.LT.0.) THEN 495 IF (QLM .GT. 0.) THEN 496 Q02=QLM 497 ELSE 498 Q02=0.0001 499 ENDIF 500 ENDIF 501 IF (Q02.GT.QSAT2)THEN 502 Q02 = QSAT2 503 ENDIF 504 505 T10=TH10*(PSHLTR*H1M5)**CAPA 506 QSAT10 = PQ0/PSHLTR*EXP(A2*(T10-A3)/(T10-A4)) 507 C Corrected to use pressure at 10m. (DJ) 508 P10=PDS*EXP(-10.0*G/(287.04*TLM)) 509 T10=TH10*(P10*H1M5)**CAPA 510 QSAT10 = PQ0/P10*EXP(A2*(T10-A3)/(T10-A4)) 511 C 512 Q10 =HLFLX/AKHS10+QZ0 513 IF (Q10.LT.0.) THEN 514 IF (QLM .GT. 0.) THEN 515 Q10=QLM 516 ELSE 517 Q10=0.0001 518 ENDIF 519 ENDIF 520 IF (Q10.GT.QSAT10)THEN 521 Q10 = QSAT10 522 ENDIF 523 524 P30=PDS*EXP(-30.0*G/(287.04*TLM)) 525 T30=TH30*(P30*H1M5)**CAPA 526 QSAT30 = PQ0/P30*EXP(A2*(T30-A3)/(T30-A4)) 527 Q30 =HLFLX/AKHS30+QZ0 528 IF (Q30.LT.0.) THEN 529 IF (QLM .GT. 0.) THEN 530 Q30=QLM 531 ELSE 532 Q30=0.0001 533 ENDIF 534 ENDIF 535 IF (Q30.GT.QSAT30)THEN 536 Q30 = QSAT30 Page 9 Source Listing SFCDIF 2014-12-17 20:47 SFCDIF.F 537 ENDIF 538 539 c new calculation of 10-m winds 540 C----------------------------------------------------- 541 U10E=U10 542 V10E=V10 543 U30E=U30 544 V30E=V30 545 C----------------------------------------------------- 546 IF(SM.LT.0.5) THEN 547 c choose the equivalent z0 here: 548 czj ZU=0.01 549 ckm zu=zu*0.1 550 C 551 zuuz=amin1(zu*0.50,0.10) 552 zu=amax1(zu*0.10,zuuz) 553 ZU10=ZU+10. 554 ZU30=ZU+30. 555 RLNU10=ALOG(ZU10/ZU) 556 RLNU30=ALOG(ZU30/ZU) 557 ZTAU=ZU*RLMP 558 ZTAU10=ZU10*RLMP 559 ZTAU30=ZU30*RLMP 560 c-------------------------------------------------------- 561 IF(RLMP.LT.0)THEN 562 XLU104=1.-16.*ZTAU10 563 XLU304=1.-16.*ZTAU30 564 XU104 =1.-16.*ZTAU 565 XLU10=SQRT(SQRT(XLU104)) 566 XLU30=SQRT(SQRT(XLU304)) 567 XU10 =SQRT(SQRT(XU104)) 568 SIMM10=PSPMU(XLU10)-PSPMU(XU10)+RLNU10 569 SIMM30=PSPMU(XLU30)-PSPMU(XU10)+RLNU30 570 ELSE 571 ZTAU10=AMIN1(ZTAU10,ZTMAX) 572 ZTAU30=AMIN1(ZTAU30,ZTMAX) 573 SIMM10=PSPMS(ZTAU10)-PSPMS(ZTAU)+RLNU10 574 SIMM30=PSPMS(ZTAU30)-PSPMS(ZTAU)+RLNU30 575 ENDIF 576 c----------------------------------------------------- 577 EKMS10=AMAX1(USTARK/SIMM10,CXCH) 578 EKMS30=AMAX1(USTARK/SIMM30,CXCH) 579 U10E=UMFLX/EKMS10+UZ0 580 V10E=VMFLX/EKMS10+VZ0 581 U30E=UMFLX/EKMS30+UZ0 582 V30E=VMFLX/EKMS30+VZ0 583 ENDIF 584 c---------------------------------------------------------------- 585 U10=U10E 586 V10=V10E 587 U30=U30E 588 V30=V30E 589 C 590 C----------------------------------------------------------------------- 591 RETURN 592 END Page 10 Source Listing SFCDIF 2014-12-17 20:47 Entry Points SFCDIF.F ENTRY POINTS Name sfcdif_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 100 Label 218 181 110 Label 220 200 Label 279 224 210 Label 281 A2 Param 46 R(4) 4 scalar 408,422,426,442 A3 Param 46 R(4) 4 scalar 408,422,426,442 A4 Param 46 R(4) 4 scalar 408,422,426,442 ABS Func 116 scalar 116,160,208,269 AKHS Dummy 3 R(4) 4 scalar ARG,INOUT 87,91,102,106,160,168,178,206,208, 209,267,269,270,299,300 AKHS02 Local 306 R(4) 4 scalar 306,384,394,409 AKHS10 Local 307 R(4) 4 scalar 307,385,395,428 AKHS30 Local 308 R(4) 4 scalar 308,386,396,443 AKMS Dummy 3 R(4) 4 scalar ARG,INOUT 81,161,167,202,205,259,266,297,298 AKMS10 Local 304 R(4) 4 scalar 304,382,390,391 AKMS30 Local 305 R(4) 4 scalar 305,383,392,393 ALOG Func 58 scalar 58,60,63,66,175,176,263,318,319,32 0,321,322,471,472 ALPHA Local 116 R(4) 4 scalar 116,117 AMAX1 Func 74 scalar 74,155,161,167,168,183,202,205,206 ,226,259,266,267,304,305,306,307,3 08,382,383,384,385,386,468,493,494 AMIN1 Func 143 scalar 143,249,250,251,252,367,368,369,37 0,371,467,487,488 APE Dummy 6 R(4) 4 1 45 ARG,INOUT 71 ATAN Func 63 scalar 63,116 BETA Param 26 R(4) 4 scalar 26 BTG Param 26 R(4) 4 scalar 27,157,159 BTGH Local 159 R(4) 4 scalar 159,160,208,269 CAPA Param 47 R(4) 4 scalar 407,421,425,441 CNV Param 27 R(4) 4 scalar 151 CT Dummy 3 R(4) 4 scalar ARG,INOUT 292 CXCH Local 146 R(4) 4 scalar 146,167,168,205,206,266,267,304,30 5,306,307,308,382,383,384,385,386, 493,494 CZIL Param 45 R(4) 4 scalar 45 CZIV Param 35 R(4) 4 scalar 40,41 DTHV Local 149 R(4) 4 scalar 149,151,157,160,178,208,209,269,27 0 DU2 Local 155 R(4) 4 scalar 155,157,161,202,259 EKMS10 Local 493 R(4) 4 scalar 493,495,496 EKMS30 Local 494 R(4) 4 scalar 494,497,498 ELFC Param 27 R(4) 4 scalar 178,209,270 EPSA Param 31 R(4) 4 scalar 116 Page 11 Source Listing SFCDIF 2014-12-17 20:47 Symbol Table SFCDIF.F Name Object Declared Type Bytes Dimen Elements Attributes References EPSIT Param 31 R(4) 4 scalar EPSU2 Param 31 R(4) 4 scalar 115,155 EPSUST Param 31 R(4) 4 scalar 161,202,259 EXCM Param 25 R(4) 4 scalar 146 EXP Func 163 scalar 163,261,406,408,422,424,426,440,44 2 FHNEU Param 24 R(4) 4 scalar 25,194,199,335,336,337,341,342,343 FZQ1 Param 40 R(4) 4 scalar 90 FZQ2 Param 41 R(4) 4 scalar 105 FZT1 Param 40 R(4) 4 scalar 86 FZT2 Param 41 R(4) 4 scalar 101 FZU1 Param 40 R(4) 4 scalar 80 G Param 23 R(4) 4 scalar 23,26,27,424,440 G2LI Scalar 55 I(4) 4 1 239 COM G2LJ Scalar 55 I(4) 4 1 389 COM GLB_TABLE Common 37 128 GLKBR Param 34 R(4) 4 scalar 34 GLKBS Param 34 R(4) 4 scalar 34,35 GRRS Param 34 R(4) 4 scalar 41 H1M5 Param 47 R(4) 4 scalar 407,421,425,441 HLFLX Local 300 R(4) 4 scalar 300,409,428,443 HPBL Dummy 3 R(4) 4 scalar ARG,INOUT 159 HSFLX Local 299 R(4) 4 scalar 299,394,395,396 HTMG Scalar 52 R(4) 4 3 4183695 COM IBROW Scalar 22 I(4) 4 scalar COM ICHUNKTAB Scalar 29 I(4) 4 1 8 COM IDIM1 Param 19 I(4) 4 scalar IDIM2 Param 19 I(4) 4 scalar IE_GLB_TABLE Scalar 38 I(4) 4 1 8 COM IE_LOC_TABLE Scalar 28 I(4) 4 1 8 COM IGSTL Param 14 I(4) 4 scalar 19 IGSTR Param 14 I(4) 4 scalar 19 ILCOL Scalar 22 I(4) 4 scalar COM ILPAD1 Scalar 23 I(4) 4 scalar COM ILPAD2 Scalar 23 I(4) 4 scalar COM ILPAD3 Scalar 23 I(4) 4 scalar COM ILPAD4 Scalar 23 I(4) 4 scalar COM ILPAD5 Scalar 23 I(4) 4 scalar COM IM Param 3 I(4) 4 scalar 16,19,43,44,45,46,47,51,52,57 INPES Param 11 I(4) 4 scalar 16,19,27,28,29,38,39 INUMQ Scalar 34 I(4) 4 1 100 COM IQUILT_GROUP Scalar 33 I(4) 4 scalar COM IRCOL Scalar 22 I(4) 4 scalar COM IRPAD1 Scalar 24 I(4) 4 scalar COM IRPAD2 Scalar 24 I(4) 4 scalar COM IRPAD3 Scalar 24 I(4) 4 scalar COM IRPAD4 Scalar 24 I(4) 4 scalar COM IRPAD5 Scalar 24 I(4) 4 scalar COM IS_GLB_TABLE Scalar 38 I(4) 4 1 8 COM IS_LOC_TABLE Scalar 27 I(4) 4 1 8 COM ITAIL Param 16 I(4) 4 scalar ITEMP Scalar 47 I(4) 4 2 91719 COM ITEMP2 Scalar 47 I(4) 4 2 91719 COM ITR Local 181 I(4) 4 scalar 181,224 ITRMX Param 28 I(4) 4 scalar 181,224 Page 12 Source Listing SFCDIF 2014-12-17 20:47 Symbol Table SFCDIF.F Name Object Declared Type Bytes Dimen Elements Attributes References ITROW Scalar 22 I(4) 4 scalar COM JBPAD1 Scalar 25 I(4) 4 scalar COM JBPAD2 Scalar 25 I(4) 4 scalar COM JBPAD3 Scalar 25 I(4) 4 scalar COM JBPAD4 Scalar 25 I(4) 4 scalar COM JBPAD5 Scalar 25 I(4) 4 scalar COM JDIM1 Param 20 I(4) 4 scalar JDIM2 Param 20 I(4) 4 scalar JE_GLB_TABLE Scalar 39 I(4) 4 1 8 COM JE_LOC_TABLE Scalar 28 I(4) 4 1 8 COM JGSTL Param 15 I(4) 4 scalar 20 JGSTR Param 15 I(4) 4 scalar 20 JM Param 3 I(4) 4 scalar 17,20,43,44,45,46,47,51,52,57 JNPES Param 11 I(4) 4 scalar 17,20,27,28,29,38,39 JS_GLB_TABLE Scalar 39 I(4) 4 1 8 COM JS_LOC_TABLE Scalar 27 I(4) 4 1 8 COM JTAIL Param 17 I(4) 4 scalar JTPAD1 Scalar 26 I(4) 4 scalar COM JTPAD2 Scalar 26 I(4) 4 scalar COM JTPAD3 Scalar 26 I(4) 4 scalar COM JTPAD4 Scalar 26 I(4) 4 scalar COM JTPAD5 Scalar 26 I(4) 4 scalar COM L2GI Scalar 57 I(4) 4 1 239 COM L2GJ Scalar 57 I(4) 4 1 389 COM LM Param 3 I(4) 4 scalar 20,50,52 LMHK Dummy 2 I(4) 4 scalar ARG,INOUT 69,71,72,141 LMHP Local 69 I(4) 4 scalar 69,141 LP1 Param 20 I(4) 4 scalar 53 LSM Param 3 I(4) 4 scalar MAPPINGS Common 56 5024 MH Local 124 I(4) 4 scalar 124,127 ML Local 118 I(4) 4 scalar 118,120,121,122,124,126,127 MOD Func 124 scalar 124 MPI_COMM_COMP Scalar 30 I(4) 4 scalar COM MPI_COMM_INTER Scalar 31 I(4) 4 scalar COM MPI_COMM_INTER_ARRAY Scalar 32 I(4) 4 1 100 COM MPPCOM Common 1 1464 MYIE Scalar 4 I(4) 4 scalar COM MYIE1 Scalar 4 I(4) 4 scalar COM MYIE1_P1 Scalar 8 I(4) 4 scalar COM MYIE1_P2 Scalar 8 I(4) 4 scalar COM MYIE1_P3 Scalar 8 I(4) 4 scalar COM MYIE1_P4 Scalar 8 I(4) 4 scalar COM MYIE2 Scalar 4 I(4) 4 scalar COM MYIE2_P1 Scalar 9 I(4) 4 scalar COM MYIE_P1 Scalar 7 I(4) 4 scalar COM MYIE_P2 Scalar 7 I(4) 4 scalar COM MYIE_P3 Scalar 7 I(4) 4 scalar COM MYIE_P4 Scalar 7 I(4) 4 scalar COM MYIE_P5 Scalar 7 I(4) 4 scalar COM MYIS Scalar 4 I(4) 4 scalar COM MYIS1 Scalar 4 I(4) 4 scalar COM MYIS1_P1 Scalar 6 I(4) 4 scalar COM MYIS1_P2 Scalar 6 I(4) 4 scalar COM MYIS1_P3 Scalar 6 I(4) 4 scalar COM Page 13 Source Listing SFCDIF 2014-12-17 20:47 Symbol Table SFCDIF.F Name Object Declared Type Bytes Dimen Elements Attributes References MYIS1_P4 Scalar 6 I(4) 4 scalar COM MYIS2 Scalar 4 I(4) 4 scalar COM MYIS_P1 Scalar 5 I(4) 4 scalar COM MYIS_P2 Scalar 5 I(4) 4 scalar COM MYIS_P3 Scalar 5 I(4) 4 scalar COM MYIS_P4 Scalar 5 I(4) 4 scalar COM MYIS_P5 Scalar 5 I(4) 4 scalar COM MYJE Scalar 15 I(4) 4 scalar COM MYJE1 Scalar 15 I(4) 4 scalar COM MYJE1_P1 Scalar 17 I(4) 4 scalar COM MYJE1_P2 Scalar 17 I(4) 4 scalar COM MYJE1_P3 Scalar 17 I(4) 4 scalar COM MYJE1_P4 Scalar 17 I(4) 4 scalar COM MYJE2 Scalar 15 I(4) 4 scalar COM MYJE2_P1 Scalar 18 I(4) 4 scalar COM MYJE2_P2 Scalar 18 I(4) 4 scalar COM MYJE2_P3 Scalar 18 I(4) 4 scalar COM MYJE2_P4 Scalar 18 I(4) 4 scalar COM MYJE3 Scalar 15 I(4) 4 scalar COM MYJE3_P4 Scalar 18 I(4) 4 scalar COM MYJE4 Scalar 15 I(4) 4 scalar COM MYJE4_P1 Scalar 19 I(4) 4 scalar COM MYJE4_P4 Scalar 19 I(4) 4 scalar COM MYJE5 Scalar 15 I(4) 4 scalar COM MYJE5_P1 Scalar 19 I(4) 4 scalar COM MYJE5_P2 Scalar 19 I(4) 4 scalar COM MYJE_P1 Scalar 16 I(4) 4 scalar COM MYJE_P2 Scalar 16 I(4) 4 scalar COM MYJE_P3 Scalar 16 I(4) 4 scalar COM MYJE_P4 Scalar 16 I(4) 4 scalar COM MYJE_P5 Scalar 16 I(4) 4 scalar COM MYJS Scalar 10 I(4) 4 scalar COM MYJS1 Scalar 10 I(4) 4 scalar COM MYJS1_P1 Scalar 12 I(4) 4 scalar COM MYJS1_P2 Scalar 12 I(4) 4 scalar COM MYJS1_P3 Scalar 12 I(4) 4 scalar COM MYJS1_P4 Scalar 12 I(4) 4 scalar COM MYJS2 Scalar 10 I(4) 4 scalar COM MYJS2_P1 Scalar 13 I(4) 4 scalar COM MYJS2_P2 Scalar 13 I(4) 4 scalar COM MYJS2_P3 Scalar 13 I(4) 4 scalar COM MYJS2_P4 Scalar 13 I(4) 4 scalar COM MYJS3 Scalar 10 I(4) 4 scalar COM MYJS3_P4 Scalar 13 I(4) 4 scalar COM MYJS4 Scalar 10 I(4) 4 scalar COM MYJS4_P1 Scalar 14 I(4) 4 scalar COM MYJS4_P4 Scalar 14 I(4) 4 scalar COM MYJS5 Scalar 10 I(4) 4 scalar COM MYJS5_P1 Scalar 14 I(4) 4 scalar COM MYJS5_P2 Scalar 14 I(4) 4 scalar COM MYJS_P1 Scalar 11 I(4) 4 scalar COM MYJS_P2 Scalar 11 I(4) 4 scalar COM MYJS_P3 Scalar 11 I(4) 4 scalar COM MYJS_P4 Scalar 11 I(4) 4 scalar COM MYJS_P5 Scalar 11 I(4) 4 scalar COM Page 14 Source Listing SFCDIF 2014-12-17 20:47 Symbol Table SFCDIF.F Name Object Declared Type Bytes Dimen Elements Attributes References MYPE Scalar 2 I(4) 4 scalar COM MY_E Scalar 20 I(4) 4 scalar COM MY_IE_GLB Scalar 2 I(4) 4 scalar COM MY_IE_LOC Scalar 3 I(4) 4 scalar COM MY_IS_GLB Scalar 2 I(4) 4 scalar COM MY_IS_LOC Scalar 3 I(4) 4 scalar COM MY_JE_GLB Scalar 2 I(4) 4 scalar COM MY_JE_LOC Scalar 3 I(4) 4 scalar COM MY_JS_GLB Scalar 2 I(4) 4 scalar COM MY_JS_LOC Scalar 3 I(4) 4 scalar COM MY_N Scalar 20 I(4) 4 scalar COM MY_NE Scalar 21 I(4) 4 scalar COM MY_NEB Scalar 21 I(4) 4 1 8 COM MY_NW Scalar 21 I(4) 4 scalar COM MY_S Scalar 20 I(4) 4 scalar COM MY_SE Scalar 21 I(4) 4 scalar COM MY_SW Scalar 21 I(4) 4 scalar COM MY_W Scalar 20 I(4) 4 scalar COM NPES Scalar 2 I(4) 4 scalar COM P10 Local 424 R(4) 4 scalar 424,425,426 P30 Local 440 R(4) 4 scalar 440,441,442 PD Dummy 6 R(4) 4 scalar ARG,INOUT 404 PDS Local 404 R(4) 4 scalar 404,406,424,440 PIFR Param 29 R(4) 4 scalar 117 PIHF Param 29 R(4) 4 scalar 64,116 PQ0 Param 46 R(4) 4 scalar 408,422,426,442 PSHLTR Local 406 R(4) 4 scalar 406,407,408,421,422 PSHZ Local 193 R(4) 4 scalar 193,194,198,199,246,247,255,256,33 5,336,337,341,342,343,363,364,365, 375,376,377 PSLHS Local 61 R(4) 4 scalar 198,199,341,342,343 PSLHU Local 60 R(4) 4 scalar 193,194,335,336,337 PSLMS Local 59 R(4) 4 scalar 196,197,339,340 PSLMU Local 58 R(4) 4 scalar 191,192,333,334 PSMZ Local 191 R(4) 4 scalar 191,192,196,197,244,245,253,254,33 3,334,339,340,361,362,373,374 PSPHS Local 67 R(4) 4 scalar 255,256,375,376,377 PSPHU Local 66 R(4) 4 scalar 246,247,363,364,365 PSPMS Local 65 R(4) 4 scalar 253,254,373,374,489,490 PSPMU Local 63 R(4) 4 scalar 244,245,361,362,484,485 PT Dummy 6 R(4) 4 scalar ARG,INOUT 404 Q Dummy 6 R(4) 4 1 45 ARG,INOUT 72 Q02 Dummy 4 R(4) 4 scalar ARG,INOUT 409,410,412,414,417,418 Q10 Dummy 4 R(4) 4 scalar ARG,INOUT 428,429,431,433,436,437 Q30 Dummy 5 R(4) 4 scalar ARG,INOUT 443,444,446,448,451,452 QLM Local 72 R(4) 4 scalar 72,92,107,149,151,300,411,412,430, 431,445,446 QS Dummy 2 R(4) 4 scalar ARG,INOUT 92,107,137 QSAT10 Local 422 R(4) 4 scalar 422,426,436,437 QSAT2 Local 408 R(4) 4 scalar 408,417,418 QSAT30 Local 442 R(4) 4 scalar 442,451,452 QVISC Param 36 R(4) 4 scalar 37,40,41 QZ0 Dummy 2 R(4) 4 scalar ARG,INOUT 92,107,137,149,151,300,409,428,443 RDZ Local 145 R(4) 4 scalar 145,146,167,168 RFAC Param 25 R(4) 4 scalar 61 Page 15 Source Listing SFCDIF 2014-12-17 20:47 Symbol Table SFCDIF.F Name Object Declared Type Bytes Dimen Elements Attributes References RFC Param 24 R(4) 4 scalar 25 RIB Local 157 R(4) 4 scalar 157,165,302 RIC Param 24 R(4) 4 scalar 25,165,302 RLMA Local 212 R(4) 4 scalar 212,216,273,277 RLMN Local 209 R(4) 4 scalar 209,212,270,273 RLMO Local 178 R(4) 4 scalar 178,183,184,185,187,188,190,211,21 2,216,226,227,228,230,231,233,272, 273,277 RLMP Local 211 R(4) 4 scalar 211,272,324,325,326,327,328,332,34 8,473,474,475,477 RLNT02 Local 320 R(4) 4 scalar 320,335,341,363,375 RLNT10 Local 321 R(4) 4 scalar 321,336,342,364,376 RLNT30 Local 322 R(4) 4 scalar 322,337,343,365,377 RLNU10 Local 318 R(4) 4 scalar 318,333,339,361,373,471,484,489 RLNU30 Local 319 R(4) 4 scalar 319,334,340,362,374,472,485,490 RLOGT Local 176 R(4) 4 scalar 176,194,199,247,256,263 RLOGU Local 175 R(4) 4 scalar 175,192,197,245,254 RQVISC Param 37 R(4) 4 scalar 91,106 RRIC Param 25 R(4) 4 scalar 59 RTVISC Param 37 R(4) 4 scalar 40,41,87,102 RVISC Param 37 R(4) 4 scalar 40,80,81,101 RWGH Local 82 R(4) 4 scalar 82,83,84 SFCDIF Subr 2 SIMH Local 194 R(4) 4 scalar 194,199,206,247,256,267 SIMH02 Local 335 R(4) 4 scalar 335,341,363,375,384 SIMH10 Local 336 R(4) 4 scalar 336,342,364,376,385 SIMH30 Local 337 R(4) 4 scalar 337,343,365,377,386 SIMM Local 192 R(4) 4 scalar 192,197,205,245,254,266 SIMM10 Local 333 R(4) 4 scalar 333,339,361,373,382,484,489,493 SIMM30 Local 334 R(4) 4 scalar 334,340,362,374,383,485,490,494 SM Dummy 2 R(4) 4 scalar ARG,INOUT 74,76,114,148,163,165,180,302,330, 462 SMALL Param 34 R(4) 4 scalar 35 SQPR Param 38 R(4) 4 scalar 38,40,41 SQRT Func 80 scalar 80,101,161,163,202,239,240,241,242 ,259,261,295,355,356,357,358,359,4 81,482,483 SQSC Param 38 R(4) 4 scalar 38 SQVISC Param 45 R(4) 4 scalar 45 T Dummy 6 R(4) 4 1 45 ARG,INOUT 71 T02 Local 407 R(4) 4 scalar 407,408 T10 Local 421 R(4) 4 scalar 421,422,425,426 T30 Local 441 R(4) 4 scalar 441,442 TEMP1 Scalar 43 R(4) 4 2 91719 COM TEMP10 Scalar 45 R(4) 4 2 91719 COM TEMP11 Scalar 45 R(4) 4 2 91719 COM TEMP12 Scalar 45 R(4) 4 2 91719 COM TEMP13 Scalar 46 R(4) 4 2 91719 COM TEMP14 Scalar 46 R(4) 4 2 91719 COM TEMP15 Scalar 46 R(4) 4 2 91719 COM TEMP16 Scalar 46 R(4) 4 2 91719 COM TEMP2 Scalar 43 R(4) 4 2 91719 COM TEMP2X Scalar 51 R(4) 4 2 92971 COM TEMP3 Scalar 43 R(4) 4 2 91719 COM TEMP4 Scalar 43 R(4) 4 2 91719 COM Page 16 Source Listing SFCDIF 2014-12-17 20:47 Symbol Table SFCDIF.F Name Object Declared Type Bytes Dimen Elements Attributes References TEMP5 Scalar 44 R(4) 4 2 91719 COM TEMP6 Scalar 44 R(4) 4 2 91719 COM TEMP7 Scalar 44 R(4) 4 2 91719 COM TEMP8 Scalar 44 R(4) 4 2 91719 COM TEMP9 Scalar 45 R(4) 4 2 91719 COM TEMPCOM Common 42 6603768 TERM1 Local 405 R(4) 4 scalar 405,406 TH02 Dummy 4 R(4) 4 scalar ARG,INOUT 394,407 TH10 Dummy 4 R(4) 4 scalar ARG,INOUT 395,421,425 TH30 Dummy 5 R(4) 4 scalar ARG,INOUT 396,441 THLM Local 71 R(4) 4 scalar 71,88,103,149,151,299 THS Dummy 2 R(4) 4 scalar ARG,INOUT 88,103,134 THZ0 Dummy 2 R(4) 4 scalar ARG,INOUT 88,103,134,149,151,299,394,395,396 TLM Dummy 6 R(4) 4 scalar ARG,INOUT 405,424,440 TOPO Common 50 17478548 TTVG Scalar 51 R(4) 4 2 92971 COM TVISC Param 36 R(4) 4 scalar 37,40,41,168,306,307,308 U10 Dummy 4 R(4) 4 scalar ARG,INOUT 390,457,501 U10E Local 457 R(4) 4 scalar 457,495,501 U30 Dummy 5 R(4) 4 scalar ARG,INOUT 392,459,503 U30E Local 459 R(4) 4 scalar 459,497,503 ULM Dummy 6 R(4) 4 scalar ARG,INOUT 83,115,116,121,155,297 UMFLX Local 297 R(4) 4 scalar 297,390,392,495,497 USTAR Dummy 3 R(4) 4 scalar ARG,INOUT 74,76,78,80,95,101,161,163,178,202 ,204,209,259,261,265,270 USTARK Local 204 R(4) 4 scalar 204,205,206,265,266,267,382,383,38 4,385,386,493,494 USTC Param 39 R(4) 4 scalar 76,95 USTFC Param 23 R(4) 4 scalar 74 USTR Param 39 R(4) 4 scalar 78,95 UZ0 Dummy 2 R(4) 4 scalar ARG,INOUT 83,98,130,155,297,390,392,495,497 V10 Dummy 4 R(4) 4 scalar ARG,INOUT 391,458,502 V10E Local 458 R(4) 4 scalar 458,496,502 V30 Dummy 5 R(4) 4 scalar ARG,INOUT 393,460,504 V30E Local 460 R(4) 4 scalar 460,498,504 VISC Param 36 R(4) 4 scalar 37,40,167,304,305 VKRM Param 24 R(4) 4 scalar 27,45,204,265 VLM Dummy 6 R(4) 4 scalar ARG,INOUT 84,116,121,155,298 VMFLX Local 298 R(4) 4 scalar 298,391,393,496,498 VZ0 Dummy 2 R(4) 4 scalar ARG,INOUT 84,99,131,155,298,391,393,496,498 WGHT Local 81 R(4) 4 scalar 81,82,87,88,91,92,102,103,106,107 WLOW Local 126 R(4) 4 scalar 126,127 WNEW Param 28 R(4) 4 scalar 212,273 WOLD Param 28 R(4) 4 scalar 28,212,273 WSTAR Dummy 3 R(4) 4 scalar ARG,INOUT 295 WSTAR2 Local 160 R(4) 4 scalar 160,161,202,208,259,269,295 WWST Param 23 R(4) 4 scalar 23,295 WWST2 Param 23 R(4) 4 scalar 160,208,269 X Local 117 R(4) 4 scalar 117,118,126 XLT Local 240 R(4) 4 scalar 240,247 XLT02 Local 357 R(4) 4 scalar 357,363 XLT024 Local 351 R(4) 4 scalar 351,357 XLT10 Local 358 R(4) 4 scalar 358,364 XLT104 Local 352 R(4) 4 scalar 352,358 XLT30 Local 359 R(4) 4 scalar 359,365 Page 17 Source Listing SFCDIF 2014-12-17 20:47 Symbol Table SFCDIF.F Name Object Declared Type Bytes Dimen Elements Attributes References XLT304 Local 353 R(4) 4 scalar 353,359 XLT4 Local 235 R(4) 4 scalar 235,240 XLU Local 239 R(4) 4 scalar 239,245 XLU10 Local 355 R(4) 4 scalar 355,361,481,484 XLU104 Local 349 R(4) 4 scalar 349,355,478,481 XLU30 Local 356 R(4) 4 scalar 356,362,482,485 XLU304 Local 350 R(4) 4 scalar 350,356,479,482 XLU4 Local 234 R(4) 4 scalar 234,239 XT Local 242 R(4) 4 scalar 242,246 XT4 Local 237 R(4) 4 scalar 237,242 XU Local 241 R(4) 4 scalar 241,244 XU10 Local 483 R(4) 4 scalar 483,484,485 XU104 Local 480 R(4) 4 scalar 480,483 XU4 Local 236 R(4) 4 scalar 236,241 Z Dummy 6 R(4) 4 1 46 ARG,INOUT 141 Z0 Dummy 3 R(4) 4 scalar ARG,INOUT 74,80,97,101,112,133,136,152,163,2 61 ZEFF Dummy 3 R(4) 4 1 4 ARG,INOUT 127 ZETALT Local 183 R(4) 4 scalar 183,184,194,199,226,227,235,252,25 6 ZETALU Local 185 R(4) 4 scalar 185,192,197,228,234,251,254 ZETAT Local 188 R(4) 4 scalar 188,193,198,231,237,250,255 ZETAU Local 187 R(4) 4 scalar 187,191,196,230,236,249,253 ZILFC Param 45 R(4) 4 scalar 163,261 ZQ Local 90 R(4) 4 scalar 90,91,105,106,136 ZQRZT Param 38 R(4) 4 scalar 40,41 ZSL Local 141 R(4) 4 scalar 141,143,145,157,172,173,262 ZSLT Local 173 R(4) 4 scalar 173,176,183,184,226,227,262,263 ZSLU Local 172 R(4) 4 scalar 172,175,185,228 ZT Local 86 R(4) 4 scalar 86,87,90,101,102,105,133,152,163,1 73,176,188,231,261,262,263,314,315 ,316,320,321,322 ZT02 Local 314 R(4) 4 scalar 314,320,326 ZT10 Local 315 R(4) 4 scalar 315,321,327 ZT30 Local 316 R(4) 4 scalar 316,322,328 ZTAT02 Local 326 R(4) 4 scalar 326,335,341,351,369,375 ZTAT10 Local 327 R(4) 4 scalar 327,336,342,352,370,376 ZTAT30 Local 328 R(4) 4 scalar 328,337,343,353,371,377 ZTAU Local 473 R(4) 4 scalar 473,480,489,490 ZTAU10 Local 324 R(4) 4 scalar 324,333,339,349,367,373,474,478,48 7,489 ZTAU30 Local 325 R(4) 4 scalar 325,334,340,350,368,374,475,479,48 8,490 ZTFC Param 43 R(4) 4 scalar 152 ZTMAX Param 32 R(4) 4 scalar 249,250,251,252,367,368,369,370,37 1,487,488 ZTMIN Param 32 R(4) 4 scalar 183,226 ZU Local 80 R(4) 4 scalar 80,81,86,97,112,127,143,172,175,18 7,230,312,313,318,319,467,468,469, 470,471,472,473 ZU10 Local 312 R(4) 4 scalar 312,318,324,469,471,474 ZU30 Local 313 R(4) 4 scalar 313,319,325,470,472,475 ZUUZ Local 467 R(4) 4 scalar 467,468 Page 18 Source Listing SFCDIF 2014-12-17 20:47 Subprograms/Common Blocks SFCDIF.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References GLB_TABLE Common 37 128 MAPPINGS Common 56 5024 MPPCOM Common 1 1464 SFCDIF Subr 2 TEMPCOM Common 42 6603768 TOPO Common 50 17478548 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__ -D __tune_pentium4__ -D __SSE2__ -D __SSE__ -D __MMX__ -double_size 64 no -d_lines no -Qdyncom -error_limit 30 no -f66 no -f77rtl Page 19 Source Listing SFCDIF 2014-12-17 20:47 SFCDIF.F no -fast -fpscomp nofilesfromcmd -fpscomp nogeneral -fpscomp noioformat -fpscomp noldio_spacing -fpscomp nologicals -fixed no -fpconstant -fpe3 -fprm nearest no -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 nono_prec_div -fp_modbits noprec_sqrt -fp_modbits nono_prec_sqrt -fp_modbits noftz -fp_modbits no_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 nono_precise_libm_functions -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 : /usrx/local/intel/composerxe/tbb/include/,/usr/include/,./,/opt/ibmhpc/pe1308/mpich2/intel/include64/, /opt/ibmhpc/pe1308/base/include64/,/usrx/local/intel/composerxe/mkl/include/,/usrx/local/intel/composerxe/tbb/include/, /gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/intel64/,/gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/, /usr/local/include/,/usr/lib/gcc/x86_64-redhat-linux/4.4.7/include/,/usr/include/,/usr/include/ -list filename : SFCDIF.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100