Changeset 186 in lmdz_wrf


Ignore:
Timestamp:
Oct 16, 2014, 11:21:06 AM (10 years ago)
Author:
lfita
Message:

Removing checking printings

Location:
trunk/WRFV3/lmdz
Files:
13 edited

Legend:

Unmodified
Added
Removed

  • trunk/WRFV3/lmdz/cltracrn.F90

    r1 r186  
    9090!AA que le pondere par la fraction de nature de sol.
    9191
    92   PRINT *,'  Lluis in cltracrn: ftsol', ftsol(550,:)
    93   PRINT *,'    UBOUND tr: ',UBOUND(tr),' trs: ',UBOUND(trs)
    94   PRINT *,'    tr: ', tr(550,:),' trs: ', trs(550)
    95  
    9692  DO i = 1,klon
    9793     zts(i) = 0.

  • trunk/WRFV3/lmdz/fisrtilp.F90

    r183 r186  
    154154! Lluis
    155155  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
    178156
    179157  if (prt_level>9)write(lunout,*)'NUAGES4 A. JAM'
    … …  
    269247  !IM : klevm1
    270248  !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 ____'
    275249  DO k = klev, 1, -1
    276250     !
    … …  
    299273           !     C        WRITE (6,*) 'cppluie ', zt(i)-(t(i,k+1)+d_t(i,k+1))
    300274        ENDDO
    301 !       PRINT *,k, zt(lp), zq(lp), zmair, zcpair, zcpeau, d_t(lp,k+1), &
    302 !         zrfl(lp)
    303275     ENDIF
    304276     !
    … …  
    665637  endif
    666638
    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 
    690639END SUBROUTINE fisrtilp

  • trunk/WRFV3/lmdz/hgardfou.F90

    r1 r186  
    9090!c
    9191      DO nsrf = 1, nbsrf
    92          PRINT *,'  Lluis nsrf: ',nsrf
    9392         DO i = 1, klon
    9493            IF (pctsrf(i,nsrf) > 0.) THEN
    … …  
    137136       &      ,jadrs(i),nsrf,zt(jadrs(i)),rlon(jadrs(i)),rlat(jadrs(i))              &
    138137       &      ,pctsrf(jadrs(i),nsrf)
    139              PRINT *,'    Lluis: nsrf k isfc tsol zt pctsrf_____'
    140              DO k=1,nbsrf
    141                PRINT *,'    ',nsrf,k,jadrs(i),tsol(jadrs(i),k),zt(jadrs(i)),pctsrf(jadrs(i),k)
    142              END DO
    143138           ENDDO
    144139         ENDIF

  • trunk/WRFV3/lmdz/limit_read_mod.F90

    r1 r186  
    119119    INTEGER :: i
    120120
    121     PRINT *,'  Lluis: knon: ',knon,' sst_out: ',UBOUND(sst_out),' : ',UBOUND(sst), 'x', UBOUND(knindex)
    122 
    123121    DO i = 1, knon
    124122       sst_out(i) = sst(knindex(i))

  • trunk/WRFV3/lmdz/phyetat0.F90

    r1 r186  
    8585
    8686      read_flag=.FALSE.
    87 
    88       PRINT *,'  Lluis in phyetat0 ftsol',ftsol(550,:),' read_flag: ',read_flag
    8987
    9088!c
    … …  
    11631161      END IF reading_file
    11641162
    1165       PRINT *,'  Lluis before init_iophy_new: ftsol: ',ftsol(550,:)
    1166 
    11671163! L. Fita, LMD. Called inside the non-used if
    11681164!!
    … …  
    11721168      CALL init_iophy_new(rlat,rlon)
    11731169       
    1174       PRINT *,'  Lluis after init_iophy_new: ftsol: ',ftsol(550,:)
    1175 
    11761170!c
    11771171!c Initialize module pbl_surface_mod
    … …  
    11911185!!      CALL fonte_neige_init(run_off_lic_0)
    11921186
    1193       PRINT *,'  Lluis leaving phyetat0 ftsol: ',ftsol(550,:)
    1194 
    11951187      RETURN
    11961188      END SUBROUTINE phyetat0

  • trunk/WRFV3/lmdz/physiq.F90

    r183 r186  
    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_____'
    21022094      DO k = 1, klev  ! re-evaporation de l'eau liquide nuageuse
    21032095      DO i = 1, klon
    … …  
    21152107         d_q_eva(i,k) = zb
    21162108      ENDDO
    2117         PRINT *,k, q_Seri(llp,k), zlvdcp, zlsdcp, zdelta, t_seri(llp,k),  &
    2118           za, zb, ql_seri(llp,k)
    21192109      ENDDO
    21202110!IM
    … …  
    21702160! non nul aux poles.
    21712161      IF (abs(solarlong0-1000.)<1.e-4) then
    2172          PRINT *,' Lluis calling zenang_an'
    21732162         call zenang_an(cycle_diurne,jH_cur,rlat,rlon,rmu0,fract)
    21742163      ELSE
    … …  
    22142203
    22152204      if (iflag_pbl/=0) then
    2216 
    2217       PRINT *,'  Lluis before pbl_surface qsol: ',qsol(llp),    &
    2218         ' rmu0: ',rmu0(llp),' jH_cur: ',jH_cur,' evap: ',evap(llp)
    22192205
    22202206      CALL pbl_surface(                                                              &
    … …  
    22412227       &     dsens,     devap,     zxsnow,                                           &
    22422228       &     zxfluxt,   zxfluxq,   q2m,     fluxq, pbl_tke )
    2243 
    2244       PRINT *,'  Lluis after pbl_surface qsol: ',qsol(llp),    &
    2245         ' rmu0: ',rmu0(llp),' jH_cur: ',jH_cur,' evap: ',evap(llp)
    22462229
    22472230!-----------------------------------------------------------------------------------------
    … …  
    24852468!cfin calcul ale et alp
    24862469!c=================================================================================================
    2487 
    2488        PRINT *,'  Lluis end calc itap: ', itap,' it_wape_prescr: ',&
    2489          it_wape_prescr,' k ale ale_wake alp alp_bl ale ale_bl _____'
    2490         DO i=1, klev
    2491          PRINT *, i, ale(i), ale_wake(i), alp(i),alp_bl(i), ale(i),  &
    2492            ale_bl(i)
    2493       END DO
    24942470
    24952471!c sb, oct02:
    … …  
    28072783!ccc fin nrlmd le 10/04/2012
    28082784
    2809          PRINT *,'  Lluis before the itap: ', itap,' it_wape_prescr: ',&
    2810            it_wape_prescr,' k ale ale_wake alp alp_bl ale ale_bl _____'
    2811           DO i=1, klev
    2812            PRINT *, i, ale(i), ale_wake(i), alp(i),alp_bl(i), ale(i),  &
    2813              ale_bl(i)
    2814          END DO
    2815 
    28162785         if (iflag_thermals>=1) then
    28172786         call calltherm(pdtphys                                                      &
    … …  
    28352804       &      ,zqla,ztva )
    28362805
    2837        PRINT *,'  Lluis after the itap: ', itap,' it_wape_prescr: ',&
    2838          it_wape_prescr,' k ale ale_wake alp alp_bl ale ale_bl _____'
    2839        DO i=1, klev
    2840          PRINT *, i, ale(i), ale_wake(i), alp(i),alp_bl(i), ale(i),  &
    2841           ale_bl(i)
    2842        END DO
    2843 
    28442806!ccc nrlmd le 10/04/2012
    28452807!c-----------Stochastic triggering-----------
    … …  
    29962958      END IF
    29972959      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
    30062960
    30072961!c-------------------------------------------------------------------------
    … …  
    30152969       &     ratqs,ratqsc)
    30162970
    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
    3025 
    30262971!c
    30272972!c Appeler le processus de condensation a grande echelle
    … …  
    30312976       print *,' ->fisrtilp '
    30322977      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
    30372978!c-------------------------------------------------------------------------
    30382979      CALL fisrtilp(dtime,paprs,pplay,                                               &
    … …  
    30442985       &           prfl, psfl, rhcl,                                                 &
    30452986       &           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
    30502987
    30512988      WHERE (rain_lsc < 0) rain_lsc = 0.
    … …  
    35623499      ENDIF
    35633500!c
    3564       PRINT *,'  Lluis before calling radlwsw rmu0: ',rmu0(llp),   &
    3565         ' dist: ',dist,' frac: ',fract(llp),' jH_cur: ',jH_cur
    35663501         CALL radlwsw                                                                &
    35673502       &        (dist, rmu0, fract,                                                  &
    … …  
    35863521       &        solsw_aero, solsw0_aero,                                             &
    35873522       &        topswcf_aero, solswcf_aero)
    3588       PRINT *,'  Lluis after calling radlwsw rmu0: ',rmu0(llp),   &
    3589         ' dist: ',dist,' frac: ',fract(llp),' jH_cur: ',jH_cur
    35903523         
    35913524!c
    … …  
    40153948      END IF
    40163949      PRINT *,'  Lluis 11 d_h_vcol: ',d_h_vcol,' d_h_vcol_phy: ',d_h_vcol_phy
    4017       PRINT *,'Lluis Reaching the SORTIES point'
    40183950
    40193951!C
    … …  
    42344166      endif
    42354167
    4236       PRINT *,'  Lluis WRITING outputs! itap: ',itap,' itau_phy: ',itau_phy,         &
    4237         ' itau_w: ',itau_w
    4238       PRINT *,'  Lluis writting outputs qsol: ',qsol(llp),    &
    4239         ' ftsol: ',ftsol(llp,:)
    4240 
    42414168#include "phys_output_write_new.h"
    42424169
    … …  
    42454172#include "write_histISCCP.h"
    42464173#endif
    4247 
    4248       PRINT *,'  Lluis WRITING histfiles!'
    42494174
    42504175#ifdef histNMC
    … …  
    43384263!      first=.false.
    43394264
    4340 ! Lluis
    4341       PRINT *,'  Lluis: ',klev,'  UBOUNDS: ',UBOUND(t_seri),      &
    4342         UBOUND(u_seri), UBOUND(d_q_con), UBOUND(d_t_con)
    4343        PRINT *,'  Lluis llp ',llp,' itap: ',itap,' zlev t_seri ' // &
    4344          'u_seri d_q_con d_t_con_____'
    4345        DO i=1,klev
    4346          PRINT *,i,t_seri(llp,i), u_seri(llp,i), d_q_con(llp,i),    &
    4347            d_t_con(llp,i)
    4348        END DO
    4349 
    4350 
    43514265      RETURN
    43524266      END SUBROUTINE physiq

  • trunk/WRFV3/lmdz/physiq_limit_variables_mod.F90

    r1 r186  
    8484    ALLOCATE(ftsoil_rst(klon, nsoilmx, nbsrf))
    8585
    86     PRINT *,'  Lluis in var_limit_init 315 qsurf_rst: ',qsurf_rst(315,:),   &
    87       'qsol_rst: ',qsol_rst(315)
    88 
    8986    RETURN
    9087

  • trunk/WRFV3/lmdz/phytrac.F90

    r1 r186  
    266266  INTEGER,SAVE  :: iflag_lscav_omp,iflag_lscav
    267267  LOGICAL,SAVE  :: convscav_omp,convscav
    268 
    269   PRINT *,'  Lluis in phytrac type_trac: ',type_trac
    270268
    271269!$OMP THREADPRIVATE(iflag_lscav,convscav)
    … …  
    374372     CASE('lmdz')
    375373        CALL traclmdz_init(pctsrf, xlat, xlon, ftsol, tr_seri, t_seri, pplay, sh, pdtphys, aerosol, lessivage)
    376         PRINT *,'  Lluis ALLOCATED trs: ',ALLOCATED(trs)
    377374
    378375     CASE('inca')

  • trunk/WRFV3/lmdz/suphel.F90

    r1 r186  
    229229      R5IES=R3IES*(RTT-R4IES)
    230230
    231       PRINT *,'Lluis exiting suphel!!!!'
    232 
    233231!C
    234232      RETURN

  • trunk/WRFV3/lmdz/surf_land_mod.F90

    r1 r186  
    149149!
    150150!****************************************************************************************
    151        PRINT *,'  Lluis before surf_land_bucket tsurf: ',tsurf(lpt), qsol(lpt)
    152151       CALL surf_land_bucket(itime, jour, knon, knindex, debut, dtime,&
    153152            tsurf, p1lay, cdragh, precip_rain, precip_snow, temp_air, &
    … …  
    157156            qsurf, z0_new, alb1_new, alb2_new, evap, &
    158157            fluxsens, fluxlat, tsurf_new, dflux_s, dflux_l)
    159        PRINT *,'  Lluis after surf_land_bucket tsurf: ',tsurf(lpt), qsol(lpt)
    160158
    161159    ENDIF ! ok_veget

  • trunk/WRFV3/lmdz/surf_ocean_mod.F90

    r1 r186  
    9393!****************************************************************************************
    9494
    95     PRINT *,'  Lluis in surf_ocean_mod lp: ',lp,' type_ocean: ',type_ocean,          &
    96       ' tsurf_in: ',tsurf_in(lp)
    9795
    9896!****************************************************************************************
    … …  
    133131    CASE('force')
    134132
    135        PRINT *,'  Lluis before ocean_forced_noice tsurf_new: ',tsurf_new(lp)
    136133       CALL ocean_forced_noice( &
    137134            itime, dtime, jour, knon, knindex, &
    … …  
    144141            qsurf, evap, fluxsens, fluxlat, flux_u1, flux_v1, &
    145142            tsurf_new, dflux_s, dflux_l)
    146        PRINT *,'  Lluis after ocean_forced_noice tsurf_new: ',tsurf_new(lp)
    147143    END SELECT
    148144

  • trunk/WRFV3/lmdz/traclmdz_mod.F90

    r1 r186  
    552552    END IF
    553553   
    554     PRINT *,'  Lluis in traclmdz_mod BOUNDS tr_seri: ',UBOUND(tr_seri),' trs: ',UBOUND(trs)
    555 
    556554    DO it=1, nbtr
    557555       IF (couchelimite .AND. pbl_flg(it) == 0 .AND. (it==id_rn .OR. it==id_pb)) THEN

  • trunk/WRFV3/lmdz/wrf_lmdz_mod.F90

    r185 r186  
    929929! Filling all 2D variables
    930930!!
    931     PRINT *,'  Lluis before get_lmdz: 5 11 ltksoil: ',pwrf_grid%ltksoil(5,:,11),     &
    932       ' ftsol: ', pftsol(315,:)
    933 
    934931    DO iy=1,ddy
    935932      DO ix=1,ddx
    … …  
    1002999    CALL vect_mat(pentr_therm, dx, dy, ddz, bdy, pwrf_grid%lentr_therm)
    10031000    CALL vect_mat(pdetr_therm, dx, dy, ddz, bdy, pwrf_grid%ldetr_therm)
    1004 
    1005     PRINT *,'  Lluis after get_lmdz: 5 11 ltksoil: ',pwrf_grid%ltksoil(5,:,11),      &
    1006       ' ftsol: ', pftsol(315,:)
    10071001
    10081002    RETURN
    … …  
    12971291        END DO
    12981292      END DO
    1299 !      PRINT *,k,'  Lluis vect: ',vect(lp,k),mat(il,k,jl),' lp: ',lp,' il jl :',il,jl
    13001293    END DO
    13011294
    … …  
    24432436    RD = 287.
    24442437    RMO3=47.9942
    2445 
    2446     PRINT *,'  Lluis in ' //TRIM(fname)//'!'
    24472438
    24482439! Remove variable assignation already done in get_lmdz? These variables are not used
Note: See TracChangeset for help on using the changeset viewer.