Changeset 183 in lmdz_wrf for trunk/WRFV3/lmdz/pbl_surface_mod.F90

Timestamp:

Oct 14, 2014, 7:40:28 PM (10 years ago)

Author:

lfita

Message:

Removing check pritings
Removing nullification of MixingRatioValues?(:,:,2) = 0 and MixingTendRatioValues?(:,:,2)
NO plul related to te absence of oliq at the initial conditions. If one introduce fake ones it works!

File:

• trunk/WRFV3/lmdz/pbl_surface_mod.F90 (modified) (24 diffs)

trunk/WRFV3/lmdz/pbl_surface_mod.F90

@@ -501 +501 @@
 !****************************************************************************************
 
-    PRINT *,'  Lluis in pbl_surface: ts ',lp,ts(lp,:)
-
     IF (first_call) THEN
        first_call=.FALSE.
@@ -606 +604 @@
        ENDDO
     ENDDO
-    PRINT *,'    Lluis after rugos: ts ',lp,ts(lp,:), ' rugos: ', rugos(lp,:),       &
-      ' alb1: ',alb1(550,:)
 ! Mean calculations of albedo
 !
@@ -628 +624 @@
        ENDDO
     ENDDO
-    PRINT *,'    Lluis after rugos: ts ',lp,ts(lp,:), ' rugos: ', rugos(lp,:),       &
-      ' alb1: ',alb1(550,:)
 ! We here suppose the fraction f1 of incoming radiance of visible radiance 
 ! as a fraction of all shortwave radiance 
@@ -652 +646 @@
        ENDDO
     ENDDO
-    PRINT *,'  Lluis ztsol ',lp,' : ',ztsol(lp)
 
 ! Linear distrubution on sub-surface of long- and shortwave net radiance
@@ -661 +654 @@
        ENDDO
     ENDDO
-    PRINT *,'  Lluis ',lp,' sollw: ',sollw(lp,:),' solsw: ',solsw(lp,:)
-
 
 ! Downwelling longwave radiation at mean surface
@@ -669 +660 @@
        lwdown_m(i) = sollw_m(i) + RSIGMA*ztsol(i)**4
     ENDDO
-    PRINT *,'  Lluis ',lp,' lwdown_m: ',lwdown_m(lp)
-
 !****************************************************************************************
 ! 4) Loop over different surfaces
@@ -683 +672 @@
        ni(:) = 0
        knon  = 0
-! Lluis
-       lpt = 1
+
        DO i = 1, klon
           IF (pctsrf(i,nsrf) > 0.) THEN
              knon = knon + 1
              ni(knon) = i
-! Lluis
-             IF (i == lp) lpt = knon
           ENDIF
        ENDDO
@@ -711 +697 @@
 !****************************************************************************************
 
-       PRINT *,'    Lluis nsrf j i rugos yts_________'
        DO j = 1, knon
           i = ni(j)
@@ -728 +713 @@
           ysollw(j)  = sollw(i,nsrf)
           yrugos(j)  = rugos(i,nsrf)
-          PRINT *,nsrf,j,i,rugos(i,nsrf),yts(i)
           yrugoro(j) = rugoro(i)
           yu1(j)     = u(i,1)
@@ -735 +719 @@
           ywindsp(j) = SQRT(u10m(i,nsrf)**2 + v10m(i,nsrf)**2 )
        END DO
-       PRINT *,'  Lluis ', lpt,' yts: ',yts(lpt)
 
        DO k = 1, klev
@@ -757 +740 @@
           END DO
        END DO
-       PRINT *,'  Lluis ', lpt,' ytsoil: ',ytsoil(lpt,:),' ftsoil: ',ftsoil(ni(lpt),:,nsrf)
 
        ! qsol(water height in soil) only for bucket continental model
@@ -766 +748 @@
           END DO
        ENDIF
-       PRINT *,'  Lluis ',lpt,' yqsol: ',yqsol(lpt)
        
 !****************************************************************************************
@@ -773 +754 @@
 !****************************************************************************************
 
-       PRINT *,'  Lluis before clcdrag: yrugos 1:', yrugos(1),' lpt: ',yrugos(lpt),  &
-         ' yts: ',yts(lpt), ' yalb: ',yalb(lpt)
        CALL clcdrag( knon, nsrf, ypaprs, ypplay, & 
             yu(:,1), yv(:,1), yt(:,1), yq(:,1), &
@@ -788 +767 @@
             ypaprs, ypplay, yu, yv, yq, yt, yts, yrugos, yqsurf, ycdragm, &
             ycoefm, ycoefh, ytke)
-      PRINT *,'  Lluis after coef_diff_turb lpt: ',lpt,' k ycoefm ycoefh ytke ____'
-      DO k=1,klev
-        PRINT *,k,ycoefm(lpt,k),ycoefh(lpt,k),ytke(lpt,k)
-      END DO
        IF (iflag_pbl>=20.AND.iflag_pbl<30) THEN
 ! In this case, coef_diff_turb is called for the Cd only
@@ -813 +788 @@
 !****************************************************************************************
 
-      PRINT *,'  Lluis before climb_hq_down knon: ',lpt,' dtime: ',dtime,      &  
-        ' AcoefH: ',AcoefH(lpt),' AcoefQ: ',AcoefQ(lpt),' BcoefH: ',BcoefH(lpt),&
-        ' BcoefQ: ',BcoefQ(lpt)
-      PRINT *,'  Lluis before k ycoefh ypaprs ypplay ydelp yt yq ____'
-      DO k=1,klev
-        PRINT *,k,ycoefh(lpt,k),ypaprs(lpt,k),ypplay(lpt,k),ydelp(lpt,k),     &
-          yt(lpt,k),yq(lpt,k)
-      END DO
 ! - Calculate the coefficients Ccoef_H, Ccoef_Q, Dcoef_H and Dcoef_Q 
        CALL climb_hq_down(knon, ycoefh, ypaprs, ypplay, &
@@ -830 +797 @@
             AcoefU, AcoefV, BcoefU, BcoefV)
       
-      PRINT *,'  Lluis after climb_hq_down knon: ',lpt,' ycoefh: ',ycoefh(lpt,1),&
-        ' ypaprs: ',ypaprs(lpt,1),' ypplay: ',ypplay(lpt,1),' ydelp: ',         &
-        ydelp(lpt,1),' yt: ',yt(lpt,1),' yq: ',yq(lpt,1),' dtime: ',dtime,      & 
-        ' AcoefH: ',AcoefH(lpt),' AcoefQ: ',AcoefQ(lpt),' BcoefH: ',BcoefH(lpt),&
-        ' BcoefQ: ',BcoefQ(lpt)
-
 !****************************************************************************************
 ! 9) Small calculations
@@ -890 +851 @@
              ylwdown(i)=lwdown_m(ni(i))
           END DO
-          PRINT *,'  Lluis before surf_land lpt=',lpt,' yts', yts(lpt),' yq: ',yq(lpt,1)
           CALL surf_land(itap, dtime, date0, jour, knon, ni,&
                rlon, rlat, &
@@ -918 +878 @@
           
        CASE(is_oce)
-          PRINT *,'  Lluis before surf_ocean lpt= ',lpt,' yts: ', yts(lpt),         &
-            ' pctsrf: ',pctsrf(lpt,:),' ytsurf_new: ',ytsurf_new(lpt)
           CALL surf_ocean(rlon, rlat, ysolsw, ysollw, yalb1, &
                yrugos, ywindsp, rmu0, yfder, yts, &
@@ -931 +889 @@
                ytsurf_new, y_dflux_t, y_dflux_q, slab_wfbils, &
                y_flux_u1, y_flux_v1)
-          PRINT *,'  Lluis after surf_ocean lpt= ',lpt,' yts: ', yts(lpt),          &
-            ' pctsrf: ',pctsrf(lpt,:),' ytsurf_new: ',ytsurf_new(lpt), 'ypct: ',    &
-            ypct(lpt)
-!          PRINT *,'     i yts  . pctsrf : ytsurf_new_________'
-!          DO i=1,knon
-!            PRINT *,'    ',i,yts(i), ' . ', pctsrf(i,:),' : ',ytsurf_new(i)
-!          END DO
           
        CASE(is_sic)
@@ -967 +918 @@
        y_d_ts(1:knon)   = ytsurf_new(1:knon) - yts(1:knon)
 
-       PRINT *,'  Lluis nsrf: ',nsrf,' lpt: ',lpt,' ytsurf_new: ',ytsurf_new(lpt),   &
-         ' yts: ',yts(lpt),' y_d_ts: ',y_d_ts(lpt)
- 
 !****************************************************************************************
 !
@@ -1114 +1062 @@
 
 !      print*,'Dans pbl OK4'
-    PRINT *,'  Lluis before diagnostics: lp= ',lp,' ts ',lp,ts(lp,:)
 
 !****************************************************************************************
@@ -1237 +1184 @@
     END DO loop_nbsrf
 
-    PRINT *,'  Lluis after loop soil types: lp= ',lp,' ts ',lp,ts(lp,:)
-
 !****************************************************************************************
 ! 16) Calculate the mean value over all sub-surfaces for som variables
@@ -1279 +1224 @@
 !   print*,'OK pbl 9'
     
-    PRINT *,'  Lluis adding tendencies: lp= ',lp,' nsrf ts d_ts_______'
     DO nsrf = 1, nbsrf
        DO i = 1, klon          
@@ -1351 +1295 @@
 !!$       ENDDO
 !!$    ENDDO
-
-    PRINT *,'  Lluis last ts: ',ts(lp,:)
 
     DO i = 1, klon