PRO PseudoObs FOR loop = 1, 6 DO BEGIN IF (loop EQ 1) THEN BEGIN ; Total analysis increments for standard system ; --------------------------------------------- PostscriptFile = 'Total_Incs_St.ps' ColType = 2 ;Blue/Red Title1 = 'Total u analysis increment' File1 = '../../PVplots_Jan07/PseudoObsTest_St_withVertRegress/PFAnal.nc' Long1 = 'longitude' Lat1 = 'latitude' Lev1 = 'hybrid_ht' Name1 = 'u' Title2 = 'Total v analysis increment' File2 = '../../PVplots_Jan07/PseudoObsTest_St_withVertRegress/PFAnal.nc' Long2 = 'longitude_1' Lat2 = 'latitude_1' Lev2 = 'hybrid_ht' Name2 = 'v' Title3 = 'Total p analysis increment' File3 = '../../PVplots_Jan07/PseudoObsTest_St_withVertRegress/PFAnal.nc' Long3 = 'longitude_1' Lat3 = 'latitude' Lev3 = 'hybrid_ht' Name3 = 'p' ENDIF IF (loop EQ 2) THEN BEGIN ; Total analysis increments for PV system ; --------------------------------------- PostscriptFile = 'Total_Incs_PV.ps' ColType = 2 ;Blue/Red Title1 = 'Total u analysis increment' File1 = 'PFAnal.nc' Long1 = 'longitude' Lat1 = 'latitude' Lev1 = 'hybrid_ht' Name1 = 'u' Title2 = 'Total v analysis increment' File2 = 'PFAnal.nc' Long2 = 'longitude_1' Lat2 = 'latitude_1' Lev2 = 'hybrid_ht' Name2 = 'v' Title3 = 'Total p analysis increment' File3 = 'PFAnal.nc' Long3 = 'longitude_1' Lat3 = 'latitude' Lev3 = 'hybrid_ht' Name3 = 'p' ENDIF IF (loop EQ 3) THEN BEGIN ; Balanced analysis increments for PV system ; ------------------------------------------ PostscriptFile = 'Balanced_Incs_PV.ps' ColType = 2 ;Blue/Red Title1 = 'Balanced u analysis increment' File1 = 'Balanced_uv.nc' Long1 = 'longitude' Lat1 = 'latitude' Lev1 = 'hybrid_ht' Name1 = 'unspecified' Title2 = 'Balanced v analysis increment' File2 = 'Balanced_uv.nc' Long2 = 'longitude_1' Lat2 = 'latitude_1' Lev2 = 'hybrid_ht' Name2 = 'unspecified_1' Title3 = 'Balanced p analysis increment' File3 = 'Pressures.nc' Long3 = 'longitude' Lat3 = 'latitude' Lev3 = 'hybrid_ht' Name3 = 'unspecified' ENDIF IF (loop EQ 4) THEN BEGIN ; Unbalanced analysis increments for PV system ; -------------------------------------------- PostscriptFile = 'Unbalanced_Incs_PV.ps' ColType = 2 ;Blue/Red Title1 = 'Unbalanced u analysis increment' File1 = 'Unbalanced_uv.nc' Long1 = 'longitude' Lat1 = 'latitude' Lev1 = 'hybrid_ht' Name1 = 'unspecified' Title2 = 'Unbalanced v analysis increment' File2 = 'Unbalanced_uv.nc' Long2 = 'longitude_1' Lat2 = 'latitude_1' Lev2 = 'hybrid_ht' Name2 = 'unspecified_1' Title3 = 'Unbalanced p analysis increment' File3 = 'Pressures.nc' Long3 = 'longitude' Lat3 = 'latitude' Lev3 = 'hybrid_ht' Name3 = 'unspecified_1' ENDIF IF (loop EQ 5) THEN BEGIN ; Standard control variables ; -------------------------- PostscriptFile = 'ControlVars_St.ps' ColType = 1 ;Multi colour Title1 = 'Streamfunction control variable' File1 = '../../PVplots_Jan07/PseudoObsTest_St_withVertRegress/ControlVars.nc' Long1 = 'longitude' Lat1 = 'latitude' Lev1 = 'hybrid_ht' Name1 = 'unspecified' Title2 = 'Unbalanced pressure control variable' File2 = '../../PVplots_Jan07/PseudoObsTest_St_withVertRegress/ControlVars.nc' Long2 = 'longitude_1' Lat2 = 'latitude_1' Lev2 = 'hybrid_ht' Name2 = 'unspecified_2' Title3 = 'Velocity potential control variable' File3 = '../../PVplots_Jan07/PseudoObsTest_St_withVertRegress/ControlVars.nc' Long3 = 'longitude_1' Lat3 = 'latitude_1' Lev3 = 'hybrid_ht' Name3 = 'unspecified_1' ENDIF IF (loop EQ 6) THEN BEGIN ; Standard control variables ; -------------------------- PostscriptFile = 'ControlVars_PV.ps' ColType = 1 ;Multi colour Title1 = 'Balanced streamfunction control variable' File1 = 'ControlVars_min.nc' Long1 = 'longitude' Lat1 = 'latitude' Lev1 = 'hybrid_ht' Name1 = 'unspecified' Title2 = 'Unbalanced pressure control variable' File2 = 'ControlVars_min.nc' Long2 = 'longitude_1' Lat2 = 'latitude_1' Lev2 = 'hybrid_ht' Name2 = 'unspecified_2' Title3 = 'Velocity potential control variable' File3 = 'ControlVars_min.nc' Long3 = 'longitude_1' Lat3 = 'latitude_1' Lev3 = 'hybrid_ht' Name3 = 'unspecified_1' ENDIF PRINT, loop ;----------------------------------------------------------- ; Set-up postscript file for standard scheme output ;----------------------------------------------------------- SET_PLOT, 'ps' DEVICE, FILE=PostscriptFile, /PORTRAIT, /COLOR, XSIZE=40.0, YSIZE=20.0, $ XOFFSET=2.0, YOFFSET=2.0, FONT_SIZE=8, /ENCAPSULATED SetColours, ColType ;----------------------------------------------------------- ; Field 1, level 20 ;----------------------------------------------------------- fileid = NCDF_OPEN(File1) varidx = NCDF_VARID(fileid, Long1) varidy = NCDF_VARID(fileid, Lat1) varidz = NCDF_VARID(fileid, Lev1) varid = NCDF_VARID(fileid, Name1) NCDF_VARGET, fileid, varidx, x_degrees NCDF_VARGET, fileid, varidy, y_degrees NCDF_VARGET, fileid, varidz, vertlevs NCDF_VARGET, fileid, varid, Data NCDF_CLOSE, fileid FOR z = 0, N_ELEMENTS(vertlevs)-1 DO BEGIN vertlevs(z) = z+1 ENDFOR BoxPlot, Data(*,*,20), 1, x_degrees, y_degrees, 0.05, 0.55, 0.25, 0.9, $ -180, Title1+' level 20', 1, 1, 1 ;----------------------------------------------------------- ; Field 1, Greenwich ;----------------------------------------------------------- DataGreenwich = FINDGEN(N_ELEMENTS(y_degrees),N_ELEMENTS(vertlevs)) ;PRINT, N_ELEMENTS(y_degrees), N_ELEMENTS(vertlevs) FOR y = 0, N_ELEMENTS(y_degrees)-1 DO BEGIN FOR z = 0, N_ELEMENTS(vertlevs)-1 DO BEGIN DataGreenwich(y,z) = Data(0,y,z) ENDFOR ENDFOR BoxPlot, DataGreenwich(*,*), 1, y_degrees, vertlevs, 0.05, 0.1, 0.25, 0.45, $ 0.0, Title1+' Greenwich', 1, 0, 1 ;----------------------------------------------------------- ; Field 2, level 20 ;----------------------------------------------------------- fileid = NCDF_OPEN(File2) varidx = NCDF_VARID(fileid, Long2) varidy = NCDF_VARID(fileid, Lat2) varidz = NCDF_VARID(fileid, Lev2) varid = NCDF_VARID(fileid, Name2) NCDF_VARGET, fileid, varidx, x_degrees NCDF_VARGET, fileid, varidy, y_degrees NCDF_VARGET, fileid, varidz, vertlevs NCDF_VARGET, fileid, varid, Data NCDF_CLOSE, fileid FOR z = 0, N_ELEMENTS(vertlevs)-1 DO BEGIN vertlevs(z) = z+1 ENDFOR BoxPlot, Data(*,*,20), 1, x_degrees, y_degrees, 0.35, 0.55, 0.55, 0.9, $ -180.0, Title2+' level 20', 1, 1, 1 ;----------------------------------------------------------- ; Total meridional wind, Greenwich ;----------------------------------------------------------- DataGreenwich = FINDGEN(N_ELEMENTS(y_degrees),N_ELEMENTS(vertlevs)) ;PRINT, N_ELEMENTS(y_degrees), N_ELEMENTS(vertlevs) FOR y = 0, N_ELEMENTS(y_degrees)-1 DO BEGIN FOR z = 0, N_ELEMENTS(vertlevs)-1 DO BEGIN DataGreenwich(y,z) = Data(0,y,z) ENDFOR ENDFOR BoxPlot, DataGreenwich(*,*), 1, y_degrees, vertlevs, 0.35, 0.1, 0.55, 0.45, $ 0.0, Title2+' Greenwich', 1, 0, 1 ;----------------------------------------------------------- ; Total pressure, level 20 ;----------------------------------------------------------- fileid = NCDF_OPEN(File3) varidx = NCDF_VARID(fileid, Long3) varidy = NCDF_VARID(fileid, Lat3) varidz = NCDF_VARID(fileid, Lev3) varid = NCDF_VARID(fileid, Name3) NCDF_VARGET, fileid, varidx, x_degrees NCDF_VARGET, fileid, varidy, y_degrees NCDF_VARGET, fileid, varidz, vertlevs NCDF_VARGET, fileid, varid, Data NCDF_CLOSE, fileid FOR z = 0, N_ELEMENTS(vertlevs)-1 DO BEGIN vertlevs(z) = z+1 ENDFOR BoxPlot, Data(*,*,20), 1, x_degrees, y_degrees, 0.65, 0.55, 0.85, 0.9, $ -180.0, Title3+' level 20', 1, 1, 1 ;----------------------------------------------------------- ; Total meridional wind, Greenwich ;----------------------------------------------------------- DataGreenwich = FINDGEN(N_ELEMENTS(y_degrees),N_ELEMENTS(vertlevs)) ;PRINT, N_ELEMENTS(y_degrees), N_ELEMENTS(vertlevs) FOR y = 0, N_ELEMENTS(y_degrees)-1 DO BEGIN FOR z = 0, N_ELEMENTS(vertlevs)-1 DO BEGIN DataGreenwich(y,z) = Data(0,y,z) ENDFOR ENDFOR BoxPlot, DataGreenwich(*,*), 1, y_degrees, vertlevs, 0.65, 0.1, 0.85, 0.45, $ 0.0, Title3+' Greenwich', 1, 0, 1 DEVICE, /CLOSE ENDFOR END