[ncl-talk] computing daily anomalies from a high resolution data
Dennis Shea
shea at ucar.edu
Thu Apr 21 07:56:50 MDT 2022
Obviously, you have exceeded available memory.
You will have to calculate 'geographical subsets' of the anomalies. Then
You must always read all the time steps.
ilat1 =
ilat2 =
ilon1 =
ilon2 =
CHLIO = data1->CHLIO(:,ilat1:ilat2,ilon1:1lon2) ; all times
printVarSummary(CHLIO)
===
See Example of function *calcDayAnomTLL*
<https://www.ncl.ucar.edu/Document/Functions/Contributed/calcDayAnomTLL.shtml>
Also, unless you are using a very old version of NCL, you do not need to
explicitly 'load' the libraries. The following can be eliminated:
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
On Thu, Apr 21, 2022 at 5:18 AM saurabh rathore via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:
> Hello NCL users
> G'day
>
> I am trying to compute the daily anomalies from a high-resolution daily
> data as shown below from the ncdump. I am just showing the dimension of the
> variable CHLIO.
>
> dimensions:
> LON161_408 = 248 ;
> LAT361_472 = 112 ;
> TIME = UNLIMITED ; // (8362 currently)
> bnds = 2 ;
> variables:
> double CHLIO(TIME, LAT361_472, LON161_408) ;
> CHLIO:missing_value = -1.e+34 ;
> CHLIO:_FillValue = -1.e+34 ;
>
> The ncl code that I am using is below but it gets killed as the memory
> gets fulled during the opration.
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>
>
> data1 = addfile("/media/srathore/Binny5/cmes-oras/chl/chlio-1998-2020.nc
> ","r")
> print(data1)
>
> ;***********************************************************
> ; Read time and create required yyyyddd
> ;***********************************************************
> time = data1->TIME ; time:units = "hours since
> 1-1-1 00:00:0.0"
> TIME = cd_calendar(time, 0) ; type float
> year = toint( TIME(:,0) )
> month = toint( TIME(:,1) )
> day = toint( TIME(:,2) )
> ddd = day_of_year(year, month, day)
> yyyyddd = year*1000 + ddd ; needed for input
> print(yyyyddd)
>
> ;***********************************************************
> ; Read data: short2flt
> ;***********************************************************
> CHLIO = data1->CHLIO
> printVarSummary( chl )
>
> So how to compute the daily anomalies in this case?
>
> Cheers, Saurabh
>
> --
>
>
> *REGARDS*
>
> *Saurabh Rathore*
> *Research Scholar (PhD.)*
> *Centre For Oceans, Rivers, Atmosphere & Land Science Technology*
> *Indian Institute Of Technology, Kharagpur*
> *contact :- 91- 8345984434*
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220421/b3b4dcd8/attachment.html>
More information about the ncl-talk
mailing list