[ncl-talk] ESMF Regridding of WRF output files - issues on variable coordinates
Mirce Morales
mirce.morales at gmail.com
Thu Dec 6 08:57:46 MST 2018
Hi Dennis,
I am sorry for the confusing email. I will do the changes you suggest me.
Thank you very much,
Mirce.
El jue., 6 dic. 2018 a las 8:15, Dennis Shea (<shea at ucar.edu>) escribió:
> Your email is a confusing.
>
> Reset the following:
> opt at CopyVarCoords = True
>
> ===
> * XLAT_SOURCE = sfile->XLAT*
> * XLONG_SOURCE = sfile->XLONG*
>
> ; writing to netCDF
>
> system(...) ; remove any preexisting file
> fout = addfile(... ,"c")
>
> fout->T2 = ... regrid variable...
>
> * fout->XLAT = XLAT_SOURCE*
> * fout->XLONG = XLONG_SOURCE*
>
>
> On Wed, Dec 5, 2018 at 10:08 PM Mirce Morales <mirce.morales at gmail.com>
> wrote:
>
>> Hello all,
>>
>> I am trying to regrid some variables in my WRF output files (wrfout_d0*),
>> to do it, I am taking as a guide the following script on the NCL web page:
>> https://www.ncl.ucar.edu/Applications/Scripts/ESMF_regrid_20.ncl
>>
>> Up to now, I was able to run the script without error, but I see that if
>> I do ncdump -h "My regridded file", there is like a matrix of numbers where
>> the XLAT and XLONG should be, so it is not possible to see the variables
>> and their attributes well....to avoid that, I found out I have to set
>> opt at CopyVarCoords = False.
>>
>> The problem is that when I do that, now I get the following kind of
>> attributes in my regridded files:
>>
>> *netcdf wrfout_d03_2010-09-01_00\:00\:00 {*
>> *dimensions:*
>> * ncl0 = 1 ;*
>> * ncl1 = 156 ;*
>> * ncl2 = 138 ;*
>> * ncl3 = 1 ;*
>> * ncl4 = 156 ;*
>> * ncl5 = 138 ;*
>> * ncl6 = 1 ;*
>> * ncl7 = 156 ;*
>> * ncl8 = 138 ;*
>> * ncl9 = 1 ;*
>> * ncl10 = 156 ;*
>> * ncl11 = 138 ;*
>> * ncl12 = 1 ;*
>> * ncl13 = 156 ;*
>> * ncl14 = 138 ;*
>> * ncl15 = 1 ;*
>> * ncl16 = 156 ;*
>> * ncl17 = 138 ;*
>> * ncl18 = 1 ;*
>> * ncl19 = 156 ;*
>> * ncl20 = 138 ;*
>> * ncl21 = 1 ;*
>> * ncl22 = 156 ;*
>> * ncl23 = 138 ;*
>> * ncl24 = 1 ;*
>> * ncl25 = 156 ;*
>> * ncl26 = 138 ;*
>> * ncl27 = 1 ;*
>> * ncl28 = 156 ;*
>> * ncl29 = 138 ;*
>> * ncl30 = 1 ;*
>> * ncl31 = 156 ;*
>> * ncl32 = 138 ;*
>> *variables:*
>> * float T2(ncl0, ncl1, ncl2) ;*
>> * T2:missing_value = 9.96921e+36f ;*
>> * T2:remap = "remapped via ESMF_regrid_with_weights:
>> Bilinear" ;*
>> * T2:_FillValue = 9.96921e+36f ;*
>> * T2:FieldType = 104 ;*
>> * T2:MemoryOrder = "XY " ;*
>> * T2:description = "TEMP at 2 M" ;*
>> * T2:units = "K" ;*
>> * T2:stagger = "" ;*
>>
>> When it should be similar to this:
>>
>> *netcdf wrfout_d03_2010-09-01_01\:00\:00 {*
>> *dimensions:*
>> * Time = UNLIMITED ; // (1 currently)*
>> * DateStrLen = 19 ;*
>> * west_east = 54 ;*
>> * south_north = 60 ;*
>> * bottom_top = 49 ;*
>> * bottom_top_stag = 50 ;*
>> * soil_layers_stag = 4 ;*
>> * west_east_stag = 55 ;*
>> * south_north_stag = 61 ;*
>> *variables:*
>> *float T2(Time, south_north, west_east) ;*
>> * T2:FieldType = 104 ;*
>> * T2:MemoryOrder = "XY " ;*
>> * T2:description = "TEMP at 2 M" ;*
>> * T2:units = "K" ;*
>> * T2:stagger = "" ;*
>> * T2:coordinates = "XLONG XLAT XTIME" ;*
>>
>> I don't know how to fix this....In case it helps, these are the
>> regridding options I am using:
>>
>> *;---Open one WRF file *
>> * sfile =
>> addfile("../IMTAProy_BMJ/wrfout_d03_2010-10-31_16:00:00","r") ; contains
>> source grid*
>>
>> *;---Read lat/lon of source grid*
>> * lat_src = sfile->XLAT(0,:,:)*
>> * lon_src = sfile->XLONG(0,:,:)*
>>
>> *;---Read lat/lon of destination grid*
>> * lat_dst = dstfile->XLAT_M(0,:,:)*
>> * lon_dst = dstfile->XLONG_M(0,:,:)*
>>
>> *;---Options to pass to ESMF_regrid*
>>
>> * interp_method = "bilinear"*
>> * opt = True*
>>
>> * opt at SrcGridLat = lat_src*
>> * opt at SrcGridLon = lon_src*
>> * opt at DstGridLat = lat_dst*
>> * opt at DstGridLon = lon_dst*
>>
>> * opt at SrcRegional = True ; Default is False*
>> * opt at DstRegional = True*
>>
>> * opt at CopyVarCoords = False ; This is the default in V6.1.0 (not
>> 6.1.0-beta)*
>> * opt at ForceOverwrite = True*
>>
>> *;---Debug information*
>> * opt at PrintTimings = True*
>> * opt at Debug = True*
>>
>> In summary, I have the following doubts:
>>
>> 1) How can I keep the coordinates of each variable in a proper way?
>> 2) How can I keep the dimensions exactly as in my source files?
>> 3) To create the weight files in a separate script, is there any template
>> I should base on for this case?
>> 4) How can I add XLAT and XLONG in my regridded files like they are in
>> the original wrfout* files?...like this:
>>
>> *variables:*
>> * char Times(Time, DateStrLen) ;*
>> * float XLAT(Time, south_north, west_east) ;*
>> * XLAT:FieldType = 104 ;*
>> * XLAT:MemoryOrder = "XY " ;*
>> * XLAT:description = "LATITUDE, SOUTH IS NEGATIVE" ;*
>> * XLAT:units = "degree_north" ;*
>> * XLAT:stagger = "" ;*
>> * XLAT:coordinates = "XLONG XLAT" ;*
>> * float XLONG(Time, south_north, west_east) ;*
>> * XLONG:FieldType = 104 ;*
>> * XLONG:MemoryOrder = "XY " ;*
>> * XLONG:description = "LONGITUDE, WEST IS NEGATIVE" ;*
>> * XLONG:units = "degree_east" ;*
>> * XLONG:stagger = "" ;*
>> * XLONG:coordinates = "XLONG XLAT" ;*
>>
>> I will really appreciate any help/clue about this.
>>
>> In case it is helpful my script and files are in the following path in
>> Cheyenne:
>>
>> */glade/scratch/mirce/LaSierra/FullDomainForcing/WRF/second_option*
>>
>> If it would be better to paste here my the script I will be happy to do
>> it!
>>
>> Thanks a lot in advance!
>> Mirce.
>>
>>
>>
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20181206/e173ff9d/attachment.html>
More information about the ncl-talk
mailing list