Changeset 4470 for LMDZ6/trunk/libf/dyn3d

Timestamp:

Mar 10, 2023, 5:55:53 PM (2 years ago)

Author:

Laurent Fairhead

Message:

Replaced STOP instructions by calls to abort_gcm for a cleaner exit

Location:

LMDZ6/trunk/libf/dyn3d

Files:

• bilan_dyn.F (modified) (1 diff)
• conf_gcm.F90 (modified) (15 diffs)
• groupe.F (modified) (1 diff)
• guide_mod.F90 (modified) (5 diffs)
• vlsplt.F (modified) (2 diffs)
• vlspltqs.F (modified) (1 diff)
• wrgrads.F (modified) (1 diff)

LMDZ6/trunk/libf/dyn3d/bilan_dyn.F

@@ -185 +185 @@
            WRITE(lunout,*)'dt_app=',dt_app
            WRITE(lunout,*)'dt_cum=',dt_cum
-           stop
+           call abort_gcm('bilan_dyn','stopped',1)
         endif
 

LMDZ6/trunk/libf/dyn3d/conf_gcm.F90

@@ -415 +415 @@
         write(lunout,*)'conf_gcm: La valeur de clat passee par run.def', &
              ' est differente de celle lue sur le fichier  start '
-        STOP
+        CALL abort_gcm("conf_gcm","stopped",1)
      ENDIF
 
@@ -429 +429 @@
         write(lunout,*)'conf_gcm: La valeur de grossismx passee par ', &
              'run.def est differente de celle lue sur le fichier  start '
-        STOP
+        CALL abort_gcm("conf_gcm","stopped",1)
      ENDIF
 
@@ -443 +443 @@
         write(lunout,*)'conf_gcm: La valeur de grossismy passee par ', &
              'run.def est differente de celle lue sur le fichier  start '
-        STOP
+        CALL abort_gcm("conf_gcm","stopped",1)
      ENDIF
 
@@ -449 +449 @@
         write(lunout,*) &
              'conf_gcm: ***  ATTENTION !! grossismx < 1 .   *** '
-        STOP
+        CALL abort_gcm("conf_gcm","stopped",1)
      ELSE
         alphax = 1. - 1./ grossismx
@@ -457 +457 @@
         write(lunout,*) &
              'conf_gcm: ***  ATTENTION !! grossismy < 1 .   *** '
-        STOP
+        CALL abort_gcm("conf_gcm","stopped",1)
      ELSE
         alphay = 1. - 1./ grossismy
@@ -479 +479 @@
            write(lunout,*)' *** fxyhypb lu sur le fichier start est ', &
                 'F alors  qu il est  T  sur  run.def  ***'
-           STOP
+           CALL abort_gcm("conf_gcm","stopped",1)
         ENDIF
      ELSE
@@ -486 +486 @@
            write(lunout,*)' ***  fxyhypb lu sur le fichier start est ', &
                 'T alors  qu il est  F  sur  run.def  ****  '
-           STOP
+           CALL abort_gcm("conf_gcm","stopped",1)
         ENDIF
      ENDIF
@@ -502 +502 @@
            write(lunout,*)'conf_gcm: La valeur de dzoomx passee par ', &
                 'run.def est differente de celle lue sur le fichier  start '
-           STOP
+           CALL abort_gcm("conf_gcm","stopped",1)
         ENDIF
      ENDIF
@@ -518 +518 @@
            write(lunout,*)'conf_gcm: La valeur de dzoomy passee par ', &
                 'run.def est differente de celle lue sur le fichier  start '
-           STOP
+           CALL abort_gcm("conf_gcm","stopped",1)
         ENDIF
      ENDIF
@@ -533 +533 @@
            write(lunout,*)'conf_gcm: La valeur de taux passee par ', &
                 'run.def est differente de celle lue sur le fichier  start '
-           STOP
+           CALL abort_gcm("conf_gcm","stopped",1)
         ENDIF
      ENDIF
@@ -548 +548 @@
            write(lunout,*)'conf_gcm: La valeur de tauy passee par ', &
                 'run.def est differente de celle lue sur le fichier  start '
-           STOP
+           CALL abort_gcm("conf_gcm","stopped",1)
         ENDIF
      ENDIF
@@ -569 +569 @@
               write(lunout,*)' *** ysinus lu sur le fichier start est F', &
                    ' alors  qu il est  T  sur  run.def  ***'
-              STOP
+              CALL abort_gcm("conf_gcm","stopped",1)
            ENDIF
         ELSE
@@ -576 +576 @@
               write(lunout,*)' *** ysinus lu sur le fichier start est T', &
                    ' alors  qu il est  F  sur  run.def  ****  '
-              STOP
+              CALL abort_gcm("conf_gcm","stopped",1)
            ENDIF
         ENDIF
@@ -701 +701 @@
         write(lunout,*) &
              'conf_gcm: ***  ATTENTION !! grossismx < 1 .   *** '
-        STOP
+        CALL abort_gcm("conf_gcm","stopped",1)
      ELSE
         alphax = 1. - 1./ grossismx
@@ -708 +708 @@
      IF( grossismy.LT.1. )  THEN
         write(lunout,*) 'conf_gcm: ***ATTENTION !! grossismy < 1 . *** '
-        STOP
+        CALL abort_gcm("conf_gcm","stopped",1)
      ELSE
         alphay = 1. - 1./ grossismy

LMDZ6/trunk/libf/dyn3d/groupe.F

@@ -51 +51 @@
       if (firstcall) then
          if (groupe_ok) then
-           if(mod(iim,2**ngroup).ne.0) stop'probleme du nombre de point'
+            if(mod(iim,2**ngroup).ne.0)         
+     &        CALL abort_gcm('groupe','probleme du nombre de point',1)
          endif
          firstcall=.false.

LMDZ6/trunk/libf/dyn3d/guide_mod.F90

@@ -359 +359 @@
     INCLUDE "dimensions.h"
     INCLUDE "paramet.h"
+    INCLUDE "iniprint.h"
+
 
     ! Variables entree
@@ -382 +384 @@
     INTEGER       :: l
     CHARACTER(LEN=20) :: modname="guide_main"
+    CHARACTER (len = 80)   :: abort_message
+
 
 !-----------------------------------------------------------------------
@@ -447 +451 @@
       IF (reste.EQ.0.) THEN
           IF (itau_test.EQ.itau) THEN
-            write(*,*)trim(modname)//' second pass in advreel at itau=',&
+            write(lunout,*)trim(modname)//' second pass in advreel at itau=',&
             itau
-            stop
+              abort_message='stopped'
+              CALL abort_gcm(modname,abort_message,1) 
           ELSE
               IF (guide_v) vnat1=vnat2
@@ -942 +947 @@
     integer i,j,ilon,ilat
     character(len=20),parameter :: modname="tau2alpha"
+    CHARACTER (len = 80)   :: abort_message
 
 
@@ -1038 +1044 @@
               if (gamma.lt.1.e-5) then
                 write(*,*)trim(modname)//' gamma =',gamma,'<1e-5'
-                stop
+                abort_message='stopped'
+                CALL abort_gcm(modname,abort_message,1)
               endif
               gamma=log(0.5)/log(gamma)

LMDZ6/trunk/libf/dyn3d/vlsplt.F

@@ -184 +184 @@
             DO ij=iip2,ip1jm-1
                dxqu(ij)=q(ij+1,l,iq)-q(ij,l,iq)
-c              IF(u_m(ij,l).lt.0.) stop'limx n admet pas les U<0'
-c              sigu(ij)=u_m(ij,l)/masse(ij,l,iq)
             ENDDO
             DO ij=iip1+iip1,ip1jm,iip1
@@ -310 +308 @@
       ENDDO
 #endif
-c       stop
 
 c       go to 9999

LMDZ6/trunk/libf/dyn3d/vlspltqs.F

@@ -248 +248 @@
             DO ij=iip2,ip1jm-1
                dxqu(ij)=q(ij+1,l,iq)-q(ij,l,iq)
-c              IF(u_m(ij,l).lt.0.) stop'limx n admet pas les U<0'
-c              sigu(ij)=u_m(ij,l)/masse(ij,l,iq)
             ENDDO
             DO ij=iip1+iip1,ip1jm,iip1

LMDZ6/trunk/libf/dyn3d/wrgrads.F

@@ -80 +80 @@
            print*,'nvar  ',nvar(if)
            print*,'vars ',(var(iv,if),iv=1,nvar(if))
-
-           stop
+           CALL abort_gcm("wrgrads","problem",1)
          endif
       endif