Skip to content

Commit 3bb5e91

Browse files
twjulianotwjuliano
authored
Adding accumulated GHI and fixing registry units for solar diagnostics package (#1610)
TYPE: enhancement, bug fix KEYWORDS: WRF-Solar, diagnostics SOURCE: Timothy Juliano and Pedro Jimenez (NCAR/RAL) DESCRIPTION OF CHANGES: We have added a new 2D variable (`ghi_accum`) and time series variable (`ts_ghi_accum`) for the WRF-Solar diagnostics package. This variable uses `SWDNB` as computed by the radiation parameterization to accumulate downwelling GHI at the surface. The new variable is set to zero at the beginning of the simulation for either `restart=.true.` or `restart=.false.`. As a separate but very minor change, we have clarified units in the registry for cloud base and cloud height WRF-Solar variables. LIST OF MODIFIED FILES: M Registry/Registry.EM_COMMON M Registry/registry.solar_fields M dyn_em/start_em.F M phys/module_diag_solar.F M phys/module_diagnostics_driver.F M phys/module_physics_init.F M share/wrf_timeseries.F M run/README.tslist TESTS CONDUCTED: 1. Conducted test on Cheyenne to ensure that GHI is summed over time and set to zero properly at initialization for both 2D and time series variables. 2. All jenkins test are passing. RELEASE NOTE: A new variable that accumulates GHI over time is added to the WRF-Solar diagnostics package. Cloud base and cloud height units (AGL) are clarified in the registry.
1 parent b186098 commit 3bb5e91

8 files changed

+60
-31
lines changed

‎Registry/Registry.EM_COMMON

+1-1
Original file line numberDiff line numberDiff line change
@@ -3226,7 +3226,7 @@ package icedepth_one seaice_thickness_opt==1 - state:icedept
32263226
#Time series options for text output
32273227
package notseries process_time_series==0 - -
32283228
package tseries process_time_series==1 - state:ts_hour,ts_u,ts_v,ts_q,ts_t,ts_psfc,ts_glw,ts_gsw,ts_hfx,ts_lh,ts_tsk,ts_tslb,ts_clw,ts_rainc,ts_rainnc,ts_u_profile,ts_v_profile,ts_gph_profile,ts_th_profile,ts_p_profile,ts_w_profile
3229-
package tseries_add_solar process_time_series==2 - state:ts_hour,ts_u,ts_v,ts_q,ts_t,ts_psfc,ts_glw,ts_gsw,ts_hfx,ts_lh,ts_tsk,ts_tslb,ts_clw,ts_rainc,ts_rainnc,ts_u_profile,ts_v_profile,ts_gph_profile,ts_th_profile,ts_cldfrac2d,ts_wvp,ts_lwp,ts_iwp,ts_swp,ts_lwp_tot,ts_iwp_tot,ts_swp_tot,ts_re_qc,ts_re_qi,ts_re_qs,ts_re_qc_tot,ts_re_qi_tot,ts_re_qs_tot,ts_tau_qc,ts_tau_qi,ts_tau_qs,ts_tau_qc_tot,ts_tau_qi_tot,ts_tau_qs_tot,ts_cbaseht,ts_ctopht,ts_cbaseht_tot,ts_ctopht_tot,ts_clrnidx,ts_p_profile,ts_w_profile,ts_swdown,ts_swddni,ts_swddif,ts_swdownc,ts_swddnic,ts_swdown2,ts_swddni2,ts_swddif2,ts_swdownc2,ts_swddnic2
3229+
package tseries_add_solar process_time_series==2 - state:ts_hour,ts_u,ts_v,ts_q,ts_t,ts_psfc,ts_glw,ts_gsw,ts_hfx,ts_lh,ts_tsk,ts_tslb,ts_clw,ts_rainc,ts_rainnc,ts_u_profile,ts_v_profile,ts_gph_profile,ts_th_profile,ts_cldfrac2d,ts_wvp,ts_lwp,ts_iwp,ts_swp,ts_lwp_tot,ts_iwp_tot,ts_swp_tot,ts_re_qc,ts_re_qi,ts_re_qs,ts_re_qc_tot,ts_re_qi_tot,ts_re_qs_tot,ts_tau_qc,ts_tau_qi,ts_tau_qs,ts_tau_qc_tot,ts_tau_qi_tot,ts_tau_qs_tot,ts_cbaseht,ts_ctopht,ts_cbaseht_tot,ts_ctopht_tot,ts_clrnidx,ts_p_profile,ts_w_profile,ts_swdown,ts_swddni,ts_swddif,ts_swdownc,ts_swddnic,ts_swdown2,ts_swddni2,ts_swddif2,ts_swdownc2,ts_swddnic2,ts_ghi_accum
32303230

32313231
# WRF-HAILCAST
32323232
state real HAILCAST_DHAIL1 ij misc 1 - r "HAILCAST_DHAIL1" "WRF-HAILCAST Hail Diameter, 1st rank order" "mm"

‎Registry/registry.solar_fields

+26-24
Original file line numberDiff line numberDiff line change
@@ -5,8 +5,8 @@
55
rconfig integer solar_diagnostics namelist,diags 1 0 - "solar diagnostics: 0=no, 1=yes"
66

77
# 3-D arrays used for diagnosed quantities
8-
state real qc_tot ikj dyn_em 1 - - "QC_TOT" "CLOUD WATER MIXING RATIO RES + UNRES" "kg kg-1"
9-
state real qi_tot ikj dyn_em 1 - - "QI_TOT" "ICE MIXING RATIO RES + UNRES" "kg kg-1"
8+
state real qc_tot ikj dyn_em 1 - - "QC_TOT" "CLOUD WATER MIXING RATIO RESOLVED + UNRESOLVED" "kg kg-1"
9+
state real qi_tot ikj dyn_em 1 - - "QI_TOT" "ICE MIXING RATIO RESOLVED + UNRESOLVED" "kg kg-1"
1010

1111
# 2-D arrays that will be filled with diagnosed quantities
1212
state real cldfrac2d ij dyn_em 1 - h "CLDFRAC2D" "2-D MAX CLOUD FRACTION" "%"
@@ -15,25 +15,26 @@ state real lwp ij dyn_em 1 - h "LWP"
1515
state real iwp ij dyn_em 1 - h "IWP" "ICE CLOUD WATER PATH" "kg m-2"
1616
state real swp ij dyn_em 1 - h "SWP" "SNOW CLOUD WATER PATH" "kg m-2"
1717
state real wp_sum ij dyn_em 1 - h "WP_SUM" "SUM OF LWP+IWP+SWP" "kg m-2"
18-
state real lwp_tot ij dyn_em 1 - h "LWP_TOT" "LIQUID CLOUD WATER PATH RES + UNRES" "kg m-2"
19-
state real iwp_tot ij dyn_em 1 - h "IWP_TOT" "ICE CLOUD WATER PATH RES + UNRES" "kg m-2"
20-
state real wp_tot_sum ij dyn_em 1 - h "WP_TOT_SUM" "SUM OF LWP+IWP+SWP RES + UNRES" "kg m-2"
18+
state real lwp_tot ij dyn_em 1 - h "LWP_TOT" "LIQUID CLOUD WATER PATH RESOLVED + UNRESOLVED" "kg m-2"
19+
state real iwp_tot ij dyn_em 1 - h "IWP_TOT" "ICE CLOUD WATER PATH RESOLVED + UNRESOLVED" "kg m-2"
20+
state real wp_tot_sum ij dyn_em 1 - h "WP_TOT_SUM" "SUM OF LWP+IWP+SWP RESOLVED + UNRESOLVED" "kg m-2"
2121
state real re_qc ij dyn_em 1 - h "RE_QC" "MASS-WEIGHTED LIQUID CLOUD EFFECTIVE RADIUS" "m"
2222
state real re_qi ij dyn_em 1 - h "RE_QI" "MASS-WEIGHTED ICE EFFECTIVE RADIUS" "m"
2323
state real re_qs ij dyn_em 1 - h "RE_QS" "MASS-WEIGHTED SNOW EFFECTIVE RADIUS" "m"
24-
state real re_qc_tot ij dyn_em 1 - h "RE_QC_TOT" "MASS-WEIGHTED LIQUID CLOUD EFFECTIVE RADIUS RES + UNRES" "m"
25-
state real re_qi_tot ij dyn_em 1 - h "RE_QI_TOT" "MASS-WEIGHTED ICE EFFECTIVE RADIUS RES + UNRES" "m"
24+
state real re_qc_tot ij dyn_em 1 - h "RE_QC_TOT" "MASS-WEIGHTED LIQUID CLOUD EFFECTIVE RADIUS RESOLVED + UNRESOLVED" "m"
25+
state real re_qi_tot ij dyn_em 1 - h "RE_QI_TOT" "MASS-WEIGHTED ICE EFFECTIVE RADIUS RESOLVED + UNRESOLVED" "m"
2626
state real tau_qc ij dyn_em 1 - h "TAU_QC" "MASS-WEIGHTED LIQUID CLOUD OPTICAL THICKNESS" " "
2727
state real tau_qi ij dyn_em 1 - h "TAU_QI" "MASS-WEIGHTED ICE OPTICAL THICKNESS" " "
2828
state real tau_qs ij dyn_em 1 - h "TAU_QS" "MASS-WEIGHTED SNOW OPTICAL THICKNESS" " "
29-
state real tau_qc_tot ij dyn_em 1 - h "TAU_QC_TOT" "MASS-WEIGHTED LIQUID CLOUD OPTICAL THICKNESS RES + UNRES" " "
30-
state real tau_qi_tot ij dyn_em 1 - h "TAU_QI_TOT" "MASS-WEIGHTED ICE OPTICAL THICKNESS RES + UNRES" " "
31-
state real cbaseht ij dyn_em 1 - h "CBASEHT" "CLOUD BASE HEIGHT" "m"
32-
state real ctopht ij dyn_em 1 - h "CTOPHT" "CLOUD TOP HEIGHT" "m"
33-
state real cbaseht_tot ij dyn_em 1 - h "CBASEHT_TOT" "CLOUD BASE HEIGHT RES + UNRES" "m"
34-
state real ctopht_tot ij dyn_em 1 - h "CTOPHT_TOT" "CLOUD TOP HEIGHT RES + UNRES" "m"
29+
state real tau_qc_tot ij dyn_em 1 - h "TAU_QC_TOT" "MASS-WEIGHTED LIQUID CLOUD OPTICAL THICKNESS RESOLVED + UNRESOLVED" " "
30+
state real tau_qi_tot ij dyn_em 1 - h "TAU_QI_TOT" "MASS-WEIGHTED ICE OPTICAL THICKNESS RESOLVED + UNRESOLVED" " "
31+
state real cbaseht ij dyn_em 1 - h "CBASEHT" "CLOUD BASE HEIGHT" "m agl"
32+
state real ctopht ij dyn_em 1 - h "CTOPHT" "CLOUD TOP HEIGHT" "m agl"
33+
state real cbaseht_tot ij dyn_em 1 - h "CBASEHT_TOT" "CLOUD BASE HEIGHT RESOLVED + UNRESOLVED" "m agl"
34+
state real ctopht_tot ij dyn_em 1 - h "CTOPHT_TOT" "CLOUD TOP HEIGHT RESOLVED + UNRESOLVED" "m agl"
3535
state real clrnidx ij dyn_em 1 - h "CLRNIDX" "CLEARNESS INDEX" " "
36-
state real sza ij dyn_em 1 - h "SZA" "SOLAR ZENITH ANGLE" "deg"
36+
state real sza ij dyn_em 1 - h "SZA" "SOLAR ZENITH ANGLE" "deg"
37+
state real ghi_accum ij dyn_em 1 - h "GHI_ACCUM" "ACCUMULATED GHI" "J m-2"
3738

3839
# Time series variables
3940
state real ts_cldfrac2d ?! misc - - - "TS_CLDFRAC2D" "2-D MAX CLOUD FRACTION"
@@ -42,25 +43,26 @@ state real ts_lwp ?! misc - - - "TS_L
4243
state real ts_iwp ?! misc - - - "TS_IWP" "ICE CLOUD WATER PATH"
4344
state real ts_swp ?! misc - - - "TS_SWP" "SNOW CLOUD WATER PATH"
4445
state real ts_wp_sum ?! misc - - - "TS_WP_SUM" "SUM OF LWP+IWP+SWP"
45-
state real ts_lwp_tot ?! misc - - - "TS_LWP_TOT" "LIQUID CLOUD WATER PATH RES + UNRES"
46-
state real ts_iwp_tot ?! misc - - - "TS_IWP_TOT" "ICE CLOUD WATER PATH RES + UNRES"
47-
state real ts_wp_tot_sum ?! misc - - - "TS_WP_TOT_SUM" "SUM OF LWP+IWP+SWP RES + UNRES"
46+
state real ts_lwp_tot ?! misc - - - "TS_LWP_TOT" "LIQUID CLOUD WATER PATH RESOLVED + UNRESOLVED"
47+
state real ts_iwp_tot ?! misc - - - "TS_IWP_TOT" "ICE CLOUD WATER PATH RESOLVED + UNRESOLVED"
48+
state real ts_wp_tot_sum ?! misc - - - "TS_WP_TOT_SUM" "SUM OF LWP+IWP+SWP RESOLVED + UNRESOLVED"
4849
state real ts_re_qc ?! misc - - - "TS_RE_QC" "MASS-WEIGHTED LIQUID CLOUD EFFECTIVE RADIUS"
4950
state real ts_re_qi ?! misc - - - "TS_RE_QI" "MASS-WEIGHTED ICE EFFECTIVE RADIUS"
5051
state real ts_re_qs ?! misc - - - "TS_RE_QS" "MASS-WEIGHTED SNOW EFFECTIVE RADIUS"
51-
state real ts_re_qc_tot ?! misc - - - "TS_RE_QC_TOT" "MASS-WEIGHTED LIQUID CLOUD EFFECTIVE RADIUS RES + UNRES"
52-
state real ts_re_qi_tot ?! misc - - - "TS_RE_QI_TOT" "MASS-WEIGHTED ICE EFFECTIVE RADIUS RES + UNRES"
52+
state real ts_re_qc_tot ?! misc - - - "TS_RE_QC_TOT" "MASS-WEIGHTED LIQUID CLOUD EFFECTIVE RADIUS RESOLVED + UNRESOLVED"
53+
state real ts_re_qi_tot ?! misc - - - "TS_RE_QI_TOT" "MASS-WEIGHTED ICE EFFECTIVE RADIUS RESOLVED + UNRESOLVED"
5354
state real ts_tau_qc ?! misc - - - "TS_TAU_QC" "MASS-WEIGHTED LIQUID CLOUD OPTICAL THICKNESS"
5455
state real ts_tau_qi ?! misc - - - "TS_TAU_QI" "MASS-WEIGHTED ICE OPTICAL THICKNESS"
5556
state real ts_tau_qs ?! misc - - - "TS_TAU_QS" "MASS-WEIGHTED SNOW OPTICAL THICKNESS"
56-
state real ts_tau_qc_tot ?! misc - - - "TS_TAU_QC_TOT" "MASS-WEIGHTED LIQUID CLOUD OPTICAL THICKNESS RES + UNRES"
57-
state real ts_tau_qi_tot ?! misc - - - "TS_TAU_QI_TOT" "MASS-WEIGHTED ICE OPTICAL THICKNESS RES + UNRES"
57+
state real ts_tau_qc_tot ?! misc - - - "TS_TAU_QC_TOT" "MASS-WEIGHTED LIQUID CLOUD OPTICAL THICKNESS RESOLVED + UNRESOLVED"
58+
state real ts_tau_qi_tot ?! misc - - - "TS_TAU_QI_TOT" "MASS-WEIGHTED ICE OPTICAL THICKNESS RESOLVED + UNRESOLVED"
5859
state real ts_cbaseht ?! misc - - - "TS_CBASEHT" "CLOUD BASE HEIGHT"
5960
state real ts_ctopht ?! misc - - - "TS_CTOPHT" "CLOUD TOP HEIGHT"
60-
state real ts_cbaseht_tot ?! misc - - - "TS_CBASEHT_TOT" "CLOUD BASE HEIGHT RES + UNRES"
61-
state real ts_ctopht_tot ?! misc - - - "TS_CTOPHT_TOT" "CLOUD TOP HEIGHT RES + UNRES"
61+
state real ts_cbaseht_tot ?! misc - - - "TS_CBASEHT_TOT" "CLOUD BASE HEIGHT RESOLVED + UNRESOLVED"
62+
state real ts_ctopht_tot ?! misc - - - "TS_CTOPHT_TOT" "CLOUD TOP HEIGHT RESOLVED + UNRESOLVED"
6263
state real ts_clrnidx ?! misc - - - "TS_CLRNIDX" "CLEARNESS INDEX"
6364
state real ts_sza ?! misc - - - "TS_SZA" "SOLAR ZENITH ANGLE"
65+
state real ts_ghi_accum ?! misc - - - "TS_GHI_ACCUM" "ACCUMULATED GHI"
6466
state real ts_swdown ?! misc - - - "TS_SWDOWN" "DOWNWARD SHORT WAVE FLUX AT GROUND SURFACE"
6567
state real ts_swddni ?! misc - - - "TS_SWDDNI" "SHORTWAVE SURFACE DOWNWARD DIRECT NORMAL IRRADIANCE"
6668
state real ts_swddif ?! misc - - - "TS_SWDDIF" "SHORTWAVE SURFACE DOWNWARD DIFFUSE IRRADIANCE"
@@ -75,4 +77,4 @@ state real ts_swddnic2 ?! misc - - - "TS_S
7577
# Package declarations
7678

7779
package skip_solar_output solar_diagnostics==0 - -
78-
package do_solar_output solar_diagnostics==1 - state:cldfrac2d,wvp,lwp,iwp,swp,wp_sum,lwp_tot,iwp_tot,wp_tot_sum,re_qc,re_qi,re_qs,re_qc_tot,re_qi_tot,tau_qc,tau_qi,tau_qs,tau_qc_tot,tau_qi_tot,cbaseht,ctopht,cbaseht_tot,ctopht_tot,clrnidx,sza
80+
package do_solar_output solar_diagnostics==1 - state:cldfrac2d,wvp,lwp,iwp,swp,wp_sum,lwp_tot,iwp_tot,wp_tot_sum,re_qc,re_qi,re_qs,re_qc_tot,re_qi_tot,tau_qc,tau_qi,tau_qs,tau_qc_tot,tau_qi_tot,cbaseht,ctopht,cbaseht_tot,ctopht_tot,clrnidx,sza,ghi_accum

‎dyn_em/start_em.F

+1
Original file line numberDiff line numberDiff line change
@@ -1255,6 +1255,7 @@ SUBROUTINE start_domain_em ( grid, allowed_to_read &
12551255
,grid%ccn_conc & ! RAS
12561256
,grid%QKE &!JOE-for mynn
12571257
,grid%pep_pbl & ! EEPS
1258+
,grid%ghi_accum & ! Solar diagnostics
12581259
,grid%landusef,grid%landusef2,grid%mosaic_cat_index & ! danli mosaic
12591260
,grid%TSK_mosaic,grid%TSLB_mosaic,grid%SMOIS_mosaic,grid%SH2O_mosaic & ! danli mosaic
12601261
,grid%CANWAT_mosaic,grid%SNOW_mosaic,grid%SNOWH_mosaic,grid%SNOWC_mosaic & ! danli mosaic

‎phys/module_diag_solar.F

+7-2
Original file line numberDiff line numberDiff line change
@@ -28,10 +28,10 @@ MODULE module_diag_solar
2828
! !
2929
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
3030

31-
SUBROUTINE solar_diag (rho, dz8w, ph, phb, cldfrac3d, coszen, swdnb, swdnt, &
31+
SUBROUTINE solar_diag (dt, rho, dz8w, ph, phb, cldfrac3d, coszen, swdnb, swdnt, &
3232
param_first_scalar, p_qc, p_qi, p_qs, qv, qc, qi, qs, &
3333
qc_tot, qi_tot, has_reqc, has_reqi, has_reqs, f_qv, f_qc, f_qi, f_qs, &
34-
re_cloud, re_ice, re_snow, clrnidx, sza, cldfrac2d, wvp2d, lwp2d, iwp2d, swp2d, &
34+
re_cloud, re_ice, re_snow, clrnidx, sza, ghi_accum, cldfrac2d, wvp2d, lwp2d, iwp2d, swp2d, &
3535
wp2d_sum, lwp2d_tot, iwp2d_tot, wp2d_tot_sum, re_cloud_path, re_ice_path, re_snow_path, &
3636
re_cloud_path_tot, re_ice_path_tot, tau_qc, tau_qi, tau_qs, tau_qc_tot, tau_qi_tot, &
3737
cbase, ctop, cbase_tot, ctop_tot, ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, &
@@ -43,6 +43,7 @@ SUBROUTINE solar_diag (rho, dz8w, ph, phb, cldfrac3d, coszen, swdnb, swdnt, &
4343

4444
IMPLICIT NONE
4545

46+
REAL, INTENT(IN) :: dt
4647
REAL, DIMENSION(ims:ime, jms:jme), INTENT(IN) :: coszen, swdnb, swdnt
4748
REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(IN) :: ph, phb, cldfrac3d, qv, qc, qi, qs, qc_tot, qi_tot, &
4849
re_cloud, re_ice, re_snow, rho, dz8w
@@ -51,6 +52,7 @@ SUBROUTINE solar_diag (rho, dz8w, ph, phb, cldfrac3d, coszen, swdnb, swdnt, &
5152
REAL, DIMENSION(ims:ime, jms:jme), INTENT(OUT) :: clrnidx, sza, cldfrac2d, wvp2d, lwp2d, iwp2d, swp2d, wp2d_sum, &
5253
lwp2d_tot, iwp2d_tot, wp2d_tot_sum, re_cloud_path, re_ice_path, re_snow_path, re_cloud_path_tot, re_ice_path_tot, &
5354
tau_qc, tau_qi, tau_qs, tau_qc_tot, tau_qi_tot, cbase, ctop, cbase_tot, ctop_tot
55+
REAL, DIMENSION(ims:ime, jms:jme), INTENT(INOUT) :: ghi_accum
5456
INTEGER, INTENT(IN) :: ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms, kme, &
5557
ips, ipe, jps, jpe, kps, kpe, kts, kte, num_tiles
5658
INTEGER, DIMENSION(num_tiles), INTENT(IN) :: i_start, i_end, j_start, j_end
@@ -88,6 +90,9 @@ SUBROUTINE solar_diag (rho, dz8w, ph, phb, cldfrac3d, coszen, swdnb, swdnt, &
8890
clrnidx(i, j) = MISSING
8991
end if
9092

93+
!!! ACCUMULATED GHI
94+
ghi_accum(i, j) = ghi_accum(i, j) + swdnb(i, j) * dt
95+
9196
!!! 2-D CLOUD FRACTION
9297
cldfrac2d(i, j) = MAXVAL( cldfrac3d(i, kts:kte-1, j) )
9398

‎phys/module_diagnostics_driver.F

+4-3
Original file line numberDiff line numberDiff line change
@@ -62,7 +62,7 @@ SUBROUTINE diagnostics_driver ( grid, config_flags, &
6262
model_config_rec
6363

6464
USE module_streams
65-
USE module_utility, ONLY : WRFU_Time
65+
USE module_utility, ONLY : WRFU_Time
6666

6767
!=============================================================
6868
! USE Association for the Diagnostic Packages
@@ -1094,7 +1094,8 @@ SUBROUTINE diagnostics_driver ( grid, config_flags, &
10941094

10951095
SOLAR_FIELDS: IF (config_flags%solar_diagnostics == do_solar_output) THEN
10961096
CALL wrf_debug (100 , '--> CALL DIAGNOSTICS PACKAGE: SOLAR_DIAG')
1097-
CALL solar_diag ( &
1097+
1098+
CALL solar_diag (dt=grid%dt, &
10981099
rho=grid%rho, dz8w=dz8w, ph=grid%ph_2, phb=grid%phb, &
10991100
cldfrac3d=grid%cldfra, coszen=grid%coszen, swdnb=grid%swdnb, &
11001101
swdnt=grid%swdnt, param_first_scalar=param_first_scalar, &
@@ -1116,7 +1117,7 @@ SUBROUTINE diagnostics_driver ( grid, config_flags, &
11161117
tau_qc_tot=grid%tau_qc_tot, tau_qi_tot=grid%tau_qi_tot, &
11171118
cbase=grid%cbaseht, ctop=grid%ctopht, &
11181119
cbase_tot=grid%cbaseht_tot, ctop_tot=grid%ctopht_tot, &
1119-
clrnidx=grid%clrnidx, sza=grid%sza, &
1120+
clrnidx=grid%clrnidx, sza=grid%sza, ghi_accum=grid%ghi_accum, &
11201121
ids=ids, ide=ide, jds=jds, jde=jde, kds=kds, kde=kde, &
11211122
ims=ims, ime=ime, jms=jms, jme=jme, kms=kms, kme=kme, &
11221123
ips=ips, ipe=ipe, jps=jps, jpe=jpe, kps=kps, kpe=kpe, &

0 commit comments

Comments
 (0)