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

 ! 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_TScst_5d_00010101_00020520_grid_U.nc"
set w 1
!  to have z down
def axis/from_data/depth/name=z1 z[g=uo]
!  to have longitude instead of index points
define axis/x/from/modulo/units=degrees xaxis=nav_lon[j=1]
!  to have latitude instead of index points
define axis/y/from/units=degrees yaxis=nav_lat[i=1]
let uo_mod = uo[gx=xaxis@asn,gy=yaxis@asn]
fill/lev=(-INF)(-1,1,0.1)(INF)/k=1/l=1/title="TScst-vel-U SURFACE at 1 timestep"  uo_mod
! "NAME OF OUTPUT FILE"
FRAME /FILE="TScst-velU-surf-initial.gif"
set w 2
fill/lev=(-INF)(-1,1,0.1)(INF)/k=1/l=100/title="TScst-vel-U SURFACE after 500 days" uo_mod
! "NAME OF OUTPUT FILE"
FRAME /FILE="TScst-velU-surf-at500days.gif"
set w 3
fill/lev=(-INF)(-1,1,0.1)(INF)/title="TScst-vel-U prof between 1S and 1N at 1 timestep"  uo_mod[y=1S:1N@ave,z=1:1000,l=1,d=1]
! "NAME OF OUTPUT FILE"
FRAME /FILE="TScst-velU-prof-1S1N-initial.gif"
set w 4
fill/lev=(-INF)(-1,1,0.1)(INF)/title="TScst-vel-U prof between 1S and 1N after 500 days" uo_mod[y=1S:1N@ave,z=1:1000,l=100,d=1]
! "NAME OF OUTPUT FILE"
FRAME /FILE="TScst-velU-prof-1S1N-at500days.gif"
!dataset=2
use "bx025_TScst_5d_00010101_00020520_grid_T.nc"
def axis/from_data/depth/name=z1 z[g=thetao]
define axis/x/from/modulo/units=degrees xaxis=nav_lon[j=1]
define axis/y/from/units=degrees yaxis=nav_lat[i=1]
let thetao_mod = thetao[gx=xaxis@asn,gy=yaxis@asn]
set w 5
fill/lev=(-INF)(3,23,1)(INF)/title="TScst-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]
! "NAME OF OUTPUT FILE"
FRAME /FILE="TScst-tempT-1S1N-initial.gif"
set w 6
fill/lev=(-INF)(3,23,1)(INF)/title="TScst-temp-T between 1S and 1N after 500 days N" thetao_mod[gz=z1@asn,y=1S:1N@ave,z=1:1000,l=100,d=2]
! "NAME OF OUTPUT FILE"
FRAME /FILE="TScst-tempT-1S1N-at500days.gif"
!