Changeset 4763 for LMDZ6/trunk

Timestamp:

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

Author:

lguez

Message:

Polish

File:

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

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

@@ -1 +1 @@
 MODULE strataer_emiss_mod
-! This module contains information about strato microphysic model emission parameters
-  
+  ! This module contains information about strato microphysic model emission parameters
+
 CONTAINS
-    
+
   SUBROUTINE strataer_emiss_init()
-    
+
     USE strataer_local_var_mod
     USE ioipsl_getin_p_mod, ONLY : getin_p
     USE print_control_mod, ONLY : lunout
     USE mod_phys_lmdz_para, ONLY : is_master
-    
+
     ! Local variables
     INTEGER                  :: ispec
-    
+
     WRITE(lunout,*) 'IN STRATAER_EMISS INIT WELCOME!'
-    
+
     IF (flag_emit.EQ.1 .OR. flag_emit.EQ.4) THEN ! Volcano
        CALL getin_p('nErupt',nErupt) !eruption nb
        CALL getin_p('injdur',injdur) !injection duration
-       
+
        IF (flag_emit==1) THEN
           CALL getin_p('nAerErupt', nAerErupt) !sulfur aer nb
@@ -25 +25 @@
           CALL getin_p('nSpeciesErupt', nSpeciesErupt) !chimical species nb
        ENDIF
-       
+
        IF (nErupt.GT.0) THEN
           ALLOCATE(year_emit_vol(nErupt),mth_emit_vol(nErupt),day_emit_vol(nErupt))
@@ -37 +37 @@
                 ALLOCATE(m_H2O_emiss_vol(nErupt))
                 ALLOCATE(m_H2O_emiss_vol_daily(nErupt))
-!                ALLOCATE(d_q_emiss(klon,klev))
+                !                ALLOCATE(d_q_emiss(klon,klev))
                 ALLOCATE(budg_emi(klon,nAerErupt+1))
                 m_H2O_emiss_vol(:)=0.
                 m_H2O_emiss_vol_daily(:)=0.
-!                d_q_emiss(:,:)=0.
+                !                d_q_emiss(:,:)=0.
              ELSE
                 ALLOCATE(budg_emi(klon,nAerErupt))
@@ -48 +48 @@
           ELSEIF (flag_emit==4) THEN
              ALLOCATE(m_Chlore_emiss_vol(nErupt))
-             ALLOCATE(id_HCl) 
+             ALLOCATE(id_HCl)
              ALLOCATE(m_Brome_emiss_vol(nErupt))
              ALLOCATE(id_HBr)
@@ -65 +65 @@
                'No eruption define in physiq_def (nErupt=0). Add at one eruption or use background condition.',1)
        ENDIF ! fin if nerupt
-       
+
        ! injection params (dates, loc, injections params)
        CALL getin_p('year_emit_vol',year_emit_vol)
@@ -78 +78 @@
        IF (flag_emit==1) THEN
           CALL getin_p('m_sulf_emiss_vol',m_sulf_emiss_vol)
-        if (ok_qemiss) then
-         CALL getin_p('m_H2O_emiss_vol',m_H2O_emiss_vol)
-        endif
+          if (ok_qemiss) then
+             CALL getin_p('m_H2O_emiss_vol',m_H2O_emiss_vol)
+          endif
        ELSEIF (flag_emit==4) THEN
           CALL getin_p('id_species',id_species)
-          CALL getin_p('m_Chlore_emiss_vol',m_Chlore_emiss_vol)  
+          CALL getin_p('m_Chlore_emiss_vol',m_Chlore_emiss_vol)
           CALL getin_p('id_HCl',id_HCl)
           CALL getin_p('m_Brome_emiss_vol',m_Brome_emiss_vol)
@@ -92 +92 @@
           CALL getin_p('id_H2O',id_H2O)
        ENDIF
-       
+
     ELSEIF (flag_emit == 2) THEN ! SAI
        CALL getin_p('m_aer_emiss_sai',m_aer_emiss_sai)
@@ -105 +105 @@
        CALL getin_p('day_emit_sai_start',day_emit_sai_start)
        CALL getin_p('day_emit_sai_end',day_emit_sai_end)
-       
+
     ELSEIF (flag_emit == 3) THEN ! SAI between latitudes
        CALL getin_p('m_aer_emiss_sai',m_aer_emiss_sai)
@@ -114 +114 @@
        CALL getin_p('xlat_min_sai',xlat_min_sai)
     ENDIF
-    
+
     IF (flag_emit == 1) THEN
        DO ispec=1,nAerErupt
@@ -124 +124 @@
              m_species_emiss_vol(:,ispec) = m_Chlore_emiss_vol(:)
           ENDIF
-          IF (id_species(ispec) == id_HBr) THEN 
+          IF (id_species(ispec) == id_HBr) THEN
              m_species_emiss_vol(:,ispec) = m_Brome_emiss_vol(:)
           ENDIF
@@ -135 +135 @@
        ENDDO
     ENDIF
-    
+
     !============= Injection ponderation =============
     IF (flag_emit > 0) THEN
@@ -141 +141 @@
        WRITE(lunout,*) 'IN STRATAER INIT : ponde_lonlat_vol',ponde_lonlat_vol
     ENDIF
-    
+
     !============= Print params =============
     IF (is_master) THEN
@@ -148 +148 @@
           WRITE(lunout,*) 'IN STRATAER injdur: ',injdur
           WRITE(lunout,*) 'IN STRATAER nAerErupt: ',nAerErupt
-          
+
           WRITE(lunout,*) 'IN STRATAER : year_emit_vol',year_emit_vol
           WRITE(lunout,*) 'IN STRATAER : mth_emit_vol',mth_emit_vol
@@ -195 +195 @@
        ENDIF
     ENDIF ! if master
-    
+
     WRITE(lunout,*) 'IN STRATAER_EMISS END'
   END SUBROUTINE strataer_emiss_init
-  
+
   ! Compute the ponderation to applicate in each grid point for all eruptions and init
   ! dlat & dlon variables
   SUBROUTINE strataer_ponde_init()
-    
+
     USE regular_lonlat_mod, ONLY: lon_reg, lat_reg
     USE dimphy, ONLY: klon
@@ -208 +208 @@
     USE print_control_mod, ONLY : lunout
     USE strataer_local_var_mod
-    
+
     INCLUDE "YOMCST.h"  !--RPI
-    
+
     ! local var
     REAL                :: lat_reg_deg,lon_reg_deg ! lat and lon of grid points in degree
     INTEGER             :: ieru, i, j
-    
+
     ALLOCATE(ponde_lonlat_vol(nErupt))
-    
+
     !Compute lon/lat ponderation for injection
     dlat=180./2./FLOAT(nbp_lat)   ! d latitude in degree
@@ -224 +224 @@
     WRITE(lunout,*) 'IN STRATAER_INIT xlat_min=',xlat_min_vol,'xlat_max=',xlat_max_vol
     WRITE(lunout,*) 'IN STRATAER_INIT xlon_min=',xlon_min_vol,'xlon_max=',xlon_max_vol
-    
+
     DO ieru=1, nErupt
        ponde_lonlat_vol(ieru) = 0
@@ -241 +241 @@
        ENDIF
     ENDDO !ieru
-    
+
     WRITE(lunout,*) 'IN STRATAER_PONDE_INIT ponde_lonlat: ', ponde_lonlat_vol
-    
+
   END SUBROUTINE strataer_ponde_init
-  
+
 END MODULE strataer_emiss_mod