<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Dear All,<br class=""><div class=""><div class=""><div class=""><div class=""><div class=""><div class=""><br class=""></div><div class="">I think this is more WRF related issue but in any case I am writing to both group.<br class=""><br class=""></div>I have downloaded sea-ice data from NOAA/NSIDC in netcdf format from this page:<br class=""><br class=""><a href="http://nsidc.org/data/G02202" class="">http://nsidc.org/data/G02202</a><br class=""><br class=""></div>I am trying to convert and write this data to a WPS intermediate file format so that I can produce required met files for WRF simulation. I am following the steps described in NCL functions page:<br class=""><br class=""><a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/wrf_wps_write_int.shtml" class="">https://www.ncl.ucar.edu/Document/Functions/Built-in/wrf_wps_write_int.shtml</a><br class=""><br class=""></div>and my script is given at the end of this email. Everything seems to work fine and WPS intermediate file for sea-ice is produced with logical values however all sea-ice values in the final met_em files are 0. My question is has anyone successfully tried to convert one of these data set to WPS intermediate file, and then produced met_em files properly?<br class=""><br class=""></div>I would be grateful if someone could help me address my problem.<br class=""><br class=""></div>Thanks in advance,<br class=""></div>Deniz<br class=""><br class=""><div class="" style="word-wrap: break-word;"><span class="" style="border-collapse: collapse; font-family: Arial;"><font color="#7c7c7c" class="">--</font></span><div class="" style="border-collapse: separate;"><font face="Arial" class="" color="#7c7c7c">Deniz Bozkurt, PhD<br class="">Postdoctoral Researcher</font></div><div class="" style="border-collapse: separate;"><font face="Arial" class="" color="#7c7c7c">Center for Climate and Resilience Research (CR2)<span class=""> <</span><span class=""><a href="http://www.cr2.cl/" class="">www.cr2.cl</a>></span></font></div></div><font color="#7c7c7c" class=""><span class="" style="border-collapse: collapse; font-family: Arial;">Department of Geophysics</span>, University of Chile<span class="" style="border-collapse: collapse; font-family: Arial;"></span></font><div class=""><font color="#7c7c7c" class=""><span class="" style="border-collapse: collapse; font-family: Arial;"></span></font><div class="" style="word-wrap: break-word;"><div class="" style="font-size: 11px; border-collapse: separate;"><span class="" style="border-collapse: collapse;"><font face="Arial" class="" color="#7c7c7c">Blanco Encalada 2002, Santiago, Chile</font></span></div><div class="" style="font-size: 11px; border-collapse: separate;"><font color="#7c7c7c" class=""><br class=""></font></div><div class="" style="font-size: 11px; border-collapse: separate;"><font color="#7c7c7c" class="">Mobil: +56 (9) 56107310</font></div><div class="" style="font-size: 11px; border-collapse: separate;"><font color="#9a9a9a" class=""><a href="https://sites.google.com/site/deniboz" class="">https://sites.google.com/site/deniboz</a></font></div><div class="" style="font-size: 11px; border-collapse: separate;"><br class=""></div><div class=""><br class=""></div></div></div><br class=""><div class=""><br class="">f=addfile("<a href="http://seaice_conc_daily_sh_n07_19810103_v03r01.nc" class="">seaice_conc_daily_sh_n07_19810103_v03r01.nc</a>","r")<br class="">dd=f->goddard_bt_seaice_conc<br class="">printVarSummary(dd)<br class="">latd = f->latitude<br class="">lond = f->longitude<br class="">printVarSummary(lond)<br class="">FIELD_T = "SEAICE"<br class="">UNITS_T = "%"<br class="">DESC_T = "Sea Ice Fraction"<br class="">DATE = "1981-01-03"<br class="">WPS_IM_root_name = "SEAICE"<br class=""><br class="">output_file_name = WPS_IM_root_name + ":" + DATE<br class="">earth_radius = 6367.470<br class=""><br class=""><br class="">opt = True<br class="">opt@level = (/200100/)<br class="">opt@projection = 0<br class="">opt@date = DATE<br class="">opt@map_source = "NSIDC"<br class="">opt@startloc = "SWCORNER" ; 8 chars exact<br class="">opt@startlon = lond(0,0);0<br class="">opt@startlat = latd(0,0);90<br class="">opt@deltalon = lond(0,1)-latd(0,0);0.25<br class="">opt@deltalat = latd(1,0)-lond(0,0);-0.25<br class="">opt@is_wind_earth_relative = False<br class=""><br class="">wrf_wps_write_int(WPS_IM_root_name,FIELD_T,UNITS_T,DESC_T,dd(0,:,:),opt)</div><div class=""><br class=""></div><div class=""><div class="" style="word-wrap: break-word;"><br class=""></div><div class=""><div class="" style="word-wrap: break-word;"><div class="" style="color: rgb(132, 132, 132); font-size: 11px; border-collapse: separate;"><br class=""></div></div></div></div><br class=""><br class="">
<br class=""></body></html>