Changeset 854

Timestamp:

Nov 27, 2012, 6:50:04 PM (12 years ago)

Author:

slebonnois

Message:

SL: correction d'un bug pour Titan (zlsm1) + quelques petits ajustements

Location:

trunk

Files:

• LMDZ.TITAN/deftank/gcm.def (modified) (1 diff)
• LMDZ.TITAN/libf/phytitan/physiq.F (modified) (2 diffs)
• LMDZ.TITAN/makegcm (deleted)
• LMDZ.VENUS/deftank/gcm.def (modified) (1 diff)
• LMDZ.VENUS/makegcm (deleted)
• UTIL/NCL/visu-utils.ncl (modified) (1 diff)

trunk/LMDZ.TITAN/deftank/gcm.def

@@ -9 +9 @@
 iperiod=5
 ## periode de la dissipation (en pas)
-idissip=5
+## dissip_period=5  PAR DEFAUT: 0 ie calcul autom. (=> 55 pour runs recents)
 ## choix de l'operateur de dissipation (star ou  non star )
 lstardis=y

trunk/LMDZ.TITAN/libf/phytitan/physiq.F

@@ -278 +278 @@
       REAL dist, rmu0(klon), fract(klon), pdecli
       REAL zday
-      REAL zls,zlsdeg
-      REAL zlsm1 ! for write_hist* (NB: it is not properly initialized!!)
+      REAL zls,zlsdeg,zlsm1
 c
       INTEGER i, k, iq, ig, j, ll, l
@@ -852 +851 @@
 
       CALL orbite(zls,dist,pdecli) 
+      IF (debut) zlsm1=zls
 
 c dans zenang, Ls en degres ; dans mucorr, Ls en radians

trunk/LMDZ.VENUS/deftank/gcm.def

@@ -9 +9 @@
 iperiod=5
 ## periode de la dissipation (en pas)
-idissip=5
+## dissip_period=5  PAR DEFAUT: 0 ie calcul autom. (=> 25 pour runs recents)
 ## choix de l'operateur de dissipation (star ou  non star )
 lstardis=y

trunk/UTIL/NCL/visu-utils.ncl

@@ -645 +645 @@
 prepared = "dummy"    ; DEFAULT
 
+if (labelvar.eq."ustar") then
+  if (isfilevar(infile,"vitu")) then
+    u = infile->vitu
+  else
+    u = infile->U
+  end if
+  myvar = dim_rmvmean_n_Wrap(u,3)
+  myvar@long_name="Deviation from zonal-mean zonal wind"
+  prepared = labelvar
+end if
+
+if (labelvar.eq."uprime") then
+  if (isfilevar(infile,"vitu")) then
+    u = infile->vitu
+  else
+    u = infile->U
+  end if
+  myvar = dim_rmvmean_n_Wrap(u,0)
+  myvar@long_name="Deviation from time-mean zonal wind"
+  prepared = labelvar
+end if
+
 if (labelvar.eq."dTdyn") then
   if (isfilevar(infile,"vitu")) then