- Timestamp:
- Aug 25, 2022, 12:14:04 PM (2 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/wake.F90
r4230 r4231 673 673 END DO 674 674 END DO 675 print*, 'ptop, pupper, ktop, kupper', ptop, pupper, ktop, kupper675 !print*, 'ptop, pupper, ktop, kupper', ptop, pupper, ktop, kupper 676 676 677 677 … … 870 870 CALL pkupper (klon, klev, ptop, ph, pupper, kupper) 871 871 872 print*, 'ptop, pupper, ktop, kupper', ptop, pupper, ktop, kupper872 !print*, 'ptop, pupper, ktop, kupper', ptop, pupper, ktop, kupper 873 873 874 874 ! wk_adv is the logical flag enabling wake evolution in the time advance … … 1151 1151 IF (wk_adv(i) .AND. k<=kupper(i)+1) THEN 1152 1152 dth(i, k) = deltatw(i, k)/ppi(i, k) 1153 1153 ! print *, 'VVVVwake k, the(i,k), dth(i,k), sigmaw(i) ', k, the(i,k), dth(i,k), sigmaw(i) 1154 1154 th1(i, k) = the(i, k) - sigmaw(i)*dth(i, k) ! undisturbed area 1155 1155 th2(i, k) = the(i, k) + (1.-sigmaw(i))*dth(i, k) ! wake
Note: See TracChangeset
for help on using the changeset viewer.