Changeset 1517 for LMDZ5/trunk/libf/phylmd/thermcell_old.F
- Timestamp:
- May 10, 2011, 3:07:35 PM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/thermcell_old.F
r1403 r1517 41 41 c ---------- 42 42 43 INTEGER ngrid,nlay,w2di,tho 43 INTEGER ngrid,nlay,w2di 44 REAL tho 44 45 real ptimestep,l_mix,r_aspect 45 46 REAL pt(ngrid,nlay),pdtadj(ngrid,nlay) … … 339 340 if(w2di.eq.2) then 340 341 entr(ig,k)=entr(ig,k)+ 341 s ptimestep*(zzz-entr(ig,k))/ REAL(tho)342 s ptimestep*(zzz-entr(ig,k))/tho 342 343 else 343 344 entr(ig,k)=zzz … … 590 591 591 592 if (w2di.eq.1) then 592 fm0=fm0+ptimestep*(fm-fm0)/ REAL(tho)593 entr0=entr0+ptimestep*(entr-entr0)/ REAL(tho)593 fm0=fm0+ptimestep*(fm-fm0)/tho 594 entr0=entr0+ptimestep*(entr-entr0)/tho 594 595 else 595 596 fm0=fm … … 846 847 c ---------- 847 848 848 INTEGER ngrid,nlay,w2di,tho 849 INTEGER ngrid,nlay,w2di 850 REAL tho 849 851 real ptimestep,l_mix,r_aspect 850 852 REAL pt(ngrid,nlay),pdtadj(ngrid,nlay) … … 2269 2271 2270 2272 if (w2di.eq.1) then 2271 fm0=fm0+ptimestep*(fm-fm0)/ REAL(tho)2272 entr0=entr0+ptimestep*(alim+entr-entr0)/ REAL(tho)2273 fm0=fm0+ptimestep*(fm-fm0)/tho 2274 entr0=entr0+ptimestep*(alim+entr-entr0)/tho 2273 2275 else 2274 2276 fm0=fm … … 2654 2656 c ---------- 2655 2657 2656 INTEGER ngrid,nlay,w2di,tho 2658 INTEGER ngrid,nlay,w2di 2659 REAL tho 2657 2660 real ptimestep,l_mix,r_aspect 2658 2661 REAL pt(ngrid,nlay),pdtadj(ngrid,nlay) … … 3492 3495 3493 3496 if (w2di.eq.1) then 3494 fm0=fm0+ptimestep*(fm-fm0)/ REAL(tho)3495 entr0=entr0+ptimestep*(entr-entr0)/ REAL(tho)3497 fm0=fm0+ptimestep*(fm-fm0)/tho 3498 entr0=entr0+ptimestep*(entr-entr0)/tho 3496 3499 else 3497 3500 fm0=fm … … 3771 3774 c ---------- 3772 3775 3773 INTEGER ngrid,nlay,w2di,tho 3776 INTEGER ngrid,nlay,w2di 3777 REAL tho 3774 3778 real ptimestep,l_mix,r_aspect 3775 3779 REAL pt(ngrid,nlay),pdtadj(ngrid,nlay) … … 4492 4496 cRC 4493 4497 if (w2di.eq.1) then 4494 fm0=fm0+ptimestep*(fm-fm0)/ REAL(tho)4495 entr0=entr0+ptimestep*(entr-entr0)/ REAL(tho)4498 fm0=fm0+ptimestep*(fm-fm0)/tho 4499 entr0=entr0+ptimestep*(entr-entr0)/tho 4496 4500 else 4497 4501 fm0=fm … … 5184 5188 c ---------- 5185 5189 5186 INTEGER ngrid,nlay,w2di,tho 5190 INTEGER ngrid,nlay,w2di 5191 REAL tho 5187 5192 real ptimestep,l_mix,r_aspect 5188 5193 REAL pt(ngrid,nlay),pdtadj(ngrid,nlay) … … 5924 5929 cRC 5925 5930 if (w2di.eq.1) then 5926 fm0=fm0+ptimestep*(fm-fm0)/ REAL(tho)5927 entr0=entr0+ptimestep*(entr-entr0)/ REAL(tho)5931 fm0=fm0+ptimestep*(fm-fm0)/tho 5932 entr0=entr0+ptimestep*(entr-entr0)/tho 5928 5933 else 5929 5934 fm0=fm
Note: See TracChangeset
for help on using the changeset viewer.