nicer way of character array constructor for GEOSIT_defs (LDAS_Forcin…

…g.F90)
GEOS-ESM · Dec 5, 2023 · 9377a7b · 9377a7b
1 parent a8227f1
commit 9377a7b
Showing 1 changed file with 18 additions and 18 deletions.
diff --git a/src/Components/GEOSldas_GridComp/GEOSmetforce_GridComp/LDAS_Forcing.F90 b/src/Components/GEOSldas_GridComp/GEOSmetforce_GridComp/LDAS_Forcing.F90
@@ -3176,6 +3176,8 @@ subroutine get_GEOS( date_time, force_dtstep, met_path, met_tag,             &
     character(40), dimension(N_MERRA2plusAerosol_vars,  N_defs_cols) :: M2INT_defs
     character(40), dimension(N_MERRA2plusAerosol_vars,  N_defs_cols) :: M2COR_defs
 
+    character(40), dimension(N_G5DAS_vars)                           :: GEOSIT_ftag                          
+
     character(40), dimension(:,:), allocatable :: GEOSgcm_defs
 
     ! NOTE: met_path, prec_path, and met_tag for current ('inst') time and fwd and bkwd 'tavg' 
@@ -3249,26 +3251,24 @@ subroutine get_GEOS( date_time, force_dtstep, met_path, met_tag,             &
     !
     ! define GEOS-IT file specs 
     !
-    ! same as G5DAS except for file tag (column 3)
-    !
-    ! GEOSIT character(40):
-    !
-    !                             1         2         3         4
-    !                    1234567890123456789012345678901234567890     
+    ! same as G5DAS except for file tags (column 3)
+
+    GEOSIT_defs = G5DAS_defs
 
-    GEOSIT_defs( 1,3) = 'lfo_tavg_1hr_glo_L576x361_slv           '
-    GEOSIT_defs( 2,3) = 'lfo_tavg_1hr_glo_L576x361_slv           '
-    GEOSIT_defs( 3,3) = 'lfo_tavg_1hr_glo_L576x361_slv           '
-    GEOSIT_defs( 4,3) = 'lfo_tavg_1hr_glo_L576x361_slv           '
-    GEOSIT_defs( 5,3) = 'lfo_tavg_1hr_glo_L576x361_slv           '
-    GEOSIT_defs( 6,3) = 'lfo_tavg_1hr_glo_L576x361_slv           '
-    GEOSIT_defs( 7,3) = 'lfo_tavg_1hr_glo_L576x361_slv           '
-    GEOSIT_defs( 8,3) = 'lfo_inst_1hr_glo_L576x361_slv           '
-    GEOSIT_defs( 9,3) = 'lfo_inst_1hr_glo_L576x361_slv           '
-    GEOSIT_defs(10,3) = 'lfo_inst_1hr_glo_L576x361_slv           '
-    GEOSIT_defs(11,3) = 'lfo_inst_1hr_glo_L576x361_slv           '
-    GEOSIT_defs(12,3) = 'lfo_inst_1hr_glo_L576x361_slv           '
+    GEOSIT_ftag( 1,3)=[character(len=40):: 'lfo_tavg_1hr_glo_L576x361_slv']
+    GEOSIT_ftag( 2,3)=[character(len=40):: 'lfo_tavg_1hr_glo_L576x361_slv']
+    GEOSIT_ftag( 3,3)=[character(len=40):: 'lfo_tavg_1hr_glo_L576x361_slv']
+    GEOSIT_ftag( 4,3)=[character(len=40):: 'lfo_tavg_1hr_glo_L576x361_slv']
+    GEOSIT_ftag( 5,3)=[character(len=40):: 'lfo_tavg_1hr_glo_L576x361_slv']
+    GEOSIT_ftag( 6,3)=[character(len=40):: 'lfo_tavg_1hr_glo_L576x361_slv']
+    GEOSIT_ftag( 7,3)=[character(len=40):: 'lfo_tavg_1hr_glo_L576x361_slv']
+    GEOSIT_ftag( 8,3)=[character(len=40):: 'lfo_inst_1hr_glo_L576x361_slv']
+    GEOSIT_ftag( 9,3)=[character(len=40):: 'lfo_inst_1hr_glo_L576x361_slv']
+    GEOSIT_ftag(10,3)=[character(len=40):: 'lfo_inst_1hr_glo_L576x361_slv']
+    GEOSIT_ftag(11,3)=[character(len=40):: 'lfo_inst_1hr_glo_L576x361_slv']
+    GEOSIT_ftag(12,3)=[character(len=40):: 'lfo_inst_1hr_glo_L576x361_slv']
 
+    GEOSIT_defs(:,3) = GEOSIT_ftag
 
     ! MERRA-2 file specs with uncorrected (AGCM) precip from the "int" Collection
     !  (ie, the precip generated by the AGCM within the MERRA-2 system)