Changeset 3893 for trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F

Timestamp:

Aug 13, 2025, 7:53:49 PM (5 months 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/dyn1d/rcm1d.F (modified) (19 diffs)

trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F

@@ -179 +179 @@
         write(*,*) 'which should contain some input parameters'
         write(*,*) ' ... might as well stop here ...'
-        stop
+        call abort_physic("rcm1d", "Unable to find rcm1d.def", 1)
       else
         close(90)
@@ -195 +195 @@
         write(*,*) 'Please remove the file and restart the run.'
         write(*,*) 'Runtime parameters are supposed to be in rcm1d.def'
-        stop
+        call abort_physic("rcm1d", 
+     &       "A run.def is already here, delete it",1)
       else
         call system('touch run.def')
@@ -241 +242 @@
                     write(*,*) "rcm1d: error reading number of tracers"
                     write(*,*) "   (first line of traceur.def) "
-                    stop
+                    call abort_physic("traceur.def", 
+     &                   "error reading number of tracers",1)
                  ENDIF
                ENDDO
@@ -325 +327 @@
                   write(*,*) "rcm1d: error reading number of tracers"
                   write(*,*) "   (first line of traceur.def) "
-                  stop
+                  call abort_physic("traceur.def", 
+     &                 "error reading number of tracers",1)
                ENDIF
              ENDDO
@@ -331 +334 @@
           nqtot=nq ! set value of nqtot (in infotrac module) as nq
           if (ierr.ne.0) then
-            write(*,*) "rcm1d: error reading number of tracers"
-            write(*,*) "   (first line of traceur.def) "
-            stop
+            call abort_physic("traceur.def", 
+     &           "Unable to read traceur.def",1)
           endif
           if (nq>0) then
@@ -345 +347 @@
             write(*,*) "rcm1d: Error. You set tracer=.true."
             write(*,*) "       but # of tracers in traceur.def is ",nq
-            stop
+            call abort_physic("rcm1d", 
+     &           "tracer=.true. but number of tracers is zero", 1)
           endif
 
@@ -353 +356 @@
             noms(iq)=tname(iq)
             if (ierr.ne.0) then
-              write(*,*) 'rcm1d: error reading tracer names...'
-              stop
+              call abort_physic("traceur.def", 
+     &             "error reading tracer names", 1)
             endif
           enddo !of do iq=1,nq
@@ -374 +377 @@
           write(*,*) ' If you want to run with tracers, I need it'
           write(*,*) ' ... might as well stop here ...'
-          stop
+          call abort_physic("traceur.def", 
+     &         "Unable to find traceur.def", 1)
         endif
         close(90)
@@ -396 +400 @@
           write(*,*) " > If you want tracers, set tracer=.true."
           write(*,*) "------------------------------"
-          stop
+          call abort_physic("rcm1d", "traceur=.false. but nq>0", 1)
        endif
       endif ! of if (tracer)
@@ -425 +429 @@
       IF (g.eq.-99999.) THEN
           PRINT *,"STOP. I NEED g IN RCM1D.DEF."
-          STOP
+          call abort_physic("rcm1d", "g is not specified", 1)
       ELSE
           PRINT *,"--> g = ",g
@@ -436 +440 @@
       IF (rad.eq.-99999. .and. rings_shadow .eqv. .true.) THEN
           PRINT *,"STOP. I NEED rad IN RCM1D.DEF."
-          STOP
+          call abort_physic("rcm1d", "rad is not specified", 1)
       ELSE
           PRINT *,"--> rad = ",rad
@@ -446 +450 @@
       IF (daysec.eq.-99999.) THEN
           PRINT *,"STOP. I NEED daysec IN RCM1D.DEF."
-          STOP
+          call abort_physic("rcm1d", "daysec is not specified", 1)
       ELSE
           PRINT *,"--> daysec = ",daysec
@@ -459 +463 @@
       IF (year_day.eq.-99999.) THEN
           PRINT *,"STOP. I NEED year_day IN RCM1D.DEF."
-          STOP
+          call abort_physic("rcm1d", "year_day is not specified", 1)
       ELSE
           PRINT *,"--> year_day = ",year_day
@@ -469 +473 @@
       IF (periastr.eq.-99999.) THEN
           PRINT *,"STOP. I NEED periastr IN RCM1D.DEF."
-          STOP
+          call abort_physic("rcm1d", "periastr is not specified", 1)
       ELSE
           PRINT *,"--> periastr = ",periastr
@@ -479 +483 @@
       IF (apoastr.eq.-99999.) THEN
           PRINT *,"STOP. I NEED apoastr IN RCM1D.DEF."
-          STOP
+          call abort_physic("rcm1d", "apoastr is not specified", 1)
       ELSE
           PRINT *,"--> apoastr = ",apoastr
@@ -489 +493 @@
       IF (peri_day.eq.-99999.) THEN
           PRINT *,"STOP. I NEED peri_day IN RCM1D.DEF."
-          STOP
+          call abort_physic("rcm1d", "peri_day is not specified", 1)
       ELSE IF (peri_day.gt.year_day) THEN
           PRINT *,"STOP. peri_day.gt.year_day"
-          STOP
+          call abort_physic("rcm1d", 
+     &         "peri_day is greater than year_day", 1)
       ELSE
           PRINT *,"--> peri_day = ", peri_day
@@ -502 +507 @@
       IF (obliquit.eq.-99999.) THEN
           PRINT *,"STOP. I NEED obliquit IN RCM1D.DEF."
-          STOP
+          call abort_physic("rcm1d", "obliquit is not specified", 1)
       ELSE
           PRINT *,"--> obliquit = ",obliquit
@@ -523 +528 @@
         IF (psurf.eq.-99999.) THEN
             PRINT *,"STOP. I NEED psurf IN RCM1D.DEF."
-            STOP
+            call abort_physic("rcm1d", "psurf is not specified", 1)
         ELSE
             PRINT *,"--> psurf = ",psurf
@@ -927 +932 @@
 
       if (autozlevs.and.hybrid) then
-         print*,'Cannot use autozlevs and hybrid together!'
-         call abort
+         call abort_physic("rcm1d", 
+     &        "Cannot use autozlevs and hybrid together", 1)
       endif