Changeset 808 for trunk/UTIL

Timestamp:

Oct 16, 2012, 12:57:35 PM (12 years ago)

Author:

slebonnois

Message:

SL: Many changes for VENUS (related to newstart) and TITAN (related to clouds). Please read DOC/chantiers/commit_importants.log (cf v808).

Location:

trunk/UTIL/NCL

Files:

• help.txt (modified) (1 diff)
• planeto.ncl (modified) (2 diffs)
• visu-ncl (modified) (5 diffs)
• visu-utils.ncl (modified) (2 diffs)

trunk/UTIL/NCL/help.txt

@@ -46 +46 @@
 
   Default: liste => list the variables in file
+
+  If varname = customvar then a user-defined variable can be plotted
+  with -v customvar,labelvar
+  This labelvar must be defined in the customVar routine in visu-utils.ncl
 
 [-v2 varname]                    : Variable name to be overplotted

trunk/UTIL/NCL/planeto.ncl

@@ -28 +28 @@
     tname=liste(i)
   end if
-  if ((liste(i).eq."lev").or.(liste(i).eq."presnivs")) then
+  if ((liste(i).eq."lev").or.(liste(i).eq."presnivs").or.(liste(i).eq."levgrnd")) then
     pname=liste(i)
   end if
@@ -86 +86 @@
 ;---- load variable(s)
 ;---------------------
-  variable=infile->$var$
-  nbdim=dimsizes(filevardimsizes(infile,var))
-
+
+  if (var.ne."custom") then
+    variable=infile->$var$
+    nbdim=dimsizes(filevardimsizes(infile,var))
+  else
+    variable=customVar(infile,labelvar)
+    nbdim=dimsizes(dimsizes(variable))
+  end if
+  
   overplot=False
   diffdims=False

trunk/UTIL/NCL/visu-ncl

@@ -1 +1 @@
 #!/bin/bash
 
-# Directory where the NCL files are located ** EDIT **
+# Directory where the NCL files are located  ** EDIT **
 
 export NCLDIR=~/LMDZ5/UTIL/NCL
@@ -11 +11 @@
 inFile="nofile"
 var="liste"
+labelvar="dummy"
 var2="dummy"
 media="x11"
@@ -56 +57 @@
           inFile=$2 ; shift ; shift ;;
       "-v")
-          var=$2 ; shift ; shift ;;
+          var=${2%%,*}
+          tmp=${2#*,}
+          labelvar=${tmp%%,*}
+          if [ $labelvar != $var ]; then
+            labelvar=${tmp#*,}
+          fi
+          shift ; shift ;;
       "-v2")
           var2=$2 ; shift ; shift ;;
@@ -190 +197 @@
 ncl inFile=\"$inFile\"  \
     var=\"$var\" \
+    labelvar=\"$labelvar\" \
     var2=\"$var2\" \
     media=\"$media\" \
@@ -220 +228 @@
     planeto.ncl
 
-if [ -f *pdf ]; then
+if [ -f *pdf ]; then 
    mv *pdf $SIMUDIR
 fi
 
-if [ -f *ps ]; then
+if [ -f *ps ]; then 
    mv *ps $SIMUDIR
 fi

trunk/UTIL/NCL/visu-utils.ncl

@@ -56 +56 @@
   if (nbdim.eq.4) then
        vartrim=var({$dimname(0)$|mindimval(0):maxdimval(0)},{$dimname(1)$|mindimval(1):maxdimval(1)},{$dimname(2)$|mindimval(2):maxdimval(2)},{$dimname(3)$|mindimval(3):maxdimval(3)})
+;       vartrim=var($dimname(0)$|0:0,{$dimname(1)$|mindimval(1):maxdimval(1)},{$dimname(2)$|mindimval(2):maxdimval(2)},{$dimname(3)$|mindimval(3):maxdimval(3)})
   end if
 
@@ -628 +629 @@
 end
 
+;;;;;;;;;;;; customVar ;;;;;;;;;;;;;;;;
+
+undef("customVar")
+function customVar(infile,labelvar)
+;============================
+
+; infile   : input file
+; labelvar : input label for variable
+
+; THIS IS TO BE CUSTOMIZED EACH TIME ACCORDING TO YOUR NEEDS 
+; because it is impossible to make it automatic...
+
+begin
+
+prepared = "dummy"    ; DEFAULT
+
+if (labelvar.eq."dTdyn") then
+  if (isfilevar(infile,"vitu")) then
+    dtdyn = infile->dtdyn
+  else
+    dtdyn = infile->DTCORE
+  end if
+  myvar = dtdyn*1e7
+  copy_VarMeta(dtdyn,myvar)
+  myvar@units    ="K/Vd"
+  myvar@long_name="Dynamics heating rate"
+  prepared = labelvar
+end if
+
+if (labelvar.eq."dTadj") then
+  if (isfilevar(infile,"vitu")) then
+    dtadj = infile->dtajs
+  else
+    dtadj = infile->DADJ
+  end if
+  myvar = dtadj*1e7
+  copy_VarMeta(dtadj,myvar)
+  myvar@units    ="K/Vd"
+  myvar@long_name="Dry convection heating rate"
+  prepared = labelvar
+end if
+
+if (labelvar.eq."dTvdf") then
+  if (isfilevar(infile,"vitu")) then
+    dtvdf = infile->dtvdf
+  else
+    dtvdf = infile->DTV
+  end if
+  myvar = dtvdf*1e7
+  copy_VarMeta(dtvdf,myvar)
+  myvar@units    ="K/Vd"
+  myvar@long_name="PBL heating rate"
+  prepared = labelvar
+end if
+
+if (labelvar.eq."dTrad") then
+  if (isfilevar(infile,"vitu")) then
+    dtlwr = infile->dtlwr
+    dtswr = infile->dtswr
+  else
+    dtlwr = infile->QRL
+    dtswr = infile->QRS
+  end if
+  myvar = (dtlwr+dtswr)*1e7
+  copy_VarMeta(dtlwr,myvar)
+  myvar@units    ="K/Vd"
+  myvar@long_name="Radiative balance"
+  prepared = labelvar
+end if
+
+if (labelvar.eq."Tbal") then
+  if (isfilevar(infile,"vitu")) then
+    dtlwr = infile->dtlwr
+    dtswr = infile->dtswr
+    dtajs = infile->dtajs
+    dtdyn = infile->dtdyn
+    dtvdf = infile->dtvdf
+  else
+    dtlwr = infile->QRL
+    dtswr = infile->QRS
+    dtajs = infile->DADJ
+    dtdyn = infile->DTCORE
+    dtvdf = infile->DTV
+  end if
+  myvar = (dtlwr+dtswr+dtajs+dtdyn+dtvdf)*1e7
+  copy_VarMeta(dtlwr,myvar)
+  myvar@units    ="K/Vd"
+  myvar@long_name="Thermal balance"
+  prepared = labelvar
+end if
+
+if (labelvar.eq."Tforc") then
+  if (isfilevar(infile,"vitu")) then
+    dtlwr = infile->dtlwr
+    dtswr = infile->dtswr
+    dtajs = infile->dtajs
+    if (isfilevar(infile,"dtvdf")) then
+          dtvdf = infile->dtvdf
+    else
+          dtvdf = dtajs*0.
+    end if
+  else
+    dtlwr = infile->QRL
+    dtswr = infile->QRS
+    dtajs = infile->DADJ
+    if (isfilevar(infile,"DTV")) then
+          dtvdf = infile->DTV
+    else
+          dtvdf = dtajs*0.
+    end if
+  end if
+  myvar = (dtlwr+dtswr+dtajs+dtvdf)*1e7
+  copy_VarMeta(dtlwr,myvar)
+  myvar@units    ="K/Vd"
+  myvar@long_name="Thermal physical forcing"
+  prepared = labelvar
+end if
+
+if (labelvar.eq."ske") then
+  if (isfilevar(infile,"vitu")) then
+    u = infile->vitu
+    v = infile->vitv
+  else
+    u = infile->U
+    v = infile->V
+  end if
+  myvar = u*u+v*v
+  copy_VarMeta(u,myvar)
+  myvar@units    ="m2/s2"
+  myvar@long_name="Specific kinetic energy"
+  prepared = labelvar
+end if
+
+if (prepared.eq."dummy") then   ; DEFAULT: LEAVE AS IS
+  print("You chose a customized variable !")
+  print("Modify the function customVar in visu-utils.ncl to prepare this variable")
+  exit
+end if
+
+  return(myvar)
+
+end
+