Changeset 2683 for LMDZ5/trunk/libf/phylmd/dyn1d/1D_read_forc_cases.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/1D_read_forc_cases.h (modified) (7 diffs)

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
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+