Changeset 1999 for LMDZ5/branches/testing/libf/phy1d/1D_decl_cases.h

Timestamp:

Mar 20, 2014, 10:57:19 AM (10 years ago)

Author:

Laurent Fairhead

Message:

Merged trunk changes r1920:1997 into testing branch

Location:

LMDZ5/branches/testing

Files:

• . (modified) (1 prop)
• libf/phy1d/1D_decl_cases.h (modified) (13 diffs)

LMDZ5/branches/testing/libf/phy1d/1D_decl_cases.h

@@ -16 +16 @@
         real    sec_print
 !!
-        integer nn
-        integer it_toga1, it_toga2
-        real time_toga1,time_toga2
-
         real ts_toga(nt_toga)
         real plev_toga(nlev_toga,nt_toga),w_toga(nlev_toga,nt_toga)
@@ -34 +30 @@
         real hq_prof(nlev_toga),vq_prof(nlev_toga)
 
-        real plev_mod(llm),w_mod(llm), t_mod(llm),q_mod(llm)
+        real w_mod(llm), t_mod(llm),q_mod(llm)
         real u_mod(llm),v_mod(llm), ht_mod(llm),vt_mod(llm)
         real hq_mod(llm),vq_mod(llm),qv_mod(llm),ql_mod(llm),qt_mod(llm)
@@ -87 +83 @@
         character*80 :: fich_amma
 ! Option du cas AMMA ou on impose la discretisation verticale (Ap,Bp)
-        logical  :: fixe_disvert=.true.
         integer nlev_amma, nt_amma
 !       parameter (nlev_amma=29, nt_amma=48)  ! Fleur, juillet 2012
@@ -104 +99 @@
 !profils initiaux:
         real plev_amma(nlev_amma)
-        real tv_amma(nlev_amma),rho_amma(nlev_amma)
-        real thv_amma(nlev_amma)
         
         real z_amma(nlev_amma)
@@ -111 +104 @@
         real u_amma(nlev_amma)
         real v_amma(nlev_amma)
-
-        real thvsurf_amma,tvsurf_amma,rhosurf_amma,thsurf
 
         real th_ammai(nlev_amma),q_ammai(nlev_amma)
@@ -130 +121 @@
 
 !champs interpoles
-        real plev_profamma(nlev_amma),vitw_profamma(nlev_amma)
+        real vitw_profamma(nlev_amma)
         real ht_profamma(nlev_amma)
         real hq_profamma(nlev_amma)
@@ -148 +139 @@
         integer year_ini_fire, day_ini_fire, mth_ini_fire
         real heure_ini_fire
-        real day_ju_ini_fire   ! Julian day of fire first day
         parameter (year_ini_fire=1987) 
         parameter (mth_ini_fire=7)
@@ -154 +144 @@
         parameter (heure_ini_fire=0.) !0h en secondes
 
-!profils initiaux:
-        real z_fire(nlev_fire)
-        real thl_fire(nlev_fire),qt_fire(nlev_fire)
-        real u_fire(nlev_fire), v_fire(nlev_fire)
-        real tke_fire(nlev_fire) 
-        
-!forcings
-        real ugeo_fire(nlev_fire),vgeo_fire(nlev_fire)
-        real wls_fire(nlev_fire),dqtdx_fire(nlev_fire)
-        real dqtdy_fire(nlev_fire)
-        real dqtdt_fire(nlev_fire),thl_rad_fire(nlev_fire)
-         
-        real ugeo_mod(llm),vgeo_mod(llm),wls_mod(llm)
-        real dqtdx_mod(llm),dqtdy_mod(llm),dqtdt_mod(llm)
-        real thl_rad_mod(llm)
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 ! Declarations specifiques au cas GCSSold
@@ -180 +155 @@
         real  Ts_gcssold
         real  dtime_frcg
-cAl1     logical :: imp_fcg_gcssold
-c        logical :: ts_fcg_gcssold 
-c        logical :: Tp_fcg_gcssold 
         logical :: Turb_fcg_gcssold
-        common /turb_forcing/ dtime_frcg,
-     $      Turb_fcg_gcssold, hthturb_gcssold, hqturb_gcssold
+
+        common /turb_forcing/
+     s  dtime_frcg,hthturb_gcssold, hqturb_gcssold,Turb_fcg_gcssold
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 ! Declarations specifiques au cas Arm_cu
@@ -206 +179 @@
         real adv_qt_armcu(nt_armcu)
         real theta_mod(llm),rv_mod(llm),play_mod(llm)
-        real d_t_dyn_ls(llm),d_q_dyn_ls(llm)
 ! profc comme "profil armcu"
-        real h_profc,play_profc,t_profc,th_profc,plev_profc
-        real u_profc,v_profc,qv_profc,rv_profc
         
 ! forcages interpoles dans le temps
@@ -229 +199 @@
         logical  :: trouve_700=.true.
         parameter (dt_sandu=6.*3600.)   ! forcages donnes ttes les 6 heures par ifa_sandu.txt
-!       parameter (tau_sandu=3600.)  ! temps de relaxation u,v,thetal,qt vers profil init et au dessus 700hPa
+        parameter (tau_sandu=3600.)  ! temps de relaxation u,v,thetal,qt vers profil init et au dessus 700hPa
 !!
-        integer it_sandu1, it_sandu2
-        real time_sandu1,time_sandu2
-
         real ts_sandu(nt_sandu)
 ! profs comme "profil sandu"
@@ -242 +209 @@
         real omega_profs(nlev_sandu),o3mmr_profs(nlev_sandu)
 
+        real, dimension(llm) :: relax_u,relax_v,relax_thl
+        real, dimension(llm,2) :: relax_q
+
         real thl_mod(llm),omega_mod(llm),o3mmr_mod(llm),tke_mod(llm)
-! pour relaxer u,v,thl et qt vers les profils initiaux au dessus de 700hPa
-        real relax_u(llm),relax_v(llm),relax_thl(llm),relax_q(llm,2)
 !vertical advection computation
         real d_t_z(llm), d_q_z(llm)
@@ -260 +228 @@
         parameter (mth_ini_astex=6)
         parameter (day_ini_astex=13)  ! 165 = 13 juin 1992
-        real dt_astex, tau_astex
+        real dt_astex
         parameter (dt_astex=3600.)    ! forcages donnes ttes les heures par ifa_astex.txt
-        integer it_astex1, it_astex2
-        real time_astex1,time_astex2
         real ts_astex(nt_astex),div_astex(nt_astex),ug_astex(nt_astex)
         real vg_astex(nt_astex),ufa_astex(nt_astex),vfa_astex(nt_astex)