Changeset 808 for trunk/UTIL/NCL/visu-utils.ncl

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).

File:

• trunk/UTIL/NCL/visu-utils.ncl (modified) (2 diffs)

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
+