Ignore:
Timestamp:
Jul 20, 2024, 12:00:23 PM (2 months ago)
Author:
abarral
Message:

remove fixed-form \s+& remaining in .f90,.F90

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/stratemit.F90

    r5082 r5087  
    5454    IF (is_mpi_root) THEN
    5555       WRITE(*,*) 'IN STRATEMIT: date from phys_cal_mod=',year_cur,'-',&
    56                       &  mth_cur,'-',day_cur,'-',hour,' flh2o=',flh2o
     56    mth_cur,'-',day_cur,'-',hour,' flh2o=',flh2o
    5757    ENDIF
    5858
     
    7575               
    7676       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 ) THEN
     77     xlat(i)<latmax+dlat_loc .AND. &
     78     xlon(i)>=lonmin-dlon .AND. &
     79     xlon(i)<lonmax+dlon ) THEN
    8080                   !
    8181          WRITE(*,*) 'coordinates of volcanic injection point=',&
    82             &    xlat(i),xlon(i),day_cur,mth_cur,year_cur
     82      xlat(i),xlon(i),day_cur,mth_cur,year_cur
    8383          WRITE(*,*) 'DD m_emiss_vol_daily=', &
    84             &    m_emiss_vol_daily
     84      m_emiss_vol_daily
    8585         
    8686          !compute altLMDz
     
    9797          IF (flag_emit==3) then
    9898             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))
    100100          ELSE
    101101             theta=1.
     
    108108            ! stretch emission over stretchlong period
    109109            emission=m_emiss_vol_daily/m_air_gridbox(i,k)*f_lay_emiss(k)/stretchlong/ &
    110                  &     (86400.-pdt)*theta
     110       (86400.-pdt)*theta
    111111           
    112112             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= ', theta
     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= ', theta
    121121             
    122122             IF(emission < 1.E-34) emission = 0.0
     
    124124             IF (flh2o==0) THEN
    125125                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*pdtphys
     126    'i= ',i,'k= ',k, 'flh2o= ',flh2o, &
     127    tr_seri(i,k,id_spec), &
     128    tr_seri(i,k,id_spec)+emission*pdtphys
    129129             
    130130                tr_seri(i,k,id_spec)=tr_seri(i,k,id_spec)+emission*pdtphys
     
    141141               
    142142                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)
    146146               
    147147                IF(d_q_emiss(i,k) > 1.E34) THEN
Note: See TracChangeset for help on using the changeset viewer.