Changeset 2693 for LMDZ5/trunk/libf/phylmd
- Timestamp:
- Nov 1, 2016, 10:08:00 AM (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/rrtm/readaerosol_optic_rrtm.F90
r2644 r2693 20 20 21 21 USE infotrac_phy 22 USE YOMCST 22 23 23 24 IMPLICIT NONE 25 24 26 include "clesphys.h" 25 include "YOMCST.h"26 27 27 28 28 ! Input arguments … … 104 104 !--convert to ug m-3 unit for consistency with offline fields 105 105 ! 106 doi=1,nbtr107 select case(trim(solsym(i)))108 case("ASBCM")106 DO i=1,nbtr 107 SELECT CASE(trim(solsym(i))) 108 CASE ("ASBCM") 109 109 id_ASBCM = i 110 case("ASPOMM")110 CASE ("ASPOMM") 111 111 id_ASPOMM = i 112 case("ASSO4M")112 CASE ("ASSO4M") 113 113 id_ASSO4M = i 114 case("ASMSAM")114 CASE ("ASMSAM") 115 115 id_ASMSAM = i 116 case("CSSO4M")116 CASE ("CSSO4M") 117 117 id_CSSO4M = i 118 case("CSMSAM")118 CASE ("CSMSAM") 119 119 id_CSMSAM = i 120 case("SSSSM")120 CASE ("SSSSM") 121 121 id_SSSSM = i 122 case("CSSSM")122 CASE ("CSSSM") 123 123 id_CSSSM = i 124 case("ASSSM")124 CASE ("ASSSM") 125 125 id_ASSSM = i 126 case("CIDUSTM")126 CASE ("CIDUSTM") 127 127 id_CIDUSTM = i 128 case("AIBCM")128 CASE ("AIBCM") 129 129 id_AIBCM = i 130 case("AIPOMM")130 CASE ("AIPOMM") 131 131 id_AIPOMM = i 132 case("ASNO3M")132 CASE ("ASNO3M") 133 133 id_ASNO3M = i 134 case("CSNO3M")134 CASE ("CSNO3M") 135 135 id_CSNO3M = i 136 case("CINO3M")136 CASE ("CINO3M") 137 137 id_CINO3M = i 138 end select 139 enddo 140 138 END SELECT 139 ENDDO 141 140 142 141 bcsol(:,:) = tr_seri(:,:,id_ASBCM) *zrho(:,:)*1.e9 ! ASBCM … … 171 170 ! 172 171 ! Read and interpolate sulfate 173 IF ( flag_aerosol .EQ. 1 .OR. & 174 flag_aerosol .EQ. 6 ) THEN 172 IF ( flag_aerosol .EQ. 1 .OR. flag_aerosol .EQ. 6 ) THEN 175 173 176 174 CALL readaerosol_interp(id_ASSO4M_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, sulfacc, sulfacc_pi,loadso4) … … 178 176 sulfacc(:,:) = 0. ; sulfacc_pi(:,:) = 0. 179 177 loadso4=0. 180 END 178 ENDIF 181 179 182 180 ! Read and interpolate bcsol and bcins 183 IF ( flag_aerosol .EQ. 2 .OR. & 184 flag_aerosol .EQ. 6 ) THEN 181 IF ( flag_aerosol .EQ. 2 .OR. flag_aerosol .EQ. 6 ) THEN 185 182 186 183 ! Get bc aerosol distribution … … 192 189 bcins(:,:) = 0. ; bcins_pi(:,:) = 0. 193 190 loadbc=0. 194 END IF 195 191 ENDIF 196 192 197 193 ! Read and interpolate pomsol and pomins 198 IF ( flag_aerosol .EQ. 3 .OR. & 199 flag_aerosol .EQ. 6 ) THEN 194 IF ( flag_aerosol .EQ. 3 .OR. flag_aerosol .EQ. 6 ) THEN 200 195 201 196 CALL readaerosol_interp(id_ASPOMM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, pomsol, pomsol_pi, load_tmp3) … … 206 201 pomins(:,:) = 0. ; pomins_pi(:,:) = 0. 207 202 loadoa=0. 208 END IF 209 203 ENDIF 210 204 211 205 ! Read and interpolate csssm, ssssm, assssm 212 IF (flag_aerosol .EQ. 4 .OR. & 213 flag_aerosol .EQ. 6 ) THEN 206 IF (flag_aerosol .EQ. 4 .OR. flag_aerosol .EQ. 6 ) THEN 214 207 215 208 CALL readaerosol_interp(id_SSSSM_phy ,itap, pdtphys,rjourvrai, & … … 228 221 229 222 ! Read and interpolate cidustm 230 IF (flag_aerosol .EQ. 5 .OR. & 231 flag_aerosol .EQ. 6 ) THEN 223 IF (flag_aerosol .EQ. 5 .OR. flag_aerosol .EQ. 6 ) THEN 232 224 233 225 CALL readaerosol_interp(id_CIDUSTM_phy, itap, pdtphys, rjourvrai, debut, pplay, paprs, t_seri, cidust, cidust_pi, loaddust) … … 299 291 DO i = 1, klon 300 292 pdel(i,k) = paprs(i,k) - paprs (i,k+1) 301 END 302 END 293 ENDDO 294 ENDDO 303 295 304 296 !--new aerosol properties
Note: See TracChangeset
for help on using the changeset viewer.