| Re: storing in array [message #60223 is a reply to message #60221] |
Mon, 12 May 2008 08:35   |
kishore1818
Messages: 24
Registered: June 2007
|
Junior Member |
|
|
On May 12, 8:18 am, Spon <christoph.b...@gmail.com> wrote:
> On May 12, 11:03 am, Spon <christoph.b...@gmail.com> wrote:
>
>> Output = { feature_type:"", $
>> feature_type_qa:"", $
>> ice_water_phase:"", $
>> ice_water_phase_qa:"", $
>> feature_subtype:"", $
>> ; etc. , $
>> }
>> case feature_type of
>> 0 : Output.FeatureType = 'invalid (bad or missing data)'
>> 1 : Output.FeatureType = 'clear air'
>> 2 : begin
>> Output.FeatureType = 'cloud'
>> case feature_subtype of
>> 0 : Output.FeatureSubType = 'low overcast, transparent'
>> 1 : Output.FeatureSubType = 'low overcast, opaque'
>
> Obviously, your variable names have to match up (including structure
> field names), which they don't in my example, sorry.
>
>> 0 : Output.FeatureType = 'invalid (bad or missing data)'
>> 1 : Output.FeatureSubType = 'low overcast, opaque'
>
> and so on
> should be
>
>> 0 : Output.Feature_Type = 'invalid (bad or missing data)'
>> 1 : Output.Feature_SubType = 'low overcast, opaque'
>
> to make my example work.
Hi,
Thanks for your suggestions, I modified according your suggestions.
Pl. go through my program
pro vfm_feature_flags1,val
; this routine demonstrates how to read and extract values from a
feature
; classification flag 16-bit integer value in CALIPSO Level 2
Vertical
; Feature Mask files
;
; INPUT:
; val - the feature classification flag value to be decoded
;
; OUTPUT:
; all information is printed into the IDL log window
print, val
feature_type = 0
feature_type_qa = 0
ice_water_phase = 0
ice_water_phase_qa = 0
feature_subtype = 0
cloud_aerosol_psc_type_qa = 0
horizontal_averaging = 0
Output = { feature_type:"", $
feature_type_qa:"", $
ice_water_phase:"", $
ice_water_phase_qa:"", $
feature_subtype:"", $
cloud_aerosol_psc_type_qa:"", $
horizontal_averaging:""}
for i=0,15 do begin
if ((val and 2L^i) NE 0) then begin
print,'Bit set: ',i+1
case i+1 of
1 : feature_type = feature_type + 1
2 : feature_type = feature_type + 2
3 : feature_type = feature_type + 4
4 : feature_type_qa = feature_type_qa + 1
5 : feature_type_qa = feature_type_qa + 2
6 : ice_water_phase = ice_water_phase + 1
7 : ice_water_phase = ice_water_phase + 2
8 : ice_water_phase_qa = ice_water_phase_qa + 1
9 : ice_water_phase_qa = ice_water_phase_qa + 2
10 : feature_subtype = feature_subtype + 1
11 : feature_subtype = feature_subtype + 2
12 : feature_subtype = feature_subtype + 4
13 : cloud_aerosol_psc_type_qa = cloud_aerosol_psc_type_qa + 1
14 : horizontal_averaging = horizontal_averaging + 1
15 : horizontal_averaging = horizontal_averaging + 2
16: horizontal_averaging = horizontal_averaging + 4
else:
endcase
endif
endfor
case feature_type of
0 : Output.Feature_Type = 'invalid (bad or missing data)'
1 : Output.Feature_Type = 'clear air"
2 : begin
Output.Feature_Type = 'cloud'
case feature_subtype of
0 : Output.Feature_Subtype = 'low overcast, transparent'
1 : Output.Feature_Subtype = 'low overcast, opaque'
2 : Output.Feature_Subtype = 'transition stratocumulus'
3 : Output.Feature_Subtype = 'low, broken cumulus'
4 : Output.Feature_Subtype = 'altocumulus (transparent)'
5 : Output.Feature_Subtype = 'altostratus (opaque)'
6 : Output.Feature_Subtype = 'cirrus (transparent)'
7 : Output.Feature_Subtype = 'deep convective (opaque)'
else : Output,'*** error getting Feature Subtype'
endcase
end
3 : begin
Output.Feature_Type = 'aerosol'
case feature_subtype of
0 : Output.Feature_Subtype = 'not determined'
1 : Output.Feature_Subtype = 'clean marine'
2 : Output.Feature_Subtype = 'dust'
3 : Output.Feature_Subtype = 'polluted continental'
4 : Output.Feature_Subtype = 'clean continental'
5 : Output.Feature_Subtype = 'polluted dust'
6 : Output.Feature_Subtype = 'smoke'
7 : Output.Feature_Subtype = 'other"
else : Output,'*** error getting Feature Subtype'
endcase
end
4 : begin
Output.Feature_Type = 'stratospheric feature--PSC or
stratospheric aerosol'
case feature_subtype of
0 : Output.Feature_Subtype = 'not determined'
1 : Output.Feature_Subtype = 'non-depolarizing PSC'
2 : Output.Feature_Subtype = 'depolarizing PSC'
3 : Output.Feature_Subtype = 'non-depolarizing aerosol'
4 : Output.Feature_Subtype = 'depolarizing aerosol'
5 : Output.Feature_Subtype = 'spare'
6 : Output.Feature_Subtype = 'spare'
7 : Output.Feature_Subtype = 'other'
else : Output,'*** error getting Feature Subtype'
endcase
end
5 : Output.Feature_Type = 'surface'
6 : Output.Feature_Type = 'subsurface'
7 : Output.Feature_Type = 'no signal (totally attenuated)'
else : Output,'*** error getting Feature Type'
endcase
case feature_type_qa of
0 : Output.Feature_Type_QA = 'none'
1 : Output.Feature_Type_QA = 'low'
2 : Output.Feature_Type_QA = 'medium'
3 : Output.Feature_Type_QA = 'high'
else : Output,'*** error getting Feature Type QA'
endcase
case ice_water_phase of
0 : Output.Ice_Water_Phase = 'unknown/not determined'
1 : Output.Ice_Water_Phase = 'ice'
2 : Output.Ice_Water_Phase = 'water'
3 : Output.Ice_Water_Phase = 'mixed phase'
else : Output,'*** error getting Ice/Water Phase'
endcase
case ice_water_phase_qa of
0 : Output.Ice_Water_Phase_QA = 'none'
1 : Output.Ice_Water_Phase_QA = 'low'
2 : Output.Ice_Water_Phase_QA = 'medium'
3 : Output.Ice_Water_Phase_QA = 'high'
else : Output,'*** error getting Ice/Water Phase QA'
endcase
if (cloud_aerosol_psc_type_qa eq 0) then begin
Output.Cloud_Aerosol_PSC_Type_QA = 'not confident'
endif else begin
Output.Cloud_Aerosol_PSC_Type_QA = 'confident'
endelse
case horizontal_averaging of
0 : Output.Horizontal_averaging = 'not applicable'
1 : Output.Horizontal_averaging = '1/3 km'
2 : Output.Horizontal_averaging = '1 km'
3 : Output.Horizontal_averaging = '5 km'
4 : Output.Horizontal_averaging = '20 km'
5 : Output.Horizontal_averaging = '80 km'
else : Output,'*** error getting Horizontal averaging'
endcase
help,Output,/structure
;Return, Output
end
Using above routine I wrote a small progam
data1=[44474, 36282, 28602, 28090]
for i=0,3 do begin
;vfm_feature_flags,data1(i)
vfm_feature_flags1,data1(i)
pause='' & read,pause
endfor
end
Output
FEATURE_SUBTYPE STRING 'cirrus (transparent)'
FEATURE_SUBTYPE STRING 'cirrus (transparent)'
FEATURE_SUBTYPE STRING 'deep convective (opaque)'
FEATURE_SUBTYPE STRING 'cirrus (transparent)'
Only, problem is how to recall output parameters.
For example, I want to recall all output.FEATURE_SUBTYPE.
These output I have to calculate how many cirrus and deep convective
system around the globe at each day at different Longitude latitude
bands.
Once again I am very much thankful to you.
Regards,
Kishore
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
