<div dir="ltr"><div><div><div>Dear NCL<br><br></div>In the given script i need to divide the variable with 1000 . I know its simple but i am getting stuck with this query in the given script<br><br></div>I also tried <br></div>var1= cloud_water_layer / 1000 <br><br><div><div>but did not happen anything<br></div><div>please suggest me <br></div><div><br><div><br>;Reading HDF4-SDS files<br>;*********** Load Libraries ************************************<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot; <br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot; <br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl&quot; <br>;**************************************************************<br>begin<br>;***************************************************************<br>; User Input<br>;***************************************************************<br>                                             ; INPUT<br>   diri    = &quot;./&quot;                            ; input directory<br>   fili    = &quot;2A12.090930.67647.6.HDF&quot;       ;       file<br><br>                                             ; OUTPUT<br>;***************************************************************<br>; Read hdf <br>;***************************************************************<br><br>    f   = addfile(&quot;/home/kunal/Pictures/NCL_Scripts/2A12.090930.67647.6.HDF&quot;,&quot;r&quot;)<br> <br><br>   lat  = f-&gt;geolocation(:,:,0)            ; (scan, pixel)<br>   lon  = f-&gt;geolocation(:,:,1)            <br><br>   lat@units = &quot;degrees_north&quot;             ; add units<br>   lon@units = &quot;degrees_east&quot;<br>   print(lat@hdf_name)<br><br>;-------------------------------<br>    var=f-&gt;cldWater(:,:,:)<br>       <br>    var@lon2d=lon<br>    var@lat2d=lat<br>    var_name=var@hdf_name<br>    var_unit=&quot;gm/m~S~3&quot;       <br>        printVarSummary(var)<br>        print(var_unit)<br>;       print(var)             <br><br><br>        cloud_water=short2flt(var)<br>    asciiwrite(&quot;cloud_water.txt&quot;,cloud_water)<br>        printVarSummary(cloud_water)<br>    <br>    layer =2<br>    cloud_water_layer=cloud_water(:,:,layer)<br>       <br>;-------------------------------<br>  printMinMax(lat, True)<br>  printMinMax(lon, True)<br>  printMinMax (var, True)<br>;-------------------------------<br>  wks_type = &quot;pdf&quot;                       <br>  wks_type@wkPaperSize = &quot;A4&quot;<br>  wks = gsn_open_wks(wks_type,&quot;2004&quot;)<br><br>  res                   = True<br>  res@tiMainOn          = True <br>  res@gsnMaximize       = True     ; maximize plot in frame<br>  res@cnFillOn          = True     ; turn on contour fill<br>  res@cnLinesOn         = False    ; turn off contour lines<br>  res@cnLineLabelsOn    = False    ; turn off line labels<br><br><br><br>res@mpMinLonF                   = 60<br>res@mpMaxLonF                   = 100<br>res@mpMinLatF                   = 0<br>res@mpMaxLatF                   = 30<br>res@mpFillOn                    = True <br>res@mpOutlineOn                 = True<br>res@mpOutlineBoundarySets       = &quot;National&quot;<br>res@mpGeophysicalLineColor      = &quot;Blue&quot;    ;Color of Conuntry Outline<br>res@mpNationalLineThicknessF    = 3.5       ;thickness of Country Outline<br>res@cnLevelSelectionMode = &quot;ExplicitLevels&quot;   ; set explicit contour levels<br>res@cnLevels             = (/0,1e+04,2e+04,3e+04,4e+04,5e+04,6e+04,7e+04,8e+04,9e+04/)<br>;res@cnLevelSelectionMode=&quot;AutomaticLevels&quot;<br>;res@cnFillPalette = &quot;BlAqGrYeOrReVi200&quot;<br>res@cnFillPalette = &quot;precip3_16lev&quot;<br><br>plot=gsn_csm_contour_map(wks,cloud_water_layer,res)<br><br>end  <br><br><br clear="all"><div><div><div class="gmail_signature"><div dir="ltr"><div><br></div><br><div><br></div><div><br></div><div><br></div><div><br></div><div><p style="margin:0px;border-collapse:collapse;font-family:Tahoma,Verdana;font-size:12px"><font color="#1F497D"><br></font></p></div></div></div></div>
</div></div></div></div></div>