Changeset 4764 for LMDZ6/trunk/libf/phylmd/StratAer/strataer_emiss_mod.F90

Timestamp:

Dec 10, 2023, 10:37:06 PM (6 months ago)

Author:

lguez

Message:

Bug fix: move assignments to strataer_emiss_init

Bug fix: move assignments from procedure strataer_init to procedure
strataer_emiss_init. The variables to which we make assignments are
allocated in strataer_emiss_init and strataer_emiss_init is called
by physiq after strataer_init.

File:

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

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

@@ -28 +28 @@
        IF (nErupt.GT.0) THEN
           ALLOCATE(year_emit_vol(nErupt),mth_emit_vol(nErupt),day_emit_vol(nErupt))
+          year_emit_vol=0 ; mth_emit_vol=0 ; day_emit_vol=0
           ALLOCATE(altemiss_vol(nErupt),sigma_alt_vol(nErupt))
           ALLOCATE(xlat_min_vol(nErupt),xlon_min_vol(nErupt))
           ALLOCATE(xlat_max_vol(nErupt),xlon_max_vol(nErupt))
+          altemiss_vol=0. ; sigma_alt_vol=0.
+          xlon_min_vol=0. ; xlon_max_vol=0.
+          xlat_min_vol=0. ; xlat_max_vol=0.
           IF (flag_emit==1) THEN
              ALLOCATE(m_sulf_emiss_vol(nErupt))
              ALLOCATE(m_aer_emiss_vol(nErupt,nAerErupt))
+             m_aer_emiss_vol=0. ; m_sulf_emiss_vol=0.
              IF (ok_qemiss) then
                 ALLOCATE(m_H2O_emiss_vol(nErupt))
@@ -52 +57 @@
              ALLOCATE(id_HBr)
              ALLOCATE(id_species(nSpeciesErupt))
+             id_species = 3
              ALLOCATE(m_species_emiss_vol(nErupt,nSpeciesErupt))
+             m_species_emiss_vol=0.
              ALLOCATE(m_NOx_emiss_vol(nErupt))
              ALLOCATE(m_H2O_emiss_vol(nErupt))
+             m_Chlore_emiss_vol=0. ; m_Brome_emiss_vol=0.
+             m_NOx_emiss_vol=0. ; m_H2O_emiss_vol=0.
              ALLOCATE(id_NOx)
              ALLOCATE(id_H2O)