All output variables are zero

[Amirhossein Nikfal]

Hi,
I’ve come across a strange issue that I haven’t seen before. There were no runtime errors, but all the output variables (BRIGHTNESS TEMPERATURES, RADIANCES, etc.) are zero.

here’s my input data:
`
! jj, ii, 1, 1
! jjSize, iiSize, 62, 62
! Specify input profiles for example_fwd.F90.
! Multiple profiles may be described: follow the same format for each one.
! Comment lines (starting with ‘!’) are optional.
!
! Gas units (must be same for all profiles)
! 0 => ppmv over dry air
! 1 => kg/kg over moist air
! 2 => ppmv over moist air
!
1
!
! — Start of profile —
! — Grid point with j and i = 1/62 and 1/62 —
!
! Pressure levels (hPa)
!
52.252583
56.960472
62.092754
67.687355
73.78628
80.4347
87.682014
95.58242
104.19444
113.582794
123.81654
134.97246
147.13368
160.39107
174.84273
190.5962
207.77011
226.49106
246.89626
269.12213
293.32138
319.6699
348.36273
379.59485
413.58786
450.5839
490.87122
534.7854
580.91876
626.928
671.906
715.59625
757.4932
796.7798
832.67065
864.6287
892.4068
916.0529
935.81696
952.04553
965.2392
975.81494
984.2149
990.85
!
! Temperature profile (K)
!
206.16922
204.01598
201.86288
199.70973
198.51907
197.94205
197.36493
196.7879
197.62273
200.01039
202.39787
205.62079
208.94731
211.60233
214.06372
218.67006
223.2643
227.71817
231.01299
234.77167
238.60544
243.40344
248.54108
252.12463
255.78265
259.69568
263.43457
267.7721
273.23196
278.90967
284.08133
288.56808
292.36502
295.5747
297.7629
299.4179
301.006
302.88007
304.68307
306.1949
307.4505
308.5171
310.20566
311.52875
!
! Water vapour profile (kg/kg)
!
0.00001
0.00001
0.00001
0.00001
0.00001
0.00001
0.00001
0.00001
0.00001
0.00001
0.00001
0.00001
0.00001
0.00001
0.00002
0.00002
0.00003
0.00007
0.00016
0.00026
0.00038
0.00062
0.00103
0.00076
0.00043
0.00035
0.00098
0.00253
0.00338
0.00331
0.00305
0.00286
0.00320
0.00402
0.00430
0.00375
0.00344
0.00343
0.00345
0.00346
0.00349
0.00359
0.00440
0.00511
!
! Near-surface variables:
! 2m T (K) 2m q (kg/kg) 2m p (hPa) 10m wind u (m/s) 10m wind v (m/s) wind fetch (m)
!
312.47952 0.005113475 990.85 -0.16482714 -5.130578 100000
!
! Skin variables:
! Skin T (K) Salinity FASTEM parameters for land surfaces
!
329.55588 35.0 3.0 5.0 15.0 0.1 0.3
!
! Surface type (0=land, 1=sea, 2=sea-ice) and water type (0=fresh, 1=ocean)
!
0 1
!
! Elevation (km), latitude and longitude (degrees)
!
0.14330047607421875 24.445251 33.040924
!
! Sat. zenith and azimuth angles, solar zenith and azimuth angles (degrees)
!
66.63 37.49 58.13 -92.77
!
! Cloud top pressure (hPa) and cloud fraction for simple cloud scheme
!
500 0.21150017
!
! — End of profile j and i = 1/62 and 1/62 —
!
`

and here’s the output:
`
—————–
Instrument oli
—————–

RTTOV options structure
General config options
apply_reg_limits F
verbose T
do_checkinput T
fix_hgpl F

General RT options
addrefrac T
plane_parallel F
switchrad F
use_q2m T
do_lambertian F
rad_down_lin_tau T
dtau_test T

VIS/IR-only RT options
ozone_data F
co2_data F
n2o_data F
co_data F
ch4_data F
so2_data F
ir_sea_emis_model 2
addsolar F
do_nlte_correction F
addaerosl F
addclouds F

PC-RTTOV options
addpc F

HTFRTC options
htfrtc F

MW-only RT options
clw_data F
fastem_version 6
supply_foam_fraction F
apply_band_correction T

Interpolation and vertical grid options
addinterp T
interp_mode 1
lgradp F
reg_limit_extrap F
spacetop T

Profile 1

RTTOV profile structure
id
date 1950/01/01
time 00:00:00
number of levels 44
number of layers 43
Viewing conditions
zenangle 66.630
azangle 37.490
sunzenangle 58.130
sunazangle -92.770
elevation 0.143
latitude 24.445
longitude 33.041
Earth magnetic field 0.000000E+00
cos of Earth magnetic field and wave propagation 0.000000E+00
Black body cloud
cloud top pressure (hPa) 0.500000E+03
cloud fraction (0 – 1) 0.211500E+00
Skin parameters
surface type 0=land, 1=sea, 2=sea-ice 0
water type 0=fresh water, 1=ocean water 1
radiative skin temperature (K) 0.329556E+03
specularity (0 – 1) 0.000000E+00
snow_fraction (0 – 1) 0.000000E+00
soil_moisture (m^3/m^3) 0.000000E+00
salinity (%o) 0.350000E+02
foam fraction (0 – 1) 0.000000E+00
land/sea-ice surface parameters for fastem
3.0000 5.0000 15.0000 0.1000 0.3000
2m air parameters
surface pressure (hPa) 0.990850E+03
temperature (K) 0.312480E+03
water vapour (kg/kg) 0.511348E-02
ozone (kg/kg) 0.000000E+00
U wind component (m/s) -0.164827E+00
V wind component (m/s) -0.513058E+01
Wind fetch (m) 0.100000E+06

Gas units: kg/kg wrt moist air
level Pressure Temp WV (kg/kg)
1 52.2526 0.206169E+03 0.100000E-04
2 56.9605 0.204016E+03 0.100000E-04
3 62.0928 0.201863E+03 0.100000E-04
4 67.6874 0.199710E+03 0.100000E-04
5 73.7863 0.198519E+03 0.100000E-04
6 80.4347 0.197942E+03 0.100000E-04
7 87.6820 0.197365E+03 0.100000E-04
8 95.5824 0.196788E+03 0.100000E-04
9 104.1944 0.197623E+03 0.100000E-04
10 113.5828 0.200010E+03 0.100000E-04
11 123.8165 0.202398E+03 0.100000E-04
12 134.9725 0.205621E+03 0.100000E-04
13 147.1337 0.208947E+03 0.100000E-04
14 160.3911 0.211602E+03 0.100000E-04
15 174.8427 0.214064E+03 0.200000E-04
16 190.5962 0.218670E+03 0.200000E-04
17 207.7701 0.223264E+03 0.300000E-04
18 226.4911 0.227718E+03 0.700000E-04
19 246.8963 0.231013E+03 0.160000E-03
20 269.1221 0.234772E+03 0.260000E-03
21 293.3214 0.238605E+03 0.380000E-03
22 319.6699 0.243403E+03 0.620000E-03
23 348.3627 0.248541E+03 0.103000E-02
24 379.5949 0.252125E+03 0.760000E-03
25 413.5879 0.255783E+03 0.430000E-03
26 450.5839 0.259696E+03 0.350000E-03
27 490.8712 0.263435E+03 0.980000E-03
28 534.7854 0.267772E+03 0.253000E-02
29 580.9188 0.273232E+03 0.338000E-02
30 626.9280 0.278910E+03 0.331000E-02
31 671.9060 0.284081E+03 0.305000E-02
32 715.5963 0.288568E+03 0.286000E-02
33 757.4932 0.292365E+03 0.320000E-02
34 796.7798 0.295575E+03 0.402000E-02
35 832.6707 0.297763E+03 0.430000E-02
36 864.6287 0.299418E+03 0.375000E-02
37 892.4068 0.301006E+03 0.344000E-02
38 916.0529 0.302880E+03 0.343000E-02
39 935.8170 0.304683E+03 0.345000E-02
40 952.0455 0.306195E+03 0.346000E-02
41 965.2392 0.307450E+03 0.349000E-02
42 975.8149 0.308517E+03 0.359000E-02
43 984.2149 0.310206E+03 0.440000E-02
44 990.8500 0.311529E+03 0.511000E-02

CHANNELS PROCESSED FOR SAT landsat 8
5 6

CALCULATED BRIGHTNESS TEMPERATURES (K):
0.00 0.00

CALCULATED RADIANCES (mW/m2/sr/cm-1):
0.00 0.00

CALCULATED OVERCAST RADIANCES:
0.00 0.00

CALCULATED SURFACE TO SPACE TRANSMITTANCE:
0.0000 0.0000

CALCULATED SURFACE EMISSIVITIES:
0.000 -0.025

Level to space transmittances for channels
5 6
1 0.0000 0.0000
2 0.0000 0.0000
3 0.0000 0.0000
4 0.0000 0.0000
5 0.0000 0.0000
6 0.0000 0.0000
7 0.0000 0.0000
8 0.0000 0.0000
9 0.0000 0.0000
10 0.0000 0.0000
11 0.0000 0.0000
12 0.0000 0.0000
13 0.0000 0.0000
14 0.0000 0.0000
15 0.0000 0.0000
16 0.0000 0.0000
17 0.0000 0.0000
18 0.0000 0.0000
19 0.0000 0.0000
20 0.0000 0.0000
21 0.0000 0.0000
22 0.0000 0.0000
23 0.0000 0.0000
24 0.0000 0.0000
25 0.0000 0.0000
26 0.0000 0.0000
27 0.0000 0.0000
28 0.0000 0.0000
29 0.0000 0.0000
30 0.0000 0.0000
31 0.0000 0.0000
32 0.0000 0.0000
33 0.0000 0.0000
34 0.0000 0.0000
35 0.0000 0.0000
36 0.0000 0.0000
37 0.0000 0.0000
38 0.0000 0.0000
39 0.0000 0.0000
40 0.0000 0.0000
41 0.0000 0.0000
42 0.0000 0.0000
43 0.0000 0.0000
44 0.0000 0.0000
5 6
`

[James Hocking]

Hi,

It looks like you are simulating the Landsat OLI instrument: all channels on this sensor are VIS/NIR (below 3 microns) which means they have no thermally emitted component in RTTOV simulations. Since you have not enabled solar radiation in the simulation (opts%rt_ir%addsolar), the outputs are zero.

Best wishes,
James

[Amirhossein Nikfal]

Thank you very much. Therefore, based on this criterion (below 3 microns), solar radiation should be enabled in the simulation for all Landsat-8 bands, including the SWIR bands.

Best regards,
Amir

[James Hocking]

Hi Amir,

Yes, to simulate any channel on Landsat OLI you need to enable solar radiation.

Best wishes,
James

[Author]

Posts