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

 ! NOAA/PMEL TMAP
 ! FERRET v6.85  
 ! Linux 2.6.32-358.23.2.el6.x86_64 64-bit - 11/12/13
 ! 24-Feb-15 14:35     

!start script :
! plot over all basin from 0 to 1000 m in z
! of velocities U and temperature T
!
!USE grid-U files
!
palette blue_darkred
!dataset=1
use "bx025_tprof_1d_00010101_00020520_grid_U.nc"
set re/z=1          
set w 1
fill/lev=(-INF)(-0.2,0.2,0.01)(INF)/title="vel-U" uo[x=1:721,y=1:241,l=1,d=1]
! "NAME OF OUTPUT FILE"
FRAME/FILE="velU-at-1day.gif"
set w 2
fill/lev=(-INF)(-0.2,0.2,0.01)(INF)/title="vel-U" uo[x=1:721,y=1:241,l=100,d=1]
! "NAME OF OUTPUT FILE"
FRAME/FILE="velU-at-500days.gif"
!!!!!!!!set w 3
!!!!!!!fill/lev=(-INF)(-0.2,0.2,0.01)(INF)/title="vel-U" uo[x=1:721,y=1:241,l=@ave,d=1]
! "NAME OF OUTPUT FILE"
!!!!!!!FRAME/FILE="velU-average-1_500days.gif"
!
set w 3
def axis/from_data/depth/name=z1 z[g=uo]
ca re
! method with indix, from list i see index corresponding to -1, 1
!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]
define axis/x/from/modulo/units=degrees xaxis=nav_lon[j=1]
define axis/y/from/units=degrees yaxis=nav_lat[i=1]
let uo_mod = uo[gx=xaxis@asn,gy=yaxis@asn]
shade/k=1/l=1/title="vel-U between 1south 1 north at 1 timestep"  uo_mod
! "NAME OF OUTPUT FILE"
FRAME /FILE="velU-1S1N-initial.gif"
set w 4
shade/k=1/l=100/title="vel-U between 1south 1 north after 500 days" uo_mod
! "NAME OF OUTPUT FILE"
FRAME /FILE="velU-1S1N-at500days.gif"
!dataset=2
ca re
use "bx025_tprof_1d_00010101_00020520_grid_T.nc"
def axis/from_data/depth/name=z1 z[g=thetao]
set re/z=0:1000
!set memory/size=20000000
set w 5
fill/lev=(-INF)(3,23,1)(INF)/title="temp-T" thetao[gz=z1@asn,y=@ave,l=100,d=2]
! "NAME OF OUTPUT FILE"
FRAME /FILE="tempT-1S1N.gif"
!