Changeset 4463 for LMDZ6/trunk/libf/phylmd

Timestamp:

Mar 9, 2023, 5:46:26 PM (16 months ago)

Author:

lguez

Message:

Replace stop by call to abort_physiq

Location:

LMDZ6/trunk/libf/phylmd

Files:

• atke_exchange_coeff_mod.F90 (modified) (1 diff)
• phyaqua_mod.F90 (modified) (1 diff)
• thermcell_down.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/atke_exchange_coeff_mod.F90

@@ -137 +137 @@
 ELSE ! TO DO
     
-    print*, 'traitement non-stationnaire de la tke pas encore prevu'
-    stop
+   call abort_physic("atke_compute_km_kh", &
+        'traitement non-stationnaire de la tke pas encore prevu', 1)
 
 END IF

LMDZ6/trunk/libf/phylmd/phyaqua_mod.F90

@@ -135 +135 @@
     if (year_len.ne.360) then
       write (*,*) year_len
-      write (*,*) 'iniaqua: 360 day calendar is required !'
-      stop
+      call abort_physic("iniaqua", 'iniaqua: 360 day calendar is required !', 1)
     endif
 

LMDZ6/trunk/libf/phylmd/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