Changeset 155
- Timestamp:
- Nov 10, 2000, 11:43:55 AM (24 years ago)
- Location:
- LMDZ.3.3/branches/rel-LF/libf/phylmd
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/branches/rel-LF/libf/phylmd/YOEGWD.h
r2 r155 3 3 C ----------------------------------------------------------------- 4 4 C 5 integer NKTOPG,NSTRA 6 real GFRCRIT,GKWAKE,GRCRIT,GVCRIT,GKDRAG,GKLIFT 7 real GHMAX,GRAHILO,GSIGCR,GSSEC,GTSEC,GVSEC 5 8 COMMON/YOEGWD/ GFRCRIT,GKWAKE,GRCRIT,GVCRIT,GKDRAG,GKLIFT 6 9 * ,GHMAX,GRAHILO,GSIGCR,NKTOPG,NSTRA,GSSEC,GTSEC,GVSEC -
LMDZ.3.3/branches/rel-LF/libf/phylmd/orografi.F
r2 r155 1 1 SUBROUTINE drag_noro (nlon,nlev,dtime,paprs,pplay, 2 2 e pmea,pstd, psig, pgam, pthe,ppic,pval, 3 e kgwd,k gwdim,kdx,ktest,3 e kgwd,kdx,ktest, 4 4 e t, u, v, 5 5 s pulow, pvlow, pustr, pvstr, … … 39 39 REAL d_t(nlon,nlev), d_u(nlon,nlev), d_v(nlon,nlev) 40 40 c 41 INTEGER i, k, kgwd, k gwdim, kdx(nlon), ktest(nlon)41 INTEGER i, k, kgwd, kdx(nlon), ktest(nlon) 42 42 c 43 43 c Variables locales: … … 96 96 c appeler la routine principale 97 97 c 98 CALL orodrag(klon,klev,kgwd,k gwdim,kdx,ktest,98 CALL orodrag(klon,klev,kgwd,kdx,ktest, 99 99 . dtime, 100 100 . papmh, papmf, zgeom, … … 119 119 END 120 120 SUBROUTINE orodrag( nlon,nlev 121 i , kgwd, k gwdim, kdx, ktest121 i , kgwd, kdx, ktest 122 122 r , ptsphy 123 123 r , paphm1,papm1,pgeom1,ptm1,pum1,pvm1 … … 126 126 r , pulow,pvlow 127 127 r , pvom,pvol,pte ) 128 129 implicit none 130 128 131 c 129 132 c … … 152 155 c -------------------- 153 156 c 154 implicit logical (l)157 c implicit logical (l) 155 158 c 156 159 c method. … … 183 186 c 184 187 c 188 integer nlon, nlev, klevm1 189 integer kgwd, jl, ilevp1, jk, ji 190 real zdelp, ztemp, zforc, ztend 191 real rover, zb, zc, zconb, zabsv 192 real zzd1, ratio, zbet, zust,zvst, zdis 185 193 real pte(nlon,nlev), 186 194 * pvol(nlon,nlev), … … 223 231 * zdedt(klon), 224 232 * zvidis(klon), 225 * ztfr(klon),226 233 * znu(klon), 227 234 * zd1(klon), 228 235 * zd2(klon), 229 236 * zdmod(klon) 237 real ztmst, ptsphy, zrtmst 230 238 c 231 239 c------------------------------------------------------------------ … … 272 280 c 273 281 call orosetup 274 * ( nlon, nlev ,ktest275 * , ikcrit, ikcrith, icrit, isect, ikhlim,ikenvh,iknu,iknu2282 * ( nlon, ktest 283 * , ikcrit, ikcrith, icrit, ikenvh,iknu,iknu2 276 284 * , paphm1, papm1 , pum1 , pvm1 , ptm1 , pgeom1, pstd 277 285 * , zrho , zri , zstab , ztau , zvph , zpsi, zzdep … … 292 300 call gwstress 293 301 * ( nlon , nlev 294 * , ikcrit, isect, ikhlim, ktest, ikcrith, icrit, ikenvh, iknu295 * , zrho , zstab, zvph , pstd, psig, pmea, ppic , pval296 * , zt fr , ztau297 * , pgeom1, pgamma,zd1,zd2,zdmod,znu)302 * , ktest , icrit, ikenvh, iknu 303 * , zrho , zstab, zvph , pstd, psig, pmea, ppic 304 * , ztau 305 * , pgeom1,zdmod) 298 306 c 299 307 c … … 306 314 call gwprofil 307 315 * ( nlon , nlev 308 * , kgwd , kdx , ktest 309 * , ikcrit, ikcrith, icrit , ikenvh, iknu 310 * ,iknu2 , paphm1, zrho , zstab , ztfr , zvph 311 * , zri , ztau , ztauf 312 c 313 * , zdmod , znu , psig , pgamma , pstd , ppic , pval) 316 * , kgwd , kdx 317 * , ikcrith, icrit 318 * , paphm1, zrho , zstab , zvph 319 * , zri , ztau 320 * , zdmod , psig , pstd) 314 321 c 315 322 c … … 401 408 end 402 409 SUBROUTINE orosetup 403 * ( nlon , nlev ,ktest404 * , kkcrit, kkcrith, kcrit , ksect , kkhlim410 * ( nlon , ktest 411 * , kkcrit, kkcrith, kcrit 405 412 * , kkenvh, kknu , kknu2 406 413 * , paphm1, papm1 , pum1 , pvm1 , ptm1 , pgeom1, pstd … … 448 455 c 449 456 c----------------------------------------------------------------------- 450 implicit logical (l)457 implicit none 451 458 c 452 459 … … 461 468 c --------- 462 469 c 463 integer kkcrit(nlon),kkcrith(nlon),kcrit(nlon),ksect(nlon), 464 * kkhlim(nlon),ktest(nlon),kkenvh(nlon) 470 integer nlon 471 integer jl, jk 472 real zdelp 473 real zb, zc 474 475 integer kkcrit(nlon),kkcrith(nlon),kcrit(nlon), 476 * ktest(nlon),kkenvh(nlon) 465 477 466 478 c … … 480 492 c 481 493 c 494 integer ilevm1, ilevm2, ilevh 495 real zcons1, zcons2,zcons3, zhgeo 496 real zu, zphi, zvt1,zvt2, zst, zvar, zdwind, zwind 497 real zstabm, zstabp, zrhom, zrhop, alpha 498 real zggeenv, zggeom1,zgvar 499 logical lo 482 500 logical ll1(klon,klev+1) 483 501 integer kknu(klon),kknu2(klon),kknub(klon),kknul(klon), … … 486 504 real zhcrit(klon,klev),zvpf(klon,klev), 487 505 * zdp(klon,klev) 488 real znorm(klon),z psi(klon),zb(klon),zc(klon),506 real znorm(klon),zb(klon),zc(klon), 489 507 * zulow(klon),zvlow(klon),znup(klon),znum(klon) 490 508 c … … 862 880 SUBROUTINE gwstress 863 881 * ( nlon , nlev 864 * , k kcrit, ksect, kkhlim, ktest, kkcrith, kcrit, kkenvh882 * , ktest, kcrit, kkenvh 865 883 * , kknu 866 884 * , prho , pstab , pvph , pstd, psig 867 * , pmea , ppic , pval , ptfr, ptau868 * , pgeom1 , p gamma , pd1 , pd2 , pdmod , pnu)885 * , pmea , ppic , ptau 886 * , pgeom1 , pdmod ) 869 887 c 870 888 c**** *gwstress* … … 906 924 c 907 925 c----------------------------------------------------------------------- 908 implicit logical (l)926 implicit none 909 927 #include "dimensions.h" 910 928 #include "dimphy.h" … … 917 935 c --------- 918 936 c 919 integer kkcrit(nlon),kkcrith(nlon),kcrit(nlon),ksect(nlon), 920 * kkhlim(nlon),ktest(nlon),kkenvh(nlon),kknu(nlon) 937 integer nlon, nlev 938 integer kcrit(nlon), 939 * ktest(nlon),kkenvh(nlon),kknu(nlon) 921 940 c 922 941 real prho(nlon,nlev+1),pstab(nlon,nlev+1),ptau(nlon,nlev+1), 923 * pvph(nlon,nlev+1), ptfr(nlon),942 * pvph(nlon,nlev+1), 924 943 * pgeom1(nlon,nlev),pstd(nlon) 925 944 c 926 real p d1(nlon),pd2(nlon),pnu(nlon),psig(nlon),pgamma(nlon)927 real pmea(nlon),ppic(nlon) ,pval(nlon)945 real psig(nlon) 946 real pmea(nlon),ppic(nlon) 928 947 real pdmod(nlon) 929 948 c … … 932 951 c* 0.2 local arrays 933 952 c ------------ 953 integer jl 954 real zblock, zvar, zeff 955 logical lo 934 956 c 935 957 c----------------------------------------------------------------------- … … 985 1007 SUBROUTINE GWPROFIL 986 1008 * ( NLON, NLEV 987 * , kgwd ,kdx , ktest988 * , KKCRIT , KKCRITH, KCRIT , kkenvh, kknu,kknu2989 * , PAPHM1, PRHO , PSTAB , PTFR, PVPH , PRI , PTAU990 * , p tauf ,pdmod , pnu , psig ,pgamma, pvar, ppic,pval)1009 * , kgwd, kdx 1010 * , KKCRITH, KCRIT 1011 * , PAPHM1, PRHO , PSTAB , PVPH , PRI , PTAU 1012 * , pdmod , psig , pvar) 991 1013 992 1014 C**** *GWPROFIL* … … 1036 1058 C PASSAGE OF THE NEW GWDRAG TO I.F.S. (F. LOTT, 22/11/93) 1037 1059 C----------------------------------------------------------------------- 1038 IMPLICIT LOGICAL (L)1060 implicit none 1039 1061 C 1040 1062 … … 1051 1073 C --------- 1052 1074 C 1053 INTEGER KKCRIT(NLON),KKCRITH(NLON),KCRIT(NLON)1054 * ,kdx(nlon),ktest(NLON)1055 * ,k kenvh(NLON),kknu(NLON),kknu2(NLON)1075 integer nlon,nlev 1076 INTEGER KKCRITH(NLON),KCRIT(NLON) 1077 * ,kdx(nlon) 1056 1078 C 1057 1079 REAL PAPHM1(NLON,NLEV+1), PSTAB(NLON,NLEV+1), 1058 1080 * PRHO (NLON,NLEV+1), PVPH (NLON,NLEV+1), 1059 * PRI (NLON,NLEV+1), PTFR (NLON), PTAU(NLON,NLEV+1), 1060 * ptauf (NLON,nlev+1) 1081 * PRI (NLON,NLEV+1), PTAU(NLON,NLEV+1) 1061 1082 1062 REAL pdmod (NLON) , p nu (NLON) , psig(NLON),1063 * p gamma(NLON) , pvar(NLON) , ppic(NLON), pval(NLON)1083 REAL pdmod (NLON) , psig(NLON), 1084 * pvar(NLON) 1064 1085 1065 1086 C----------------------------------------------------------------------- … … 1068 1089 C ------------ 1069 1090 C 1091 integer ilevh, ji, kgwd, jl, jk 1092 real zsqr, zalfa, zriw, zdel, zb, zalpha,zdz2n 1093 real zdelp, zdelpt 1070 1094 REAL ZDZ2 (KLON,KLEV) , ZNORM(KLON) , zoro(KLON) 1071 1095 REAL ZTAU (KLON,KLEV+1) … … 1221 1245 END 1222 1246 SUBROUTINE lift_noro (nlon,nlev,dtime,paprs,pplay, 1223 e plat,pmea,pstd, p sig, pgam, pthe, ppic,pval,1224 e k gwd,kgwdim,kdx,ktest,1247 e plat,pmea,pstd, ppic, 1248 e ktest, 1225 1249 e t, u, v, 1226 1250 s pulow, pvlow, pustr, pvstr, … … 1255 1279 REAL pplay(klon,klev) 1256 1280 REAL plat(nlon),pmea(nlon) 1257 REAL pstd(nlon) ,psig(nlon),pgam(nlon),pthe(nlon)1258 REAL ppic(nlon) ,pval(nlon)1281 REAL pstd(nlon) 1282 REAL ppic(nlon) 1259 1283 REAL pulow(nlon),pvlow(nlon),pustr(nlon),pvstr(nlon) 1260 1284 REAL t(nlon,nlev), u(nlon,nlev), v(nlon,nlev) 1261 1285 REAL d_t(nlon,nlev), d_u(nlon,nlev), d_v(nlon,nlev) 1262 1286 c 1263 INTEGER i, k, k gwd, kgwdim, kdx(nlon), ktest(nlon)1287 INTEGER i, k, ktest(nlon) 1264 1288 c 1265 1289 c Variables locales: … … 1318 1342 c appeler la routine principale 1319 1343 c 1320 CALL OROLIFT(klon,klev,k gwd,kgwdim,kdx,ktest,1344 CALL OROLIFT(klon,klev,ktest, 1321 1345 . dtime, 1322 . papmh, papmf,zgeom,1346 . papmh, zgeom, 1323 1347 . pt, pu, pv, 1324 . plat,pmea, pstd, p sig, pgam, pthe, ppic,pval,1348 . plat,pmea, pstd, ppic, 1325 1349 . pulow,pvlow, 1326 1350 . pdudt,pdvdt,pdtdt) … … 1341 1365 END 1342 1366 SUBROUTINE OROLIFT( NLON,NLEV 1343 I , K GWD, KGWDIM, KDX, KTEST1367 I , KTEST 1344 1368 R , PTSPHY 1345 R , PAPHM1,P APM1,PGEOM1,PTM1,PUM1,PVM11369 R , PAPHM1,PGEOM1,PTM1,PUM1,PVM1 1346 1370 R , PLAT 1347 R , PMEA, PVAROR, PSIG, PGAMMA, PTHETA,ppic,pval1371 R , PMEA, PVAROR, ppic 1348 1372 C OUTPUTS 1349 1373 R , PULOW,PVLOW … … 1365 1389 C F.LOTT LMD 22/11/95 1366 1390 C 1367 IMPLICIT LOGICAL (L)1391 implicit none 1368 1392 C 1369 1393 C … … 1378 1402 C 1379 1403 C 1404 integer nlon, nlev 1380 1405 REAL PTE(NLON,NLEV), 1381 1406 * PVOL(NLON,NLEV), … … 1387 1412 * PTM1(NLON,NLEV), 1388 1413 * PLAT(NLON),PMEA(NLON), 1389 * PVAROR(NLON), PSIG(NLON),PGAMMA(NLON),1390 * PTHETA(NLON),ppic(NLON),pval(NLON),1414 * PVAROR(NLON), 1415 * ppic(NLON), 1391 1416 * PGEOM1(NLON,NLEV), 1392 * PAPM1(NLON,NLEV),1393 1417 * PAPHM1(NLON,NLEV+1) 1394 1418 C 1395 INTEGER KDX(NLON),KTEST(NLON) 1419 INTEGER KTEST(NLON) 1420 real ptsphy 1396 1421 C----------------------------------------------------------------------- 1397 1422 C 1398 1423 C* 0.2 LOCAL ARRAYS 1399 1424 C ------------ 1400 INTEGER ISECT(KLON), 1401 * ICRIT(KLON), 1402 * IKCRITH(KLON), 1403 * IKenvh(klon), 1425 logical lifthigh, ll1 1426 integer klevm1, jl, ilevh, jk 1427 real zcons1, ztmst, zrtmst,zpi, zhgeo 1428 real zdelp, zslow, zsqua, zscav, zbet 1429 INTEGER 1404 1430 * IKNUB(klon), 1405 * IKNUL(klon), 1406 * IKCRIT(KLON), 1407 * IKHLIM(KLON) 1431 * IKNUL(klon) 1408 1432 LOGICAL LL1(KLON,KLEV+1) 1409 1433 C 1410 1434 REAL ZTAU(KLON,KLEV+1), 1411 1435 * ZTAV(KLON,KLEV+1), 1412 * ZSTAB(KLON,KLEV+1), 1413 * ZVPH(KLON,KLEV+1), 1414 * ZRHO(KLON,KLEV+1), 1415 * ZRI(KLON,KLEV+1), 1416 * ZpsI(KLON,KLEV+1), 1417 * Zzdep(KLON,KLEV) 1436 * ZRHO(KLON,KLEV+1) 1418 1437 REAL ZDUDT(KLON), 1419 * ZDVDT(KLON), 1420 * ZDTDT(KLON), 1421 * ZDEDT(KLON), 1422 * ZVIDIS(KLON), 1423 * ZTFR(KLON), 1424 * Znu(KLON), 1425 * Zd1(KLON), 1426 * Zd2(KLON), 1427 * Zdmod(KLON) 1438 * ZDVDT(KLON) 1428 1439 REAL ZHCRIT(KLON,KLEV) 1429 1440 C----------------------------------------------------------------------- … … 1680 1691 C ORIGINAL : 90-01-01 1681 1692 C ------------------------------------------------------------------ 1682 IMPLICIT LOGICAL(L)1693 implicit none 1683 1694 C 1684 1695 C ----------------------------------------------------------------- … … 1686 1697 C ---------------------------------------------------------------- 1687 1698 C 1688 integer nlon,nlev 1699 integer nlon,nlev, jk 1689 1700 REAL paprs(nlon,nlev+1) 1690 1701 REAL pplay(nlon,nlev) 1702 real zpr,zstra,zsigt,zpm1r 1691 1703 C 1692 1704 C* 1. SET THE VALUES OF THE PARAMETERS … … 1724 1736 C 1725 1737 GKDRAG=0.2 1726 GRAHILO=1 1738 GRAHILO=1. 1727 1739 GRCRIT=0.01 1728 1740 GFRCRIT=1.0
Note: See TracChangeset
for help on using the changeset viewer.