- Timestamp:
- Jul 19, 2024, 5:41:58 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/stratemit.F90
r4601 r5082 74 74 theta_max = min(xlat(i)+dlat_loc,latmax) 75 75 76 IF ( xlat(i) .GE.latmin-dlat_loc .AND. &77 & xlat(i) .LT.latmax+dlat_loc .AND. &78 & xlon(i) .GE.lonmin-dlon .AND. &79 & xlon(i) .LT.lonmax+dlon ) THEN76 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 80 80 ! 81 81 WRITE(*,*) 'coordinates of volcanic injection point=',& … … 95 95 CALL STRATDISTRIB(altLMDz,altemiss,sigma_alt,f_lay_emiss) 96 96 97 IF (flag_emit .EQ.3) then97 IF (flag_emit==3) then 98 98 theta=(sin(theta_max/180.*RPI)-sin(theta_min/180.*RPI))/ & 99 99 & (sin(xlat_max_sai/180.*RPI)-sin(xlat_min_sai/180.*RPI)) … … 122 122 IF(emission < 1.E-34) emission = 0.0 123 123 124 IF (flh2o .EQ.0) THEN124 IF (flh2o==0) THEN 125 125 IF(flag_verbose_strataer) WRITE(*,*) 'IN STRATEMIT: tr_ser avant/apres',& 126 126 & 'i= ',i,'k= ',k, 'flh2o= ',flh2o, & … … 129 129 130 130 tr_seri(i,k,id_spec)=tr_seri(i,k,id_spec)+emission*pdtphys 131 IF (id_species_total .NE.0) THEN131 IF (id_species_total/=0) THEN 132 132 tr_seri(i,k,id_species_total)=tr_seri(i,k,id_species_total)+emission*pdtphys 133 133 ENDIF 134 ELSE IF(flh2o .EQ.1) THEN134 ELSE IF(flh2o==1) THEN 135 135 d_q_emiss(i,k)=emission*pdtphys 136 136 IF(d_q_emiss(i,k) > 1.E34) THEN
Note: See TracChangeset
for help on using the changeset viewer.