[ncl-talk] ESMF Regridding of WRF output files - issues on variable coordinates
Mirce Morales
mirce.morales at gmail.com
Thu Dec 6 09:03:51 MST 2018
Hi Laura,
Thank you very much for the advice, it really helps me a lot to understand
what is happening.
I will look at the description of the filedimdef as you tell me.
Thanks again!.
Cheers,
Mirce.
El jue., 6 dic. 2018 a las 8:19, Laura Fowler (<laura at ucar.edu>) escribió:
> Hi Mirce:
> I think that if you look at the description of filedimdef, you will
> understand how to specify the dimensions that you would like to see in your
> regridded netcdf file. When using filedimdef, then the actual named
> dimensions of your variables will be written instead of all the ncl01,
> ncl02, etc dimensions. Hope this helps.
>
> Cheers,
> Laura
>
>
>
>
>
>
> 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
>>
>
>
> --
>
> !-------------------------------------------------------------------------------------------------------------
> Laura D. Fowler
>
> Mesoscale and Microscale Meteorology Division (MMM)
> National Center for Atmospheric Research
> P.O. Box 3000, Boulder CO 80307-3000
>
> e-mail: laura at ucar.edu
> phone: 303-497-1628
>
>
> !-------------------------------------------------------------------------------------------------------------
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20181206/b5c49bf8/attachment.html>
More information about the ncl-talk
mailing list