Changeset 3893 for trunk/LMDZ.GENERIC/libf/phystd/soil_settings.F

Timestamp:

Aug 13, 2025, 7:53:49 PM (7 days ago)

Author:

gmilcareck

Message:

Remove all "call abort" and "stop" statement in LMDZ.GENERIC and replacing them by call abort_physic().

File:

• trunk/LMDZ.GENERIC/libf/phystd/soil_settings.F (modified) (7 diffs)

trunk/LMDZ.GENERIC/libf/phystd/soil_settings.F

@@ -92 +92 @@
         allocate(oldmlayer(dimlen),stat=ierr)
         if (ierr.ne.0) then
-          write(*,*) 'soil_settings: failed allocation of oldmlayer!'
-          stop
+          call abort_physic("soil_settings",
+     &         "failed allocation of oldmlayer",1)
         endif
 
@@ -128 +128 @@
           allocate(oldmlayer(dimlen),stat=ierr)
           if (ierr.ne.0) then
-            write(*,*) 'soil_settings: failed allocation of oldmlayer!'
-            stop
+            call abort_physic("soil_settings",
+     &           "failed allocation of oldmlayer",1)
           endif
         endif
@@ -167 +167 @@
          write(*,*) "soil_settings error: mlay1_soil=",mlay1_soil
          write(*,*) " does not match comsoil_h lay1_soil=",lay1_soil
-         stop
+         call abort_physic("soil_settings",
+     &        "mlay1_soil does not match lay1_soil",1)
       endif
       if ((abs(malpha-alpha_soil)/alpha_soil).gt.1.e-6) then
          write(*,*) "soil_settings error: malpha=",malpha
          write(*,*) " does not match comsoil_h alpha_soil=",alpha_soil
-         stop
+         call abort_physic("soil_settings",
+     &        "malpha does not match alpha_soil",1)
       endif
       do iloop=1,nsoil
@@ -207 +209 @@
        call get_field(nid,"inertiedat",surfinertia,found)
        if (.not.found) then
-         write(*,*) "soil_settings: Failed loading <inertiedat>"
-         call abort
+         call abort_physic("soil_settings",
+     &        "Failed loading <inertiedat>", 1)
        endif
        
@@ -223 +225 @@
            allocate(oldinertiedat(ngrid,dimlen),stat=ierr)
            if (ierr.ne.0) then
-            write(*,*) 'soil_settings: failed allocation of ',
-     &                 'oldinertiedat!'
-            stop
+            call abort_physic("soil_settings",
+     &           "failed allocation of oldinertiedat!",1)
            endif
          endif ! of if (.not.allocated(oldinertiedat))
         call get_field(nid,"inertiedat",oldinertiedat,found)
         if (.not.found) then
-          write(*,*) "soil_settings: Failed loading <inertiedat>"
-          call abort
+          call abort_physic("soil_settings",
+     &         "Failed loading <inertiedat>", 1)
         endif
        else ! put values in therm_i
          call get_field(nid,"inertiedat",inertiedat,found)
          if (.not.found) then
-           write(*,*) "soil_settings: Failed loading <inertiedat>"
-           call abort
+           call abort_physic("soil_settings", 
+     &          "Failed loading <inertiedat>", 1)
          endif
 !        endif
@@ -259 +260 @@
            allocate(oldtsoil(ngrid,dimlen),stat=ierr)
            if (ierr.ne.0) then
-             write(*,*) 'soil_settings: failed allocation of ',
-     &                  'oldtsoil!'
-             stop
+             call abort_physic("soil_settings",
+     &            "failed allocation of oldtsoil",1)
            endif
          endif
          call get_field(nid,"tsoil",oldtsoil,found)
          if (.not.found) then
-           write(*,*) "soil_settings: Failed loading <tsoil>"
-           call abort
+           call abort_physic("soil_settings", 
+     &          "Failed loading <tsoil>", 1)
          endif
        else ! put values in tsoil
@@ -273 +273 @@
      &                  timeindex=indextime)
          if (.not.found) then
-           write(*,*) "soil_settings: Failed loading <tsoil>"
-           call abort
+           call abort_physic("soil_settings", 
+     &          "Failed loading <tsoil>", 1)
          endif
        endif ! of if (interpol)