Changeset 4625 for LMDZ6/trunk/libf/phylmd/StratAer

Timestamp:

Jul 12, 2023, 11:31:50 PM (17 months ago)

Author:

dcugnet

Message:

StratAer? commit, N. Lebas

Location:

LMDZ6/trunk/libf/phylmd/StratAer

Files:

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

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

@@ -43 +43 @@
   REAL paprs_glo(klon_glo,klev+1)
 
-  REAL, POINTER:: latitude(:)
+  REAL, ALLOCATABLE :: latitude(:)
 ! (of input data sorted in strictly ascending order)
 
-  REAL, POINTER:: longitude(:)
+  REAL, ALLOCATABLE :: longitude(:)
 ! (of input data sorted in strictly ascending order)
 
-  REAL, POINTER:: time(:)
+  REAL, ALLOCATABLE :: time(:)
 ! (of input data sorted in strictly ascending order)
 
-  REAL, POINTER:: lev(:)
+  REAL, ALLOCATABLE :: lev(:)
 ! levels of input data
 

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

@@ -99 +99 @@
        CALL getin_p('xlat_sai',xlat_sai)
        CALL getin_p('xlon_sai',xlon_sai)
+       CALL getin_p('year_emit_sai_start',year_emit_sai_start)
+       CALL getin_p('year_emit_sai_end',year_emit_sai_end)
+       CALL getin_p('mth_emit_sai_start',mth_emit_sai_start)
+       CALL getin_p('mth_emit_sai_end',mth_emit_sai_end)
+       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
@@ -169 +175 @@
           WRITE(lunout,*) 'IN STRATAER : altemiss_sai',altemiss_sai
           WRITE(lunout,*) 'IN STRATAER : sigma_alt_sai',sigma_alt_sai
+          WRITE(lunout,*) 'IN STRATAER : year_emit_sai start/end',year_emit_sai_start,year_emit_sai_end
+          WRITE(lunout,*) 'IN STRATAER : mth_emit_sai start/end',mth_emit_sai_start,mth_emit_sai_end
+          WRITE(lunout,*) 'IN STRATAER : day_emit_sai start/end',day_emit_sai_start,day_emit_sai_end
           WRITE(lunout,*) 'IN STRATAER : xlat_min_sai',xlat_min_sai
           WRITE(lunout,*) 'IN STRATAER : xlat_max_sai',xlat_max_sai

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

@@ -99 +99 @@
   REAL,SAVE    :: sigma_alt_sai          ! standard deviation of emission altitude in m
   !$OMP THREADPRIVATE(m_aer_emiss_sai, altemiss_sai, sigma_alt_sai)
+  INTEGER,SAVE    :: year_emit_sai_start
+  INTEGER,SAVE    :: year_emit_sai_end
+  INTEGER,SAVE    :: mth_emit_sai_start
+  INTEGER,SAVE    :: mth_emit_sai_end
+  INTEGER,SAVE    :: day_emit_sai_start
+  INTEGER,SAVE    :: day_emit_sai_end
+  !$OMP THREADPRIVATE(year_emit_sai_start, year_emit_sai_end)
+  !$OMP THREADPRIVATE(mth_emit_sai_start, mth_emit_sai_end)
+  !$OMP THREADPRIVATE(day_emit_sai_start, day_emit_sai_end)
   REAL,SAVE    :: xlat_sai               ! latitude of SAI in degree
   REAL,SAVE    :: xlon_sai               ! longitude of SAI in degree