Ignore:
Timestamp:
May 30, 2019, 4:21:47 PM (5 years ago)
Author:
oboucher
Message:

Removing some WRITEs as some could print a lot of useless info in the output

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/IPSLCM6.0.15/libf/phylmd/StratAer/traccoag_mod.F90

    r3525 r3529  
    5959!----------------
    6060    REAL                                   :: m_aer_emiss_vol_daily ! daily injection mass emission
    61     REAL                                   :: sum_emi_so2         ! Test sum of all LON for budg_emi_so2
    6261    INTEGER                                :: it, k, i, ilon, ilev, itime, i_int, ieru
    6362    LOGICAL,DIMENSION(klon,klev)           :: is_strato           ! true = above tropopause, false = below
     
    8483       WRITE(lunout,*) 'in traccoag: date from phys_cal_mod =',year_cur,'-',mth_cur,'-',day_cur,'-',hour
    8584       WRITE(lunout,*) 'IN traccoag flag_sulf_emit: ',flag_sulf_emit
    86        IF (flag_sulf_emit == 1) THEN
    87           WRITE(lunout,*) 'IN traccoag nErupt: ',nErupt
    88           WRITE(lunout,*) 'IN traccoag injdur: ',injdur
    89           WRITE(lunout,*) 'IN traccoag : year_emit_vol',year_emit_vol
    90           WRITE(lunout,*) 'IN traccoag : mth_emit_vol',mth_emit_vol
    91           WRITE(lunout,*) 'IN traccoag : day_emit_vol',day_emit_vol
    92           WRITE(lunout,*) 'IN traccoag : m_aer_emiss_vol',m_aer_emiss_vol
    93           WRITE(lunout,*) 'IN traccoag : altemiss_vol',altemiss_vol
    94           WRITE(lunout,*) 'IN traccoag : sigma_alt_vol',sigma_alt_vol
    95           WRITE(lunout,*) 'IN traccoag : ponde_lonlat_vol',ponde_lonlat_vol
    96           WRITE(lunout,*) 'IN traccoag : xlat_min_vol',xlat_min_vol
    97           WRITE(lunout,*) 'IN traccoag : xlat_max_vol',xlat_max_vol
    98           WRITE(lunout,*) 'IN traccoag : xlon_min_vol',xlon_min_vol
    99           WRITE(lunout,*) 'IN traccoag : xlon_max_vol',xlon_max_vol
    100        ELSEIF (flag_sulf_emit == 2) THEN
    101           WRITE(lunout,*) 'IN traccoag : m_aer_emiss_sai',m_aer_emiss_sai
    102           WRITE(lunout,*) 'IN traccoag : altemiss_sai',altemiss_sai
    103           WRITE(lunout,*) 'IN traccoag : sigma_alt_sai',sigma_alt_sai
    104           WRITE(lunout,*) 'IN traccoag : xlat_sai',xlat_sai
    105           WRITE(lunout,*) 'IN traccoag : xlon_sai',xlon_sai
    106        ELSEIF (flag_sulf_emit == 3) THEN
    107           WRITE(lunout,*) 'IN traccoag : m_aer_emiss_sai',m_aer_emiss_sai
    108           WRITE(lunout,*) 'IN traccoag : altemiss_sai',altemiss_sai
    109           WRITE(lunout,*) 'IN traccoag : sigma_alt_sai',sigma_alt_sai
    110           WRITE(lunout,*) 'IN traccoag : xlat_min_sai',xlat_min_sai
    111           WRITE(lunout,*) 'IN traccoag : xlat_max_sai',xlat_max_sai
    112           WRITE(lunout,*) 'IN traccoag : xlon_sai',xlon_sai
    113        ENDIF
    114        WRITE(lunout,*) 'IN traccoag : flag_nuc_rate_box = ',flag_nuc_rate_box
    115        IF (flag_nuc_rate_box) THEN
    116           WRITE(lunout,*) 'IN traccoag : nuclat_min = ',nuclat_min,', nuclat_max = ',nuclat_max
    117           WRITE(lunout,*) 'IN traccoag : nucpres_min = ',nucpres_min,', nucpres_max = ',nucpres_max
    118        ENDIF
    11985    ENDIF
    12086   
     
    190156      ! stretch emission over one day of Pinatubo eruption
    191157       DO ieru=1, nErupt
    192           IF (is_mpi_root) THEN
    193              sum_emi_so2 = 0.0 ! Init sum
    194           ENDIF
    195158          IF (year_cur==year_emit_vol(ieru).AND.mth_cur==mth_emit_vol(ieru).AND.&
    196159               day_cur>=day_emit_vol(ieru).AND.day_cur<(day_emit_vol(ieru)+injdur)) THEN
     
    201164                  'ponde_lonlat_vol(ieru)=',ponde_lonlat_vol(ieru),'(injdur*ponde_lonlat_vol(ieru))', &
    202165                  (injdur*ponde_lonlat_vol(ieru)),'m_aer_emiss_vol_daily=',m_aer_emiss_vol_daily,'ieru=',ieru
    203              WRITE(lunout,*) 'IN traccoag, dlon=',dlon
    204166             DO i=1,klon
    205167                !Pinatubo eruption at 15.14N, 120.35E
    206168                dlat_loc=180./RPI/2.*(boundslat(i,1)-boundslat(i,3)) ! dlat = half difference of boundary latitudes
    207                 WRITE(lunout,*) 'IN traccoag, dlat=',dlat_loc
    208169                IF ( xlat(i).GE.xlat_min_vol(ieru)-dlat_loc .AND. xlat(i).LT.xlat_max_vol(ieru)+dlat_loc .AND. &
    209170                     xlon(i).GE.xlon_min_vol(ieru)-dlon .AND. xlon(i).LT.xlon_max_vol(ieru)+dlon ) THEN
     
    259220                      budg_emi_so2(i)=budg_emi_so2(i)+emission*zdm(k)*mSatom/mSO2mol
    260221                   ENDDO
    261                    sum_emi_so2 = sum_emi_so2 + budg_emi_so2(i) ! Sum all LON
    262222                ENDIF ! emission grid cell
    263223             ENDDO ! klon loop
    264              WRITE(lunout,*) "IN traccoag (ieru=",ieru,") global sum_emi_so2=",sum_emi_so2
    265224             WRITE(lunout,*) "IN traccoag (ieru=",ieru,") m_aer_emiss_vol_daily=",m_aer_emiss_vol_daily
    266225          ENDIF ! emission period
     
    274233!       SAI continuous emission o
    275234        dlat_loc=180./RPI/2.*(boundslat(i,1)-boundslat(i,3)) ! dlat = half difference of boundary latitudes
    276         WRITE(lunout,*) "IN traccoag, dlon=",dlon
    277         WRITE(lunout,*) "IN traccoag, dlat=",dlat_loc
    278235        IF  ( xlat(i).GE.xlat_sai-dlat_loc .AND. xlat(i).LT.xlat_sai+dlat_loc .AND. &
    279236          &   xlon(i).GE.xlon_sai-dlon .AND. xlon(i).LT.xlon_sai+dlon ) THEN
    280237!
    281           WRITE(lunout,*) 'coordinates of SAI point=',xlat(i), xlon(i), day_cur, mth_cur, year_cur
    282238!         compute altLMDz
    283239          altLMDz(:)=0.0
     
    346302!       SAI scenario with continuous emission
    347303        dlat_loc=180./RPI/2.*(boundslat(i,1)-boundslat(i,3)) ! dlat = half difference of boundary latitudes
    348         WRITE(lunout,*) 'IN traccoag, dlat = ',dlat_loc
    349304        theta_min = max(xlat(i)-dlat_loc,xlat_min_sai)
    350305        theta_max = min(xlat(i)+dlat_loc,xlat_max_sai)
Note: See TracChangeset for help on using the changeset viewer.