Changeset 2683

Timestamp:

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

Author:

fhourdin

Message:

Corrections pour le cas Dice (Marie-Pierre Lefebvre)

Location:

LMDZ5/trunk/libf/phylmd/dyn1d

Files:

• 1DUTILS.h (modified) (13 diffs)
• 1D_decl_cases.h (modified) (3 diffs)
• 1D_interp_cases.h (modified) (2 diffs)
• 1D_read_forc_cases.h (modified) (7 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
 
+

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

@@ -146 +146 @@
         
         real zz_dice(nlev_dice)
-        real th_dice(nlev_dice),qv_dice(nlev_dice)
+        real t_dice(nlev_dice),qv_dice(nlev_dice)
         real u_dice(nlev_dice), v_dice(nlev_dice),o3_dice(nlev_dice)
         real ht_dice(nlev_dice,nt_dice)
@@ -153 +153 @@
         real w_dice(nlev_dice,nt_dice),omega_dice(nlev_dice,nt_dice)
         real o3_mod(llm),hu_mod(llm),hv_mod(llm)
-        real th_dicei(nlev_dice),qv_dicei(nlev_dice)
+        real t_dicei(nlev_dice),qv_dicei(nlev_dice)
         real u_dicei(nlev_dice), v_dicei(nlev_dice),o3_dicei(nlev_dice)
         real ht_dicei(nlev_dice)
@@ -287 +287 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
+

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

@@ -118 +118 @@
 ! vertical interpolation:
       CALL interp_dice_vertical(play,nlev_dice,nt_dice,plev_dice        &
-     &         ,th_dice,qv_dice,u_dice,v_dice,o3_dice                   &
+     &         ,t_dice,qv_dice,u_dice,v_dice,o3_dice                   &
      &         ,ht_profd,hq_profd,hu_profd,hv_profd,w_profd,omega_profd &
-     &         ,th_mod,qv_mod,u_mod,v_mod,o3_mod                        &
+     &         ,t_mod,qv_mod,u_mod,v_mod,o3_mod                        &
      &         ,ht_mod,hq_mod,hu_mod,hv_mod,w_mod,omega_mod,mxcalc)
 !     do l = 1, llm
@@ -810 +810 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
+

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

@@ -367 +367 @@
       fich_dice='dice_driver.nc'
       call read_dice(fich_dice,nlev_dice,nt_dice                    &
-     &     ,zz_dice,plev_dice,th_dice,qv_dice,u_dice,v_dice,o3_dice &
+     &     ,zz_dice,plev_dice,t_dice,qv_dice,u_dice,v_dice,o3_dice &
      &     ,shf_dice,lhf_dice,lwup_dice,swup_dice,tg_dice,ustar_dice& 
      &     ,psurf_dice,ug_dice,vg_dice,ht_dice,hq_dice              &
@@ -376 +376 @@
 !champs initiaux:
       do k=1,nlev_dice
-         th_dicei(k)=th_dice(k)
+         t_dicei(k)=t_dice(k)
          qv_dicei(k)=qv_dice(k)
          u_dicei(k)=u_dice(k)
@@ -405 +405 @@
 
       CALL interp_dice_vertical(play,nlev_dice,nt_dice,plev_dice       &
-     &         ,th_dicei,qv_dicei,u_dicei,v_dicei,o3_dicei             &
+     &         ,t_dicei,qv_dicei,u_dicei,v_dicei,o3_dicei             &
      &         ,ht_dicei,hq_dicei,hu_dicei,hv_dicei,w_dicei,omega_dicei&
-     &         ,th_mod,qv_mod,u_mod,v_mod,o3_mod                       &
+     &         ,t_mod,qv_mod,u_mod,v_mod,o3_mod                       &
      &         ,ht_mod,hq_mod,hu_mod,hv_mod,w_mod,omega_mod,mxcalc)
 
@@ -425 +425 @@
       do l = 1, llm
 ! Ligne du dessous ?? decommenter si on lit theta au lieu de temp
-       temp(l) = th_mod(l)*(play(l)/pzero)**rkappa 
-!      temp(l) = t_mod(l) 
+!      temp(l) = th_mod(l)*(play(l)/pzero)**rkappa 
+       temp(l) = t_mod(l) 
        q(l,1) = qv_mod(l)
        q(l,2) = 0.0
@@ -476 +476 @@
 !---------------------------------------------------------------------
 
-!!!! Si la temperature de surface n'est pas imposée:
+!!!! Si la temperature de surface n'est pas impos??e:
   
       if (forcing_gabls4) then
@@ -531 +531 @@
       write(*,*) 'SST initiale mxcalc: ',tsurf,mxcalc
       do l = 1, llm
-! Ligne du dessous à decommenter si on lit theta au lieu de temp
+! Ligne du dessous ?? decommenter si on lit theta au lieu de temp
 !      temp(l) = th_mod(l)*(play(l)/pzero)**rkappa 
        temp(l) = t_mod(l) 
@@ -909 +909 @@
       endif !forcing_case
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+