Changeset 315 for codes

Timestamp:

12/23/14 19:11:36 (9 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.

Location:

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

Files:

• callcorrk.F90 (modified) (1 diff)
• physiq.F90 (modified) (1 diff)

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

@@ -781 +781 @@
 
          if(fluxtop_dn(ig).lt.0.0)then
-            print*,'Achtung! fluxtop_dn has lost the plot!'
+            print*,'Achtung! fluxtop_dn has lost the plot!',ig
             print*,'fluxtop_dn=',fluxtop_dn(ig)
-            print*,'acosz=',acosz
-            print*,'aerosol=',aerosol(ig,:,:)
+            !print*,'acosz=',acosz
+            !print*,'aerosol=',aerosol(ig,:,:)
             print*,'temp=   ',pt(ig,:)
-            print*,'pplay=  ',pplay(ig,:)
-            call abort_physiq
+            !print*,'pplay=  ',pplay(ig,:)
+            !call abort_physiq
+            print*,'I am putting results from another point',ig-1
+            fluxtop_dn(ig)=fluxtop_dn(ig-1)
+            fluxtop_lw(ig)=fluxtop_lw(ig-1)
+            fluxabs_sw(ig)=fluxabs_sw(ig-1)
+            fluxsurf_lw(ig)=fluxsurf_lw(ig-1)
+            fluxsurf_sw(ig)=fluxsurf_sw(ig-1)
+            nfluxtopi = fluxtop_lw(ig)
+            fluxtopvdn = fluxtop_dn(ig)
+            nfluxtopv = -fluxabs_sw(ig)
          endif
 

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