New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
kelv-current-temp_1S1N.jnl in branches/2015/dev_r5102_CNRS11_TestCase/NEMOGCM/TOOLS/TEST_CASES/CREATE_PLOTS – NEMO

source: branches/2015/dev_r5102_CNRS11_TestCase/NEMOGCM/TOOLS/TEST_CASES/CREATE_PLOTS/kelv-current-temp_1S1N.jnl @ 5435

Last change on this file since 5435 was 5435, checked in by flavoni, 9 years ago

add ferret scripts creating plot for TEST CASE

  • Property svn:executable set to *
File size: 2.1 KB
Line 
1 ! NOAA/PMEL TMAP
2 ! FERRET v6.85 
3 ! Linux 2.6.32-358.23.2.el6.x86_64 64-bit - 11/12/13
4 ! 24-Feb-15 14:35     
5
6!start script :
7! plot over all basin from 0 to 1000 m in z
8! of velocities U and temperature T
9!
10!USE grid-U files
11!
12palette blue_darkred
13!dataset=1
14use "bx025_kelvin_1d_00010911_00020520_grid_U.nc"
15set w 1
16!  to have z down
17def axis/from_data/depth/name=z1 z[g=uo]
18!  to have longitude instead of index points
19define axis/x/from/modulo/units=degrees xaxis=nav_lon[j=1]
20!  to have latitude instead of index points
21define axis/y/from/units=degrees yaxis=nav_lat[i=1]
22let uo_mod = uo[gx=xaxis@asn,gy=yaxis@asn]
23fill/lev=(-INF)(-1,1,0.1)(INF)/k=1/l=1/title="kelv-vel-U SURFACE at 1 timestep"  uo_mod
24! "NAME OF OUTPUT FILE"
25FRAME /FILE="Kelv-velU-initial.gif"
26set w 2
27fill/lev=(-INF)(-1,1,0.1)(INF)/k=1/l=250/title="kelv-vel-U between 1S and 1N after 500 days" uo_mod
28! "NAME OF OUTPUT FILE"
29FRAME /FILE="Kelv-velU-at500days.gif"
30set re/z=0:1000
31set w 3
32fill/lev=(-INF)(-1,1,0.1)(INF)/title="kelv-vel-U prof between 1S and 1N at 1 timestep"  uo_mod[gz=z1@asn,gx=xaxis@asn,y=1S:1N@ave,z=1:1000,l=1,d=1]
33! "NAME OF OUTPUT FILE"
34FRAME /FILE="Kelv-velU-prof-1S1N-initial.gif"
35set w 4
36fill/lev=(-INF)(-1,1,0.1)(INF)/title="kelv-vel-U prof between 1S and 1N after 500 days" uo_mod[gz=z1@asn,gx=xaxis@asn,y=1S:1N@ave,z=1:1000,l=250,d=1]
37! "NAME OF OUTPUT FILE"
38FRAME /FILE="Kelv-velU-prof-1S1N-at500days.gif"
39!dataset=2
40use "bx025_kelvin_1d_00010911_00020520_grid_T.nc"
41def axis/from_data/depth/name=z1 z[g=thetao]
42define axis/x/from/modulo/units=degrees xaxis=nav_lon[j=1]
43define axis/y/from/units=degrees yaxis=nav_lat[i=1]
44let thetao_mod = thetao[gx=xaxis@asn,gy=yaxis@asn]
45set w 5
46fill/lev=(-INF)(3,23,1)(INF)/title="kelv-temp-T between 1S and 1N at 1 timestep" thetao_mod[gz=z1@asn,y=1S:1N@ave,z=1:1000,l=1,d=2]
47! "NAME OF OUTPUT FILE"
48FRAME /FILE="Kelv-tempT-1S1N-initial.gif"
49set w 6
50fill/lev=(-INF)(3,23,1)(INF)/title="kelv-temp-T between 1S and 1N after 500 days N" thetao_mod[gz=z1@asn,y=1S:1N@ave,z=1:1000,l=250,d=2]
51! "NAME OF OUTPUT FILE"
52FRAME /FILE="Kelv-tempT-1S1N-at500days.gif"
53!
Note: See TracBrowser for help on using the repository browser.