[gradsusr] converting rotate pole grid
Muhammad Rahiz
muhammad.rahiz at ouce.ox.ac.uk
Sat Oct 16 04:32:38 EDT 2010
Hi Robert,
I have the same problem and got it resolved. Regrid the data by running
the following CDO commands
1 - make a mygrid file containing the following (change values as
per your data). you can obtain most from the following command;
cdo griddes nhqhga.pdg6dec.pnw.nc
Notes:
xnpole/ynpole = position of rotated north pole (obtainable from model
output / technical manual)
pay attention to "-" sign because it matters
#
# gridID 0
#
gridtype = lonlat
xsize = 272
ysize = 214
xfirst = -35.14
xinc = 0.22
yfirst = -23.54
yinc = -0.22
xnpole = -162.
ynpole = 39.25
2 - add the rotated north pole position to file
cdo regrid,mygrid nhqhga.pdg6dec.pnw.nc out.nc
3 - regrid to 0.22deg resolution world map
cdo remapbil,r1800x600 out.nc new.nc
4- open grads and select your domain.
It should work.
HTH,
Muhammad Rahiz
On Fri, 15 Oct 2010, Roberto Mera wrote:
> Hi, I've been trying to tackle this problem for a while now. This
> particular file is supposed to be over the Pacific Northwest. The file
> plots over Africa and I need to shift it so that it maps correctly. Here
> is the ncdump -c of the file:
> [rjmera at climlab02 CPDN]$ ncdump -c prev/nhqhga.pdg6dec.pnw.nc
> netcdf nhqhga.pdg6dec.pnw {
> dimensions:
> X = 110 ;
> Y = 117 ;
> Z = 1 ;
> variables:
> int time(Z) ;
> time:long_name = "time" ;
> time:units = "month" ;
> float latitude(Y) ;
> latitude:long_name = "latitude" ;
> latitude:units = "degrees_north" ;
> float longitude(X) ;
> longitude:long_name = "longitude" ;
> longitude:units = "degrees_east" ;
> float Tmin(Z, Y, X) ;
> Tmin:long_name = "Monthly average of daily Tmin" ;
> Tmin:units = "C" ;
> Tmin:Fill_Value = -9999.f ;
> float Frost(Z, Y, X) ;
> Frost:long_name = "Number of frost days per month" ;
> Frost:units = "days" ;
> Frost:Fill_Value = -9999.f ;
> float Tmin_1(Z, Y, X) ;
> Tmin_1:long_name = "First monthly minimum of Tmin" ;
> Tmin_1:units = "C" ;
> Tmin_1:Fill_Value = -9999.f ;
> float Tmin_2(Z, Y, X) ;
> Tmin_2:long_name = "Second montlhy minimum of Tmin" ;
> Tmin_2:units = "C" ;
> Tmin_2:Fill_Value = -9999.f ;
> float Tmin_3(Z, Y, X) ;
> Tmin_3:long_name = "Third monthly minimum of Tmin" ;
> Tmin_3:units = "C" ;
> Tmin_3:Fill_Value = -9999.f ;
> float Tmax(Z, Y, X) ;
> Tmax:long_name = "Monthly average of daily Tmax" ;
> Tmax:units = "C" ;
> Tmax:Fill_Value = -9999.f ;
> float TT30(Z, Y, X) ;
> TT30:long_name = "Number of days with Tmax > 30 degC" ;
> TT30:units = "days" ;
> TT30:Fill_Value = -9999.f ;
> float TT35(Z, Y, X) ;
> TT35:long_name = "Number of days with Tmax > 35 degC" ;
> TT35:units = "days" ;
> TT35:Fill_Value = -9999.f ;
> float TT40(Z, Y, X) ;
> TT40:long_name = "Number of days with Tmax > 40 degC" ;
> TT40:units = "days" ;
> TT40:Fill_Value = -9999.f ;
> float Tmax_1(Z, Y, X) ;
> Tmax_1:long_name = "First monthly maximum of Tmax" ;
> Tmax_1:units = "C" ;
> Tmax_1:Fill_Value = -9999.f ;
> float Tmax_2(Z, Y, X) ;
> Tmax_2:long_name = "Second monthly maximum of Tmax" ;
> Tmax_2:units = "C" ;
> Tmax_2:Fill_Value = -9999.f ;
> float Tmax_3(Z, Y, X) ;
> Tmax_3:long_name = "Third monthly maximum of Tmax" ;
> Tmax_3:units = "C" ;
> Tmax_3:Fill_Value = -9999.f ;
> float U-std(Z, Y, X) ;
> U-std:long_name = "Monthly standard deviation of U surface wind" ;
> U-std:units = "m s-1" ;
> U-std:Fill_Value = -9999.f ;
> float V-std(Z, Y, X) ;
> V-std:long_name = "Monthly standard deviation of V surface wind" ;
> V-std:units = "m s-1" ;
> V-std:Fill_Value = -9999.f ;
> float Max_wind_speed(Z, Y, X) ;
> Max_wind_speed:long_name = "Monthly maximum of daily wind speed" ;
> Max_wind_speed:units = "m s-1" ;
> Max_wind_speed:Fill_Value = -9999.f ;
> float Dry_max_1(Z, Y, X) ;
> Dry_max_1:long_name = "Max number of consecutive dry days at
> start of the month" ;
> Dry_max_1:units = "days" ;
> Dry_max_1:Fill_Value = -9999.f ;
> float Dry_max_2(Z, Y, X) ;
> Dry_max_2:long_name = "Max number of consecutive dry days at
> end of the month" ;
> Dry_max_2:units = "days" ;
> Dry_max_2:Fill_Value = -9999.f ;
> float Dry_max_3(Z, Y, X) ;
> Dry_max_3:long_name = "Max number of consecutive dry days
> during the month" ;
> Dry_max_3:units = "days" ;
> Dry_max_3:Fill_Value = -9999.f ;
> float Wet_days(Z, Y, X) ;
> Wet_days:long_name = "Number of wet days during the month" ;
> Wet_days:units = "days" ;
> Wet_days:Fill_Value = -9999.f ;
> float Precip_1(Z, Y, X) ;
> Precip_1:long_name = "Total precipitation within the 1st day of
> the month" ;
> Precip_1:units = "mm" ;
> Precip_1:Fill_Value = -9999.f ;
> float Precip_2(Z, Y, X) ;
> Precip_2:long_name = "Total precipitation on the 21st days of
> the month" ;
> Precip_2:units = "mm" ;
> Precip_2:Fill_Value = -9999.f ;
> float Precip_3(Z, Y, X) ;
> Precip_3:long_name = "Total precipitation within 3 days" ;
> Precip_3:units = "mm" ;
> Precip_3:Fill_Value = -9999.f ;
> float Precip_4(Z, Y, X) ;
> Precip_4:long_name = "Total precipitation within the 2 last
> days of the month" ;
> Precip_4:units = "mm" ;
> Precip_4:Fill_Value = -9999.f ;
> float Precip_5(Z, Y, X) ;
> Precip_5:long_name = "Total precipitation within the last day
> of the month" ;
> Precip_5:units = "mm" ;
> Precip_5:Fill_Value = -9999.f ;
> float Pge40(Z, Y, X) ;
> Pge40:long_name = "Number of days with daily P > 40 mm" ;
> Pge40:units = "days" ;
> Pge40:Fill_Value = -9999.f ;
> float Pge80(Z, Y, X) ;
> Pge80:long_name = "Number of days with daily P > 80 mm" ;
> Pge80:units = "days" ;
> Pge80:Fill_Value = -9999.f ;
> float Pge120(Z, Y, X) ;
> Pge120:long_name = "Number of days with daily P > 120 mm" ;
> Pge120:units = "days" ;
> Pge120:Fill_Value = -9999.f ;
> float Pmax_1(Z, Y, X) ;
> Pmax_1:long_name = "First max daily total precipitation rate
> over the month" ;
> Pmax_1:units = "mm" ;
> Pmax_1:Fill_Value = -9999.f ;
> float Pmax_2(Z, Y, X) ;
> Pmax_2:long_name = "Second max daily total precipitation rate
> over the month" ;
> Pmax_2:units = "mm" ;
> Pmax_2:Fill_Value = -9999.f ;
> float Pmax_3(Z, Y, X) ;
> Pmax_3:long_name = "Third max daily total precipitation rate
> over the month" ;
> Pmax_3:units = "mm" ;
> Pmax_3:Fill_Value = -9999.f ;
> float Z(Z, Y, X) ;
> Z:long_name = "Monthly mean geopotential height at 500mb" ;
> Z:units = "m" ;
> Z:Fill_Value = -9999.f ;
> float Z-std(Z, Y, X) ;
> Z-std:long_name = "Monthly standard deviation of geopotential
> height at 500mb" ;
> Z-std:units = "m" ;
> Z-std:Fill_Value = -9999.f ;
> float Snow_mass(Z, Y, X) ;
> Snow_mass:long_name = "Snow mass on the first day of the month" ;
> Snow_mass:units = "kg m-2" ;
> Snow_mass:Fill_Value = -9999.f ;
> float DlastSM(Z, Y, X) ;
> DlastSM:long_name = "Date of last snow mass" ;
> DlastSM:units = "day" ;
> DlastSM:Fill_Value = -9999.f ;
> float Max_runoff(Z, Y, X) ;
> Max_runoff:long_name = "Maximum daily runoff" ;
> Max_runoff:units = "kg m-2 s-1" ;
> Max_runoff:Fill_Value = -9999.f ;
> float Date_max_runoff(Z, Y, X) ;
> Date_max_runoff:long_name = "Date of maximum daily runoff" ;
> Date_max_runoff:units = "day" ;
> Date_max_runoff:Fill_Value = -9999.f ;
> data:
> }
> ncdump -v latitude gives:
> latitude = 13.31, 13.09, 12.87, 12.65, 12.43, 12.21, 11.99, 11.77, 11.55,
> 11.33, 11.11, 10.89, 10.67, 10.45, 10.23, 10.01, 9.79, 9.57, 9.349999,
> 9.13, 8.91, 8.69, 8.469999, 8.25, 8.03, 7.81, 7.59, 7.37, 7.15, 6.93,
> 6.71, 6.49, 6.27, 6.05, 5.83, 5.61, 5.39, 5.17, 4.95, 4.73, 4.51,
> 4.29,
> 4.07, 3.85, 3.63, 3.41, 3.19, 2.97, 2.75, 2.53, 2.31, 2.09, 1.87,
> 1.65,
> 1.43, 1.21, 0.9899998, 0.7699998, 0.5499998, 0.3299998, 0.1099998,
> -0.1100002, -0.3300002, -0.5500002, -0.7700002, -0.9900002, -1.21,
> -1.43,
> -1.65, -1.87, -2.09, -2.31, -2.53, -2.75, -2.97, -3.19, -3.41, -3.63,
> -3.85, -4.07, -4.29, -4.51, -4.73, -4.95, -5.17, -5.39, -5.61, -5.83,
> -6.05, -6.27, -6.49, -6.71, -6.93, -7.15, -7.37, -7.59, -7.81, -8.03,
> -8.25, -8.47, -8.690001, -8.91, -9.13, -9.35, -9.57, -9.79, -10.01,
> -10.23, -10.45, -10.67, -10.89, -11.11, -11.33, -11.55, -11.77,
> -11.99, 0 ;
> }
> and longitude:
>
> longitude = 347.79, 348.01, 348.23, 348.45, 348.67, 348.89, 349.11,
> 349.33,
> 349.55, 349.77, 349.99, 350.21, 350.43, 350.65, 350.87, 351.09,
> 351.31,
> 351.53, 351.75, 351.97, 352.19, 352.41, 352.63, 352.85, 353.07,
> 353.29,
> 353.51, 353.73, 353.95, 354.17, 354.39, 354.61, 354.83, 355.05,
> 355.27,
> 355.49, 355.71, 355.93, 356.15, 356.37, 356.59, 356.81, 357.03,
> 357.25,
> 357.47, 357.69, 357.91, 358.13, 358.35, 358.57, 358.79, 359.01,
> 359.23,
> 359.45, 359.67, 359.89, 360.11, 360.33, 360.55, 360.77, 360.99,
> 361.21,
> 361.43, 361.65, 361.87, 362.09, 362.31, 362.53, 362.75, 362.97,
> 363.19,
> 363.41, 363.63, 363.85, 364.07, 364.29, 364.51, 364.73, 364.95,
> 365.17,
> 365.39, 365.61, 365.83, 366.05, 366.27, 366.49, 366.71, 366.93,
> 367.15,
> 367.37, 367.59, 367.81, 368.03, 368.25, 368.47, 368.69, 368.91,
> 369.13,
> 369.35, 369.57, 369.79, 370.01, 370.23, 370.45, 370.67, 370.89,
> 371.11,
> 371.33, 371.55, 371.77 ;
> }
>
> I set up a ctl file that reads the variables fine but plots it over Africa.
> DSET ^nhqhga.pdg6dec.pnw.nc
> dtype netcdf
> TITLE Pacific northwest
> UNDEF -99999.f
> XDEF 110 LEVELS
> 347.79
> 348.01
> 348.23
> 348.45
> 348.67
> 348.89
> 349.11
> 349.33
> 349.55
> 349.77
> 349.99
> 350.21
> 350.43
> 350.65
> 350.87
> 351.09
> 351.31
> 351.53
> 351.75
> 351.97
> 352.19
> 352.41
> 352.63
> 352.85
> 353.07
> 353.29
> 353.51
> 353.73
> 353.95
> 354.17
> 354.39
> 354.61
> 354.83
> 355.05
> 355.27
> 355.49
> 355.71
> 355.93
> 356.15
> 356.37
> 356.59
> 356.81
> 357.03
> 357.25
> 357.47
> 357.69
> 357.91
> 358.13
> 358.35
> 358.57
> 358.79
> 359.01
> 359.23
> 359.45
> 359.67
> 359.89
> 360.11
> 360.33
> 360.55
> 360.77
> 360.99
> 361.21
> 361.43
> 361.65
> 361.87
> 362.09
> 362.31
> 362.53
> 362.75
> 362.97
> 363.19
> 363.41
> 363.63
> 363.85
> 364.07
> 364.29
> 364.51
> 364.73
> 364.95
> 365.17
> 365.39
> 365.61
> 365.83
> 366.05
> 366.27
> 366.49
> 366.71
> 366.93
> 367.15
> 367.37
> 367.59
> 367.81
> 368.03
> 368.25
> 368.47
> 368.69
> 368.91
> 369.13
> 369.35
> 369.57
> 369.79
> 370.01
> 370.23
> 370.45
> 370.67
> 370.89
> 371.11
> 371.33
> 371.55
> 371.77
> YDEF 117 LEVELS
> -11.99
> -11.77
> -11.55
> -11.33
> -11.11
> -10.89
> -10.67
> -10.45
> -10.23
> -10.01
> -9.79
> -9.57
> -9.35
> -9.13
> -8.91
> -8.69
> -8.47
> -8.25
> -8.03
> -7.81
> -7.59
> -7.37
> -7.15
> -6.93
> -6.71
> -6.49
> -6.27
> -6.05
> -5.83
> -5.61
> -5.39
> -5.17
> -4.95
> -4.73
> -4.51
> -4.29
> -4.07
> -3.85
> -3.63
> -3.41
> -3.19
> -2.97
> -2.75
> -2.53
> -2.31
> -2.09
> -1.87
> -1.65
> -1.43
> -1.21
> -0.99
> -0.77
> -0.55
> -0.33
> -0.11
> 0.00
> 0.11
> 0.33
> 0.55
> 0.77
> 0.99
> 1.21
> 1.43
> 1.65
> 1.87
> 2.09
> 2.31
> 2.53
> 2.75
> 2.97
> 3.19
> 3.41
> 3.63
> 3.85
> 4.07
> 4.29
> 4.51
> 4.73
> 4.95
> 5.17
> 5.39
> 5.61
> 5.83
> 6.05
> 6.27
> 6.49
> 6.71
> 6.93
> 7.15
> 7.37
> 7.59
> 7.81
> 8.03
> 8.25
> 8.47
> 8.69
> 8.91
> 9.13
> 9.35
> 9.57
> 9.79
> 10.01
> 10.23
> 10.45
> 10.67
> 10.89
> 11.11
> 11.33
> 11.55
> 11.77
> 11.99
> 12.21
> 12.43
> 12.65
> 12.87
> 13.09
> 13.31
> ZDEF 1 LINEAR 1 1
> TDEF 1 LINEAR 0Z06DEC1967 1MO
> VARS 36
> Tmin=>Tmin 1 Z,Y,X Tmin (C)
> Frost=>frost 1 Z,Y,X Frost (days)
> Tmin_1=>tmin1 1 Z,Y,X Min1 (C)
> Tmin_2=>tmin2 1 Z,Y,X Min2 (C)
> Tmin_3=>tmin3 1 Z,Y,X Min3 (C)
> Tmax=>tmax 1 Z,Y,X Tmax (C)
> TT30=>tt30 1 Z,Y,X 30Ctmax (days)
> TT35=>tt35 1 Z,Y,X 35Ctmax (days)
> TT40=>tt40 1 Z,Y,X 40Ctmax (days)
> Tmax_1=>tmax1 1 Z,Y,X Max1 (C)
> Tmax_2=>tmax2 1 Z,Y,X Max2 (C)
> Tmax_3=>tmax3 1 Z,Y,X Max3 (C)
> U-std=>ustd 1 Z,Y,X monstdU (m s-1)
> V-std=>vstd 1 Z,Y,X monstdU (m s-1)
> Max_wind_speed=>mxwspeed 1 Z,Y,X maxwindsped (m s-1)
> Dry_max_1=>drymax1 1 Z,Y,X daysmaxdrystart (days)
> Dry_max_2=>drymax2 1 Z,Y,X daysmaxdrymid (days)
> Dry_max_3=>drymax3 1 Z,Y,X daysmaxdryend (days)
> Wet_days=>wetdays 1 Z,Y,X numwetdays (days)
> Precip_1=>p1 1 Z,Y,X 1stdayprec (mm)
> Precip_2=>p2 1 Z,Y,X 21stdayprec (mm)
> Precip_3=>p3 1 Z,Y,X 3daysprec (mm)
> Precip_4=>p4 1 Z,Y,X preclast2days (mm)
> Precip_5=>p5 1 Z,Y,X preclastday (mm)
> Pge40=>pge40 1 Z,Y,X 40mmprec (days)
> Pge80=>pge80 1 Z,Y,X 80mmprec (days)
> Pge120=>pge120 1 Z,Y,X 120mmprec (days)
> Pmax_1=>pmax1 1 Z,Y,X 1stmaxprec (C)
> Pmax_2=>pmax2 1 Z,Y,X 2ndmaxprec (C)
> Pmax_3=>pmax3 1 Z,Y,X 3rdmaxprec (C)
> Z=>z 1 Z,Y,X hgtmean500 (m)
> Z-std=>zstd 1 Z,Y,X hgtstdv (m)
> Snow_mass=>snowmass 1 Z,Y,X snowmass (kg m-2)
> DlastSM=>lastsnow 1 Z,Y,X daylastsnw (day)
> Max_runoff=>mxrunoff 1 Z,Y,X maxrunoff (kg m-2 s-1)
> Date_max_runoff=>datemxrunoff 1 Z,Y,X maxrunoffdate (day)
> ENDVARS
>
> Grads version:
> Starting "/opt/opengrads/Linux/Versions/2.0.a7.oga.3/i686/grads -HC 1 " ...
> Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3
> Copyright (c) 1988-2008 by Brian Doty and the
> Institute for Global Environment and Society (IGES)
> GrADS comes with ABSOLUTELY NO WARRANTY
> See file COPYRIGHT for more information
> Config: v2.0.a7.oga.3 little-endian readline printim grib2 netcdf
> hdf4-sds hdf5 opendap-grids,stn athena geotiff
> Issue 'q config' command for more information.
> Loading User Defined Extensions table
> </opt/opengrads/Linux/Versions/2.0.a7.oga.3/i686/gex/udxt> ... ok.
> Landscape mode? ('n' for portrait):
> GX Package Initialization: Size = 11 8.5
> Command line history in /home/rjmera/.grads.log
>
> A different file shows the projection used. here is the ncdump -c for
> that file:
> ncdump -c v588ga.peg8oct.nc
> netcdf v588ga.peg8oct {
> dimensions:
> time0 = UNLIMITED ; // (6 currently)
> time1 = 6 ;
> z0 = 1 ;
> z1 = 1 ;
> z2 = 1 ;
> z3 = 1 ;
> z4 = 1 ;
> latitude0 = 117 ;
> longitude0 = 110 ;
> latitude1 = 116 ;
> longitude1 = 110 ;
> variables:
> float time0(time0) ;
> time0:standard_name = "time" ;
> time0:axis = "T" ;
> time0:calendar = "360_day" ;
> time0:units = "days since 1959-12-01 00:00:00" ;
> time0:meaning_period = "120 hours" ;
> time0:meaning_sampling_frequency = "3 hours" ;
> float time1(time1) ;
> time1:standard_name = "time" ;
> time1:axis = "T" ;
> time1:calendar = "360_day" ;
> time1:units = "days since 1959-12-01 00:00:00" ;
> time1:meaning_period = "120 hours" ;
> time1:meaning_sampling_frequency = "1 hours" ;
> float z0(z0) ;
> z0:standard_name = "height" ;
> z0:units = "m" ;
> z0:direction = "up" ;
> float z1(z1) ;
> z1:standard_name = "tropopause_altitude" ;
> z1:units = "m" ;
> float z2(z2) ;
> z2:standard_name = "height" ;
> z2:units = "m" ;
> z2:direction = "up" ;
> float z3(z3) ;
> z3:standard_name = "air_pressure" ;
> z3:units = "hPa" ;
> z3:direction = "down" ;
> float z4(z4) ;
> z4:standard_name = "height" ;
> z4:units = "m" ;
> z4:direction = "up" ;
> float latitude0(latitude0) ;
> latitude0:standard_name = "grid_latitude" ;
> latitude0:units = "degrees" ;
> latitude0:axis = "Y" ;
> float longitude0(longitude0) ;
> longitude0:standard_name = "grid_longitude" ;
> longitude0:units = "degrees" ;
> longitude0:axis = "X" ;
> float global_latitude0(latitude0, longitude0) ;
> global_latitude0:standard_name = "latitude" ;
> global_latitude0:units = "degrees_north" ;
> float global_longitude0(latitude0, longitude0) ;
> global_longitude0:standard_name = "longitude" ;
> global_longitude0:units = "degrees_east" ;
> char rotated_pole0 ;
> rotated_pole0:grid_mapping_name = "rotated_latitude_longitude" ;
> rotated_pole0:grid_north_pole_latitude = 49.03f ;
> rotated_pole0:grid_north_pole_longitude = 58.4f ;
> float latitude1(latitude1) ;
> latitude1:standard_name = "grid_latitude" ;
> latitude1:units = "degrees" ;
> latitude1:axis = "Y" ;
> float longitude1(longitude1) ;
> longitude1:standard_name = "grid_longitude" ;
> longitude1:units = "degrees" ;
> longitude1:axis = "X" ;
> float global_latitude1(latitude1, longitude1) ;
> global_latitude1:standard_name = "latitude" ;
> global_latitude1:units = "degrees_north" ;
> float global_longitude1(latitude1, longitude1) ;
> global_longitude1:standard_name = "longitude" ;
> global_longitude1:units = "degrees_east" ;
> char rotated_pole1 ;
> rotated_pole1:grid_mapping_name = "rotated_latitude_longitude" ;
> rotated_pole1:grid_north_pole_latitude = 49.03f ;
> rotated_pole1:grid_north_pole_longitude = 58.4f ;
> float field186(time0, z0, latitude0, longitude0) ;
> field186:_FillValue = -1.073742e+09f ;
> field186:grid_mapping = "rotated_pole0" ;
> field186:coordinates = "global_longitude0 global_latitude0" ;
> field186:stash_item = "201" ;
> field186:stash_section = "1" ;
> field186:field_code = "186" ;
> field186:standard_name = "surface_net_downward_shortwave_flux" ;
> field186:long_name = "NET DOWN SURFACE SW FLUX: SW TS ONLY" ;
> field186:units = "W m-2" ;
> field186:cell_method = "time: mean " ;
> float field203(time0, z0, latitude0, longitude0) ;
> field203:_FillValue = -1.073742e+09f ;
> field203:grid_mapping = "rotated_pole0" ;
> field203:coordinates = "global_longitude0 global_latitude0" ;
> field203:stash_item = "235" ;
> field203:stash_section = "1" ;
> field203:field_code = "203" ;
> field203:standard_name = "surface_downwelling_shortwave_flux" ;
> field203:long_name = "TOTAL DOWNWARD SURFACE SW FLUX" ;
> field203:units = "W m-2" ;
> field203:cell_method = "time: mean " ;
> float field187(time0, z0, latitude0, longitude0) ;
> field187:_FillValue = -1.073742e+09f ;
> field187:grid_mapping = "rotated_pole0" ;
> field187:coordinates = "global_longitude0 global_latitude0" ;
> field187:stash_item = "201" ;
> field187:stash_section = "2" ;
> field187:field_code = "187" ;
> field187:standard_name = "surface_net_downward_longwave_flux" ;
> field187:long_name = "NET DOWN SURFACE LW RAD FLUX" ;
> field187:units = "W m-2" ;
> field187:cell_method = "time: mean " ;
> float field205(time0, z1, latitude0, longitude0) ;
> field205:_FillValue = -1.073742e+09f ;
> field205:grid_mapping = "rotated_pole0" ;
> field205:coordinates = "global_longitude0 global_latitude0" ;
> field205:stash_item = "238" ;
> field205:stash_section = "2" ;
> field205:field_code = "205" ;
> field205:standard_name = "tropopause_downwelling_longwave_flux" ;
> field205:long_name = "TOTAL DOWNWARD LW FLUX AT THE TROP." ;
> field205:units = "W m-2" ;
> field205:cell_method = "time: mean " ;
> float field180(time1, z0, latitude0, longitude0) ;
> field180:_FillValue = -1.073742e+09f ;
> field180:grid_mapping = "rotated_pole0" ;
> field180:coordinates = "global_longitude0 global_latitude0" ;
> field180:stash_item = "234" ;
> field180:stash_section = "3" ;
> field180:field_code = "180" ;
> field180:standard_name = "surface_upward_latent_heat_flux" ;
> field180:long_name = "SURFACE LATENT HEAT FLUX W/M2" ;
> field180:units = "W m-2" ;
> field180:cell_method = "time: mean " ;
> float field88(time1, z2, latitude0, longitude0) ;
> field88:_FillValue = -1.073742e+09f ;
> field88:grid_mapping = "rotated_pole0" ;
> field88:coordinates = "global_longitude0 global_latitude0" ;
> field88:stash_item = "245" ;
> field88:stash_section = "3" ;
> field88:field_code = "88" ;
> field88:standard_name = "relative_humidity" ;
> field88:long_name = "RELATIVE HUMIDITY AT 1.5M" ;
> field88:units = "1" ;
> field88:cell_method = "time: mean " ;
> float field106(time1, z0, latitude0, longitude0) ;
> field106:_FillValue = -1.073742e+09f ;
> field106:grid_mapping = "rotated_pole0" ;
> field106:coordinates = "global_longitude0 global_latitude0" ;
> field106:stash_item = "208" ;
> field106:stash_section = "8" ;
> field106:field_code = "106" ;
> field106:standard_name = "liquid_water_content_of_soil_layer" ;
> field106:long_name = "SOIL MOISTURE CONTENT" ;
> field106:units = "kg m-2" ;
> field106:cell_method = "time: mean " ;
> float field56(time1, z3, latitude1, longitude1) ;
> field56:_FillValue = -1.073742e+09f ;
> field56:grid_mapping = "rotated_pole1" ;
> field56:coordinates = "global_longitude1 global_latitude1" ;
> field56:stash_item = "201" ;
> field56:stash_section = "15" ;
> field56:field_code = "56" ;
> field56:standard_name = "eastward_wind" ;
> field56:long_name = "U COMPNT OF WIND ON PRESSURE LEVELS" ;
> field56:units = "m s-1" ;
> field56:cell_method = "time: mean " ;
> float field57(time1, z3, latitude1, longitude1) ;
> field57:_FillValue = -1.073742e+09f ;
> field57:grid_mapping = "rotated_pole1" ;
> field57:coordinates = "global_longitude1 global_latitude1" ;
> field57:stash_item = "202" ;
> field57:stash_section = "15" ;
> field57:field_code = "57" ;
> field57:standard_name = "northward_wind" ;
> field57:long_name = "V COMPNT OF WIND ON PRESSURE LEVELS" ;
> field57:units = "m s-1" ;
> field57:cell_method = "time: mean " ;
> float field8(time1, z4, latitude0, longitude0) ;
> field8:_FillValue = -1.073742e+09f ;
> field8:grid_mapping = "rotated_pole0" ;
> field8:coordinates = "global_longitude0 global_latitude0" ;
> field8:stash_item = "222" ;
> field8:stash_section = "16" ;
> field8:field_code = "8" ;
> field8:standard_name = "air_pressure_at_sea_level" ;
> field8:long_name = "PRESSURE AT MEAN SEA LEVEL" ;
> field8:units = "Pa" ;
> field8:cell_method = "time: mean " ;
>
> // global attributes:
> :Conventions = "CF-1.3" ;
> :alpham = 0.5f ;
> data:
>
> time0 = 3182.5, 3187.5, 3192.5, 3197.5, 3202.5, 3207.5 ;
>
> time1 = 3182.5, 3187.5, 3192.5, 3197.5, 3202.5, 3207.5 ;
>
> z0 = 0 ;
>
> z1 = 0 ;
>
> z2 = -1 ;
>
> z3 = 200 ;
>
> z4 = 0 ;
>
> latitude0 = 13.42, 13.2, 12.98, 12.76, 12.54, 12.32, 12.1, 11.88, 11.66,
> 11.44, 11.22, 11, 10.78, 10.56, 10.34, 10.12, 9.900001, 9.68, 9.46,
> 9.240001, 9.02, 8.8, 8.58, 8.360001, 8.14, 7.920001, 7.7, 7.48, 7.26,
> 7.04, 6.82, 6.6, 6.380001, 6.16, 5.940001, 5.72, 5.5, 5.28, 5.06,
> 4.84,
> 4.62, 4.400001, 4.18, 3.960001, 3.74, 3.52, 3.3, 3.08, 2.86, 2.64,
> 2.42,
> 2.2, 1.98, 1.76, 1.54, 1.32, 1.1, 0.8800004, 0.6600004, 0.4400004,
> 0.2200004, 4.172325e-07, -0.2199996, -0.4399996, -0.6599996,
> -0.8799996,
> -1.1, -1.32, -1.54, -1.76, -1.98, -2.2, -2.42, -2.64, -2.86,
> -3.079999,
> -3.299999, -3.52, -3.74, -3.96, -4.179999, -4.4, -4.619999, -4.84,
> -5.059999, -5.28, -5.5, -5.719999, -5.94, -6.159999, -6.38, -6.599999,
> -6.82, -7.039999, -7.26, -7.48, -7.699999, -7.92, -8.139999, -8.36,
> -8.58, -8.799999, -9.02, -9.24, -9.459999, -9.679999, -9.9, -10.12,
> -10.34, -10.56, -10.78, -11, -11.22, -11.44, -11.66, -11.88, -12.1 ;
>
> longitude0 = -12.32001, -12.10001, -11.88, -11.66, -11.44, -11.22, -11,
> -10.78, -10.56, -10.34, -10.12, -9.899994, -9.680023, -9.460022,
> -9.240021, -9.02002, -8.800018, -8.580017, -8.360016, -8.140015,
> -7.920013, -7.700012, -7.480011, -7.26001, -7.040009, -6.820007,
> -6.600006, -6.380005, -6.160004, -5.940002, -5.720001, -5.5,
> -5.279999,
> -5.059998, -4.839996, -4.619995, -4.399994, -4.180023, -3.960022,
> -3.740021, -3.52002, -3.300018, -3.080017, -2.860016, -2.640015,
> -2.420013, -2.200012, -1.980011, -1.76001, -1.540009, -1.320007,
> -1.100006, -0.8800049, -0.6600037, -0.4400024, -0.2200012, 0,
> 0.2200012,
> 0.4400024, 0.6600037, 0.8800049, 1.100006, 1.319977, 1.539978,
> 1.759979,
> 1.97998, 2.199982, 2.419983, 2.639984, 2.859985, 3.079987, 3.299988,
> 3.519989, 3.73999, 3.959991, 4.179993, 4.399994, 4.619995, 4.839996,
> 5.059998, 5.279999, 5.5, 5.720001, 5.940002, 6.160004, 6.380005,
> 6.600006, 6.819977, 7.039978, 7.259979, 7.47998, 7.699982, 7.919983,
> 8.139984, 8.359985, 8.579987, 8.799988, 9.019989, 9.23999, 9.459991,
> 9.679993, 9.899994, 10.12, 10.34, 10.56, 10.78, 11, 11.22, 11.44,
> 11.66 ;
>
> latitude1 = 13.31, 13.09, 12.87, 12.65, 12.43, 12.21, 11.99, 11.77,
> 11.55,
> 11.33, 11.11, 10.89, 10.67, 10.45, 10.23, 10.01, 9.79, 9.57, 9.349999,
> 9.13, 8.91, 8.69, 8.469999, 8.25, 8.03, 7.81, 7.59, 7.37, 7.15, 6.93,
> 6.71, 6.49, 6.27, 6.05, 5.83, 5.61, 5.39, 5.17, 4.95, 4.73, 4.51,
> 4.29,
> 4.07, 3.85, 3.63, 3.41, 3.19, 2.97, 2.75, 2.53, 2.31, 2.09, 1.87,
> 1.65,
> 1.43, 1.21, 0.9899998, 0.7699998, 0.5499998, 0.3299998, 0.1099998,
> -0.1100002, -0.3300002, -0.5500002, -0.7700002, -0.9900002, -1.21,
> -1.43,
> -1.65, -1.87, -2.09, -2.31, -2.53, -2.75, -2.97, -3.19, -3.41, -3.63,
> -3.85, -4.07, -4.29, -4.51, -4.73, -4.95, -5.17, -5.39, -5.61, -5.83,
> -6.05, -6.27, -6.49, -6.71, -6.93, -7.15, -7.37, -7.59, -7.81, -8.03,
> -8.25, -8.47, -8.690001, -8.91, -9.13, -9.35, -9.57, -9.79, -10.01,
> -10.23, -10.45, -10.67, -10.89, -11.11, -11.33, -11.55, -11.77,
> -11.99 ;
>
> longitude1 = -12.20999, -11.98999, -11.76999, -11.54999, -11.32999,
> -11.10999, -10.88998, -10.66998, -10.44998, -10.22998, -10.00998,
> -9.789978, -9.570007, -9.350006, -9.130005, -8.910004, -8.690002,
> -8.470001, -8.25, -8.029999, -7.809998, -7.589996, -7.369995,
> -7.149994,
> -6.929993, -6.709991, -6.48999, -6.269989, -6.049988, -5.829987,
> -5.609985, -5.389984, -5.169983, -4.949982, -4.72998, -4.509979,
> -4.289978, -4.070007, -3.850006, -3.630005, -3.410004, -3.190002,
> -2.970001, -2.75, -2.529999, -2.309998, -2.089996, -1.869995,
> -1.649994,
> -1.429993, -1.209991, -0.9899902, -0.769989, -0.5499878, -0.3299866,
> -0.1099854, 0.1100159, 0.3300171, 0.5500183, 0.7700195, 0.9900208,
> 1.210022, 1.429993, 1.649994, 1.869995, 2.089996, 2.309998, 2.529999,
> 2.75, 2.970001, 3.190002, 3.410004, 3.630005, 3.850006, 4.070007,
> 4.290009, 4.51001, 4.730011, 4.950012, 5.170013, 5.390015, 5.610016,
> 5.830017, 6.050018, 6.27002, 6.490021, 6.710022, 6.929993, 7.149994,
> 7.369995, 7.589996, 7.809998, 8.029999, 8.25, 8.470001, 8.690002,
> 8.910004, 9.130005, 9.350006, 9.570007, 9.790009, 10.01001, 10.23001,
> 10.45001, 10.67001, 10.89001, 11.11002, 11.33002, 11.55002, 11.77002 ;
> }
>
> I'm not sure what to put in the PDEF function and/or if that works or if
> I should regrid.
>
> Please help
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
More information about the gradsusr
mailing list