Skip to content

Commit

Permalink
Add reading GEOSIT forcing data
Browse files Browse the repository at this point in the history
  • Loading branch information
@gmao-qliu
gmao-qliu committed Dec 4, 2023
1 parent 5669de0 commit 1e55616
Showing 1 changed file with 88 additions and 4 deletions.
92 changes: 88 additions & 4 deletions src/Components/GEOSldas_GridComp/GEOSmetforce_GridComp/LDAS_Forcing.F90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -3170,6 +3170,7 @@ subroutine get_GEOS( date_time, force_dtstep, met_path, met_tag, &
real, parameter :: nodata_GEOSgcm = 1.e15

character(40), dimension(N_G5DAS_vars, N_defs_cols) :: G5DAS_defs
character(40), dimension(N_G5DAS_vars, N_defs_cols) :: GEOSIT_defs
character(40), dimension(N_MERRA_vars, N_defs_cols) :: MERRA_defs
character(40), dimension(N_MERRA2plusAerosol_vars, N_defs_cols) :: M2INT_defs
character(40), dimension(N_MERRA2plusAerosol_vars, N_defs_cols) :: M2COR_defs
Expand Down Expand Up @@ -3243,6 +3244,33 @@ subroutine get_GEOS( date_time, force_dtstep, met_path, met_tag, &
G5DAS_defs(11,:)=[character(len=40):: 'QLML ','inst','inst1_2d_lfo_Nx','diag','S']
G5DAS_defs(12,:)=[character(len=40):: 'SPEEDLML','inst','inst1_2d_lfo_Nx','diag','S']

! -----------------------------------------------------------------------
!
! define GEOSIT file specs
!
! columns of GEOSgcm_defs are as follows:
! 1 - short variable name in gfio file
! 2 - averaging mode in G5DAS/MERRA file ('tavg' or 'inst')
! 3 - file tag (eg. 'bkg.sfc' or 'diag_sfc', or 'tavg1_2d_rad_Nx')
! 4 - file dir ('ana' or 'diag')
! 5 - treated as S="state" or F="flux" in subroutine interpolate_to_timestep()

! GEOSIT file specs (default, unless otherwise specified via "met_tag")


GEOSIT_defs( 1,:)=[character(len=40):: 'SWGDN ','tavg','lfo_tavg_1hr_glo_L576x361_slv','diag','F']
GEOSIT_defs( 2,:)=[character(len=40):: 'LWGAB ','tavg','lfo_tavg_1hr_glo_L576x361_slv','diag','F']
GEOSIT_defs( 3,:)=[character(len=40):: 'PARDR ','tavg','lfo_tavg_1hr_glo_L576x361_slv','diag','F']
GEOSIT_defs( 4,:)=[character(len=40):: 'PARDF ','tavg','lfo_tavg_1hr_glo_L576x361_slv','diag','F']
GEOSIT_defs( 5,:)=[character(len=40):: 'PRECCU ','tavg','lfo_tavg_1hr_glo_L576x361_slv','diag','F']
GEOSIT_defs( 6,:)=[character(len=40):: 'PRECLS ','tavg','lfo_tavg_1hr_glo_L576x361_slv','diag','F']
GEOSIT_defs( 7,:)=[character(len=40):: 'PRECSNO ','tavg','lfo_tavg_1hr_glo_L576x361_slv','diag','F']
GEOSIT_defs( 8,:)=[character(len=40):: 'PS ','inst','lfo_inst_1hr_glo_L576x361_slv','diag','S']
GEOSIT_defs( 9,:)=[character(len=40):: 'HLML ','inst','lfo_inst_1hr_glo_L576x361_slv','diag','S']
GEOSIT_defs(10,:)=[character(len=40):: 'TLML ','inst','lfo_inst_1hr_glo_L576x361_slv','diag','S']
GEOSIT_defs(11,:)=[character(len=40):: 'QLML ','inst','lfo_inst_1hr_glo_L576x361_slv','diag','S']
GEOSIT_defs(12,:)=[character(len=40):: 'SPEEDLML','inst','lfo_inst_1hr_glo_L576x361_slv','diag','S']


! MERRA-2 file specs with uncorrected (AGCM) precip from the "int" Collection
! (ie, the precip generated by the AGCM within the MERRA-2 system)
Expand Down Expand Up @@ -3576,8 +3604,12 @@ subroutine get_GEOS( date_time, force_dtstep, met_path, met_tag, &

allocate(GEOSgcm_defs(N_GEOSgcm_vars,N_defs_cols))

GEOSgcm_defs(1:N_G5DAS_vars, :) = G5DAS_defs

if (index(met_tag, 'GEOSIT') > 0 .or. index(met_tag, 'geosit') > 0) then
GEOSgcm_defs(1:N_G5DAS_vars, :) = GEOSIT_defs
else
GEOSgcm_defs(1:N_G5DAS_vars, :) = G5DAS_defs
end if

call parse_G5DAS_met_tag( met_path, met_tag, date_time_inst, &
met_path_inst, prec_path_inst, met_tag_inst, use_prec_corr, &
use_Predictor )
Expand Down Expand Up @@ -4658,6 +4690,10 @@ subroutine parse_G5DAS_met_tag( met_path_in, met_tag_in, date_time, &
type(date_time_type) :: dt_end_d5124_rpit2
type(date_time_type) :: dt_end_d5124_rpit3

type(date_time_type) :: dt_end_d5294_geosit1
type(date_time_type) :: dt_end_d5294_geosit2
type(date_time_type) :: dt_end_d5294_geosit3

type(date_time_type) :: dt_end_e5110_fp
type(date_time_type) :: dt_end_e5130_fp
type(date_time_type) :: dt_end_e5131_fp
Expand Down Expand Up @@ -4738,6 +4774,33 @@ subroutine parse_G5DAS_met_tag( met_path_in, met_tag_in, date_time, &
dt_end_d5124_rpit3%min = 0
dt_end_d5124_rpit3%sec = 0

! | stream start | stream end
! ----------------------------------------
! d5294_geosit1 | 1 Jan 1998 | 1 Jan 2008
! d5294_geosit2 | 1 Jan 2008 | 1 Jan 2018
! d5294_geosit3 | 1 Jan 2018 | (present)

dt_end_d5294_geosit1%year = 2008
dt_end_d5294_geosit1%month = 1
dt_end_d5294_geosit1%day = 1
dt_end_d5294_geosit1%hour = 0
dt_end_d5294_geosit1%min = 0
dt_end_d5294_geosit1%sec = 0

dt_end_d5294_geosit2%year = 2018
dt_end_d5294_geosit2%month = 1
dt_end_d5294_geosit2%day = 1
dt_end_d5294_geosit2%hour = 0
dt_end_d5294_geosit2%min = 0
dt_end_d5294_geosit2%sec = 0

dt_end_d5294_geosit3%year = 9999
dt_end_d5294_geosit3%month = 1
dt_end_d5294_geosit3%day = 1
dt_end_d5294_geosit3%hour = 0
dt_end_d5294_geosit3%min = 0
dt_end_d5294_geosit3%sec = 0

! ---------------------------------
!
! FP streams
Expand Down Expand Up @@ -4957,6 +5020,23 @@ subroutine parse_G5DAS_met_tag( met_path_in, met_tag_in, date_time, &

end if

elseif (met_tag_out(1:18)=='cross_GEOSIT') then

if (datetime_lt_refdatetime( date_time, dt_end_d5294_geosit1 )) then

stream = 'd5294_geosit_jan98' ! use d5294_geosit stream 1

elseif (datetime_lt_refdatetime( date_time, dt_end_d5294_geosit2 )) then

stream = 'd5294_geosit_jan08' ! use d5294_geosit stream 2

else

stream = 'd5294_geosit_jan18' ! use d5294_geosit stream 2

end if


elseif (met_tag_out(1:8)=='cross_FP') then

if (datetime_lt_refdatetime( date_time, dt_end_e5110_fp )) then
Expand Down Expand Up @@ -5103,7 +5183,7 @@ subroutine get_GEOS_forcing_filename(fname_full,file_exists, date_time, daily_fi
! local variables

character(300) :: fname, fname_full_tmp1, fname_full_tmp2
character( 14) :: time_stamp
character( 20) :: time_stamp
character( 4) :: YYYY, HHMM, day_dir
character( 2) :: MM, DD

Expand Down Expand Up @@ -5131,8 +5211,12 @@ subroutine get_GEOS_forcing_filename(fname_full,file_exists, date_time, daily_fi

time_stamp(1:8) = YYYY // MM // DD

else
elseif (index(met_tag,'GEOSIT') > 0 .or. index(met_tag,'geosit') > 0) then

time_stamp = YYYY //'-'// MM //'-'// DD // 'T' // trim(HHMM) // 'Z'

else

time_stamp = YYYY // MM // DD // '_' // trim(HHMM) // 'z'

end if
Expand Down

0 comments on commit 1e55616

Please sign in to comment.