[gradsusr] converting rotate pole grid

Roberto Mera rjmera at coas.oregonstate.edu
Fri Oct 15 15:58:30 EDT 2010


  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



More information about the gradsusr mailing list