Changeset 3893 for trunk/LMDZ.GENERIC/libf/aeronostd/photochemistry_asis.F90

Timestamp:

Aug 13, 2025, 7:53:49 PM (4 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/aeronostd/photochemistry_asis.F90 (modified) (19 diffs)

trunk/LMDZ.GENERIC/libf/aeronostd/photochemistry_asis.F90

@@ -132 +132 @@
        print*,'set jonline=.false. in callphys.def'
        print*,'or set photolysis reactions in monoreact file'
-       call abort
+       call abort_physic("photochemistry","Jonline cannot run without photolysis",1)
      end if
      call init_photolysis(nlayer,nreact)     ! on-line photolysis
@@ -296 +296 @@
    call dgesv(nesp,1,mat,nesp,indx,cnew,nesp,code)
 #else
-   write(*,*) "photochemistry_asis error, missing LAPACK routine dgesv"
-   call abort
+   call abort_physic("photochemistry","missing LAPACK routine dgesv",1)
 #endif
 
@@ -412 +411 @@
       call dgesv(nesp,1,mat,nesp,indx,cnew,nesp,code)
 #else
-   write(*,*) "photochemistry_asis error, missing LAPACK routine dgesv"
-   call abort
+   call abort_physic("photochemistry","missing LAPACK routine dgesv",1)
 #endif
 
@@ -568 +566 @@
     print*, 'Please check the reaction ',ireact
     if (ireact>nreact) print*, 'Please check the reaction count into the code beacause ireact > nreact is not possible'
-    call abort
+    call abort_physic("photochemistry","ireact>nreact is not possible",1)
   end if
 
@@ -599 +597 @@
     print*, 'Please check the reaction ',ireact
     if (ireact>nreact) print*, 'Please check the reaction count into the code beacause ireact > nreact is not possible'
-    call abort
+    call abort_physic("photochemistry","ireact>nreact is not possible",1)
   end if
 
@@ -625 +623 @@
     print*, 'n4_hard_coding    = ', n4_hard_coding
     print*, 'n3_hard_coding    = ', n3_hard_coding
-    call abort
+    call abort_physic("photochemistry","wrong dimensions in reactionrates",1)
  end if
  if ((nb_hv /= nb_hv_max)         .or.  &
@@ -640 +638 @@
     print*, 'nb_pfunc2_max = ', nb_pfunc2_max
     print*, 'nb_pfunc3_max = ', nb_pfunc3_max
-    call abort
+    call abort_physic("photochemistry","wrong dimensions in reaction rates",1)
  end if
 firstcall = .false.
@@ -913 +911 @@
       print*, 'WARNING: ', chemnoms(iq),' is only destroy'
     else
-      print*, 'Error in indice: bad value in specheckr or specheckp'
-      call abort
+      call abort_physic("photochemistry", "Error in indice: bad value in specheckr or specheckp",1)
     end if
   end if
@@ -950 +947 @@
    print*, 'nb_phot_hv_max    = ', nb_phot_hv_max
    print*, 'nb_hv_max         = ', nb_hv_max
-   call abort
+   call abort_physic("photochemistry","wrong dimensions in indice",1)
 end if  
 
@@ -1079 +1076 @@
             n = n + 1
             if (n>nmax) then
-              print*,'Error in split_str: to much words'
+              print*,'Error in split_str: too much words'
               print*,'limit = ',nmax
               print*,'change it, if you want, with increasing nmax'
               print*,'line is:',line
-              call abort
+              call abort_physic("photochemistry","too much words in split_str",1)
             end if
             read(line,*,iostat=ios) buf(1:n)  ! use list-directed input
@@ -1141 +1138 @@
                 vtype = "v4"
             else
-                print*,'Error in photochemistry_asis (find_vtype):'
-                print*,'3 different reactants not implemented'
-                call abort
+                call abort_physic("photochemistry","(find_vtype): 3 different reactants not implemented",1)
             end if
         else
-            print*,'Error in photochemistry_asis (find_vtype):'
-            print*,'more than 2 different reactants not implemented'
-            call abort
+            call abort_physic("photochemistry","(find_vtype): more than 2 different reactants not implemented",1)
         end if
 
@@ -1192 +1185 @@
            write(*,*)'28/03/2025: New default value is network.def (in simulation folder)'
            write(*,*)'Former default is chemnetwork/reactfile'
-           call abort
+           call abort_physic("photochemistry","Unable to open chemical reaction file",1)
         end if
     
@@ -1227 +1220 @@
                           .and. trim(words(2))/=trim(words(3))) nb_reaction_3 = nb_reaction_3 + 1
             else
-                print*,'Error in photochemistry_asis (count_react):'
-                print*,'vtype not found'
-                call abort
+                call abort_physic("photochemistry","(count_react): vtype not found",1)
             end if
 
@@ -1250 +1241 @@
               print*, 'It should be 0 for photolysis reations and 1 or 2 for the others'
               print*, 'And not : ', typerate
-              call abort
+              call abort_physic("photochemistry","wrong index coefficient rate line",1)
             end if
 
@@ -1314 +1305 @@
            write(*,*)'28/03/2025: New default value is network.def (in simulation folder)'
            write(*,*)'Former default is chemnetwork/reactfile'
-           call abort
+           call abort_physic("photochemistry","Unable to open chemical reaction file",1)
         end if
 
@@ -1337 +1328 @@
                 nb_phot = nb_phot + 1
             else
-                print*,'Error in photochemistry_asis (get_react):'
-                print*,'vtype not found'
-                call abort
+                call abort_physic("photochemistry","(get_react): vtype not found",1)
             end if
 
@@ -1345 +1334 @@
             if (trim(paramline)=='') then
                 print*, 'Error in reactfile: where are the parameters - line',ilines
-                call abort
+                call abort_physic("photochemistry","Error in reactfile, the parameters - line not found",1)
             else
                 read(paramline,*) reactions(ireact)%rtype
@@ -1393 +1382 @@
                     else
                         print*, 'Error in reactfile: just only one specie can be after M corresponding to the third body - line',ilines
-                        call abort
+                        call abort_physic("Error in reactfile","just only one specie can be after M corresponding to the third body - line",1)
                     end if
                 end if
@@ -1509 +1498 @@
                 print*,'Error in photochemistry_asis (get_react):'
                 print*,'vtype',reactions(ireact)%vtype,' not implemented'
-                call abort
+                call abort_physic("Error in photochemistry_asis (get_react)", "vtype',reactions(ireact)%vtype,' not implemented",1)
             end if