Skip to content

Commit 807bd87

Browse files
MennoVeermanMennoVeerman
committed
change all references to ice effective radius to ice effective diameter. Not working yet, still needs to be linked to most recent rrtmgp-data
1 parent ce3e86e commit 807bd87

20 files changed

+181
-181
lines changed

‎allsky/allsky_init.py‎

Lines changed: 8 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -31,7 +31,7 @@ def calc_p_q_T(z):
3131

3232
i_above_zt = np.where(z > z_trop)
3333
q[i_above_zt] = q_t
34-
34+
3535
T_0 = 300.
3636
Tv_0 = (1. + 0.608*q_0)*T_0
3737

@@ -51,16 +51,16 @@ def calc_p_q_T(z):
5151
Tv = T * (1. + 0.608*q)
5252

5353
# T = Tv / (1. + 0.608*q)
54-
54+
5555
g = 9.79764
5656
Rd = 287.04
5757
p0 = 101480.
5858

5959
p = p0 * (Tv / Tv_0)**(g/(Rd*gamma))
60-
60+
6161
p_tmp = p0 * (Tv/Tv_0)**(g/(Rd*gamma)) \
6262
* np.exp( -( (g*(z-z_trop)) / (Rd*Tv) ) )
63-
63+
6464
p[i_above_zt] = p_tmp[i_above_zt]
6565

6666
return p, q, T
@@ -159,20 +159,20 @@ def calc_p_q_T(z):
159159
nc_lwp = nc_file.createVariable('lwp', float_type, ('lay', 'y', 'x'))
160160
nc_iwp = nc_file.createVariable('iwp', float_type, ('lay', 'y', 'x'))
161161
nc_rel = nc_file.createVariable('rel', float_type, ('lay', 'y', 'x'))
162-
nc_rei = nc_file.createVariable('rei', float_type, ('lay', 'y', 'x'))
162+
nc_dei = nc_file.createVariable('dei', float_type, ('lay', 'y', 'x'))
163163

164164
min_rel, max_rel = 2.5, 21.5
165-
min_rei, max_rei = 10., 180.
165+
min_dei, max_dei = 10., 180.
166166

167167
rel_val = 0.5*(min_rel + max_rel)
168-
rei_val = 0.5*(min_rei + max_rei)
168+
dei_val = 0.5*(min_dei + max_dei)
169169

170170
cloud_flag = (np.arange(1, n_col_x+1)%3 > 0)
171171
cloud_mask = np.where((nc_p_lay[:,:,:] > 1.e4) & (nc_p_lay[:,:,:] < 9.e4) & cloud_flag[None, None,:], True, False)
172172

173173
nc_lwp[:,:,:] = np.where(cloud_mask & (nc_T_lay[:,:,:] > 263.), 10., 0.)
174174
nc_iwp[:,:,:] = np.where(cloud_mask & (nc_T_lay[:,:,:] < 273.), 10., 0.)
175175
nc_rel[:,:,:] = np.where(nc_lwp[:,:,:] > 0., rel_val, 0.)
176-
nc_rei[:,:,:] = np.where(nc_iwp[:,:,:] > 0., rei_val, 0.)
176+
nc_dei[:,:,:] = np.where(nc_iwp[:,:,:] > 0., dei_val, 0.)
177177

178178
nc_file.close()

‎include/Cloud_optics.h‎

Lines changed: 10 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -41,18 +41,18 @@ class Cloud_optics : public Optical_props
4141
Cloud_optics(
4242
const Array<Float,2>& band_lims_wvn,
4343
const Float radliq_lwr, const Float radliq_upr, const Float radliq_fac,
44-
const Float radice_lwr, const Float radice_upr, const Float radice_fac,
44+
const Float diamice_lwr, const Float diamice_upr, const Float diamice_fac,
4545
const Array<Float,2>& lut_extliq, const Array<Float,2>& lut_ssaliq, const Array<Float,2>& lut_asyliq,
4646
const Array<Float,3>& lut_extice, const Array<Float,3>& lut_ssaice, const Array<Float,3>& lut_asyice);
4747

4848
void cloud_optics(
4949
const Array<Float,2>& clwp, const Array<Float,2>& ciwp,
50-
const Array<Float,2>& reliq, const Array<Float,2>& reice,
50+
const Array<Float,2>& reliq, const Array<Float,2>& deice,
5151
Optical_props_1scl& optical_props);
5252

5353
void cloud_optics(
5454
const Array<Float,2>& clwp, const Array<Float,2>& ciwp,
55-
const Array<Float,2>& reliq, const Array<Float,2>& reice,
55+
const Array<Float,2>& reliq, const Array<Float,2>& deice,
5656
Optical_props_2str& optical_props);
5757

5858
private:
@@ -64,8 +64,8 @@ class Cloud_optics : public Optical_props
6464
// Lookup table constants.
6565
Float radliq_lwr;
6666
Float radliq_upr;
67-
Float radice_lwr;
68-
Float radice_upr;
67+
Float diamice_lwr;
68+
Float diamice_upr;
6969

7070
// Lookup table coefficients.
7171
Array<Float,2> lut_extliq;
@@ -84,18 +84,18 @@ class Cloud_optics_gpu : public Optical_props_gpu
8484
Cloud_optics_gpu(
8585
const Array<Float,2>& band_lims_wvn,
8686
const Float radliq_lwr, const Float radliq_upr, const Float radliq_fac,
87-
const Float radice_lwr, const Float radice_upr, const Float radice_fac,
87+
const Float diamice_lwr, const Float diamice_upr, const Float diamice_fac,
8888
const Array<Float,2>& lut_extliq, const Array<Float,2>& lut_ssaliq, const Array<Float,2>& lut_asyliq,
8989
const Array<Float,3>& lut_extice, const Array<Float,3>& lut_ssaice, const Array<Float,3>& lut_asyice);
9090

9191
void cloud_optics(
9292
const Array_gpu<Float,2>& clwp, const Array_gpu<Float,2>& ciwp,
93-
const Array_gpu<Float,2>& reliq, const Array_gpu<Float,2>& reice,
93+
const Array_gpu<Float,2>& reliq, const Array_gpu<Float,2>& deice,
9494
Optical_props_1scl_gpu& optical_props);
9595

9696
void cloud_optics(
9797
const Array_gpu<Float,2>& clwp, const Array_gpu<Float,2>& ciwp,
98-
const Array_gpu<Float,2>& reliq, const Array_gpu<Float,2>& reice,
98+
const Array_gpu<Float,2>& reliq, const Array_gpu<Float,2>& deice,
9999
Optical_props_2str_gpu& optical_props);
100100

101101
private:
@@ -107,8 +107,8 @@ class Cloud_optics_gpu : public Optical_props_gpu
107107
// Lookup table constants.
108108
Float radliq_lwr;
109109
Float radliq_upr;
110-
Float radice_lwr;
111-
Float radice_upr;
110+
Float diamice_lwr;
111+
Float diamice_upr;
112112

113113
// Lookup table coefficients.
114114
Array<Float,2> lut_extliq;

‎include_rt/Cloud_optics_rt.h‎

Lines changed: 6 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -32,27 +32,27 @@
3232
// Forward declarations.
3333
class Optical_props_rt;
3434

35-
#ifdef USECUDA
35+
#ifdef USECUDA
3636
class Cloud_optics_rt : public Optical_props_rt
3737
{
3838
public:
3939
Cloud_optics_rt(
4040
const Array<Float,2>& band_lims_wvn,
4141
const Float radliq_lwr, const Float radliq_upr, const Float radliq_fac,
42-
const Float radice_lwr, const Float radice_upr, const Float radice_fac,
42+
const Float diamice_lwr, const Float diamice_upr, const Float diamice_fac,
4343
const Array<Float,2>& lut_extliq, const Array<Float,2>& lut_ssaliq, const Array<Float,2>& lut_asyliq,
4444
const Array<Float,3>& lut_extice, const Array<Float,3>& lut_ssaice, const Array<Float,3>& lut_asyice);
4545

4646
void cloud_optics(
4747
const int ibnd,
4848
const Array_gpu<Float,2>& clwp, const Array_gpu<Float,2>& ciwp,
49-
const Array_gpu<Float,2>& reliq, const Array_gpu<Float,2>& reice,
49+
const Array_gpu<Float,2>& reliq, const Array_gpu<Float,2>& deice,
5050
Optical_props_1scl_rt& optical_props);
5151

5252
void cloud_optics(
5353
const int ibnd,
5454
const Array_gpu<Float,2>& clwp, const Array_gpu<Float,2>& ciwp,
55-
const Array_gpu<Float,2>& reliq, const Array_gpu<Float,2>& reice,
55+
const Array_gpu<Float,2>& reliq, const Array_gpu<Float,2>& deice,
5656
Optical_props_2str_rt& optical_props);
5757

5858
private:
@@ -64,8 +64,8 @@ class Cloud_optics_rt : public Optical_props_rt
6464
// Lookup table constants.
6565
Float radliq_lwr;
6666
Float radliq_upr;
67-
Float radice_lwr;
68-
Float radice_upr;
67+
Float diamice_lwr;
68+
Float diamice_upr;
6969

7070
// Lookup table coefficients.
7171
Array<Float,2> lut_extliq;

‎include_test/Radiation_solver.h‎

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -54,7 +54,7 @@ class Radiation_solver_longwave
5454
const Array<Float,2>& col_dry,
5555
const Array<Float,1>& t_sfc, const Array<Float,2>& emis_sfc,
5656
const Array<Float,2>& lwp, const Array<Float,2>& iwp,
57-
const Array<Float,2>& rel, const Array<Float,2>& rei,
57+
const Array<Float,2>& rel, const Array<Float,2>& dei,
5858
Array<Float,3>& tau, Array<Float,3>& lay_source,
5959
Array<Float,3>& lev_source, Array<Float,2>& sfc_source,
6060
Array<Float,2>& lw_flux_up, Array<Float,2>& lw_flux_dn, Array<Float,2>& lw_flux_net,
@@ -81,7 +81,7 @@ class Radiation_solver_longwave
8181
const Array_gpu<Float,2>& col_dry,
8282
const Array_gpu<Float,1>& t_sfc, const Array_gpu<Float,2>& emis_sfc,
8383
const Array_gpu<Float,2>& lwp, const Array_gpu<Float,2>& iwp,
84-
const Array_gpu<Float,2>& rel, const Array_gpu<Float,2>& rei,
84+
const Array_gpu<Float,2>& rel, const Array_gpu<Float,2>& dei,
8585
Array_gpu<Float,3>& tau, Array_gpu<Float,3>& lay_source,
8686
Array_gpu<Float,3>& lev_source, Array_gpu<Float,2>& sfc_source,
8787
Array_gpu<Float,2>& lw_flux_up, Array_gpu<Float,2>& lw_flux_dn, Array_gpu<Float,2>& lw_flux_net,
@@ -151,7 +151,7 @@ class Radiation_solver_shortwave
151151
const Array<Float,2>& sfc_alb_dir, const Array<Float,2>& sfc_alb_dif,
152152
const Array<Float,1>& tsi_scaling, const Array<Float,1>& mu0,
153153
const Array<Float,2>& lwp, const Array<Float,2>& iwp,
154-
const Array<Float,2>& rel, const Array<Float,2>& rei,
154+
const Array<Float,2>& rel, const Array<Float,2>& dei,
155155
const Array<Float,2>& rh,
156156
const Aerosol_concs& aerosol_concs,
157157
Array<Float,3>& tau, Array<Float,3>& ssa, Array<Float,3>& g,
@@ -188,7 +188,7 @@ class Radiation_solver_shortwave
188188
const Array_gpu<Float,2>& sfc_alb_dir, const Array_gpu<Float,2>& sfc_alb_dif,
189189
const Array_gpu<Float,1>& tsi_scaling, const Array_gpu<Float,1>& mu0,
190190
const Array_gpu<Float,2>& lwp, const Array_gpu<Float,2>& iwp,
191-
const Array_gpu<Float,2>& rel, const Array_gpu<Float,2>& rei,
191+
const Array_gpu<Float,2>& rel, const Array_gpu<Float,2>& dei,
192192
const Array_gpu<Float,2>& rh,
193193
const Aerosol_concs_gpu& aerosol_concs,
194194
Array_gpu<Float,3>& tau, Array_gpu<Float,3>& ssa, Array_gpu<Float,3>& g,

‎include_test/Radiation_solver_bw.h‎

Lines changed: 7 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -58,7 +58,7 @@ class Radiation_solver_longwave
5858
const Array<Float,2>& col_dry,
5959
const Array<Float,1>& t_sfc, const Array<Float,2>& emis_sfc,
6060
const Array<Float,2>& lwp, const Array<Float,2>& iwp,
61-
const Array<Float,2>& rel, const Array<Float,2>& rei,
61+
const Array<Float,2>& rel, const Array<Float,2>& dei,
6262
Array<Float,3>& tau, Array<Float,3>& lay_source,
6363
Array<Float,3>& lev_source_inc, Array<Float,3>& lev_source_dec, Array<Float,2>& sfc_source,
6464
Array<Float,2>& lw_flux_up, Array<Float,2>& lw_flux_dn, Array<Float,2>& lw_flux_net,
@@ -76,7 +76,7 @@ class Radiation_solver_longwave
7676
Array_gpu<Float,2>& col_dry,
7777
const Array_gpu<Float,1>& t_sfc, const Array_gpu<Float,2>& emis_sfc,
7878
const Array_gpu<Float,2>& lwp, const Array_gpu<Float,2>& iwp,
79-
const Array_gpu<Float,2>& rel, const Array_gpu<Float,2>& rei,
79+
const Array_gpu<Float,2>& rel, const Array_gpu<Float,2>& dei,
8080
Array_gpu<Float,3>& tau, Array_gpu<Float,3>& lay_source,
8181
Array_gpu<Float,3>& lev_source_inc, Array_gpu<Float,3>& lev_source_dec, Array_gpu<Float,2>& sfc_source,
8282
Array_gpu<Float,2>& lw_flux_up, Array_gpu<Float,2>& lw_flux_dn, Array_gpu<Float,2>& lw_flux_net,
@@ -136,7 +136,7 @@ class Radiation_solver_shortwave
136136
const Array<Float,2>& sfc_alb,
137137
const Array<Float,1>& tsi_scaling, const Array<Float,1>& mu0,
138138
const Array<Float,2>& lwp, const Array<Float,2>& iwp,
139-
const Array<Float,2>& rel, const Array<Float,2>& rei,
139+
const Array<Float,2>& rel, const Array<Float,2>& dei,
140140
Array<Float,3>& tau, Array<Float,3>& ssa, Array<Float,3>& g,
141141
Array<Float,2>& toa_src,
142142
Array<Float,2>& sw_flux_up, Array<Float,2>& sw_flux_dn,
@@ -174,7 +174,7 @@ class Radiation_solver_shortwave
174174
const Array_gpu<Float,1>& tsi_scaling,
175175
const Array_gpu<Float,1>& mu0, const Array_gpu<Float,1>& azi,
176176
const Array_gpu<Float,2>& lwp, const Array_gpu<Float,2>& iwp,
177-
const Array_gpu<Float,2>& rel, const Array_gpu<Float,2>& rei,
177+
const Array_gpu<Float,2>& rel, const Array_gpu<Float,2>& dei,
178178
const Array_gpu<Float,1>& land_use_map,
179179
const Array_gpu<Float,2>& rh,
180180
const Aerosol_concs_gpu& aerosol_concs,
@@ -209,7 +209,7 @@ class Radiation_solver_shortwave
209209
const Array_gpu<Float,1>& tsi_scaling,
210210
const Array_gpu<Float,1>& mu0, const Array_gpu<Float,1>& azi,
211211
const Array_gpu<Float,2>& lwp, const Array_gpu<Float,2>& iwp,
212-
const Array_gpu<Float,2>& rel, const Array_gpu<Float,2>& rei,
212+
const Array_gpu<Float,2>& rel, const Array_gpu<Float,2>& dei,
213213
const Array_gpu<Float,1>& land_use_map,
214214
const Array_gpu<Float,2>& rh,
215215
const Aerosol_concs_gpu& aerosol_concs,
@@ -250,12 +250,12 @@ class Radiation_solver_shortwave
250250
Array_gpu<Float,3> mie_cdfs_bb;
251251
Array_gpu<Float,4> mie_angs_bb;
252252
Array_gpu<Float,4> mie_phase_bb;
253-
253+
254254
Array_gpu<Float,1> mie_phase_angs;
255255
Array_gpu<Float,3> mie_cdfs;
256256
Array_gpu<Float,4> mie_angs;
257257
Array_gpu<Float,4> mie_phase;
258-
258+
259259
#endif
260260
};
261261
#endif

‎include_test/Radiation_solver_rt.h‎

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -54,7 +54,7 @@ class Radiation_solver_longwave
5454
Array_gpu<Float,2>& col_dry,
5555
const Array_gpu<Float,1>& t_sfc, const Array_gpu<Float,2>& emis_sfc,
5656
const Array_gpu<Float,2>& lwp, const Array_gpu<Float,2>& iwp,
57-
const Array_gpu<Float,2>& rel, const Array_gpu<Float,2>& rei,
57+
const Array_gpu<Float,2>& rel, const Array_gpu<Float,2>& dei,
5858
Array_gpu<Float,2>& tau, Array_gpu<Float,2>& lay_source,
5959
Array_gpu<Float,2>& lev_source_inc, Array_gpu<Float,2>& lev_source_dec, Array_gpu<Float,1>& sfc_source,
6060
Array_gpu<Float,2>& lw_flux_up, Array_gpu<Float,2>& lw_flux_dn, Array_gpu<Float,2>& lw_flux_net,
@@ -121,11 +121,11 @@ class Radiation_solver_shortwave
121121
const Array_gpu<Float,1>& tsi_scaling,
122122
const Array_gpu<Float,1>& mu0, const Array_gpu<Float,1>& azi,
123123
const Array_gpu<Float,2>& lwp, const Array_gpu<Float,2>& iwp,
124-
const Array_gpu<Float,2>& rel, const Array_gpu<Float,2>& rei,
124+
const Array_gpu<Float,2>& rel, const Array_gpu<Float,2>& dei,
125125
const Array_gpu<Float,2>& rh,
126126
const Aerosol_concs_gpu& aerosol_concs,
127-
Array_gpu<Float,2>& tot_tau_out, Array_gpu<Float,2>& tot_ssa_out,
128-
Array_gpu<Float,2>& cld_tau_out, Array_gpu<Float,2>& cld_ssa_out, Array_gpu<Float,2>& cld_asy_out,
127+
Array_gpu<Float,2>& tot_tau_out, Array_gpu<Float,2>& tot_ssa_out,
128+
Array_gpu<Float,2>& cld_tau_out, Array_gpu<Float,2>& cld_ssa_out, Array_gpu<Float,2>& cld_asy_out,
129129
Array_gpu<Float,2>& aer_tau_out, Array_gpu<Float,2>& aer_ssa_out, Array_gpu<Float,2>& aer_asy_out,
130130
Array_gpu<Float,2>& sw_flux_up, Array_gpu<Float,2>& sw_flux_dn,
131131
Array_gpu<Float,2>& sw_flux_dn_dir, Array_gpu<Float,2>& sw_flux_net,

‎rfmip/rfmip_init.py‎

Lines changed: 14 additions & 14 deletions
Original file line numberDiff line numberDiff line change
@@ -12,20 +12,20 @@
1212
# Save all the input data to NetCDF
1313
nc_file = nc.Dataset('rte_rrtmgp_input_expt_{:02d}.nc'.format(expt), mode='w', datamodel='NETCDF4', clobber=True)
1414
nc_file_rfmip = nc.Dataset('multiple_input4MIPs_radiation_RFMIP_UColorado-RFMIP-1-2_none.nc', mode='r', datamodel='NETCDF4', clobber=False)
15-
15+
1616
# Create a group for the radiation and set up the values.
1717
nc_file.createDimension('lay', nc_file_rfmip.dimensions['layer'].size)
1818
nc_file.createDimension('lev', nc_file_rfmip.dimensions['level'].size)
1919
nc_file.createDimension('y', 1)
2020
nc_file.createDimension('x', nc_file_rfmip.dimensions['site'].size)
2121
nc_file.createDimension('band_lw', band_lw)
2222
nc_file.createDimension('band_sw', band_sw)
23-
23+
2424
nc_pres_layer = nc_file.createVariable('p_lay', float_type, ('lay', 'y', 'x'))
2525
nc_pres_level = nc_file.createVariable('p_lev', float_type, ('lev', 'y', 'x'))
2626
nc_temp_layer = nc_file.createVariable('t_lay', float_type, ('lay', 'y', 'x'))
2727
nc_temp_level = nc_file.createVariable('t_lev', float_type, ('lev', 'y', 'x'))
28-
28+
2929
nc_pres_layer[:,:,:] = nc_file_rfmip.variables['pres_layer'][:,:].transpose()
3030
nc_pres_level[:,:,:] = nc_file_rfmip.variables['pres_level'][:,:].transpose()
3131

@@ -35,25 +35,25 @@
3535

3636
nc_temp_layer[:,:,:] = (nc_file_rfmip.variables['temp_layer'][expt,:,:]).transpose()
3737
nc_temp_level[:,:,:] = (nc_file_rfmip.variables['temp_level'][expt,:,:]).transpose()
38-
38+
3939
nc_surface_emissivity = nc_file.createVariable('emis_sfc', float_type, ('y', 'x', 'band_lw'))
4040
nc_surface_temperature = nc_file.createVariable('t_sfc', float_type, ('y', 'x'))
41-
41+
4242
nc_surface_emissivity[:,:] = np.tile( (nc_file_rfmip.variables['surface_emissivity'][:]) [:,None], (1, band_lw) )
4343
nc_surface_temperature[:,:] = nc_file_rfmip.variables['surface_temperature'][expt,:]
44-
44+
4545
nc_surface_albedo_dir = nc_file.createVariable('sfc_alb_dir', float_type, ('y', 'x', 'band_sw'))
4646
nc_surface_albedo_dif = nc_file.createVariable('sfc_alb_dif', float_type, ('y', 'x', 'band_sw'))
47-
47+
4848
nc_surface_albedo_dir[:,:,:] = np.tile( (nc_file_rfmip.variables['surface_albedo'][:]) [:,None], (1, band_sw) )
4949
nc_surface_albedo_dif[:,:,:] = np.tile( (nc_file_rfmip.variables['surface_albedo'][:]) [:,None], (1, band_sw) )
50-
50+
5151
nc_mu0 = nc_file.createVariable('mu0', float_type, ('y', 'x'))
5252
nc_mu0[:,:] = np.maximum(0., np.cos( np.deg2rad(nc_file_rfmip.variables['solar_zenith_angle'][:]) ) )
53-
53+
5454
nc_tsi = nc_file.createVariable('tsi', float_type, ('y', 'x'))
5555
nc_tsi[:,:] = nc_file_rfmip.variables['total_solar_irradiance'][:]
56-
56+
5757
nc_h2o = nc_file.createVariable('vmr_h2o' , float_type, ('lay', 'y', 'x'))
5858
nc_o3 = nc_file.createVariable('vmr_o3' , float_type, ('lay', 'y', 'x'))
5959
nc_co2 = nc_file.createVariable('vmr_co2' , float_type)
@@ -77,7 +77,7 @@
7777
# Profiles
7878
nc_h2o[:,:,:] = nc_file_rfmip.variables['water_vapor'][expt,:,:].transpose() * float(nc_file_rfmip.variables['water_vapor'].units)
7979
nc_o3 [:,:,:] = nc_file_rfmip.variables['ozone'][expt,:,:].transpose() * float(nc_file_rfmip.variables['ozone'].units)
80-
80+
8181
# Constants
8282
nc_co2 [:] = nc_file_rfmip.variables['carbon_dioxide_GM'][expt] * float(nc_file_rfmip.variables['carbon_dioxide_GM'].units)
8383
nc_n2o [:] = nc_file_rfmip.variables['nitrous_oxide_GM'][expt] * float(nc_file_rfmip.variables['nitrous_oxide_GM'].units)
@@ -100,13 +100,13 @@
100100
# CvH: To be removed if settings can be set.
101101
nc_lwp = nc_file.createVariable('lwp', float_type, ('lay', 'y', 'x'))
102102
nc_iwp = nc_file.createVariable('iwp', float_type, ('lay', 'y', 'x'))
103-
nc_rel = nc_file.createVariable('rei', float_type, ('lay', 'y', 'x'))
104-
nc_rei = nc_file.createVariable('rel', float_type, ('lay', 'y', 'x'))
103+
nc_rel = nc_file.createVariable('dei', float_type, ('lay', 'y', 'x'))
104+
nc_dei = nc_file.createVariable('rel', float_type, ('lay', 'y', 'x'))
105105

106106
nc_lwp[:,:,:] = 0.
107107
nc_iwp[:,:,:] = 0.
108108
nc_rel[:,:,:] = 0.
109-
nc_rei[:,:,:] = 0.
109+
nc_dei[:,:,:] = 0.
110110
# CvH end.
111111

112112
nc_file_rfmip.close()

‎rrtmgp-data‎

‎rte-rrtmgp‎

0 commit comments

Comments
 (0)