Ignore:
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().

Location:
trunk/LMDZ.GENERIC/libf/aeronostd/chemistry_early_earth
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.GENERIC/libf/aeronostd/chemistry_early_earth/calchim_asis_early_earth.F90

    r1810 r3893  
    206206         i_co2 = igcm_co2
    207207         if (i_co2 == 0) then
    208             write(*,*) "calchim: Error; no CO2 tracer !!!"
    209             stop
     208            call abort_physic("calchim_early_earth", "No CO2 tracer",1)
    210209         else
    211210            nbq = nbq + 1
    … …  
    214213         i_co = igcm_co
    215214         if (i_co == 0) then
    216             write(*,*) "calchim: Error; no CO tracer !!!"
    217             stop
     215            call abort_physic("calchim_early_earth", "No CO tracer",1)
    218216         else
    219217            nbq = nbq + 1
    … …  
    222220         i_o = igcm_o
    223221         if (i_o == 0) then
    224             write(*,*) "calchim: Error; no O tracer !!!"
    225             stop
     222            call abort_physic("calchim_early_earth", "No O tracer",1)
    226223         else
    227224            nbq = nbq + 1
    … …  
    230227         i_o1d = igcm_o1d
    231228         if (i_o1d == 0) then
    232             write(*,*) "calchim: Error; no O1D tracer !!!"
    233             stop
     229            call abort_physic("calchim_early_earth", "No O1D tracer",1)
    234230         else
    235231            nbq = nbq + 1
    … …  
    238234         i_o2 = igcm_o2
    239235         if (i_o2 == 0) then
    240             write(*,*) "calchim: Error; no O2 tracer !!!"
    241             stop
     236            call abort_physic("calchim_early_earth", "No O2 tracer",1)
    242237         else
    243238            nbq = nbq + 1
    … …  
    246241         i_o3 = igcm_o3
    247242         if (i_o3 == 0) then
    248             write(*,*) "calchim: Error; no O3 tracer !!!"
    249             stop
     243            call abort_physic("calchim_early_earth", "No O3 tracer",1)
    250244         else
    251245            nbq = nbq + 1
    … …  
    254248         i_h = igcm_h
    255249         if (i_h == 0) then
    256             write(*,*) "calchim: Error; no H tracer !!!"
    257             stop
     250            call abort_physic("calchim_early_earth", "No H tracer",1)
    258251         else
    259252            nbq = nbq + 1
    … …  
    262255         i_h2 = igcm_h2
    263256         if (i_h2 == 0) then
    264             write(*,*) "calchim: Error; no H2 tracer !!!"
    265             stop
     257            call abort_physic("calchim_early_earth", "No H2 tracer",1)
    266258         else
    267259            nbq = nbq + 1
    … …  
    270262         i_oh = igcm_oh
    271263         if (i_oh == 0) then
    272             write(*,*) "calchim: Error; no OH tracer !!!"
    273             stop
     264            call abort_physic("calchim_early_earth", "No OH tracer",1)
    274265         else
    275266            nbq = nbq + 1
    … …  
    278269         i_ho2 = igcm_ho2
    279270         if (i_ho2 == 0) then
    280             write(*,*) "calchim: Error; no HO2 tracer !!!"
    281             stop
     271            call abort_physic("calchim_early_earth", "No HO2 tracer",1)
    282272         else
    283273            nbq = nbq + 1
    … …  
    286276         i_h2o2 = igcm_h2o2
    287277         if (i_h2o2 == 0) then
    288             write(*,*) "calchim: Error; no H2O2 tracer !!!"
    289             stop
     278            call abort_physic("calchim_early_earth", "No H2O2 tracer",1)
    290279         else
    291280            nbq = nbq + 1
    … …  
    294283         i_h2o = igcm_h2o_vap
    295284         if (i_h2o == 0) then
    296             write(*,*) "calchim: Error; no water vapor tracer !!!"
    297             stop
     285            call abort_physic("calchim_early_earth", "No water vapor tracer",1)
    298286         else
    299287            nbq = nbq + 1
    … …  
    303291         if (i_n == 0) then
    304292            if (photochem) then
    305                write(*,*) "calchim: Error; no N tracer !!!"
    306293               write(*,*) "N will be ignored in the chemistry"
    307                stop
     294               call abort_physic("calchim_early_earth", "No N tracer",1)
    308295            end if
    309296         else
    … …  
    314301         if (i_n2d == 0) then
    315302            if (photochem) then
    316                write(*,*) "calchim: Error; no N2D tracer !!!"
    317303               write(*,*) "N2d will be ignored in the chemistry"
    318                stop
     304               call abort_physic("calchim_early_earth", "No N2D tracer",1)
    319305            end if
    320306         else
    … …  
    325311         if (i_no == 0) then
    326312            if (photochem) then
    327                write(*,*) "calchim: Error; no NO tracer !!!"
    328313               write(*,*) "NO will be ignored in the chemistry"
    329                stop
     314               call abort_physic("calchim_early_earth", "No NO tracer",1)
    330315            end if
    331316         else
    … …  
    336321         if (i_no2 == 0) then
    337322            if (photochem) then
    338                write(*,*) "calchim: Error; no NO2 tracer !!!"
    339323               write(*,*) "NO2 will be ignored in the chemistry"
    340                stop
     324               call abort_physic("calchim_early_earth", "No NO2 tracer",1)
    341325            end if
    342326         else
    … …  
    346330         i_n2 = igcm_n2
    347331         if (i_n2 == 0) then
    348             write(*,*) "calchim: Error; no N2 tracer !!!"
    349             stop
     332            call abort_physic("calchim_early_earth", "No N2 tracer",1)
    350333         else
    351334            nbq = nbq + 1
    … …  
    354337         i_ch4 = igcm_ch4
    355338         if (i_ch4 == 0) then
    356             write(*,*) "calchim: Error; no CH4 tracer !!!"
    357             stop
     339            call abort_physic("calchim_early_earth", "No CH4 tracer",1)
    358340         else
    359341            nbq = nbq + 1
    … …  
    362344         i_ch3 = igcm_ch3
    363345         if (i_ch3 == 0) then
    364             write(*,*) "calchim: Error; no ch3 tracer !!!"
    365             stop
     346            call abort_physic("calchim_early_earth", "No CH3 tracer",1)
    366347         else
    367348            nbq = nbq + 1
    … …  
    370351         i_ch = igcm_ch
    371352         if (i_ch == 0) then
    372             write(*,*) "calchim: Error; no ch3 tracer !!!"
    373             stop
     353            call abort_physic("calchim_early_earth", "No CH tracer",1)
    374354         else
    375355            nbq = nbq + 1
    … …  
    378358         i_3ch2 = igcm_3ch2
    379359         if (i_3ch2 == 0) then
    380             write(*,*) "calchim: Error; no 3ch2 tracer !!!"
    381             stop
     360            call abort_physic("calchim_early_earth", "No 3CH2 tracer",1)
    382361         else
    383362            nbq = nbq + 1
    … …  
    386365         i_1ch2 = igcm_1ch2
    387366         if (i_1ch2 == 0) then
    388             write(*,*) "calchim: Error; no 1ch2 tracer !!!"
    389             stop
     367            call abort_physic("calchim_early_earth", "No 1CH2 tracer",1)
    390368         else
    391369            nbq = nbq + 1
    … …  
    394372         i_cho = igcm_cho
    395373         if (i_cho == 0) then
    396             write(*,*) "calchim: Error; no cho tracer !!!"
    397             stop
     374            call abort_physic("calchim_early_earth", "No CHO tracer",1)
    398375         else
    399376            nbq = nbq + 1
    … …  
    402379         i_ch2o = igcm_ch2o
    403380         if (i_ch2o == 0) then
    404             write(*,*) "calchim: Error; no ch2o tracer !!!"
    405             stop
     381            call abort_physic("calchim_early_earth", "No CH2O tracer",1)
    406382         else
    407383            nbq = nbq + 1
    … …  
    410386         i_ch3o = igcm_ch3o
    411387         if (i_ch3o == 0) then
    412             write(*,*) "calchim: Error; no ch3o tracer !!!"
    413             stop
     388            call abort_physic("calchim_early_earth", "No CH3O tracer",1)
    414389         else
    415390            nbq = nbq + 1
    … …  
    418393         i_c = igcm_c
    419394         if (i_c == 0) then
    420             write(*,*) "calchim: Error; no c tracer !!!"
    421             stop
     395            call abort_physic("calchim_early_earth", "No C tracer",1)
    422396         else
    423397            nbq = nbq + 1
    … …  
    426400         i_c2 = igcm_c2
    427401         if (i_c2 == 0) then
    428             write(*,*) "calchim: Error; no c2 tracer !!!"
    429             stop
     402            call abort_physic("calchim_early_earth", "No C2 tracer",1)
    430403         else
    431404            nbq = nbq + 1
    … …  
    434407         i_c2h = igcm_c2h
    435408         if (i_c2h == 0) then
    436             write(*,*) "calchim: Error; no c2h tracer !!!"
    437             stop
     409            call abort_physic("calchim_early_earth", "No C2H tracer",1)
    438410         else
    439411            nbq = nbq + 1
    … …  
    442414         i_c2h2 = igcm_c2h2
    443415         if (i_c2h2 == 0) then
    444             write(*,*) "calchim: Error; no c2h2 tracer !!!"
    445             stop
     416            call abort_physic("calchim_early_earth", "No C2H2 tracer",1)
    446417         else
    447418            nbq = nbq + 1
    … …  
    450421         i_c2h3 = igcm_c2h3
    451422         if (i_c2h3 == 0) then
    452             write(*,*) "calchim: Error; no c2h3 tracer !!!"
    453             stop
     423            call abort_physic("calchim_early_earth", "No C2H3 tracer",1)
    454424         else
    455425            nbq = nbq + 1
    … …  
    458428         i_c2h4 = igcm_c2h4
    459429         if (i_c2h4 == 0) then
    460             write(*,*) "calchim: Error; no c2h4 tracer !!!"
    461             stop
     430            call abort_physic("calchim_early_earth", "No C2H4 tracer",1)
    462431         else
    463432            nbq = nbq + 1
    … …  
    466435         i_c2h6 = igcm_c2h6
    467436         if (i_c2h6 == 0) then
    468             write(*,*) "calchim: Error; no c2h6 tracer !!!"
    469             stop
     437            call abort_physic("calchim_early_earth", "No C2H6 tracer",1)
    470438         else
    471439            nbq = nbq + 1
    … …  
    474442         i_ch2co = igcm_ch2co
    475443         if (i_ch2co == 0) then
    476             write(*,*) "calchim: Error; no ch2co tracer !!!"
    477             stop
     444            call abort_physic("calchim_early_earth", "No CH2CO tracer",1)
    478445         else
    479446            nbq = nbq + 1
    … …  
    482449         i_ch3co = igcm_ch3co
    483450         if (i_ch3co == 0) then
    484             write(*,*) "calchim: Error; no ch3co tracer !!!"
    485             stop
     451            call abort_physic("calchim_early_earth", "No CH3CO tracer",1)
    486452         else
    487453            nbq = nbq + 1
    … …  
    490456         i_hcaer = igcm_hcaer
    491457         if (i_hcaer == 0) then
    492             write(*,*) "calchim: Error; no hcaer tracer !!!"
    493             stop
     458            call abort_physic("calchim_early_earth", "No hcaer tracer",1)
    494459         else
    495460            nbq = nbq + 1

  • trunk/LMDZ.GENERIC/libf/aeronostd/chemistry_early_earth/photochemistry_asis_early_earth.F90

    r1810 r3893  
    275275#else
    276276#   write(*,*) "photochemistry_asis error, missing LAPACK routine dgesv"
    277 #   stop
     277#   call abort_physic("photochemistry_early_earth", "missing LAPACK routine dgesv",1)
    278278#endif
    279279
    … …  
    399399      call dgesv(nesp,1,mat,nesp,indx,cnew,nesp,code)
    400400#else
    401    write(*,*) "photochemistry_asis error, missing LAPACK routine dgesv"
    402    stop
     401   call abort_physic("photochemistry_early_earth", "missing LAPACK routine dgesv",1)
    403402#endif
    404403
    … …  
    33203319   print*, 'nb_reaction_4_max = ', nb_reaction_4_max
    33213320   print*, 'nb_reaction_3_max = ', nb_reaction_3_max
    3322    stop
     3321   call abort_physic("photochemistry_early_earth", "wrong dimensions in indice",1)
    33233322end if 
    33243323
    … …  
    33953394
    33963395         if (igcm_co2 == 0) then
    3397             write(*,*) "gcmtochim: Error; no CO2 tracer !!!"
    3398             stop
     3396            call abort_physic("phochemistry_early_earth", "no CO2 tracer", 1)
    33993397         endif
    34003398         if (igcm_co == 0) then
    3401             write(*,*) "gcmtochim: Error; no CO tracer !!!"
    3402             stop
     3399            call abort_physic("phochemistry_early_earth", "no CO tracer", 1)
    34033400         end if
    34043401         if (igcm_o == 0) then
    3405             write(*,*) "gcmtochim: Error; no O tracer !!!"
    3406             stop
     3402            call abort_physic("phochemistry_early_earth", "no O tracer", 1)
    34073403         end if
    34083404         if (igcm_o1d == 0) then
    3409             write(*,*) "gcmtochim: Error; no O1D tracer !!!"
    3410             stop
     3405            call abort_physic("phochemistry_early_earth", "no O1D tracer", 1)
    34113406         end if
    34123407         if (igcm_o2 == 0) then
    3413             write(*,*) "gcmtochim: Error; no O2 tracer !!!"
    3414             stop
     3408            call abort_physic("phochemistry_early_earth", "no O2 tracer", 1)
    34153409         end if
    34163410         if (igcm_o3 == 0) then
    3417             write(*,*) "gcmtochim: Error; no O3 tracer !!!"
    3418             stop
     3411            call abort_physic("phochemistry_early_earth", "no O3 tracer", 1)
    34193412         end if
    34203413         if (igcm_h == 0) then
    3421             write(*,*) "gcmtochim: Error; no H tracer !!!"
    3422             stop
     3414            call abort_physic("phochemistry_early_earth", "no H tracer", 1)
    34233415         end if
    34243416         if (igcm_h2 == 0) then
    3425             write(*,*) "gcmtochim: Error; no H2 tracer !!!"
    3426             stop
     3417            call abort_physic("phochemistry_early_earth", "no H2 tracer", 1)
    34273418         end if
    34283419         if (igcm_oh == 0) then
    3429             write(*,*) "gcmtochim: Error; no OH tracer !!!"
    3430             stop
     3420            call abort_physic("phochemistry_early_earth", "no OH tracer", 1)
    34313421         end if
    34323422         if (igcm_ho2 == 0) then
    3433             write(*,*) "gcmtochim: Error; no HO2 tracer !!!"
    3434             stop
     3423            call abort_physic("phochemistry_early_earth", "no HO2 tracer", 1)
    34353424         end if
    34363425         if (igcm_h2o2 == 0) then
    3437             write(*,*) "gcmtochim: Error; no H2O2 tracer !!!"
    3438             stop
     3426            call abort_physic("phochemistry_early_earth", "no H2O2 tracer", 1)
    34393427         end if
    34403428         if (igcm_n == 0) then
    3441             write(*,*) "gcmtochim: Error; no N tracer !!!"
    3442             stop
     3429            call abort_physic("phochemistry_early_earth", "no N tracer", 1)
    34433430         end if
    34443431         if (igcm_n2d == 0) then
    3445             write(*,*) "gcmtochim: Error; no N2D tracer !!!"
    3446             stop
     3432            call abort_physic("phochemistry_early_earth", "no N2D tracer", 1)
    34473433         end if
    34483434         if (igcm_no == 0) then
    3449             write(*,*) "gcmtochim: Error; no NO tracer !!!"
    3450             stop
     3435            call abort_physic("phochemistry_early_earth", "no NO tracer", 1)
    34513436         end if
    34523437         if (igcm_no2 == 0) then
    3453             write(*,*) "gcmtochim: Error; no NO2 tracer !!!"
    3454             stop
     3438            call abort_physic("phochemistry_early_earth", "no NO2 tracer", 1)
    34553439         end if
    34563440         if (igcm_n2 == 0) then
    3457             write(*,*) "gcmtochim: Error; no N2 tracer !!!"
    3458             stop
     3441            call abort_physic("phochemistry_early_earth", "no N2 tracer", 1)
    34593442         end if
    34603443         if (igcm_h2o_vap == 0) then
    3461             write(*,*) "gcmtochim: Error; no water vapor tracer !!!"
    3462             stop
     3444            call abort_physic("phochemistry_early_earth", "no water vapor tracer", 1)
    34633445         end if
    34643446
    34653447         if (igcm_ch4 == 0) then
    3466             write(*,*) "gcmtochim: Error; no CH4 tracer !!!"
    3467             stop
     3448            call abort_physic("phochemistry_early_earth", "no CH4 tracer", 1)
    34683449         end if
    34693450         if (igcm_ch3 == 0) then
    3470             write(*,*) "gcmtochim: Error; no CH3 tracer !!!"
    3471             stop
     3451            call abort_physic("phochemistry_early_earth", "no CH3 tracer", 1)
    34723452         end if
    34733453         if (igcm_ch == 0) then
    3474             write(*,*) "gcmtochim: Error; no CH tracer !!!"
    3475             stop
     3454            call abort_physic("phochemistry_early_earth", "no CH tracer", 1)
    34763455         end if
    34773456         if (igcm_3ch2 == 0) then
    3478             write(*,*) "gcmtochim: Error; no 3CH2 tracer !!!"
    3479             stop
     3457            call abort_physic("phochemistry_early_earth", "no 3CH2 tracer", 1)
    34803458         end if
    34813459         if (igcm_1ch2 == 0) then
    3482             write(*,*) "gcmtochim: Error; no 1CH2 tracer !!!"
    3483             stop
     3460            call abort_physic("phochemistry_early_earth", "no 1CH2 tracer", 1)
    34843461         end if
    34853462         if (igcm_cho == 0) then
    3486             write(*,*) "gcmtochim: Error; no CHO tracer !!!"
    3487             stop
     3463            call abort_physic("phochemistry_early_earth", "no CHO tracer", 1)
    34883464         end if
    34893465         if (igcm_ch2o == 0) then
    3490             write(*,*) "gcmtochim: Error; no CH2O tracer !!!"
    3491             stop
     3466            call abort_physic("phochemistry_early_earth", "no CH2O tracer", 1)
    34923467         end if
    34933468         if (igcm_ch3o == 0) then
    3494             write(*,*) "gcmtochim: Error; no CH3O tracer !!!"
    3495             stop
     3469            call abort_physic("phochemistry_early_earth", "no CH3O tracer", 1)
    34963470         end if
    34973471         if (igcm_c == 0) then
    3498             write(*,*) "gcmtochim: Error; no C tracer !!!"
    3499             stop
     3472            call abort_physic("phochemistry_early_earth", "no C tracer", 1)
    35003473         end if
    35013474         if (igcm_c2 == 0) then
    3502             write(*,*) "gcmtochim: Error; no C2 tracer !!!"
    3503             stop
     3475            call abort_physic("phochemistry_early_earth", "no C2 tracer", 1)
    35043476         end if
    35053477         if (igcm_c2h == 0) then
    3506             write(*,*) "gcmtochim: Error; no C2H tracer !!!"
    3507             stop
     3478            call abort_physic("phochemistry_early_earth", "no C2H tracer", 1)
    35083479         end if
    35093480         if (igcm_c2h2 == 0) then
    3510             write(*,*) "gcmtochim: Error; no C2H2 tracer !!!"
    3511             stop
     3481            call abort_physic("phochemistry_early_earth", "no C2H2 tracer", 1)
    35123482         end if
    35133483         if (igcm_c2h3 == 0) then
    3514             write(*,*) "gcmtochim: Error; no C2H3 tracer !!!"
    3515             stop
     3484            call abort_physic("phochemistry_early_earth", "no C2H3 tracer", 1)
    35163485         end if
    35173486         if (igcm_c2h4 == 0) then
    3518             write(*,*) "gcmtochim: Error; no C2H4 tracer !!!"
    3519             stop
     3487            call abort_physic("phochemistry_early_earth", "no C2H4 tracer", 1)
    35203488         end if
    35213489         if (igcm_c2h6 == 0) then
    3522             write(*,*) "gcmtochim: Error; no C2H6 tracer !!!"
    3523             stop
     3490            call abort_physic("phochemistry_early_earth", "no C2H6 tracer", 1)
    35243491         end if
    35253492         if (igcm_ch2co == 0) then
    3526             write(*,*) "gcmtochim: Error; no CH2CO tracer !!!"
    3527             stop
     3493            call abort_physic("phochemistry_early_earth", "no CH2CO tracer", 1)
    35283494         end if
    35293495         if (igcm_ch3co == 0) then
    3530             write(*,*) "gcmtochim: Error; no CH3CO tracer !!!"
    3531             stop
     3496            call abort_physic("phochemistry_early_earth", "no CH3CO tracer", 1)
    35323497         end if
    35333498         if (igcm_hcaer == 0) then
    3534             write(*,*) "gcmtochim: Error; no HCAER tracer !!!"
    3535             stop
     3499            call abort_physic("phochemistry_early_earth", "no hcaer tracer", 1)
    35363500         end if
    35373501

  • trunk/LMDZ.GENERIC/libf/aeronostd/chemistry_early_earth/read_phototable_early_earth.F90

    r1796 r3893  
    6464        write(*,*)'   callphys.def with:'
    6565        write(*,*)'   phototable=filename'
    66         stop
     66        call abort_physic("read_phototable_early_earth", "cannot open photolysis lookup table",1)
    6767      end if
    6868
Note: See TracChangeset for help on using the changeset viewer.