Changeset 183 in lmdz_wrf


Ignore:
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!

Location:
trunk/WRFV3
Files:
8 edited

Legend:

Unmodified
Added
Removed

  • trunk/WRFV3/lmdz/calcul_fluxs_mod.F90

    r1 r183  
    9696    INTEGER                                              :: lpt
    9797    lpt = 550
    98 
    99     IF (nisurf == 3) THEN
    100       PRINT *,'  Lluis in calcul_fluxs________________________'
    101       PRINT *,'    knon: ',knon,' nisurf:' , nisurf, ' dtime: ',dtime
    102       PRINT *,'    tsurf: ',tsurf(lpt),' p1lay: ', p1lay(lpt),' cal: ',cal(lpt),' beta: ', beta(lpt)
    103       PRINT *,'    coef1lay: ', coef1lay(lpt),' ps: ', ps(lpt)
    104       PRINT *,'    p_rain: ',precip_rain(lpt),' p_snow: ', precip_snow(lpt),' snow: ', snow(lpt)
    105       PRINT *,'    qsurf: ', qsurf(lpt)
    106       PRINT *,'    radsol: ',radsol(lpt),' dif_grnd: ', dif_grnd(lpt),' t1lay: ', t1lay(lpt)
    107       PRINT *,'    q1lay: ', q1lay(lpt),' u1lay: ', u1lay(lpt),' v1lay: ', v1lay(lpt)
    108       PRINT *,'    ptAcoef: ', petAcoef(lpt),' peqAcoef: ', peqAcoef(lpt)
    109       PRINT *,'    pteBcoef: ', petBcoef(lpt),' peqBcoef: ', peqBcoef(lpt)
    110       PRINT *,'    tsurf_new: ', tsurf_new(lpt), ' evap: ',evap(lpt),' fluxlat: ',fluxlat(lpt)
    111       PRINT *,'    fluxsens: ', fluxsens(lpt),' dflux_s: ', dflux_s(lpt),' dflux_l: ', dflux_l(lpt)
    112     END IF
    11398
    11499    IF (check) WRITE(*,*)'Entree ', modname,' surface = ',nisurf

  • trunk/WRFV3/lmdz/climb_hq_mod.F90

    r1 r183  
    186186    Bcoef_Q_out = Bcoef_Q
    187187
    188     PRINT *,'    Lluis in climb_hq_down 1 k dplay dpaprs Kcoefhq gamaq gamah '// &
    189       'delp q Ccoef_Q Dcoef_Q Acoef_Q Bcoef_Q local_H ________'
    190 
    191     DO k=2,klev
    192       PRINT *,k,pplay(1,k-1)-pplay(1,k),paprs(1,k-1)-paprs(1,k),Kcoefhq(1,k), &
    193         gamaq(1,k),gamah(1,k),delp(1,k),q(1,k),Ccoef_Q(1,k),Dcoef_Q(1,k),     &
    194         Acoef_Q(1),Bcoef_Q(1), local_H(1,k)
    195     END DO
    196 
    197188  END SUBROUTINE climb_hq_down
    198189!

  • trunk/WRFV3/lmdz/coefkzmin.F90

    r7 r183  
    114114         pblhmin(ig)=0.07*ustar(ig)/max(abs(coriol(ig)),2.546e-5)
    115115      enddo
    116       PRINT *,'    Fin of coefkzmin k ustar zlev pblhmin qmin kap kn km _____'
    117116         
    118117      do k=2,klev
    … …  
    126125            kn(ig,k)=kmin
    127126            km(ig,k)=kmin
    128             IF (ig == 1) PRINT *,'    ',k,ustar(ig),zlev(ig,k),pblhmin(ig),&
    129               qmin,kap,kn(ig,k),km(ig,k)
    130127         enddo
    131128      enddo

  • trunk/WRFV3/lmdz/fisrtilp.F90

    r1 r183  
    139139  REAL zzz
    140140
     141! Lluis
     142  INTEGER lp
     143
     144
    141145  include "YOETHF.h"
    142146  include "FCTTRE.h"
    … …  
    147151  !ym
    148152  zdelq=0.0
     153
     154! Lluis
     155  lp = 449
     156  PRINT *,'  Lluis beginning fisrtilp 15,15: 485 ___________________'
     157  PRINT *,'    dtime:',dtime,' rain: ',rain(lp), ' snow: ', snow(lp),     &
     158   ' iflag_cldcon: ', iflag_cldcon
     159  PRINT *,'    k paprs paprs prfl psfl ____'
     160  DO i=1, klev+1
     161    PRINT *,'   ',i,paprs(lp,i), prfl(lp,i), psfl(lp,i)
     162  END DO
     163  PRINT *,'    k pplay t q ptconv ratqs d_t d_q d_ql rneb radliq '
     164  DO i=1, klev
     165    PRINT *,i, pplay(lp,i), t(lp,i), q(lp,i), ptconv(lp,i), ratqs(lp,i),   &
     166     d_t(lp,i), d_q(lp,i), d_ql(lp,i), rneb(lp,i), radliq(lp,i)
     167  END DO
     168  PRINT *,'    k pfrac_impa pfrac_nucl pfrac_1nucl frac_impa ____'
     169  DO i=1, klev
     170     PRINT *, i, pfrac_impa(lp,i), pfrac_nucl(lp,i), pfrac_1nucl(lp,i),    &
     171     frac_impa(lp,i)
     172  END DO
     173  PRINT *, 'k frac_nucl beta rhcl zqta fraca ztv zpspsk ztla zthl ____'
     174  DO i=1, klev
     175   PRINT *, i, frac_nucl(lp,i), beta(lp,i), rhcl(lp,i), zqta(lp,i),        &
     176     fraca(lp,i), ztv(lp,i), zpspsk(lp,i), ztla(lp,i), zthl(lp,i)
     177  END DO
    149178
    150179  if (prt_level>9)write(lunout,*)'NUAGES4 A. JAM'
    … …  
    240269  !IM : klevm1
    241270  !ym      klevm1=klev-1
     271! Lluis
     272!  PRINT *,'  RCPD: ', RCPD,' RVTMP2: ', RVTMP2
     273!  PRINT *,'  Lluis lp: ',lp,' k zt zq zmair zcpair zcpeau d_t zrfl rhcl ____'
     274  PRINT *,'  Lluis lp: ',lp,' k zt zq zmair zcpair zcpeau d_t zrfl rhcl ____'
    242275  DO k = klev, 1, -1
    243276     !
    … …  
    266299           !     C        WRITE (6,*) 'cppluie ', zt(i)-(t(i,k+1)+d_t(i,k+1))
    267300        ENDDO
     301!       PRINT *,k, zt(lp), zq(lp), zmair, zcpair, zcpeau, d_t(lp,k+1), &
     302!         zrfl(lp)
    268303     ENDIF
    269304     !
    … …  
    630665  endif
    631666
     667  PRINT *,'  Lluis end fisrtilp 15,15: 485 ___________________'
     668  PRINT *,'    dtime:',dtime,' rain: ',rain(lp), ' snow: ', snow(lp),     &
     669   ' iflag_cldcon: ', iflag_cldcon
     670  PRINT *,'    k paprs paprs prfl psfl ____'
     671  DO i=1, klev+1
     672    PRINT *,'   ',i,paprs(lp,i), prfl(lp,i), psfl(lp,i)
     673  END DO
     674  PRINT *,'    k pplay t q ptconv ratqs d_t d_q d_ql rneb radliq '
     675  DO i=1, klev
     676    PRINT *,i, pplay(lp,i), t(lp,i), q(lp,i), ptconv(lp,i), ratqs(lp,i),   &
     677     d_t(lp,i), d_q(lp,i), d_ql(lp,i), rneb(lp,i), radliq(lp,i)
     678  END DO
     679  PRINT *,'    k pfrac_impa pfrac_nucl pfrac_1nucl frac_impa ____'
     680  DO i=1, klev
     681     PRINT *, i, pfrac_impa(lp,i), pfrac_nucl(lp,i), pfrac_1nucl(lp,i),    &
     682     frac_impa(lp,i)
     683  END DO
     684  PRINT *, 'k frac_nucl beta rhcl zqta fraca ztv zpspsk ztla zthl ____'
     685  DO i=1, klev
     686   PRINT *, i, frac_nucl(lp,i), beta(lp,i), rhcl(lp,i), zqta(lp,i),        &
     687     fraca(lp,i), ztv(lp,i), zpspsk(lp,i), ztla(lp,i), zthl(lp,i)
     688  END DO
     689
    632690END SUBROUTINE fisrtilp

  • trunk/WRFV3/lmdz/ocean_forced_mod.F90

    r1 r183  
    100100    END DO
    101101
    102     PRINT *,'   Lluis pt tsurf_limit_check tsurf_lim_____'
    103     DO i=-10,10
    104       PRINT *,lp+i,tsurf_limit(lp+i),tsurf_lim(lpt+i)
    105     END DO
    106 
    107     PRINT *,'Lluis in ocean_surface_mod lp= ',lp,' tsurf_limit: ',tsurf_limit(lp),   &
    108       ' lpt: ',lpt,' tsurf_lim: ',tsurf_lim(lpt),' knindex: ',knindex(lpt)
    109 
    110102!****************************************************************************************
    111103! 2)
    … …  
    125117    v1_lay(:) = v1(:) - v0(:)
    126118
    127     PRINT *,'  Lluis before calcul_fluxs tsurf_lim:', tsurf_lim(lpt),                &
    128       ' tsurf_new: ', tsurf_new(lpt)
    129119! Calcul de tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l and qsurf
    130120    CALL calcul_fluxs(knon, is_oce, dtime, &
    … …  
    134124         AcoefH, AcoefQ, BcoefH, BcoefQ, &
    135125         tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l)
    136     PRINT *,'  Lluis after calcul_fluxs tsurf_lim: ', tsurf_lim(lp),' tsurf_new: ',  &
    137       tsurf_new(lpt)
    138     PRINT *,'    Lluis i knindex tsurf_limit tsurf_lim tsurf_new ___________'
    139     DO i=1,knon
    140       PRINT *,i,knindex(i),tsurf_limit(knindex(i)),tsurf_lim(i),tsurf_new(i)
    141     END DO
    142126
    143127! - Flux calculation at first modele level for U and V

  • trunk/WRFV3/lmdz/pbl_surface_mod.F90

    r1 r183  
    501501!****************************************************************************************
    502502
    503     PRINT *,'  Lluis in pbl_surface: ts ',lp,ts(lp,:)
    504 
    505503    IF (first_call) THEN
    506504       first_call=.FALSE.
    … …  
    606604       ENDDO
    607605    ENDDO
    608     PRINT *,'    Lluis after rugos: ts ',lp,ts(lp,:), ' rugos: ', rugos(lp,:),       &
    609       ' alb1: ',alb1(550,:)
    610606! Mean calculations of albedo
    611607!
    … …  
    628624       ENDDO
    629625    ENDDO
    630     PRINT *,'    Lluis after rugos: ts ',lp,ts(lp,:), ' rugos: ', rugos(lp,:),       &
    631       ' alb1: ',alb1(550,:)
    632626! We here suppose the fraction f1 of incoming radiance of visible radiance
    633627! as a fraction of all shortwave radiance
    … …  
    652646       ENDDO
    653647    ENDDO
    654     PRINT *,'  Lluis ztsol ',lp,' : ',ztsol(lp)
    655648
    656649! Linear distrubution on sub-surface of long- and shortwave net radiance
    … …  
    661654       ENDDO
    662655    ENDDO
    663     PRINT *,'  Lluis ',lp,' sollw: ',sollw(lp,:),' solsw: ',solsw(lp,:)
    664 
    665656
    666657! Downwelling longwave radiation at mean surface
    … …  
    669660       lwdown_m(i) = sollw_m(i) + RSIGMA*ztsol(i)**4
    670661    ENDDO
    671     PRINT *,'  Lluis ',lp,' lwdown_m: ',lwdown_m(lp)
    672 
    673662!****************************************************************************************
    674663! 4) Loop over different surfaces
    … …  
    683672       ni(:) = 0
    684673       knon  = 0
    685 ! Lluis
    686        lpt = 1
     674
    687675       DO i = 1, klon
    688676          IF (pctsrf(i,nsrf) > 0.) THEN
    689677             knon = knon + 1
    690678             ni(knon) = i
    691 ! Lluis
    692              IF (i == lp) lpt = knon
    693679          ENDIF
    694680       ENDDO
    … …  
    711697!****************************************************************************************
    712698
    713        PRINT *,'    Lluis nsrf j i rugos yts_________'
    714699       DO j = 1, knon
    715700          i = ni(j)
    … …  
    728713          ysollw(j)  = sollw(i,nsrf)
    729714          yrugos(j)  = rugos(i,nsrf)
    730           PRINT *,nsrf,j,i,rugos(i,nsrf),yts(i)
    731715          yrugoro(j) = rugoro(i)
    732716          yu1(j)     = u(i,1)
    … …  
    735719          ywindsp(j) = SQRT(u10m(i,nsrf)**2 + v10m(i,nsrf)**2 )
    736720       END DO
    737        PRINT *,'  Lluis ', lpt,' yts: ',yts(lpt)
    738721
    739722       DO k = 1, klev
    … …  
    757740          END DO
    758741       END DO
    759        PRINT *,'  Lluis ', lpt,' ytsoil: ',ytsoil(lpt,:),' ftsoil: ',ftsoil(ni(lpt),:,nsrf)
    760742
    761743       ! qsol(water height in soil) only for bucket continental model
    … …  
    766748          END DO
    767749       ENDIF
    768        PRINT *,'  Lluis ',lpt,' yqsol: ',yqsol(lpt)
    769750       
    770751!****************************************************************************************
    … …  
    773754!****************************************************************************************
    774755
    775        PRINT *,'  Lluis before clcdrag: yrugos 1:', yrugos(1),' lpt: ',yrugos(lpt),  &
    776          ' yts: ',yts(lpt), ' yalb: ',yalb(lpt)
    777756       CALL clcdrag( knon, nsrf, ypaprs, ypplay, &
    778757            yu(:,1), yv(:,1), yt(:,1), yq(:,1), &
    … …  
    788767            ypaprs, ypplay, yu, yv, yq, yt, yts, yrugos, yqsurf, ycdragm, &
    789768            ycoefm, ycoefh, ytke)
    790       PRINT *,'  Lluis after coef_diff_turb lpt: ',lpt,' k ycoefm ycoefh ytke ____'
    791       DO k=1,klev
    792         PRINT *,k,ycoefm(lpt,k),ycoefh(lpt,k),ytke(lpt,k)
    793       END DO
    794769       IF (iflag_pbl>=20.AND.iflag_pbl<30) THEN
    795770! In this case, coef_diff_turb is called for the Cd only
    … …  
    813788!****************************************************************************************
    814789
    815       PRINT *,'  Lluis before climb_hq_down knon: ',lpt,' dtime: ',dtime,      & 
    816         ' AcoefH: ',AcoefH(lpt),' AcoefQ: ',AcoefQ(lpt),' BcoefH: ',BcoefH(lpt),&
    817         ' BcoefQ: ',BcoefQ(lpt)
    818       PRINT *,'  Lluis before k ycoefh ypaprs ypplay ydelp yt yq ____'
    819       DO k=1,klev
    820         PRINT *,k,ycoefh(lpt,k),ypaprs(lpt,k),ypplay(lpt,k),ydelp(lpt,k),     &
    821           yt(lpt,k),yq(lpt,k)
    822       END DO
    823790! - Calculate the coefficients Ccoef_H, Ccoef_Q, Dcoef_H and Dcoef_Q
    824791       CALL climb_hq_down(knon, ycoefh, ypaprs, ypplay, &
    … …  
    830797            AcoefU, AcoefV, BcoefU, BcoefV)
    831798     
    832       PRINT *,'  Lluis after climb_hq_down knon: ',lpt,' ycoefh: ',ycoefh(lpt,1),&
    833         ' ypaprs: ',ypaprs(lpt,1),' ypplay: ',ypplay(lpt,1),' ydelp: ',         &
    834         ydelp(lpt,1),' yt: ',yt(lpt,1),' yq: ',yq(lpt,1),' dtime: ',dtime,      &
    835         ' AcoefH: ',AcoefH(lpt),' AcoefQ: ',AcoefQ(lpt),' BcoefH: ',BcoefH(lpt),&
    836         ' BcoefQ: ',BcoefQ(lpt)
    837 
    838799!****************************************************************************************
    839800! 9) Small calculations
    … …  
    890851             ylwdown(i)=lwdown_m(ni(i))
    891852          END DO
    892           PRINT *,'  Lluis before surf_land lpt=',lpt,' yts', yts(lpt),' yq: ',yq(lpt,1)
    893853          CALL surf_land(itap, dtime, date0, jour, knon, ni,&
    894854               rlon, rlat, &
    … …  
    918878         
    919879       CASE(is_oce)
    920           PRINT *,'  Lluis before surf_ocean lpt= ',lpt,' yts: ', yts(lpt),         &
    921             ' pctsrf: ',pctsrf(lpt,:),' ytsurf_new: ',ytsurf_new(lpt)
    922880          CALL surf_ocean(rlon, rlat, ysolsw, ysollw, yalb1, &
    923881               yrugos, ywindsp, rmu0, yfder, yts, &
    … …  
    931889               ytsurf_new, y_dflux_t, y_dflux_q, slab_wfbils, &
    932890               y_flux_u1, y_flux_v1)
    933           PRINT *,'  Lluis after surf_ocean lpt= ',lpt,' yts: ', yts(lpt),          &
    934             ' pctsrf: ',pctsrf(lpt,:),' ytsurf_new: ',ytsurf_new(lpt), 'ypct: ',    &
    935             ypct(lpt)
    936 !          PRINT *,'     i yts  . pctsrf : ytsurf_new_________'
    937 !          DO i=1,knon
    938 !            PRINT *,'    ',i,yts(i), ' . ', pctsrf(i,:),' : ',ytsurf_new(i)
    939 !          END DO
    940891         
    941892       CASE(is_sic)
    … …  
    967918       y_d_ts(1:knon)   = ytsurf_new(1:knon) - yts(1:knon)
    968919
    969        PRINT *,'  Lluis nsrf: ',nsrf,' lpt: ',lpt,' ytsurf_new: ',ytsurf_new(lpt),   &
    970          ' yts: ',yts(lpt),' y_d_ts: ',y_d_ts(lpt)
    971  
    972920!****************************************************************************************
    973921!
    … …  
    11141062
    11151063!      print*,'Dans pbl OK4'
    1116     PRINT *,'  Lluis before diagnostics: lp= ',lp,' ts ',lp,ts(lp,:)
    11171064
    11181065!****************************************************************************************
    … …  
    12371184    END DO loop_nbsrf
    12381185
    1239     PRINT *,'  Lluis after loop soil types: lp= ',lp,' ts ',lp,ts(lp,:)
    1240 
    12411186!****************************************************************************************
    12421187! 16) Calculate the mean value over all sub-surfaces for som variables
    … …  
    12791224!   print*,'OK pbl 9'
    12801225   
    1281     PRINT *,'  Lluis adding tendencies: lp= ',lp,' nsrf ts d_ts_______'
    12821226    DO nsrf = 1, nbsrf
    12831227       DO i = 1, klon         
    … …  
    13511295!!$       ENDDO
    13521296!!$    ENDDO
    1353 
    1354     PRINT *,'  Lluis last ts: ',ts(lp,:)
    13551297
    13561298    DO i = 1, klon

  • trunk/WRFV3/lmdz/physiq.F90

    r178 r183  
    13011301! L. Fita, LMD. Point for checkings...
    13021302      INTEGER                                            :: llp
    1303       llp = 644
     1303      llp = 449
    13041304
    13051305!c======================================================================
    … …  
    20922092!c Re-evaporer l'eau liquide nuageuse
    20932093!c
     2094      PRINT *,'    Lluis: lp: ',llp,' RLVTT: ', RLVTT,' RCPD: ', RCPD,  &
     2095        ' RVTMP2: ',RVTMP2, 'RTT: ', RTT
     2096      PRINT *,'    Lluis k qx_1 qx_2 ql_seri _____'
     2097      DO k=1,klev
     2098        PRINT *,k, qx(llp,k,:),ql_seri(llp,k)
     2099      END DO
     2100      PRINT *,'    Lluis re-evaporation: k q_seri zlvdcp '//         &
     2101        ' zlsdc zdelta t_seri za zb ql_seri_____'
    20942102      DO k = 1, klev  ! re-evaporation de l'eau liquide nuageuse
    20952103      DO i = 1, klon
    … …  
    21072115         d_q_eva(i,k) = zb
    21082116      ENDDO
     2117        PRINT *,k, q_Seri(llp,k), zlvdcp, zlsdcp, zdelta, t_seri(llp,k),  &
     2118          za, zb, ql_seri(llp,k)
    21092119      ENDDO
    21102120!IM
    … …  
    29862996      END IF
    29872997      PRINT *,'  Lluis 6 d_h_vcol: ',d_h_vcol,' d_h_vcol_phy: ',d_h_vcol_phy
     2998      PRINT *,'  Lluis before lp: ',llp,' iflag_ratqs: ', iflag_ratqs, &
     2999        ' iflag_con: ', iflag_con,' iflag_cldcon: ',iflag_cldcon,      &
     3000        ' ratqsbas: ', ratqsbas,' ratqshaut: ',ratqshaut,              &
     3001        ' tau_ratqs: ', tau_ratqs
     3002      PRINT *,'    ratqs ratqsc _______________'
     3003      DO k=1,klev
     3004        PRINT *, k, ratqs(llp,k), ratqsc(llp,k)
     3005      END DO
    29883006
    29893007!c-------------------------------------------------------------------------
    … …  
    29973015       &     ratqs,ratqsc)
    29983016
     3017      PRINT *,'  Lluis after lp: ',llp,' iflag_ratqs: ', iflag_ratqs,  &
     3018        ' iflag_con: ', iflag_con,' iflag_cldcon: ',iflag_cldcon,      &
     3019        ' ratqsbas: ', ratqsbas,' ratqshaut: ',ratqshaut,              &
     3020        ' tau_ratqs: ',tau_ratqs
     3021      PRINT *,'    ratqs ratqsc _______________'
     3022      DO k=1,klev
     3023        PRINT *, k, ratqs(llp,k), ratqsc(llp,k)
     3024      END DO
    29993025
    30003026!c
    … …  
    30053031       print *,' ->fisrtilp '
    30063032      ENDIF
     3033      PRINT *,'    Lluis before fisrtilp k t_seri q_seri ql_seri ___'
     3034      DO k=1,klev
     3035        PRINT *,k, t_seri(llp,k), q_seri(llp,k), ql_seri(llp,k)
     3036      END DO
    30073037!c-------------------------------------------------------------------------
    30083038      CALL fisrtilp(dtime,paprs,pplay,                                               &
    … …  
    30143044       &           prfl, psfl, rhcl,                                                 &
    30153045       &           zqasc, fraca,ztv,zpspsk,ztla,zthl,iflag_cldcon )
     3046      PRINT *,'    Lluis after fisrtilp k t_seri q_seri ql_seri ___'
     3047      DO k=1,klev
     3048        PRINT *,k, t_seri(llp,k), q_seri(llp,k), ql_seri(llp,k)
     3049      END DO
    30163050
    30173051      WHERE (rain_lsc < 0) rain_lsc = 0.

  • trunk/WRFV3/phys/module_lmdz_phys.F

    r179 r183  
    10941094
    10951095! Lluis Testing water species
    1096     MixingRatioValues(:,:,2)=0.
     1096!    MixingRatioValues(:,:,2)=MixingRatioValues(:,:,1)*0.2
    10971097
    10981098    IF (wrf_dbg >= 200) THEN
    … …  
    11071107!      lmdzp=ddimx*(jlmdz-1)+ilmdz
    11081108
    1109       PRINT *,'  WRF+LMDZ lmdzp: ',lmdzp,' : ',ilmdz,jlmdz
     1109      PRINT *,'  WRF+LMDZ before LMDZ physiq lmdzp: ',lmdzp,' : ',ilmdz,jlmdz
    11101110      PRINT *,'k wrf_T wrf_Temp TValues_______________'
    11111111      DO iz=1,dimz
    … …  
    11791179!!    MixingRatiotendValues(:,:,1)=1000.*MixingRatiotendValues(:,:,1)
    11801180!!    MixingRatiotendValues(:,:,1)=MixingRatiotendValues(:,:,1)+MixingRatiotendValues(:,:,2)
    1181     MixingRatiotendValues(:,:,2)=0.
     1181!    MixingRatiotendValues(:,:,2)=0.
    11821182
    11831183!!    PRINT *,'  Lluis: outgoing from physiq'
    … …  
    12581258
    12591259    IF (wrf_dbg >= 200) THEN
    1260       PRINT *,'  WRF+LMDZ lmdzp: ',lmdzp
     1260      PRINT *,'  WRF+LMDZ after LMDZ physiq lmdzp: ',lmdzp
    12611261      PRINT *,'k T tendencies _______________'
    12621262      DO iz=1,dimz
Note: See TracChangeset for help on using the changeset viewer.