Changeset 2393 for LMDZ5/trunk/libf/phylmd/dyn1d

Timestamp:

Nov 18, 2015, 12:25:20 PM (9 years ago)

Author:

jyg

Message:

Add various intializations of arrays in lmdz1d.F90
and in the convection scheme. Add output variables
for boundary layer splitting.

File:

• LMDZ5/trunk/libf/phylmd/dyn1d/lmdz1d.F90 (modified) (9 diffs)

LMDZ5/trunk/libf/phylmd/dyn1d/lmdz1d.F90

@@ -8 +8 @@
       PROGRAM lmdz1d
 
-      USE ioipsl, only: ju2ymds, ymds2ju, ioconf_calendar
-      use phys_state_var_mod
-      use dimphy
-      use surface_data, only : type_ocean,ok_veget
-      use pbl_surface_mod, only : ftsoil, pbl_surface_init,                     &
+   USE ioipsl, only: ju2ymds, ymds2ju, ioconf_calendar
+   USE phys_state_var_mod, ONLY : phys_state_var_init, phys_state_var_end, &
+       clwcon, detr_therm, &
+       qsol, fevap, z0m, z0h, agesno, &
+       du_gwd_rando, du_gwd_front, entr_therm, f0, fm_therm, &
+       falb_dir, falb_dif, &
+       ftsol, pbl_tke, pctsrf, radsol, rain_fall, snow_fall, ratqs, &
+       rlat, rlon, rnebcon, rugoro, sig1, w01, solaire_etat0, sollw, sollwdown, &
+       solsw, t_ancien, q_ancien, u_ancien, v_ancien, wake_cstar, wake_deltaq, &
+       wake_deltat, wake_delta_pbl_TKE, delta_tsurf, wake_fip, wake_pe, &
+       wake_s, zgam, &
+       zmax0, zmea, zpic, zsig, &
+       zstd, zthe, zval, ale_bl, ale_bl_trig, alp_bl
+   use dimphy
+   use surface_data, only : type_ocean,ok_veget
+   use pbl_surface_mod, only : ftsoil, pbl_surface_init,                     &
      &                            pbl_surface_final
       use fonte_neige_mod, only : fonte_neige_init, fonte_neige_final 
 
-      use infotrac ! new
-      use control_mod
-      USE indice_sol_mod
-      USE phyaqua_mod
-      USE mod_1D_cases_read
-      USE mod_1D_amma_read
-      USE print_control_mod, ONLY: prt_level
-      USE iniphysiq_mod, ONLY: iniphysiq
-      USE mod_const_mpi, ONLY: comm_lmdz
+   use infotrac ! new
+   use control_mod
+   USE indice_sol_mod
+   USE phyaqua_mod
+   USE mod_1D_cases_read
+   USE mod_1D_amma_read
+   USE print_control_mod, ONLY: prt_level
+   USE iniphysiq_mod, ONLY: iniphysiq
+   USE mod_const_mpi, ONLY: comm_lmdz
 
       implicit none
@@ -127 +138 @@
 !vertical advection computation
 !       real d_t_z(llm), d_q_z(llm)
-!       real d_t_dyn_z(llm), d_q_dyn_z(llm)
+!       real d_t_dyn_z(llm), dq_dyn_z(llm)
 !       real zz(llm)
 !       real zfact
@@ -516 +527 @@
 !      allocate(d_th_adv(llm)) 
 
+      q(:,:) = 0.
+      dq(:,:) = 0.
+      dq_dyn(:,:) = 0.
+      d_q_adv(:,:) = 0.
+      d_q_nudge(:,:) = 0.
+
 !
 !   No ozone climatology need be read in this pre-initialization
@@ -730 +747 @@
         solsw=0.
         sollw=0.
+        sollwdown=rsigma*tsurf**4
         radsol=0.
         rnebcon=0.
         ratqs=0.
         clwcon=0.
+        zmax0 = 0.
         zmea=0.
         zstd=0.
@@ -742 +761 @@
         sig1=0.
         w01=0.
+        wake_cstar = 0.
+        wake_deltaq = 0.
+        wake_deltat = 0.
+        wake_delta_pbl_TKE = 0.
+        delta_tsurf = 0.
+        wake_fip = 0.
+        wake_pe = 0.
+        wake_s = 0.
+        ale_bl = 0.
+        ale_bl_trig = 0.
+        alp_bl = 0.
+        IF (ALLOCATED(du_gwd_rando)) du_gwd_rando = 0.
+        IF (ALLOCATED(du_gwd_front)) du_gwd_front = 0.
+        entr_therm = 0.
+        detr_therm = 0.
+        f0 = 0.
+        fm_therm = 0.
         u_ancien(1,:)=u(:)
         v_ancien(1,:)=v(:)
@@ -752 +788 @@
 ! pctsrf(:,is_sic),ftsol(:,nsrf),tsoil(:,isoil,nsrf),qsurf(:,nsrf)
 ! qsol,falb_dir(:,nsrf),falb_dif(:,nsrf),evap(:,nsrf),snow(:,nsrf)
-! radsol,solsw,sollw,fder,rain_fall,snow_fall,frugs(:,nsrf)
+! radsol,solsw,sollw, sollwdown,fder,rain_fall,snow_fall,frugs(:,nsrf)
 ! agesno(:,nsrf),zmea,zstd,zsig,zgam,zthe,zpic,zval,rugoro
 ! t_ancien,q_ancien,,frugs(:,is_oce),clwcon(:,1),rnebcon(:,1),ratqs(:,1)
@@ -939 +975 @@
          d_q_adv=0.0
        endif
-      print*, 'calcul de fcoriolis ', fcoriolis
+!      print*, 'calcul de fcoriolis ', fcoriolis
 
        if (forcing_toga .or. forcing_GCSSold .or. forcing_twpice            &
@@ -949 +985 @@
         endif
 
-      print*, 'fcoriolis ', fcoriolis, xlat,mxcalc
+      IF (prt_level >= 1) print*, 'fcoriolis, xlat,mxcalc ', &
+                                   fcoriolis, xlat,mxcalc
 
        du_age(1:mxcalc)=fcoriolis*(v(1:mxcalc)-vg(1:mxcalc))
        dv_age(1:mxcalc)=-fcoriolis*(u(1:mxcalc)-ug(1:mxcalc))
-       print *,'u-ug=',u-ug
+!       print *,'u-ug=',u-ug
 
 !!!!!!!!!!!!!!!!!!!!!!!!
@@ -960 +997 @@
        sfdt = sin(0.5*fcoriolis*timestep)
        cfdt = cos(0.5*fcoriolis*timestep)
-       print *,'fcoriolis,sfdt,cfdt,timestep',fcoriolis,sfdt,cfdt,timestep
+!       print *,'fcoriolis,sfdt,cfdt,timestep',fcoriolis,sfdt,cfdt,timestep
 !
         du_age(1:mxcalc)= -2.*sfdt/timestep*                                &