Changeset 4764

Timestamp:

Dec 10, 2023, 10:37:06 PM (5 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.

Location:

LMDZ6/trunk/libf/phylmd/StratAer

Files:

• strataer_emiss_mod.F90 (modified) (2 diffs)
• strataer_local_var_mod.F90 (modified) (1 diff)

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)

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

@@ -192 +192 @@
     nErupt = 0 ! eruption number
     injdur = 0 ! init injection duration
-    nAerErupt = 1 ; nSpeciesErupt = 1 ; id_species = 3
-    year_emit_vol=0 ; mth_emit_vol=0 ; day_emit_vol=0
-    m_species_emiss_vol=0. ; m_aer_emiss_vol=0. ; m_sulf_emiss_vol=0.
-    altemiss_vol=0. ; sigma_alt_vol=0.
-    xlon_min_vol=0. ; xlon_max_vol=0.
-    xlat_min_vol=0. ; xlat_max_vol=0.
-    m_Chlore_emiss_vol=0. ; m_Brome_emiss_vol=0.
-    m_NOx_emiss_vol=0. ; m_H2O_emiss_vol=0.
     day_emit_roc=0 ; ifreqroc=2 ; flh2o=0
+    nAerErupt = 1 ; nSpeciesErupt = 1
     
     !============= Read params =============