Changeset 2502 for LMDZ5/trunk/libf/phylmd
- Timestamp:
- May 4, 2016, 11:51:01 AM (9 years ago)
- Location:
- LMDZ5/trunk/libf/phylmd
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/cv3p1_closure.F90
r2458 r2502 113 113 IF (prt_level>=20) PRINT *, 'cv3p1_param nloc ncum nd icb inb nl', nloc, & 114 114 ncum, nd, icb(nloc), inb(nloc), nl 115 DO k = 1, n l115 DO k = 1, nd !jyg: initialization up to nd 116 116 DO il = 1, ncum 117 117 m(il, k) = 0.0 … … 120 120 121 121 !CR: initializations for erosion of adiabatic ascent 122 DO k = 1,n l122 DO k = 1,nd !jyg: initialization up to nd 123 123 DO il = 1, ncum 124 124 mad(il,k)=0. … … 628 628 629 629 !Verification sum(me)=sum(m) 630 DO k = 1,n l+1630 DO k = 1,nd !jyg: initialization up to nd 631 631 DO il = 1, ncum 632 632 md(il,k)=0. -
LMDZ5/trunk/libf/phylmd/cv3p2_closure.F90
r2458 r2502 126 126 DO k = 1, nl 127 127 DO il = 1, ncum 128 m(il, k) = 0.0129 128 rhodp(il,k) = 0.007*p(il, k)*(ph(il,k)-ph(il,k+1))/tv(il, k) 130 129 END DO 131 130 END DO 132 131 133 !CR : initializations for erosion of adiabatic ascent134 DO k = 1,n l132 !CR+jyg: initializations (up to nd) for erosion of adiabatic ascent and of m and wlim 133 DO k = 1,nd 135 134 DO il = 1, ncum 136 135 mad(il,k)=0. … … 138 137 betalim(il,k)=1. 139 138 wlim(il,k)=0. 139 m(il, k) = 0.0 140 140 ENDDO 141 141 ENDDO … … 711 711 712 712 !Verification sum(me)=sum(m) 713 DO k = 1,n l+1713 DO k = 1,nd 714 714 DO il = 1, ncum 715 715 md(il,k)=0.
Note: See TracChangeset
for help on using the changeset viewer.