Changeset 2683 for LMDZ5/trunk/libf/phylmd/dyn1d/1DUTILS.h

Timestamp:

Oct 24, 2016, 6:46:30 PM (8 years ago)

Author:

fhourdin

Message:

Corrections pour le cas Dice (Marie-Pierre Lefebvre)

File:

• LMDZ5/trunk/libf/phylmd/dyn1d/1DUTILS.h (modified) (13 diffs)

LMDZ5/trunk/libf/phylmd/dyn1d/1DUTILS.h

@@ -55 +55 @@
 
 !Config  Key  = prt_level
-!Config  Desc = niveau d'impressions de débogage
+!Config  Desc = niveau d'impressions de d?bogage
 !Config  Def  = 0
-!Config  Help = Niveau d'impression pour le débogage
+!Config  Help = Niveau d'impression pour le d?bogage
 !Config         (0 = minimum d'impression)
 !      prt_level = 0
@@ -134 +134 @@
         ENDIF
 
-!Paramètres de forçage
+!Param?tres de for?age
 !Config  Key  = tend_t
 !Config  Desc = forcage ou non par advection de T
@@ -1123 +1123 @@
 !----------------------------------------------------------------------
 !   Calcul de l'advection verticale (ascendance et subsidence) de 
-!   température et d'humidité. Hypothèse : ce qui rentre de l'extérieur
-!   a les mêmes caractéristiques que l'air de la colonne 1D (WTG) ou 
+!   temp?rature et d'humidit?. Hypoth?se : ce qui rentre de l'ext?rieur
+!   a les m?mes caract?ristiques que l'air de la colonne 1D (WTG) ou 
 !   sans WTG rajouter une advection horizontale 
 !----------------------------------------------------------------------  
@@ -1197 +1197 @@
 !----------------------------------------------------------------------
 !   Calcul de l'advection verticale (ascendance et subsidence) de 
-!   température et d'humidité. Hypothèse : ce qui rentre de l'extérieur
-!   a les mêmes caractéristiques que l'air de la colonne 1D (WTG) ou 
+!   temp?rature et d'humidit?. Hypoth?se : ce qui rentre de l'ext?rieur
+!   a les m?mes caract?ristiques que l'air de la colonne 1D (WTG) ou 
 !   sans WTG rajouter une advection horizontale 
 !----------------------------------------------------------------------  
@@ -2951 +2951 @@
        endif
        if (annee_ref.eq.1992 .and. day1.lt.day_ini_toga) then
-        print*,'TOGA-COARE a débuté le 1er Nov 1992 (jour julien=306)'
+        print*,'TOGA-COARE a d?but? le 1er Nov 1992 (jour julien=306)'
         print*,'Changer dayref dans run.def'
         stop
@@ -3783 +3783 @@
 !=====================================================================
       subroutine read_dice(fich_dice,nlevel,ntime                         &
-     &     ,zz,pres,th,qv,u,v,o3                                          &
+     &     ,zz,pres,t,qv,u,v,o3                                          &
      &     ,shf,lhf,lwup,swup,tg,ustar,psurf,ug,vg                        &
      &     ,hadvt,hadvq,hadvu,hadvv,w,omega)
@@ -3793 +3793 @@
 
 #include "netcdf.inc"
+#include "YOMCST.h"
 
       integer ntime,nlevel
@@ -3800 +3801 @@
       real*8 zz(nlevel)
 
-      real*8 th(nlevel),pres(nlevel)
+      real*8 th(nlevel),pres(nlevel),t(nlevel)
       real*8 qv(nlevel),u(nlevel),v(nlevel),o3(nlevel)
       real*8 shf(ntime),lhf(ntime),lwup(ntime),swup(ntime),tg(ntime)
@@ -3806 +3807 @@
       real*8 hadvt(nlevel,ntime),hadvq(nlevel,ntime),hadvu(nlevel,ntime)
       real*8 hadvv(nlevel,ntime),w(nlevel,ntime),omega(nlevel,ntime)
+      real*8 pzero
 
       integer nid, ierr
@@ -3812 +3814 @@
       integer var3didin(nbvar3d)
 
+      pzero=100000.
       ierr = NF_OPEN(fich_dice,NF_NOWRITE,nid)
       if (ierr.NE.NF_NOERR) then
@@ -3986 +3989 @@
          endif
 !          write(*,*)'lecture th ok',th
+           do k=1,nlevel
+             t(k)=th(k)*(pres(k)/pzero)**rkappa
+           enddo
 
 #ifdef NC_DOUBLE
@@ -4798 +4804 @@
 !
 !  Cette formule remplace d_q = (1/tau) [rh_targ - rh] qsat(T_new)
-!   qui n'était pas correcte.
+!   qui n'?tait pas correcte.
 !
             IF (tnew.LT.RTT) THEN
@@ -4873 +4879 @@
       END
 
+