[ncl-talk] phase difference in plots

Tue Jan 23 21:53:52 MST 2018

Sajede,

First, when you have continuing questions about NCL, please don't
explicitly email somebody who has helped you before, because then this
makes it their problem. Please just email ncl-talk at ucar.edu.

Your script is hard to follow because you have a lot of code commented out,
and I can't tell what code is actually getting executed. You must provide a
clean script that we can follow, otherwise it's very time-consuming to look
it over.

What it looks like is happening is you are opening the series of files in
"cat" mode, which causes the variables to be concatenated. This is because
it's assuming that you are opening a series of files in time, where each
file represents a different time or series of times. This is what your 2D
latitude array is now 3200 by 360, because it's the leftmost dimension that
gets concatenated every time.

In your case, I think you want to open the files in "join" mode, which will
cause a new leftmost dimension to be prepended to all the variables, that
will be the size of the number of files you have. So, if you have 16 files,
your lat/lon variables will now be 16 x 200 x 360.

I'm not sure what you mean by "phase plot", or exactly what kind of plot
you are trying to create from the 16 files, because again, the code is
confusing to follow. If you simply want to average SST across the leftmost
dimension, then:

SSTmean=dim_avg_n_Wrap(SST,0)

Your lat/lon variables will be 16 x 200 x 360, but assuming that the
lat/lon doesn't change on the 16 files, then you can simply use a subset
for plotting:

SSTmean at lat2d = lat2d(0,:,:)
SSTmean at lon2d = lon2d(0,:,:)

As for your second question about finding the index of a region for a
curvilinear grid, please see the lat/lon subsetting examples page, and in
particular, examples latlon_subset_2.ncl latlon_subset_3.ncl and
latlon_subset_compare_3.ncl at:

http://www.ncl.ucar.edu/Applications/latlon_subset.shtml

--Mary

On Tue, Jan 23, 2018 at 2:17 AM, sjd.marjani via ncl-talk <ncl-talk at ucar.edu
> wrote:

> On 2018-01-03 06:14, Alan Brammer wrote:
>
> anomexel=exel
> anomel=el - meanClimatology
> anomexla=exla-meanClimatology
>
> These lines do not copy the metadata.  You should do a
> copy_VarCoords(meanClimatology, anomel)  for each respective variable.
>
> As such, you will have seen a warning about a lack of coordinates when
> plotting. NCL will have tried to guess at the location data probably 0:360
> but your data is likely -180:180 thus the resulting plot is 180° off.
>
> Pay attention to any warnings that NCL produces, it's pretty easy to have
> scripts run without any undesired print outs (don't just ignore them).
>
> Also printVarSummary(anomel) will show a complete lack of any metadata.
> Always a good starting point when something isn't quite right.
>
> Good luck,
>
> Alan
>
>
> On Wed, Jan 3, 2018 at 4:32 AM, sjd.marjani via ncl-talk <
> ncl-talk at ucar.edu> wrote:
>
>> Hi and HAPPY new YEAR .
>>
>> my name is Sajede
>>
>> I Need your help
>>
>> I send you my script in Word Document
>>
>> in that ,i highlightet 2 parts
>>
>> => part1 is "YELLOW" which contains my problem
>>
>> =>part2 is "GREEN" which has no problem in ploting
>>
>> and also i send you the related plots in png format
>>
>> **************************************************
>>
>> my problem is why the PART.png has PHASE differences?
>>
>> Thank YOU in advance
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>> Hello again
>
> I have use your advise and its help me so much ,but I have 3 important
> question and all of them are about curvilinear Oceanic Data.
>
> I really need HElp.
>
> ============================================
>
> first Question:
>
> these are my data files
>
> tos_Omon_GFDL-CM3_historical_r1i1p1_192501-192912.nc
>
> tos_Omon_GFDL-CM3_historical_r1i1p1_193001-193412.nc
>
> tos_Omon_GFDL-CM3_historical_r1i1p1_193501-193912.nc
>
> tos_Omon_GFDL-CM3_historical_r1i1p1_194001-194412.nc
>
> .
>
> .
>
> .
>
> tos_Omon_GFDL-CM3_historical_r1i1p1_199501-199912.nc
>
> tos_Omon_GFDL-CM3_historical_r1i1p1_200001-200412.nc
>
> *****************************************************
>
> and it is the attribute of one of these files :
>
> netcdf tos_Omon_GFDL-CM3_historical_r1i1p1_192501-192912 {
> dimensions:
> time = UNLIMITED ; // (60 currently)
> bnds = 2 ;
> rlat = 200 ;
> rlon = 360 ;
> vertices = 4 ;
> variables:
> double average_DT(time) ;
> average_DT:long_name = "Length of average period" ;
> average_DT:units = "days" ;
> double average_T1(time) ;
> average_T1:long_name = "Start time for average period" ;
> average_T1:units = "days since 1860-01-01 00:00:00" ;
> double average_T2(time) ;
> average_T2:long_name = "End time for average period" ;
> average_T2:units = "days since 1860-01-01 00:00:00" ;
> double bnds(bnds) ;
> bnds:long_name = "vertex number" ;
> bnds:cartesian_axis = "N" ;
> float tos(time, rlat, rlon) ;
> tos:missing_value = 1.e+20f ;
> tos:original_name = "sst + 273.15" ;
> tos:cell_methods = "time: mean" ;
> tos:valid_range = 263.f, 373.f ;
> tos:units = "K" ;
> tos:long_name = "Sea Surface Temperature" ;
> tos:coordinates = "lat lon" ;
> tos:_FillValue = 1.e+20f ;
> tos:standard_name = "surface_temperature" ;
> tos:cell_measures = "area: areacello" ;
> tos:comments = "Note: This model output is presented on the model\'s
> tripolar grid. The North Pole singularity is avoided by using this
> nonregular grid north of 65N. More information about the ocean and sea ice
> model\'s grid can be found at http://nomads.gfdl.noaa.gov/." ;
> tos:associated_files = "baseURL: http://cmip-pcmdi.llnl.gov/
> CMIP5/dataLocation gridspecFile: gridspec_ocean_fx_GFDL-CM3_
> historical_r0i0p0.nc areacello: areacello_fx_GFDL-CM3_historical_r0i0p0.nc"
> ;
> tos:original_units = "K" ;
> double time(time) ;
> time:long_name = "time" ;
> time:units = "days since 1860-01-01 00:00:00" ;
> time:cartesian_axis = "T" ;
> time:calendar_type = "noleap" ;
> time:calendar = "noleap" ;
> time:bounds = "time_bnds" ;
> time:standard_name = "time" ;
> time:axis = "T" ;
> double time_bnds(time, bnds) ;
> time_bnds:long_name = "time axis boundaries" ;
> time_bnds:units = "days since 1860-01-01 00:00:00" ;
> double rlon(rlon) ;
> rlon:long_name = "longitude in rotated pole grid" ;
> rlon:units = "degrees" ;
> rlon:axis = "X" ;
> rlon:standard_name = "grid_longitude" ;
> double rlat(rlat) ;
> rlat:long_name = "latitude in rotated pole grid" ;
> rlat:units = "degrees" ;
> rlat:axis = "Y" ;
> rlat:standard_name = "grid_latitude" ;
> float lat(rlat, rlon) ;
> lat:long_name = "latitude coordinate" ;
> lat:units = "degrees_north" ;
> lat:standard_name = "latitude" ;
> lat:bounds = "lat_vertices" ;
> double lat_vertices(rlat, rlon, vertices) ;
> lat_vertices:units = "degrees_north" ;
> float lon(rlat, rlon) ;
> lon:long_name = "longitude coordinate" ;
> lon:units = "degrees_east" ;
> lon:standard_name = "longitude" ;
> lon:bounds = "lon_vertices" ;
> double lon_vertices(rlat, rlon, vertices) ;
> lon_vertices:units = "degrees_east" ;
> int vertices(vertices) ;
>
> // global attributes:
> :title = "NOAA GFDL GFDL-CM3, historical (run 1) experiment output for
> CMIP5 AR5" ;
> :institute_id = "NOAA GFDL" ;
> :source = "GFDL-CM3 2010 atmosphere: AM3 (AM3p9,C48L48); sea ice: SIS
> (SISp2,Tripolar360x200); land: LM3 (LM3p7_cCM3,C48); ocean: MOM4
> (MOM4p1_x1_Z50_cCM3,Tripolar360x200L50)" ;
> :contact = "gfdl.climate.model.info at noaa.gov" ;
> :project_id = "CMIP5" ;
> :table_id = "Table Omon (31 Jan 2011)" ;
> :experiment_id = "historical" ;
> :realization = 1 ;
> :modeling_realm = "ocean" ;
> :tracking_id = "b062901f-bdec-46e5-ba81-f6dd63777987" ;
> :Conventions = "CF-1.4" ;
> :references = "The GFDL Data Portal (http://nomads.gfdl.noaa.gov/)
> provides access to NOAA/GFDL\'s publicly available model input and output
> data sets. From this web site one can view and download data sets and
> documentation, including those related to the GFDL coupled models
> experiments run for the IPCC\'s 5th Assessment Report and the US CCSP." ;
> :comment = "GFDL experiment name: CM3Z_D1_1860-2005_AllForc_H1\n",
> "GFDL parent experiment name: CM3Z_Control-1860_D1\n",
> "CMIP5 experiment ID: historical\n",
> "Initial conditions from this experiment (CMIP5 expt. 3.2 Core) were taken
> from 01Jan0001 of the 1860 control model experiment
> CM3Z_Control-1860_D1.\n",
> "Forcing agents (7 WMGGs, tropospheric and stratospheric ozone, land use,
> volcanoes, solar, sulfate, black carbon,\n",
> "organic carbon, dust, and sea salt) vary historically from 1860-2005
> following prescribed emission or concentration\n",
> "timeseries or internally calculated emissions. This model includes direct
> and indirect aerosol effects." ;
> :gfdl_experiment_name = "CM3Z_D1_1860-2005_AllForc_H1" ;
> :creation_date = "2011-08-22T06:05:17Z" ;
> :model_id = "GFDL-CM3" ;
> :branch_time = "0.0" ;
> :experiment = "historical" ;
> :forcing = "GHG,SA,Oz,LU,Sl,Vl,SS,BC,MD,OC (GHG includes CO2, CH4, N2O,
> CFC11, CFC12, HCFC22, CFC113)" ;
> :frequency = "mon" ;
> :initialization_method = 1 ;
> :parent_experiment_id = "piControl" ;
> :physics_version = 1 ;
> :product = "output1" ;
> :institution = "NOAA GFDL(201 Forrestal Rd, Princeton, NJ, 08540)" ;
> :history = "File was processed by fremetar (GFDL analog of CMOR).
> TripleID: [exper_id_1XhP9uft45,realiz_id_4e0zLckGg3,run_id_HhaU8NYVBh]" ;
> :parent_experiment_rip = "r1i1p1" ;
>
> ************************************************************
> ***********************************************************
>
> when I read just one of these files ,my script can plot the Sea Surface
> Temperature
>
> but
>
> when I read all the 16 files, the ncl plot it with "phase difference".
>
> the result of printVarSummary(lat2d) for the single file and multiple
> files ,show that the ncl multiple the number of data files by rlat
>
>   3200=16*200
>
> =======================
>
> Variable: lat2d
> Type: float
> Total Size: 288000 bytes
> 72000 values
> Number of Dimensions: 2
> Dimensions and sizes: [rlat | 200] x [rlon | 360]
> Coordinates:
> rlat: [-81.5..89.50000080000001]
> rlon: [-279.5..79.5]
> Number Of Attributes: 4
> long_name : latitude coordinate
> units : degrees_north
> standard_name : latitude
> bounds : lat_vertices
>
> ====================
>
> Variable: lat2dH
> Type: float
> Total Size: 4608000 bytes
> 1152000 values
> Number of Dimensions: 2
> Dimensions and sizes: [rlat | 3200] x [rlon | 360]
> Coordinates:
> rlat: [-81.5..89.50000080000001]
> rlon: [-279.5..79.5]
> Number Of Attributes: 4
> long_name : latitude coordinate
> units : degrees_north
> standard_name : latitude
> bounds : lat_vertices
>
> ==============
>
> my script is in below:
>
> ==============================================>>>
>
> ;***********************************************
> begin
> ;************************************************
>
> diri = "/home1/omid/Marjani/DATA/CMIP5/GFDL-CM3_curv2/hisGFDL-CM3/" ;
> input directory
> fili = "tos_Omon_GFDL-CM3_historical_r1i1p1_192501-192912.nc" ; your file
> tos(1872*300*360)
> pthi = diri+fili
> f = addfile( pthi, "r")
> SST = f->tos
> lat2d = f->lat
> lon2d =f->lon
>
>
> ; filess = systemfunc ("ls /home1/omid/Marjani/DATA/
> CMIP5/GFDL-CM3_curv2/hisGFDL-CM3/tos_Omon_GFDL-CM3_historical_r1i1p1*.nc")
> ;H = addfiles(filess,"r")
> ;ListSetType (H, "cat") ; concatenate (=default)
> ;SSTH= H[:]->tos
> ; lat2dH = H[:]->lat
> ; lon2dH=H[:]->lon
>
> ;***********************************************************
> **********************************
>
> ; printVarSummary(SST)
> ;printVarSummary(SSTH)
>
> ;printVarSummary(lat2d)
> ;printVarSummary(lat2dH)
>
> ;===========================================================
> ==================
> ;   First Part=> PLOT SST
> ;===========================================================
> ====================
> wks = gsn_open_wks("x11","marjani")
> gsn_define_colormap(wks,"BlAqGrYeOrReVi200") ; choose colormap
> res = True
> res at cnFillOn = True ; turn on color
> res at cnFillMode = "RasterFill" ; turn on raster mode ... optional
> res at cnLinesOn = False ; turn off contour lines
> res at cnLineLabelsOn = False
>
> res at cnLineLabelsOn = False
> res at cnFillPalette = "BlAqGrYeOrReVi200" ; your choice
> res at gsnAddCyclic = True ;-- lon < 360 degrees
> res at tiMainString = "SST_I" ; title
>
> SST at lat2d = lat2d
> SST at lon2d = lon2d
>
> ; ******************************************************
> res at mpMinLonF = 0
> res at mpMaxLonF = 360
> res at mpCenterLonF = 180 ; This is necessary to get the correct ma
>
>
> ;===========================================================
> ===============
> plot = gsn_csm_contour_map_ce(wks,SSTH(4,:,:),res) ; Can not plot it by
> the use of systemfunc and just can do it when i use  1data file
>
> ;.................................................................
> ;SSTmean=dim_avg_n_Wrap(SST(0:59,:,:),0)
> ;copy_VarCoords(SST, SSTmean)
> ;plot = gsn_csm_contour_map_ce(wks,SSTmean,res)
> ;......................................................
> ;s=SST(10,:,:)
> ;anom = s - SSTmean
> ;copy_VarCoords(SSTmean, anom)
> ;copy_VarAtts(SSTmean, anom)
> ;plot = gsn_csm_contour_map_ce(wks,anom,res) ;nemitoone
> ;.............................................................
> ;anomSST = dim_rmvmean_n_Wrap(SST(11:611:12,:,:),0)
> ;plot = gsn_csm_contour_map_ce(wks,anomSST(1,:,:),res)
>
>
> ;===========================================================
> =============================
>
> ;       Second Part
>
> ;========================================================================
>
> ....................................................................
> ...................................................................
> ;print("--------print(lat2d(:,160))----------------")
> ;print(lat2d(:,160))
> ;print("--------print(lon2d(:,160))----------------")
> ;print(lon2d(:,160))
>
> ;print("--------print(lat2d(:,190)----------------")
> ;print(lat2d(:,190))
> ;print("--------print(lon2d(:,190)----------------")
>
> ;==============================
>
> ;print("--------print(lat2d(:,64))----------------")
> ;print(lat2d(:,64))
> ;print("--------print(lon2d(:,64))----------------")
> ;print(lon2d(:,64))
>
> ;print("--------print(lat2d(:,87))----------------")
> ;print(lat2d(:,87))
> ;print("--------print(lon2d(:,87))----------------")
> ;print(lon2d(:,87))
> ;==============================================================
> end
>
>
> ======================
>
> I have 2 other question,
>
> the second question is a bout the Second part of my script.
>
> the second part is a way that i use to find for example Nino3.4 region .
>
> would you please tell me the right way to read the lat and lon of
> curvilinear grid,,and is there any method to find an exact region?
>
> ;=============================
>
> I will ask the third QUestion later.
>
> I wait for you Help .
>
> THank YOu in ADvance
>
>
>    Sajede
>
>
>
>
>
>
>
>
> _______________________________________________
> 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/20180123/65edc3af/attachment.html>