Changeset 1923
- Timestamp:
- Jan 6, 2014, 5:07:02 PM (11 years ago)
- File:
-
- 1 edited
-
LMDZ5/trunk/libf/phylmd/conflx.F (modified) (29 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/conflx.F
r1907 r1923 18 18 c et lessivage des traceurs passifs. 19 19 c====================================================================== 20 c ym#include "dimensions.h"21 c ym#include "dimphy.h"22 #include "YOMCST.h"23 #include "YOETHF.h"20 c yminclude "dimensions.h" 21 c yminclude "dimphy.h" 22 include "YOMCST.h" 23 include "YOETHF.h" 24 24 c Entree: 25 25 REAL dtime ! pas d'integration (s) … … 77 77 REAL zdelta, zqsat 78 78 c 79 #include "FCTTRE.h"79 include "FCTTRE.h" 80 80 c 81 81 c initialiser les variables de sortie (pour securite) … … 207 207 208 208 RETURN 209 END 209 END subroutine conflx 210 210 c-------------------------------------------------------------------- 211 211 SUBROUTINE flxmain(pdtime, pten, pqen, pqsen, pqhfl, pap, paph, … … 218 218 IMPLICIT none 219 219 C ------------------------------------------------------------------ 220 c ym#include "dimensions.h"221 c ym#include "dimphy.h"222 #include "YOMCST.h"223 #include "YOETHF.h"224 #include "YOECUMF.h"220 c yminclude "dimensions.h" 221 c yminclude "dimphy.h" 222 include "YOMCST.h" 223 include "YOETHF.h" 224 include "YOECUMF.h" 225 225 C ---------------------------------------------------------------- 226 226 REAL pten(klon,klev), pqen(klon,klev), pqsen(klon,klev) … … 502 502 1000 CONTINUE 503 503 RETURN 504 END 504 END SUBROUTINE flxmain 505 505 SUBROUTINE flxini(pten, pqen, pqsen, pgeo, paph, pgeoh, ptenh, 506 506 . pqenh, pqsenh, ptu, pqu, ptd, pqd, pmfd, pmfds, pmfdq, … … 514 514 C AND INITIALIZES VALUES FOR UPDRAFTS 515 515 C---------------------------------------------------------------------- 516 c ym#include "dimensions.h"517 c ym#include "dimphy.h"518 #include "YOMCST.h"519 #include "YOETHF.h"516 c yminclude "dimensions.h" 517 c yminclude "dimphy.h" 518 include "YOMCST.h" 519 include "YOETHF.h" 520 520 C 521 521 REAL pten(klon,klev) ! temperature (environnement) … … 630 630 C 631 631 RETURN 632 END 632 END SUBROUTINE flxini 633 633 SUBROUTINE flxbase(ptenh, pqenh, pgeoh, paph, 634 634 * ptu, pqu, plu, ldcum, kcbot, klab) … … 646 646 C (NON ENTRAINING PLUME,I.E.CONSTANT MASSFLUX) 647 647 C---------------------------------------------------------------------- 648 c ym#include "dimensions.h"649 c ym#include "dimphy.h"650 #include "YOMCST.h"651 #include "YOETHF.h"648 c yminclude "dimensions.h" 649 c yminclude "dimphy.h" 650 include "YOMCST.h" 651 include "YOETHF.h" 652 652 C ---------------------------------------------------------------- 653 653 REAL ptenh(klon,klev), pqenh(klon,klev) … … 712 712 c 713 713 RETURN 714 END 714 END SUBROUTINE flxbase 715 715 SUBROUTINE flxasc(pdtime, ptenh, pqenh, pten, pqen, pqsen, 716 716 . pgeo, pgeoh, pap, paph, pqte, pvervel, … … 725 725 C FOR CUMULUS PARAMETERIZATION 726 726 C---------------------------------------------------------------------- 727 c ym#include "dimensions.h"728 c ym#include "dimphy.h"729 #include "YOMCST.h"730 #include "YOETHF.h"731 #include "YOECUMF.h"727 c yminclude "dimensions.h" 728 c yminclude "dimphy.h" 729 include "YOMCST.h" 730 include "YOETHF.h" 731 include "YOECUMF.h" 732 732 C 733 733 REAL pdtime … … 763 763 REAL zwmax(klon), zzzmb 764 764 INTEGER klwmin(klon) ! level of maximum vertical velocity 765 real fact 765 766 C---------------------------------------------------------------------- 766 767 ztglace = RTT - 13. … … 835 836 DO 480 k = klev-1,3,-1 836 837 c 837 IF (LMFMID .AND. k.LT.klev-1 .AND. k.GT.klev/2) THEN838 IF (LMFMID .AND. k.LT.klev-1) THEN 838 839 DO i = 1, klon 839 840 IF (.NOT.ldcum(i) .AND. klab(i,k+1).EQ.0 .AND. 840 . pqen(i,k).GT.0.9*pqsen(i,k)) THEN 841 . pqen(i,k).GT.0.9*pqsen(i,k) 842 $ .and. pap(i, k) / paph(i, klev + 1)> 0.4) THEN 841 843 ptu(i,k+1) = pten(i,k) +(pgeo(i,k)-pgeoh(i,k+1))/RCPD 842 844 pqu(i,k+1) = pqen(i,k) … … 893 895 llo1=pen_u(i,k).GT.0..AND.(ktype(i).EQ.1.OR.ktype(i).EQ.2) 894 896 IF(llo1) THEN 895 zentr=zentr*(1.+3.*(1.-MIN(1.,(zpbot(i)-pap(i,k))/1.5E4))) 896 pen_u(i,k)=pen_u(i,k)*(1.+3.*(1.-MIN(1., 897 . (zpbot(i)-pap(i,k))/1.5E4))) 898 pde_u(i,k)=pde_u(i,k)*(1.+3.*(1.-MIN(1., 899 . (zpbot(i)-pap(i,k))/1.5E4))) 897 fact=1.+3.*(1.-MIN(1.,(zpbot(i)-pap(i,k))/1.5E4)) 898 zentr=zentr*fact 899 pen_u(i,k)=pen_u(i,k)*fact 900 pde_u(i,k)=pde_u(i,k)*fact 900 901 ENDIF 901 902 IF(llo2.AND.pqenh(i,k+1).GT.1.E-5) … … 1017 1018 800 CONTINUE 1018 1019 RETURN 1019 END 1020 END SUBROUTINE flxasc 1020 1021 SUBROUTINE flxflux(pdtime, pqen, pqsen, ptenh, pqenh, pap 1021 1022 . , paph, ldland, pgeoh, kcbot, kctop, lddraf, kdtop … … 1030 1031 C FLUXES IN THE CLOUD LAYER AND IN THE SUBCLOUD LAYER 1031 1032 C---------------------------------------------------------------------- 1032 cym#include "dimensions.h" 1033 cym#include "dimphy.h" 1034 #include "YOMCST.h" 1035 #include "YOETHF.h" 1036 #include "YOECUMF.h" 1033 c yminclude "dimensions.h" 1034 c yminclude "dimphy.h" 1035 include "YOMCST.h" 1036 include "YOETHF.h" 1037 include "YOECUMF.h" 1038 include "iniprint.h" 1037 1039 C 1038 1040 REAL cevapcu(klon,klev) … … 1068 1070 INTEGER kdtop(klon) 1069 1071 c 1070 #include "FCTTRE.h"1072 include "FCTTRE.h" 1071 1073 c 1072 1074 DO 101 k=1,klev … … 1255 1257 pdpmel(i,k) = 0.0 1256 1258 ENDIF 1257 if (pmflxr(i,k) + pmflxs(i,k) .lt.-1.e-26)1259 if (pmflxr(i,k) + pmflxs(i,k) < -1.e-26 .and. prt_level >= 1) 1258 1260 . write(*,*) 'precip. < 1e-16 ',pmflxr(i,k) + pmflxs(i,k) 1259 1261 ENDIF … … 1267 1269 c 1268 1270 RETURN 1269 END 1271 END SUBROUTINE flxflux 1270 1272 SUBROUTINE flxdtdq(pdtime, ktopm2, paph, ldcum, pten 1271 1273 . , pmfus, pmfds, pmfuq, pmfdq, pmful, pdmfup, pdmfdp … … 1276 1278 c calculer les tendances T et Q 1277 1279 c---------------------------------------------------------------------- 1278 c ym#include "dimensions.h"1279 c ym#include "dimphy.h"1280 #include "YOMCST.h"1281 #include "YOETHF.h"1282 #include "YOECUMF.h"1280 c yminclude "dimensions.h" 1281 c yminclude "dimphy.h" 1282 include "YOMCST.h" 1283 include "YOETHF.h" 1284 include "YOECUMF.h" 1283 1285 C ----------------------------------------------------------------- 1284 1286 LOGICAL llo1 … … 1339 1341 C 1340 1342 RETURN 1341 END 1343 END SUBROUTINE flxdtdq 1342 1344 SUBROUTINE flxdlfs(ptenh, pqenh, pgeoh, paph, ptu, pqu, 1343 1345 . ldcum, kcbot, kctop, pmfub, prfl, ptd, pqd, … … 1361 1363 C MOIST ENVIRONMENTAL AIR AND CLOUD AIR. 1362 1364 C---------------------------------------------------------------------- 1363 c ym#include "dimensions.h"1364 c ym#include "dimphy.h"1365 #include "YOMCST.h"1366 #include "YOETHF.h"1367 #include "YOECUMF.h"1365 c yminclude "dimensions.h" 1366 c yminclude "dimphy.h" 1367 include "YOMCST.h" 1368 include "YOETHF.h" 1369 include "YOECUMF.h" 1368 1370 C 1369 1371 REAL ptenh(klon,klev) … … 1451 1453 C 1452 1454 RETURN 1453 END 1455 END SUBROUTINE flxdlfs 1454 1456 SUBROUTINE flxddraf(ptenh, pqenh, pgeoh, paph, prfl, 1455 1457 . ptd, pqd, pmfd, pmfds, pmfdq, pdmfdp, … … 1473 1475 C 1474 1476 C---------------------------------------------------------------------- 1475 c ym#include "dimensions.h"1476 c ym#include "dimphy.h"1477 #include "YOMCST.h"1478 #include "YOETHF.h"1479 #include "YOECUMF.h"1477 c yminclude "dimensions.h" 1478 c yminclude "dimphy.h" 1479 include "YOMCST.h" 1480 include "YOETHF.h" 1481 include "YOECUMF.h" 1480 1482 C 1481 1483 REAL ptenh(klon,klev), pqenh(klon,klev) … … 1569 1571 180 CONTINUE 1570 1572 RETURN 1571 END 1573 END SUBROUTINE flxddraf 1572 1574 SUBROUTINE flxadjtq(pp, pt, pq, ldflag, kcall) 1573 1575 USE dimphy … … 1581 1583 C kcall=2 EVAPORATION IN DOWNDRAFTS (E.G. CUDLFS,CUDDRAF) 1582 1584 C 1583 c ym#include "dimensions.h"1584 c ym#include "dimphy.h"1585 #include "YOMCST.h"1585 c yminclude "dimensions.h" 1586 c yminclude "dimphy.h" 1587 include "YOMCST.h" 1586 1588 C 1587 1589 REAL pt(klon), pq(klon), pp(klon) … … 1593 1595 REAL zdelta, zcvm5, zldcp, zqsat, zcor 1594 1596 INTEGER is, i 1595 #include "YOETHF.h"1596 #include "FCTTRE.h"1597 include "YOETHF.h" 1598 include "FCTTRE.h" 1597 1599 C 1598 1600 z5alvcp = r5les*RLVTT/RCPD … … 1648 1650 230 CONTINUE 1649 1651 RETURN 1650 END 1652 END SUBROUTINE flxadjtq 1651 1653 SUBROUTINE flxsetup 1652 1654 IMPLICIT none … … 1654 1656 C THIS ROUTINE DEFINES DISPOSABLE PARAMETERS FOR MASSFLUX SCHEME 1655 1657 C 1656 #include "YOECUMF.h"1658 include "YOECUMF.h" 1657 1659 C 1658 1660 ENTRPEN=1.0E-4 ! ENTRAINMENT RATE FOR PENETRATIVE CONVECTION … … 1674 1676 c 1675 1677 RETURN 1676 END 1678 END SUBROUTINE flxsetup
Note: See TracChangeset
for help on using the changeset viewer.
