Changeset 5099 for LMDZ6/branches/Amaury_dev/libf/phylmdiso/cva_driver.F90
- Timestamp:
- Jul 22, 2024, 9:29:09 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmdiso/cva_driver.F90
r5087 r5099 431 431 ! 14 No moist convection; atmosphere is very 432 432 ! stable (=> no computation) 433 !434 433 435 434 ! ft: Array of temperature tendency (K/s) of dimension ND, defined at same … … 1066 1065 ! If the fraction of convective points is larger than comp_threshold, then compression 1067 1066 ! is assumed useless. 1068 ! 1067 1069 1068 compress = ncum < len*comp_threshold 1070 ! 1069 1071 1070 IF (.not. compress) THEN 1072 1071 DO i = 1,len … … 1074 1073 ENDDO 1075 1074 ENDIF 1076 ! 1075 1077 1076 #ifdef ISO 1078 1077 #ifdef ISOVERIF … … 1480 1479 #endif 1481 1480 ) 1482 ! 1481 1483 1482 IF (prt_level >= 10) THEN 1484 1483 Print *, 'cva_driver after cv3_unsat:mp , water, ice, evap, fondue ' … … 1494 1493 ENDDO 1495 1494 ENDIF 1496 ! 1495 1497 1496 END IF !(iflag_con==3) 1498 1497 … … 1594 1593 #endif 1595 1594 ) 1596 ! 1595 1597 1596 ! Test conseravtion de l'eau 1598 !1599 1597 1600 1598 #ifdef ISOVERIF … … 1628 1626 PRINT *, ' cv3_yield -> fqd(1) = ', fqd(igout, 1) 1629 1627 END IF !(debut) THEN 1630 ! 1628 1631 1629 IF (prt_level >= 10) THEN 1632 1630 Print *, 'cva_driver after cv3_yield:ft(1) , ftd(1) ', & … … 1635 1633 fq(igout,1), fqd(igout,1) 1636 1634 ENDIF 1637 ! 1635 1638 1636 END IF 1639 1637 … … 1747 1745 #endif 1748 1746 1749 !1750 1747 IF (prt_level >= 10) THEN 1751 1748 Print *, 'cva_driver after cv3_uncompress:ft1(1) , ftd1(1) ', & … … 1754 1751 fq1(igout,1), fqd1(igout,1) 1755 1752 ENDIF 1756 ! 1753 1757 1754 END IF 1758 1755 … … 1772 1769 1773 1770 END IF ! ncum>0 1774 ! 1775 ! 1771 1772 1776 1773 DO i = 1,len 1777 1774 IF (iflag1(i) == 14) THEN … … 1781 1778 ENDDO 1782 1779 1783 !1784 1780 ! In order take into account the possibility of changing the compression, 1785 1781 ! reset m, sig and w0 to zero for non-convective points.
Note: See TracChangeset
for help on using the changeset viewer.