Changeset 1625 for trunk/LMDZ.VENUS/libf/phyvenus
- Timestamp:
- Nov 4, 2016, 12:53:36 PM (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.VENUS/libf/phyvenus/cpdet_phy_mod.F90
r1621 r1625 80 80 !---------------------- 81 81 82 ! Ehouan debug 83 write(*,*) "t2tpot: cpp=",cpp," nu_venus=",nu_venus," t0_venus=",t0_venus 84 !write(*,*) "t2tpot: log(ypk/cpp)=",log(ypk/cpp) 85 do k=1,npoints 86 write(*,*) "t2tpot: k=",k," yt(k)=",yt(k) 87 write(*,*) "t2tpot: k=",k," yt(k)**nu_venus=",yt(k)**nu_venus," ypk(k)=",ypk(k) 88 yteta(k) = yt(k)**nu_venus & 89 - nu_venus * t0_venus**nu_venus * log(ypk(k)/cpp) 90 write(*,*) "t2tpot: k=",k," yteta(k)=",yteta(k) 91 enddo 82 yteta = yt**nu_venus & 83 & - nu_venus * t0_venus**nu_venus * log(ypk/cpp) 92 84 yteta = yteta**(1./nu_venus) 93 85 !----------------------
Note: See TracChangeset
for help on using the changeset viewer.