Changeset 315 for codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/physiq.F90

Timestamp:

12/23/14 19:11:36 (10 years ago)

Author:

aslmd

Message:

added a dirty fix to continue integrations even if fluxtop lost the plot. just copy the neighboring point and keep on rocking in the free world.

File:

• codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/physiq.F90 (modified) (1 diff)

codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/physiq.F90

@@ -1786 +1786 @@
          call planetwide_sumval(area(:)*fluxgrd(:)/totarea_planet,GND)
          call planetwide_sumval(area(:)*fluxdyn(:)/totarea_planet,DYN)
-         do ig=1,ngrid
-            if(fluxtop_dn(ig).lt.0.0)then
-               print*,'fluxtop_dn has gone crazy'
-               print*,'fluxtop_dn=',fluxtop_dn(ig)
-               print*,'tau_col=',tau_col(ig)
-               print*,'aerosol=',aerosol(ig,:,:)
-               print*,'temp=   ',pt(ig,:)
-               print*,'pplay=  ',pplay(ig,:)
-               call abort_physiq
-            endif
-         end do
+         !do ig=1,ngrid
+         !   if(fluxtop_dn(ig).lt.0.0)then
+         !      print*,'fluxtop_dn has gone crazy'
+         !      print*,'fluxtop_dn=',fluxtop_dn(ig)
+         !      print*,'tau_col=',tau_col(ig)
+         !      print*,'aerosol=',aerosol(ig,:,:)
+         !      print*,'temp=   ',pt(ig,:)
+         !      print*,'pplay=  ',pplay(ig,:)
+         !      call abort_physiq
+         !   endif
+         !end do
                      
          if(ngrid.eq.1)then