Skip to: [content] [navigation]

Description


Ocean Physical Diagnostics

Plotting currents:

vector/k=1/len=0.5/xskip=1/yskip=1/x=115e:245w u[l=@ave],v[l=@ave];go land

Contouring Thermohaline circulation in depth space (gives Sverdrups):

sha/palette=no_green_centered/lev=(-30,30,3) ty_trans[x=@sum,k=@rsum,l=@ave]+ty_trans_gm[x=@sum,l=@ave]
contour/ov/lev=(-30,30,3) ty_trans[x=@sum,k=@rsum,l=@ave]+ty_trans_gm[x=@sum,l=@ave]

Contouring Thermohaline circulation in density space (gives Sverdrups):

sha/palette=no_green_centered/lev=(-30,30,3) ty_trans_rho[x=@sum,k=@rsum,l=@ave]
contour/ov/lev=(-30,30,3) ty_trans_rho[x=@sum,k=@rsum,l=@ave]

Listing ACC volume transport in Sverdrups (should be ~134 Sv):

list tx_trans[x=70w,y=-90:-40@sum,l=@ave,k=@sum]

Listing North Atlantic Deep Water volume transport at 25S in Sverdrups (should be ~13 Sv):

let nadw= ty_trans[x=70w:20E@sum,y=-25,k=@rsum]+ty_trans_gm[x=70w:20E@sum,y=-25]
list nadw[k=@min,l=@ave]

Listing Indonesian Throughflow volume transport in Sverdrups (should be ~-10 Sv):

list ty_trans[x=100e:135e@sum,y=0,l=@ave,k=@sum]

Listing Gulf Stream volume transport in Sverdrups (should be ~30 Sv):

list ty_trans[x=82w:70w@sum,y=27,l=@ave,z=0:1000@sum]

Listing Kuroshio volume transport in Sverdrups (should be ~50 Sv):

list tx_trans[x=144e,y=27:45@sum,l=@ave,z=0:1000@sum]

Listing Bering Strait volume transport in Sverdrups (should be 1 Sv):

list ty_trans[x=170w:160w@sum,y=65,l=@ave,k=@sum]

Listing deep (1000 m) equatorial upwelling (should be less than 2-3 SV)

list wt[z=1000,x=@din,y=-30:30@din,l=@ave]/1e6

List meridional heat flux (should be ~2.0 PW at 20N and ~-1.7 PW at 15-20S):

let htrans=(temp_adv_flux_y_int_z+temp_yflux_neutral_int_z+temp_yflux_sigma+temp_yflux_overflow_int_z)/1e15
plot htrans[l=@ave,x=@sum]

Salinity in Southern ocean, after Toggweiler and Samuels (1995), near the date line, deep salinity minumum should extend below higher salinity tongue to 65 S around 1000 m:

use "/archive/jpd/jakarta/preprocessing/om1new/woa01_salt_annual_om1p5.nc"
let salti=salt[d=2]*salt[d=1,l=240]/salt[d=1,l=240]
set win 1
sha/lev=(33.8),(34.2,34.8,0.01),(35.8)/palette=modulo salti[x=170w:178w@ave,y=-80:-30]
set win 2
sha/lev=(33.8),(34.2,34.8,0.01),(35.8)/palette=modulo salt[x=170w:178w@ave,y=-80:-30,l=240,d=1]

and at the prime meridian, the S=34.8 max should extend only as far as 55S. 1000m:

set win 1
sha/lev=(33.8,35,0.1),(35.8) salti[x=0,y=-80:-30]
set win 2
sha/lev=(33.8,35,0.1),(35.8) salt[x=0,y=-80:-30,l=240,d=1]

Ocean Biogeochemical Diagnostics

Calculating global integral of particle export:

list fpon[z=75,x=@din,y=@din,l=@ave]*86400*365*117/16*12/1e15

Calculating global integral of DOC export:

list (jdon[z=0:75@din,x=@din,y=@din,l=@ave]*117/16+jldoc[z=0:75@din,x=@din,y=@din,l=@ave])*86400*365*12/1e15

Creating an animation:

Create the .gifs in Ferret:
REPEAT/l=1:46 (sha/lev=(-2,2,0.1) log(CHL[k=1);go land; FRAME/FILE=whirl-`l`.gif)
create a file (whirlgif_infile in this case) that lists the images in order:
whirl-1.gif
whirl-2.gif
whirl-3.gif
...
run whirlfile to create the movie:
whirlgif -o chl_movie.gif -i whirlgif_infile -time 50

To remove labels in Ferret, use the pplus commands:
sha/lev=(-2,2,0.1)/set_up log(CHL[k=1) ! the /set_up initializes pplus
ppl list labels ! find the numeric value of the label
go unlabel 3 ! remove the third label
ppl shade ! make the plot

Creating a timeseries: Un-cpio all of the files containing the data of interest

cpio -iv < 19380101.nc.cpio
cpio -iv < 19380501.nc.cpio
cpio -iv < 19380901.nc.cpio
...

Concatenate all of the instances into a single file:

ncrcat -O -h -v chl *.ocean_bgc_phyto_sc.nc chl_timeseries.nc

• About us
  • Welcome
  • Visiting GFDL
  • Overview
  • Tour
  • Gallery
  • People
  • Employment
  • Help
  • Notices
• Research
  • Atmospheric Physics and Chemistry
  • Biospheric Processes
  • Climate Diagnostics
  • Climate Dynamics and Prediction
  • Oceans and Climate
  • Weather and Atmospheric Dynamics
• Products and Services
  • Flexible Modeling System
  • Our Data Portal
  • Our GForge
  • Technical Services
  • Visualization
• Reference
  • Our Online Bibliography
  • Our Meetings and Seminars
  • Our Library
• GFDL Only
  • External forms and links
  • External web mail
  • Help Desk
  • Internal Administrative Services
  • Internal Technical Services
  • NOAA's Hot Items
  • webTA Login

smaller bigger reset