Changeset 4585 for LMDZ6/branches/Portage_acc/libf/phylmdiso/thermcell_down.F90

Timestamp:

Jun 23, 2023, 6:03:57 PM (15 months ago)

Author:

Laurent Fairhead

Message:

Adding some files that didn't make it on the last commit ("Node filename has unexpectedly changed kind" error)

File:

• LMDZ6/branches/Portage_acc/libf/phylmdiso/thermcell_down.F90 (modified) (5 diffs)

LMDZ6/branches/Portage_acc/libf/phylmdiso/thermcell_down.F90

@@ -58 +58 @@
 
    if ( iflag_thermals_down < 10 ) then
-        stop 'thermcell_down_dq = 0 or >= 10'
+      call abort_physic("thermcell_updown_dq", &
+           'thermcell_down_dq = 0 or >= 10', 1)
    else
         iflag_impl=iflag_thermals_down-10
@@ -128 +129 @@
          !!!! tentative de prise en compte d'un flux compensatoire montant  !!!!
          if (fup(ig,ilay)-fdn(ig,ilay) .lt. 0.) then
-            write(*,*) 'flux compensatoire montant, cas non traite par thermcell_updown_dq'
-            stop          
+            call abort_physic("thermcell_updown_dq", 'flux compensatoire '&
+                 // 'montant, cas non traite par thermcell_updown_dq', 1)
             !fthe(ig,ilay)=(fup(ig,ilay)-fdn(ig,ilay))*trac(ig,ilay-1)
          else
@@ -161 +162 @@
        do ig=1,ngrid
          if (fup(ig,ilay)-fdn(ig,ilay) .lt. 0.) then
-            write(*,*) 'flux compensatoire montant, cas non traite par thermcell_updown_dq'
-            stop          
+            call abort_physic("thermcell_updown_dq", 'flux compensatoire ' &
+                 // 'montant, cas non traite par thermcell_updown_dq', 1)
          else
             fthe(ig,ilay)=-(fup(ig,ilay)-fdn(ig,ilay))*trac(ig,ilay)
@@ -191 +192 @@
          if((fup(ig,ilay)-fdn(ig,ilay)) .lt. 0) then
             write(*,*) 'flux compensatoire montant, cas non traite par thermcell_updown_dq dans le cas d une resolution implicite, ilay : ', ilay
-            stop  
+            call abort_physic("thermcell_updown_dq", "", 1)
          else
            mstar_inv=ptimestep/masse(ig,ilay)
@@ -205 +206 @@
 
    else
-      write(*,*) 'valeur de iflag_impl non prevue'
-      stop
-
+      call abort_physic("thermcell_updown_dq", &
+           'valeur de iflag_impl non prevue', 1)
    endif