Index: /LMDZ6/trunk/libf/phylmd/StratAer/micphy_tstep.F90 =================================================================== --- /LMDZ6/trunk/libf/phylmd/StratAer/micphy_tstep.F90 (revision 3526) +++ /LMDZ6/trunk/libf/phylmd/StratAer/micphy_tstep.F90 (revision 3527) @@ -98,6 +98,6 @@ !NL - add nucleation box (if flag on) IF (flag_nuc_rate_box) THEN - IF (latitude_deg(ilon).LE.(nuclat_min) .OR. latitude_deg(ilon).GE.(nuclat_max) & - .OR. pplay(ilon,ilev).GE.nucpres_max .AND. pplay(ilon,ilev) .LE. nucpres_min ) THEN + IF (latitude_deg(ilon).LE.nuclat_min .OR. latitude_deg(ilon).GE.nuclat_max & + .OR. pplay(ilon,ilev).GE.nucpres_max .AND. pplay(ilon,ilev).LE.nucpres_min) THEN nucl_rate=0.0 ENDIF Index: /LMDZ6/trunk/libf/phylmd/StratAer/nucleation_tstep_mod.F90 =================================================================== --- /LMDZ6/trunk/libf/phylmd/StratAer/nucleation_tstep_mod.F90 (revision 3526) +++ /LMDZ6/trunk/libf/phylmd/StratAer/nucleation_tstep_mod.F90 (revision 3527) @@ -112,5 +112,4 @@ SUBROUTINE binapara(pt,prh,rhoa_in,jnuc,x,ntot,rc) - ! Fortran 90 subroutine binapara Index: /LMDZ6/trunk/libf/phylmd/StratAer/strataer_mod.F90 =================================================================== --- /LMDZ6/trunk/libf/phylmd/StratAer/strataer_mod.F90 (revision 3526) +++ /LMDZ6/trunk/libf/phylmd/StratAer/strataer_mod.F90 (revision 3527) @@ -109,8 +109,10 @@ ALLOCATE(xlat_max_vol(nErupt),xlon_max_vol(nErupt)) - year_emit_vol=0 ; mth_emit_vol=0 ; day_emit_vol=0 - m_aer_emiss_vol=0. ; altemiss_vol=0. ; sigma_alt_vol=0. - xlon_min_vol=0. ; xlon_max_vol=0. - xlat_min_vol=0. ; xlat_max_vol=0. + IF (nErupt.GT.0) THEN + year_emit_vol=0 ; mth_emit_vol=0 ; day_emit_vol=0 + m_aer_emiss_vol=0. ; altemiss_vol=0. ; sigma_alt_vol=0. + xlon_min_vol=0. ; xlon_max_vol=0. + xlat_min_vol=0. ; xlat_max_vol=0. + ENDIF CALL getin_p('year_emit_vol',year_emit_vol) @@ -140,45 +142,41 @@ CALL getin_p('nucpres_max',nucpres_max) - WRITE(lunout,*) 'IN STRATAER INIT2 year_emit_vol = ',year_emit_vol - WRITE(lunout,*) 'IN STRATAER INIT2 mth_emit_vol = ',mth_emit_vol - WRITE(lunout,*) 'IN STRATAER INIT2 day_emit_vol = ',day_emit_vol - !IF (is_master) THEN - WRITE(lunout,*) 'IN STRATAER INIT2 year_emit_vol = ',year_emit_vol - WRITE(lunout,*) 'IN STRATAER INIT2 mth_emit_vol=',mth_emit_vol - WRITE(lunout,*) 'IN STRATAER INIT2 day_emit_vol=',day_emit_vol - WRITE(lunout,*) 'IN STRATAER INIT2 =m_aer_emiss_vol',m_aer_emiss_vol - WRITE(lunout,*) 'IN STRATAER INIT2 =altemiss_vol',altemiss_vol - WRITE(lunout,*) 'IN STRATAER INIT2 =sigma_alt_vol',sigma_alt_vol - WRITE(lunout,*) 'IN STRATAER INIT2 xlon_min_vol=',xlon_min_vol - WRITE(lunout,*) 'IN STRATAER INIT2 xlon_max_vol=',xlon_max_vol - WRITE(lunout,*) 'IN STRATAER INIT2 xlat_min_vol=',xlat_min_vol - WRITE(lunout,*) 'IN STRATAER INIT2 xlat_max_vol=',xlat_max_vol - WRITE(lunout,*) 'flag_nuc_rate_box = ',flag_nuc_rate_box - WRITE(lunout,*) 'nuclat_min = ',nuclat_min - WRITE(lunout,*) 'nuclat_max = ',nuclat_max - WRITE(lunout,*) 'nucpres_min = ',nucpres_min - WRITE(lunout,*) 'nucpres_max = ',nucpres_max - WRITE(lunout,*) 'flag_sulf_emit = ',flag_sulf_emit - WRITE(lunout,*) 'injdur = ',injdur + WRITE(lunout,*) 'flag_sulf_emit = ',flag_sulf_emit + IF (flag_sulf_emit == 1) THEN + WRITE(lunout,*) 'IN STRATAER nErupt: ',nErupt + WRITE(lunout,*) 'IN STRATAER injdur: ',injdur + WRITE(lunout,*) 'IN STRATAER : year_emit_vol',year_emit_vol + WRITE(lunout,*) 'IN STRATAER : mth_emit_vol',mth_emit_vol + WRITE(lunout,*) 'IN STRATAER : day_emit_vol',day_emit_vol + WRITE(lunout,*) 'IN STRATAER : m_aer_emiss_vol',m_aer_emiss_vol + WRITE(lunout,*) 'IN STRATAER : altemiss_vol',altemiss_vol + WRITE(lunout,*) 'IN STRATAER : sigma_alt_vol',sigma_alt_vol + WRITE(lunout,*) 'IN STRATAER : ponde_lonlat_vol',ponde_lonlat_vol + WRITE(lunout,*) 'IN STRATAER : xlat_min_vol',xlat_min_vol + WRITE(lunout,*) 'IN STRATAER : xlat_max_vol',xlat_max_vol + WRITE(lunout,*) 'IN STRATAER : xlon_min_vol',xlon_min_vol + WRITE(lunout,*) 'IN STRATAER : xlon_max_vol',xlon_max_vol + ELSEIF (flag_sulf_emit == 2) THEN + WRITE(lunout,*) 'IN STRATAER : m_aer_emiss_sai',m_aer_emiss_sai + WRITE(lunout,*) 'IN STRATAER : altemiss_sai',altemiss_sai + WRITE(lunout,*) 'IN STRATAER : sigma_alt_sai',sigma_alt_sai + WRITE(lunout,*) 'IN STRATAER : xlat_sai',xlat_sai + WRITE(lunout,*) 'IN STRATAER : xlon_sai',xlon_sai WRITE(lunout,*) 'flag_sulf_emit_distrib = ',flag_sulf_emit_distrib - WRITE(lunout,*) 'nErupt = ',nErupt - WRITE(lunout,*) 'year_emit_vol = ',year_emit_vol - WRITE(lunout,*) 'mth_emit_vol = ',mth_emit_vol - WRITE(lunout,*) 'day_emit_vol = ',day_emit_vol - WRITE(lunout,*) 'm_aer_emiss_vol = ',m_aer_emiss_vol - WRITE(lunout,*) 'altemiss_vol = ',altemiss_vol - WRITE(lunout,*) 'sigma_alt_vol = ',sigma_alt_vol - WRITE(lunout,*) 'xlat_min_vol = ',xlat_min_vol - WRITE(lunout,*) 'xlat_max_vol = ',xlat_max_vol - WRITE(lunout,*) 'xlon_min_vol = ',xlon_min_vol - WRITE(lunout,*) 'xlon_max_vol = ',xlon_max_vol - WRITE(lunout,*) 'm_aer_emiss_sai = ',m_aer_emiss_sai - WRITE(lunout,*) 'altemiss_sai = ',altemiss_sai - WRITE(lunout,*) 'sigma_alt_sai = ',sigma_alt_sai - WRITE(lunout,*) 'xlat_sai = ',xlat_sai - WRITE(lunout,*) 'xlon_sai = ',xlon_sai - WRITE(lunout,*) 'xlat_min_sai = ',xlat_min_sai - WRITE(lunout,*) 'xlat_max_sai = ',xlat_max_sai + ELSEIF (flag_sulf_emit == 3) THEN + WRITE(lunout,*) 'IN STRATAER : m_aer_emiss_sai',m_aer_emiss_sai + WRITE(lunout,*) 'IN STRATAER : altemiss_sai',altemiss_sai + WRITE(lunout,*) 'IN STRATAER : sigma_alt_sai',sigma_alt_sai + WRITE(lunout,*) 'IN STRATAER : xlat_min_sai',xlat_min_sai + WRITE(lunout,*) 'IN STRATAER : xlat_max_sai',xlat_max_sai + WRITE(lunout,*) 'IN STRATAER : xlon_sai',xlon_sai + WRITE(lunout,*) 'flag_sulf_emit_distrib = ',flag_sulf_emit_distrib + ENDIF + WRITE(lunout,*) 'IN STRATAER : flag_nuc_rate_box = ',flag_nuc_rate_box + IF (flag_nuc_rate_box) THEN + WRITE(lunout,*) 'IN STRATAER : nuclat_min = ',nuclat_min,', nuclat_max = ',nuclat_max + WRITE(lunout,*) 'IN STRATAER : nucpres_min = ',nucpres_min,', nucpres_max = ',nucpres_max + ENDIF !ENDIF Index: /LMDZ6/trunk/libf/phylmd/StratAer/traccoag_mod.F90 =================================================================== --- /LMDZ6/trunk/libf/phylmd/StratAer/traccoag_mod.F90 (revision 3526) +++ /LMDZ6/trunk/libf/phylmd/StratAer/traccoag_mod.F90 (revision 3527) @@ -59,5 +59,4 @@ !---------------- REAL :: m_aer_emiss_vol_daily ! daily injection mass emission - REAL :: sum_emi_so2 ! Test sum of all LON for budg_emi_so2 INTEGER :: it, k, i, ilon, ilev, itime, i_int, ieru LOGICAL,DIMENSION(klon,klev) :: is_strato ! true = above tropopause, false = below @@ -84,37 +83,4 @@ WRITE(lunout,*) 'in traccoag: date from phys_cal_mod =',year_cur,'-',mth_cur,'-',day_cur,'-',hour WRITE(lunout,*) 'IN traccoag flag_sulf_emit: ',flag_sulf_emit - IF (flag_sulf_emit == 1) THEN - WRITE(lunout,*) 'IN traccoag nErupt: ',nErupt - WRITE(lunout,*) 'IN traccoag injdur: ',injdur - WRITE(lunout,*) 'IN traccoag : year_emit_vol',year_emit_vol - WRITE(lunout,*) 'IN traccoag : mth_emit_vol',mth_emit_vol - WRITE(lunout,*) 'IN traccoag : day_emit_vol',day_emit_vol - WRITE(lunout,*) 'IN traccoag : m_aer_emiss_vol',m_aer_emiss_vol - WRITE(lunout,*) 'IN traccoag : altemiss_vol',altemiss_vol - WRITE(lunout,*) 'IN traccoag : sigma_alt_vol',sigma_alt_vol - WRITE(lunout,*) 'IN traccoag : ponde_lonlat_vol',ponde_lonlat_vol - WRITE(lunout,*) 'IN traccoag : xlat_min_vol',xlat_min_vol - WRITE(lunout,*) 'IN traccoag : xlat_max_vol',xlat_max_vol - WRITE(lunout,*) 'IN traccoag : xlon_min_vol',xlon_min_vol - WRITE(lunout,*) 'IN traccoag : xlon_max_vol',xlon_max_vol - ELSEIF (flag_sulf_emit == 2) THEN - WRITE(lunout,*) 'IN traccoag : m_aer_emiss_sai',m_aer_emiss_sai - WRITE(lunout,*) 'IN traccoag : altemiss_sai',altemiss_sai - WRITE(lunout,*) 'IN traccoag : sigma_alt_sai',sigma_alt_sai - WRITE(lunout,*) 'IN traccoag : xlat_sai',xlat_sai - WRITE(lunout,*) 'IN traccoag : xlon_sai',xlon_sai - ELSEIF (flag_sulf_emit == 3) THEN - WRITE(lunout,*) 'IN traccoag : m_aer_emiss_sai',m_aer_emiss_sai - WRITE(lunout,*) 'IN traccoag : altemiss_sai',altemiss_sai - WRITE(lunout,*) 'IN traccoag : sigma_alt_sai',sigma_alt_sai - WRITE(lunout,*) 'IN traccoag : xlat_min_sai',xlat_min_sai - WRITE(lunout,*) 'IN traccoag : xlat_max_sai',xlat_max_sai - WRITE(lunout,*) 'IN traccoag : xlon_sai',xlon_sai - ENDIF - WRITE(lunout,*) 'IN traccoag : flag_nuc_rate_box = ',flag_nuc_rate_box - IF (flag_nuc_rate_box) THEN - WRITE(lunout,*) 'IN traccoag : nuclat_min = ',nuclat_min,', nuclat_max = ',nuclat_max - WRITE(lunout,*) 'IN traccoag : nucpres_min = ',nucpres_min,', nucpres_max = ',nucpres_max - ENDIF ENDIF @@ -190,7 +156,4 @@ ! stretch emission over one day of Pinatubo eruption DO ieru=1, nErupt - IF (is_mpi_root) THEN - sum_emi_so2 = 0.0 ! Init sum - ENDIF IF (year_cur==year_emit_vol(ieru).AND.mth_cur==mth_emit_vol(ieru).AND.& day_cur>=day_emit_vol(ieru).AND.day_cur<(day_emit_vol(ieru)+injdur)) THEN @@ -259,8 +222,6 @@ budg_emi_so2(i)=budg_emi_so2(i)+emission*zdm(k)*mSatom/mSO2mol ENDDO - sum_emi_so2 = sum_emi_so2 + budg_emi_so2(i) ! Sum all LON ENDIF ! emission grid cell ENDDO ! klon loop - WRITE(lunout,*) "IN traccoag (ieru=",ieru,") global sum_emi_so2=",sum_emi_so2 WRITE(lunout,*) "IN traccoag (ieru=",ieru,") m_aer_emiss_vol_daily=",m_aer_emiss_vol_daily ENDIF ! emission period @@ -323,5 +284,5 @@ DO k=1, klev ! stretch emission over whole year (360d) - emission=m_aer_emiss_sai*(mSO2mol/mSatom)/m_air_gridbox(i,k)*f_lay_emiss(k)/year_len/86400. + emission=m_aer_emiss_sai*(mSO2mol/mSatom)/m_air_gridbox(i,k)*f_lay_emiss(k)/FLOAT(year_len)/86400. tr_seri(i,k,id_SO2_strat)=tr_seri(i,k,id_SO2_strat)+emission*pdtphys budg_emi_so2(i)=budg_emi_so2(i)+emission*zdm(k)*mSatom/mSO2mol @@ -332,5 +293,5 @@ ! DO k=1, klev ! ! stretch emission over whole year (360d) -! emission=m_aer_emiss*(mH2SO4mol/mSatom)/m_part_dry(21)/m_air_gridbox(i,k)*f_lay_emiss(k)/year_len/86400 +! emission=m_aer_emiss*(mH2SO4mol/mSatom)/m_part_dry(21)/m_air_gridbox(i,k)*f_lay_emiss(k)/FLOAT(year_len)/86400. ! tr_seri(i,k,id_BIN01_strat+20)=tr_seri(i,k,id_BIN01_strat+20)+emission*pdtphys ! budg_emi_part(i)=budg_emi_part(i)+emission*zdm(k)*mSatom/mH2SO4mol @@ -346,5 +307,4 @@ ! SAI scenario with continuous emission dlat_loc=180./RPI/2.*(boundslat(i,1)-boundslat(i,3)) ! dlat = half difference of boundary latitudes - WRITE(lunout,*) 'IN traccoag, dlat = ',dlat_loc theta_min = max(xlat(i)-dlat_loc,xlat_min_sai) theta_max = min(xlat(i)+dlat_loc,xlat_max_sai) @@ -397,5 +357,5 @@ ! stretch emission over whole year (360d) emission=m_aer_emiss_sai*(mSO2mol/mSatom)/m_air_gridbox(i,k)*f_lay_emiss(k)/ & - & year_len/86400.*(sin(theta_max/180.*RPI)-sin(theta_min/180.*RPI))/ & + & FLOAT(year_len)/86400.*(sin(theta_max/180.*RPI)-sin(theta_min/180.*RPI))/ & & (sin(xlat_max_sai/180.*RPI)-sin(xlat_min_sai/180.*RPI)) tr_seri(i,k,id_SO2_strat)=tr_seri(i,k,id_SO2_strat)+emission*pdtphys