(no subject)

From: mahmoud elkarim <mgaberelkarim_at_nyahnyahspammersnyahnyah>
Date: Fri Feb 01 2013 - 14:39:59 MST

dear all
find attached ncl scripits
I'm beginner user in ncl scripts
when i was trying to find precipitation b**ias which is the difference
between regcm output and *GPCP observed data i found this error
(fatal:Minus: Dimension size, for dimension number 0, of operands does not
match, can't continue)
and i don't know how i can regrid gpcp to regcm

rgcm file dump is
     x = 129 ;
    y = 127 ;
    m10 = 1 ;
    m2 = 1 ;
    lev = 2 ;
    time = UNLIMITED ; // (1 currently)
    nb2 = 2 ;
variables:
    float xlon(y, x) ;
        xlon:standard_name = "longitude" ;
        xlon:long_name = "Longitude at cross points" ;
        xlon:units = "degrees_east" ;
        xlon:_CoordinateAxisType = "Lon" ;
    float xlat(y, x) ;
        xlat:standard_name = "latitude" ; smw:cell_methods = "time:
point" ;
   float tpr (time, y, x) ;
        tpr:standard_name = "precipitation_flux" ;
        tpr:long_name = "Total precipitation" ;
        tpr:units = "kg m-2 day-1" ;
        tpr:coordinates = "xlon xlat" ;
        tpr:cell_methods = "time: mean (interval: 3 hours)" ;

        xlat:long_name = "Latitude at cross points" ;
        xlat:units = "degrees_north" ;**
        xlat:_CoordinateAxisType = "Lat" ;netcdf gpccDJFmeanbox {

*
*and gpgp dump is
dimensions:
    lon = 130 ;
    lat = 111 ;
    time = UNLIMITED ; // (1 currently)
variables:
    double lon(lon) ;
        lon:standard_name = "longitude" ;
        lon:long_name = "longitude" ;
        lon:units = "degrees_east" ;
        lon:axis = "X" ;
    double lat(lat) ;
        lat:standard_name = "latitude" ;
        lat:long_name = "latitude" ;
        lat:units = "degrees_north" ;
        lat:axis = "Y" ;
    double time(time) ;
        time:standard_name = "time" ;
        time:units = "days since 1901-01-01 00:00:00" ;
        time:calendar = "proleptic_gregorian" ;
    float p(time, lat, lon) ;
        p:long_name = "full data precipitation version6" ;
        p:units = "mm/month" ;
        p:code = 16 ;
        p:_FillValue = -99999.99f ;

i ask if any one could help me
thanks for attention

*
*> Mahmoud Gaber*
*> Nile Forecasting center (NFC) *
*> minister of water resources and irrigation (MWRI) * * *

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk

Received on Fri Feb 1 14:40:08 2013

This archive was generated by hypermail 2.1.8 : Fri Feb 01 2013 - 15:31:04 MST