Changeset 178 in lmdz_wrf for trunk/WRFV3/lmdz

Timestamp:

Oct 9, 2014, 5:25:55 PM (10 years ago)

Author:

lfita

Message:

Removing Fita initialization of variables in 'lmdz/thermcell_plume.F90'
Checking why ale/alep are 0 all the time. Why zw2 = 0, it seems correct! ?

Location:

trunk/WRFV3/lmdz

Files:

• physiq.F90 (modified) (4 diffs)
• thermcell_main.F90 (modified) (8 diffs)
• thermcell_plume.F90 (modified) (2 diffs)

trunk/WRFV3/lmdz/physiq.F90

@@ -2476 +2476 @@
 !c=================================================================================================
 
+       PRINT *,'  Lluis end calc itap: ', itap,' it_wape_prescr: ',&
+         it_wape_prescr,' k ale ale_wake alp alp_bl ale ale_bl _____'
+        DO i=1, klev
+         PRINT *, i, ale(i), ale_wake(i), alp(i),alp_bl(i), ale(i),  &
+           ale_bl(i)
+      END DO
 
 !c sb, oct02:
@@ -2791 +2797 @@
 !ccc fin nrlmd le 10/04/2012
 
+         PRINT *,'  Lluis before the itap: ', itap,' it_wape_prescr: ',&
+           it_wape_prescr,' k ale ale_wake alp alp_bl ale ale_bl _____'
+          DO i=1, klev
+           PRINT *, i, ale(i), ale_wake(i), alp(i),alp_bl(i), ale(i),  &
+             ale_bl(i)
+         END DO
+
          if (iflag_thermals>=1) then
          call calltherm(pdtphys                                                      &
@@ -2811 +2824 @@
 !ccc fin nrlmd le 10/04/2012                                                         &
        &      ,zqla,ztva )
+
+       PRINT *,'  Lluis after the itap: ', itap,' it_wape_prescr: ',&
+         it_wape_prescr,' k ale ale_wake alp alp_bl ale ale_bl _____'
+       DO i=1, klev
+         PRINT *, i, ale(i), ale_wake(i), alp(i),alp_bl(i), ale(i),  &
+          ale_bl(i)
+       END DO
 
 !ccc nrlmd le 10/04/2012
@@ -4285 +4305 @@
 
 ! Lluis
-  PRINT *,'  Lluis: ',klev,'  UBOUNDS: ',UBOUND(t_seri), UBOUND(u_seri),             &
-    UBOUND(d_q_con), UBOUND(d_t_con)
-  PRINT *,'  Lluis llp ',llp,' itap: ',itap,' zlev t_seri u_seri d_q_con d_t_con_____'
-  DO i=1,klev
-     PRINT *,i,t_seri(llp,i), u_seri(llp,i), d_q_con(llp,i), d_t_con(llp,i)
-  END DO 
+      PRINT *,'  Lluis: ',klev,'  UBOUNDS: ',UBOUND(t_seri),      &
+        UBOUND(u_seri), UBOUND(d_q_con), UBOUND(d_t_con)
+       PRINT *,'  Lluis llp ',llp,' itap: ',itap,' zlev t_seri ' // &
+         'u_seri d_q_con d_t_con_____'
+       DO i=1,klev
+         PRINT *,i,t_seri(llp,i), u_seri(llp,i), d_q_con(llp,i),    &
+           d_t_con(llp,i)
+       END DO 
 
 

trunk/WRFV3/lmdz/thermcell_main.F90

@@ -253 +253 @@
 ! L. Fita, LMD July 2014. Initializing variables. 
 !   Some not initializated according to values: iflag_trig_bl, iflag_clos_bl
-     zdthladj = 0.
-     pbl_tke_max0 = 0.
-     fraca0 = 0.
-     w_conv = 0.
-     w0 = 0.
-     therm_tke_max0 = 0.
-     env_tke_max0 = 0.
-     alp_bl_det = 0.
-     alp_bl_fluct_m = 0.
-     alp_bl_fluct_tke = 0.
-     alp_bl_conv = 0.
-     alp_bl_stat = 0.
-     interp = 0.
-     klcl = 0
+      PRINT *,'  Lluis debut thermcell_main i Alp_bl Ale_bl __' 
+      DO l=1, ngrid
+        PRINT *, l, Alp_bl(l), Ale_bl(l)
+      END DO
+
+!!     zdthladj = 0.
+!!     pbl_tke_max0 = 0.
+!!     fraca0 = 0.
+!!     w_conv = 0.
+!!     w0 = 0.
+!!     therm_tke_max0 = 0.
+!!     env_tke_max0 = 0.
+!!     alp_bl_det = 0.
+!!     alp_bl_fluct_m = 0.
+!!     alp_bl_fluct_tke = 0.
+!!     alp_bl_conv = 0.
+!!     alp_bl_stat = 0.
+!!     interp = 0.
+!!     klcl = 0
 
 !-----------------------------------------------------------------------
@@ -487 +492 @@
 !      print*,'THERM thermcell_main iflag_thermals_ed=',iflag_thermals_ed
       if (iflag_thermals_ed<=9) then
+         PRINT *,'  Lluis going to thermcell_plume!!!!'
 !         print*,'THERM NOUVELLE/NOUVELLE Arnaud'
          CALL thermcell_plume(itap,ngrid,nlay,ptimestep,ztv,zthl,po,zl,rhobarz,&
@@ -495 +501 @@
 
       elseif (iflag_thermals_ed>9) then
+         PRINT *,'  Lluis going to thermcellV1_plume!!!!'
 !        print*,'THERM RIO et al 2010, version d Arnaud'
          CALL thermcellV1_plume(itap,ngrid,nlay,ptimestep,ztv,zthl,po,zl,rhobarz,&
@@ -503 +510 @@
 
       endif
+      PRINT *,'  Lluis in thermcell_main iflag_thermals_ed: ',       &
+        iflag_thermals_ed, ' ig=528 ',zw2(528,:)
 
       if (prt_level.ge.1) print*,'apres thermcell_plume ',lev_out
@@ -601 +610 @@
      &       detr,zqla,lev_out,lunout1,igout)
 !IM 060508    &       detr,zqla,zmax,lev_out,lunout,igout)
+      PRINT *,'  Lluis after thermcell_flux2 ig=528 zw2: ',zw2(528,:)
 
       if (prt_level.ge.1) print*,'thermcell_main apres thermcell_flux'
@@ -745 +755 @@
       if (prt_level.ge.10)write(lunout,*)                                &
     &     'WARNING thermcell_main wth2=0. si zw2 > 1.e-10'
+      PRINT *,'  Lluis in thermcell_main ig=528 l wth3 zf2 fraca zw2 ___' 
       do l=1,nlay
          do ig=1,ngrid
@@ -758 +769 @@
             wth3(ig,l)=zf2*(1-2.*fraca(ig,l))/(1-fraca(ig,l))  &
      &                *zw2(ig,l)*zw2(ig,l)*zw2(ig,l)
+            IF (ig == 528) PRINT *, l, wth3(ig,l), zf2,           &
+               fraca(ig,l), zw2(ig,l)
+!             PRINT *, ig, l, wth3(ig,l), zf2,           &
+!               fraca(ig,l), zw2(ig,l)
             q2(ig,l)=zf2*(zqta(ig,l)*1000.-po(ig,l)*1000.)**2
 !test: on calcul q2/po=ratqsc
@@ -1001 +1016 @@
       enddo
       enddo
+      PRINT *,'  Lluis in thermcell_main i Alp_bl Srhobarz Swth3 ' //&
+        'Ale_bl Szw2 ______', alp_bl_k
+      DO l=1, ngrid
+        PRINT *, l, Alp_bl(l), SUM(rhobarz(l,:)), SUM(wth3(l,:)),    &
+           Ale_bl(l), SUM(zw2(l,:)) 
+      END DO
 
 !test:calcul de la ponderation des couches pour KE

trunk/WRFV3/lmdz/thermcell_plume.F90

@@ -191 +191 @@
     endif
 enddo
+  PRINT *,'  Lluis in themcell_plume ig actvie l zw2 ztv zlev pphi ___'
+  DO ig=1,ngrid
+    PRINT *,ig, active(ig), 1, zw2(ig,1), ztv(ig,1), zlev(ig,1), pphi(ig,1), &
+      2, zw2(ig,2), ztv(ig,2), zlev(ig,2), pphi(ig,2)
+  END DO
 !
 !==============================================================================
@@ -645 +650 @@
     endif
 enddo
+  PRINT *,'  Lluis in themcellV1_plume ig actvie l zw2 ztv zlev pphi ___'
+  DO ig=1,ngrid
+    PRINT *,ig, active(ig), 1, zw2(ig,1), ztv(ig,1), zlev(ig,1), pphi(ig,1), &
+      2, zw2(ig,2), ztv(ig,2), zlev(ig,2), pphi(ig,2)
+  END DO
 !