;iasnfs_forecast_2010010112_2010033112 ; float sea_surface_height_above_sea_level ( time, latitude, longitude ) ; float sea_water_potential_temperature ( time, zlev, latitude, longitude ) ; and for future use... ; float eastward_sea_water_velocity ( time, zlev, latitude, longitude ) ; float northward_sea_water_velocity ( time, zlev, latitude, longitude ) ; ; float longitude ( longitude ) units : degree_east ; float latitude ( latitude ) units : degree_north ; float zlev ( zlev ) units : m positive : down ; float time ( time ) units : day start : 2010-01-01 12:00:00 UTC load "/opt/NCL/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/calendar_decode2.ncl" ; ### ; ### This is the first prototype of the model intercomparisons. ; ### This version uses Step 2 data instead of Step 3 data but that's an easy fix. ; ### ; ### Info on file and some variable names follows: ; ### iasnfs_forecast_2010010112_2010033112 ; ### float sea_surface_height_above_sea_level ( time, latitude, longitude ) ; ### float sea_water_potential_temperature ( time, zlev, latitude, longitude ) ; ### and for future use... ; ### float eastward_sea_water_velocity ( time, zlev, latitude, longitude ) ; ### float northward_sea_water_velocity ( time, zlev, latitude, longitude ) ; ### ; ### float longitude ( longitude ) units : degree_east ; ### float latitude ( latitude ) units : degree_north ; ### float zlev ( zlev ) units : m positive : down ; ### float time ( time ) units : day start : 2010-01-01 12:00:00 UTC fnameModel1="/home/gerry/GOMEX-PPP/Step2/iasnfs_forecast_2010010112_2010033112.nc" model1=addfile(fnameModel1, "r") ; type = "pdf" ; type = "ps" ; type = "ncgm" type = "png" ; type="x11" Time=model1->time ; Elapsed time in DAYS since 2010-01-01 12:00:00 UTC model1Lat=model1->latitude ; Degrees NORTH model1Lon=model1->longitude ; Degrees EAST model1Depth=model1->zlev ; Positive down, meters (22 levels) ;model1Ssh=100 * model1->sea_surface_height_above_sea_level ;Derived from ncl_filedump variable names (and interactively tested) model1Ssh=model1->sea_surface_height_above_sea_level model1Ssh = model1Ssh*100 model1Sst=model1->sea_water_potential_temperature model1CurrentU=model1->eastward_sea_water_velocity model1CurrentV=model1->northward_sea_water_velocity wks = gsn_open_wks(type,"Model1-SSH") gsn_define_colormap(wks,"BlAqGrYeOrRe") ; choose colormap ntimes=dimsizes(Time) res = True ; plot mods desired res@cnFillOn = True ; turn on color for contours res@cnLinesOn = False ; turn off contour lines res@cnLineLabelsOn = False ; turn off contour line labels ; res@gsnScalarContour = True ; contours desired res@gsnSpreadColors = True ; use full color map res@gsnSpreadColorStart = 4 ; start at color 17 res@gsnSpreadColorEnd = 200 ; end at color 200 i = NhlNewColor(wks,0.7,0.7,0.7) ; add gray to colormap res@tiMainString = "Model 1" res@gsnRightString = "Sea Surface Height (cm ref. to MSL)" res@tiMainOn = True res@mpLandFillColor = "gray" ; set land to be gray res@mpMinLonF = -98. ; Select a subregion res@mpMaxLonF = -78. res@mpMinLatF = 18. res@mpMaxLatF = 31. res@lbOrientation = "Vertical" ; vertical label bar res@pmLabelBarOrthogonalPosF = -0.01 ; move label bar closer res@lbLabelStride = 1 res@cnExplicitLabelBarLabelsOn = True ; Explicit label-bar values res@cnMinLevelValF = -100 ; Minimum SSH variation from MSL res@cnMaxLevelValF = 80 ; Maximum SSH above MSL res@cnLevelSelectionMode = "ManualLevels" ; Manually define levels ; res@mpFillOn = True ; gray landmass ; res@mpGeophysicalLineColor = "Navy" ; GeoLine color navy blue ; gsn_define_colormap(wks,"BlAqGrYeOrRevi200") ; choose color map do timeCounter=0,ntimes-1,6 ; res@txString = "2010 DOY:" + timeCounter ; 2010 Day-Of-Year res@gsnLeftString = "MJD = 2010:" + timeCounter ; plotSsh=gsn_csm_contour_map_ce(wks,model1Ssh(timeCounter,:,:),res) plotSsh=gsn_csm_contour_map(wks,model1Ssh(timeCounter,:,:),res) end do