- Timestamp:
- Jul 20, 2024, 12:00:23 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/stratemit.F90
r5082 r5087 54 54 IF (is_mpi_root) THEN 55 55 WRITE(*,*) 'IN STRATEMIT: date from phys_cal_mod=',year_cur,'-',& 56 &mth_cur,'-',day_cur,'-',hour,' flh2o=',flh2o56 mth_cur,'-',day_cur,'-',hour,' flh2o=',flh2o 57 57 ENDIF 58 58 … … 75 75 76 76 IF ( xlat(i)>=latmin-dlat_loc .AND. & 77 &xlat(i)<latmax+dlat_loc .AND. &78 &xlon(i)>=lonmin-dlon .AND. &79 &xlon(i)<lonmax+dlon ) THEN77 xlat(i)<latmax+dlat_loc .AND. & 78 xlon(i)>=lonmin-dlon .AND. & 79 xlon(i)<lonmax+dlon ) THEN 80 80 ! 81 81 WRITE(*,*) 'coordinates of volcanic injection point=',& 82 &xlat(i),xlon(i),day_cur,mth_cur,year_cur82 xlat(i),xlon(i),day_cur,mth_cur,year_cur 83 83 WRITE(*,*) 'DD m_emiss_vol_daily=', & 84 &m_emiss_vol_daily84 m_emiss_vol_daily 85 85 86 86 !compute altLMDz … … 97 97 IF (flag_emit==3) then 98 98 theta=(sin(theta_max/180.*RPI)-sin(theta_min/180.*RPI))/ & 99 &(sin(xlat_max_sai/180.*RPI)-sin(xlat_min_sai/180.*RPI))99 (sin(xlat_max_sai/180.*RPI)-sin(xlat_min_sai/180.*RPI)) 100 100 ELSE 101 101 theta=1. … … 108 108 ! stretch emission over stretchlong period 109 109 emission=m_emiss_vol_daily/m_air_gridbox(i,k)*f_lay_emiss(k)/stretchlong/ & 110 &(86400.-pdt)*theta110 (86400.-pdt)*theta 111 111 112 112 IF(flag_verbose_strataer) WRITE(*,*) 'IN STRATEMIT: emission avant/apres', & 113 &'alt= ', altLMDz(k), &114 &'flh2o= ',flh2o, &115 &'id_speac= ',id_spec,id_species_total, &116 &'emission= ',emission, &117 &'pdtphys= ',pdtphys, &118 &'rapport m_emiss/m_air*f_lay= ', m_emiss_vol_daily/m_air_gridbox(i,k)*f_lay_emiss(k), &119 &'stretchlong= ', stretchlong, &120 &'theta= ', theta113 'alt= ', altLMDz(k), & 114 'flh2o= ',flh2o, & 115 'id_speac= ',id_spec,id_species_total, & 116 'emission= ',emission, & 117 'pdtphys= ',pdtphys, & 118 'rapport m_emiss/m_air*f_lay= ', m_emiss_vol_daily/m_air_gridbox(i,k)*f_lay_emiss(k), & 119 'stretchlong= ', stretchlong, & 120 'theta= ', theta 121 121 122 122 IF(emission < 1.E-34) emission = 0.0 … … 124 124 IF (flh2o==0) THEN 125 125 IF(flag_verbose_strataer) WRITE(*,*) 'IN STRATEMIT: tr_ser avant/apres',& 126 &'i= ',i,'k= ',k, 'flh2o= ',flh2o, &127 &tr_seri(i,k,id_spec), &128 &tr_seri(i,k,id_spec)+emission*pdtphys126 'i= ',i,'k= ',k, 'flh2o= ',flh2o, & 127 tr_seri(i,k,id_spec), & 128 tr_seri(i,k,id_spec)+emission*pdtphys 129 129 130 130 tr_seri(i,k,id_spec)=tr_seri(i,k,id_spec)+emission*pdtphys … … 141 141 142 142 IF(flag_verbose_strataer) WRITE(*,*) 'IN STRATEMIT: ',& 143 &'i= ',i,'k= ',k, 'flh2o= ',flh2o, &144 &'emission= ',emission, &145 &'d_q_emiss(i,k)= ',d_q_emiss(i,k)143 'i= ',i,'k= ',k, 'flh2o= ',flh2o, & 144 'emission= ',emission, & 145 'd_q_emiss(i,k)= ',d_q_emiss(i,k) 146 146 147 147 IF(d_q_emiss(i,k) > 1.E34) THEN
Note: See TracChangeset
for help on using the changeset viewer.