Changeset 4765

Timestamp:

Dec 10, 2023, 10:37:07 PM (5 months ago)

Author:

lguez

Message:

Bug fix: only call getin_p if nErupt > 0

The variables for which we call getin_p here should be allocated and
they are allocated only if nErupt > 0.

File:

• LMDZ6/trunk/libf/phylmd/StratAer/strataer_emiss_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/StratAer/strataer_emiss_mod.F90

@@ -35 +35 @@
           xlon_min_vol=0. ; xlon_max_vol=0.
           xlat_min_vol=0. ; xlat_max_vol=0.
+          ! injection params (dates, loc, injections params)
+          CALL getin_p('year_emit_vol',year_emit_vol)
+          CALL getin_p('mth_emit_vol',mth_emit_vol)
+          CALL getin_p('day_emit_vol',day_emit_vol)
+          CALL getin_p('altemiss_vol',altemiss_vol)
+          CALL getin_p('sigma_alt_vol',sigma_alt_vol)
+          CALL getin_p('xlon_min_vol',xlon_min_vol)
+          CALL getin_p('xlon_max_vol',xlon_max_vol)
+          CALL getin_p('xlat_min_vol',xlat_min_vol)
+          CALL getin_p('xlat_max_vol',xlat_max_vol)
           IF (flag_emit==1) THEN
              ALLOCATE(m_sulf_emiss_vol(nErupt))
@@ -75 +85 @@
        ENDIF ! fin if nerupt
 
-       ! injection params (dates, loc, injections params)
-       CALL getin_p('year_emit_vol',year_emit_vol)
-       CALL getin_p('mth_emit_vol',mth_emit_vol)
-       CALL getin_p('day_emit_vol',day_emit_vol)
-       CALL getin_p('altemiss_vol',altemiss_vol)
-       CALL getin_p('sigma_alt_vol',sigma_alt_vol)
-       CALL getin_p('xlon_min_vol',xlon_min_vol)
-       CALL getin_p('xlon_max_vol',xlon_max_vol)
-       CALL getin_p('xlat_min_vol',xlat_min_vol)
-       CALL getin_p('xlat_max_vol',xlat_max_vol)
        IF (flag_emit==1) THEN
           CALL getin_p('m_sulf_emiss_vol',m_sulf_emiss_vol)