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.
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/current-temp_1S1N.jnl @ 5387

Last change on this file since 5387 was 5387, checked in by flavoni, 10 years ago

start adding ferret scripts to generate plots for test cases, related to ticket #1469

  • Property svn:executable set to *
File size: 1.7 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_tprof_1d_00010101_00020520_grid_U.nc"
15set re/z=1         
16set w 1
17fill/lev=(-INF)(-0.2,0.2,0.01)(INF)/title="vel-U" uo[x=1:721,y=1:241,l=1,d=1]
18! "NAME OF OUTPUT FILE"
19FRAME/FILE="velU-at-1day.gif"
20set w 2
21fill/lev=(-INF)(-0.2,0.2,0.01)(INF)/title="vel-U" uo[x=1:721,y=1:241,l=100,d=1]
22! "NAME OF OUTPUT FILE"
23FRAME/FILE="velU-at-500days.gif"
24!!!!!!!!set w 3
25!!!!!!!fill/lev=(-INF)(-0.2,0.2,0.01)(INF)/title="vel-U" uo[x=1:721,y=1:241,l=@ave,d=1]
26! "NAME OF OUTPUT FILE"
27!!!!!!!FRAME/FILE="velU-average-1_500days.gif"
28!
29set w 3
30def axis/from_data/depth/name=z1 z[g=uo]
31ca re
32! method with indix, from list i see index corresponding to -1, 1
33!fill/lev=(-INF)(-0.0002,0.0002,0.00001)(INF)/title="vel-U" uo[gz=z1@asn,j=117:125@ave,l=100,d=1]
34define axis/x/from/modulo/units=degrees xaxis=nav_lon[j=1]
35define axis/y/from/units=degrees yaxis=nav_lat[i=1]
36let uo_mod = uo[gx=xaxis@asn,gy=yaxis@asn]
37shade/k=1/l=1/title="vel-U between 1south 1 north at 1 timestep"  uo_mod
38! "NAME OF OUTPUT FILE"
39FRAME /FILE="velU-1S1N-initial.gif"
40set w 4
41shade/k=1/l=100/title="vel-U between 1south 1 north after 500 days" uo_mod
42! "NAME OF OUTPUT FILE"
43FRAME /FILE="velU-1S1N-at500days.gif"
44!dataset=2
45ca re
46use "bx025_tprof_1d_00010101_00020520_grid_T.nc"
47def axis/from_data/depth/name=z1 z[g=thetao]
48set re/z=0:1000
49!set memory/size=20000000
50set w 5
51fill/lev=(-INF)(3,23,1)(INF)/title="temp-T" thetao[gz=z1@asn,y=@ave,l=100,d=2]
52! "NAME OF OUTPUT FILE"
53FRAME /FILE="tempT-1S1N.gif"
54!
Note: See TracBrowser for help on using the repository browser.