- Timestamp:
- 12/23/14 19:11:36 (10 years ago)
- Location:
- codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/callcorrk.F90
r313 r315 781 781 782 782 if(fluxtop_dn(ig).lt.0.0)then 783 print*,'Achtung! fluxtop_dn has lost the plot!' 783 print*,'Achtung! fluxtop_dn has lost the plot!',ig 784 784 print*,'fluxtop_dn=',fluxtop_dn(ig) 785 print*,'acosz=',acosz786 print*,'aerosol=',aerosol(ig,:,:)785 !print*,'acosz=',acosz 786 !print*,'aerosol=',aerosol(ig,:,:) 787 787 print*,'temp= ',pt(ig,:) 788 print*,'pplay= ',pplay(ig,:) 789 call abort_physiq 788 !print*,'pplay= ',pplay(ig,:) 789 !call abort_physiq 790 print*,'I am putting results from another point',ig-1 791 fluxtop_dn(ig)=fluxtop_dn(ig-1) 792 fluxtop_lw(ig)=fluxtop_lw(ig-1) 793 fluxabs_sw(ig)=fluxabs_sw(ig-1) 794 fluxsurf_lw(ig)=fluxsurf_lw(ig-1) 795 fluxsurf_sw(ig)=fluxsurf_sw(ig-1) 796 nfluxtopi = fluxtop_lw(ig) 797 fluxtopvdn = fluxtop_dn(ig) 798 nfluxtopv = -fluxabs_sw(ig) 790 799 endif 791 800 -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/physiq.F90
r313 r315 1786 1786 call planetwide_sumval(area(:)*fluxgrd(:)/totarea_planet,GND) 1787 1787 call planetwide_sumval(area(:)*fluxdyn(:)/totarea_planet,DYN) 1788 do ig=1,ngrid1789 if(fluxtop_dn(ig).lt.0.0)then1790 print*,'fluxtop_dn has gone crazy'1791 print*,'fluxtop_dn=',fluxtop_dn(ig)1792 print*,'tau_col=',tau_col(ig)1793 print*,'aerosol=',aerosol(ig,:,:)1794 print*,'temp= ',pt(ig,:)1795 print*,'pplay= ',pplay(ig,:)1796 call abort_physiq1797 endif1798 end do1788 !do ig=1,ngrid 1789 ! if(fluxtop_dn(ig).lt.0.0)then 1790 ! print*,'fluxtop_dn has gone crazy' 1791 ! print*,'fluxtop_dn=',fluxtop_dn(ig) 1792 ! print*,'tau_col=',tau_col(ig) 1793 ! print*,'aerosol=',aerosol(ig,:,:) 1794 ! print*,'temp= ',pt(ig,:) 1795 ! print*,'pplay= ',pplay(ig,:) 1796 ! call abort_physiq 1797 ! endif 1798 !end do 1799 1799 1800 1800 if(ngrid.eq.1)then
Note: See TracChangeset
for help on using the changeset viewer.