- Timestamp:
- Nov 28, 2014, 4:36:29 PM (10 years ago)
- Location:
- LMDZ5/branches/testing
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/branches/testing
- Property svn:mergeinfo changed
/LMDZ5/trunk merged: 2072,2075-2115,2117-2126,2128-2158
- Property svn:mergeinfo changed
-
LMDZ5/branches/testing/libf/phylmd/rrtm/aeropt_6bands_rrtm.F90
r2073 r2160 6 6 tau_allaer, piz_allaer, & 7 7 cg_allaer, m_allaer_pi, & 8 flag_aerosol, pplay, t_seri)8 flag_aerosol, zrho ) 9 9 10 10 USE dimphy 11 11 USE aero_mod 12 12 USE phys_local_var_mod, only: absvisaer 13 USE YOMCST , only: RD , RG14 13 15 14 ! Yves Balkanski le 12 avril 2006 … … 22 21 IMPLICIT NONE 23 22 23 INCLUDE "YOMCST.h" 24 24 INCLUDE "iniprint.h" 25 25 INCLUDE "clesphys.h" … … 30 30 REAL, DIMENSION(klon,klev), INTENT(in) :: pdel 31 31 REAL, INTENT(in) :: delt 32 REAL, DIMENSION(klon,klev,naero_ spc), INTENT(in) :: m_allaer33 REAL, DIMENSION(klon,klev,naero_ spc), INTENT(in) :: m_allaer_pi32 REAL, DIMENSION(klon,klev,naero_tot), INTENT(in) :: m_allaer 33 REAL, DIMENSION(klon,klev,naero_tot), INTENT(in) :: m_allaer_pi 34 34 REAL, DIMENSION(klon,klev), INTENT(in) :: RHcl ! humidite relative ciel clair 35 35 INTEGER, INTENT(in) :: flag_aerosol 36 REAL, DIMENSION(klon,klev), INTENT(in) :: pplay 37 REAL, DIMENSION(klon,klev), INTENT(in) :: t_seri 36 REAL, DIMENSION(klon,klev), INTENT(in) :: zrho 38 37 ! 39 38 ! Output arguments: … … 41 40 ! 2= natural aerosols 42 41 ! 43 REAL, DIMENSION(klon,klev,2,nbands_ rrtm), INTENT(out) :: tau_allaer ! epaisseur optique aerosol44 REAL, DIMENSION(klon,klev,2,nbands_ rrtm), INTENT(out) :: piz_allaer ! single scattering albedo aerosol45 REAL, DIMENSION(klon,klev,2,nbands_ rrtm), INTENT(out) :: cg_allaer ! asymmetry parameter aerosol42 REAL, DIMENSION(klon,klev,2,nbands_sw_rrtm), INTENT(out) :: tau_allaer ! epaisseur optique aerosol 43 REAL, DIMENSION(klon,klev,2,nbands_sw_rrtm), INTENT(out) :: piz_allaer ! single scattering albedo aerosol 44 REAL, DIMENSION(klon,klev,2,nbands_sw_rrtm), INTENT(out) :: cg_allaer ! asymmetry parameter aerosol 46 45 47 46 ! … … 67 66 REAL :: cg_ae2b_int ! Intermediate computation of Assymetry parameter 68 67 REAL :: Fact_RH(nbre_RH) 69 REAL :: zrho70 68 REAL :: fac 71 69 REAL :: zdp1(klon,klev) … … 73 71 INTEGER :: nb_aer 74 72 75 REAL, DIMENSION(klon,klev,naero_ spc) :: mass_temp76 REAL, DIMENSION(klon,klev,naero_ spc) :: mass_temp_pi77 REAL, DIMENSION(klon,klev,naero_tot,nbands_ rrtm) :: tau_ae78 REAL, DIMENSION(klon,klev,naero_tot,nbands_ rrtm) :: tau_ae_pi79 REAL, DIMENSION(klon,klev,naero_tot,nbands_ rrtm) :: piz_ae80 REAL, DIMENSION(klon,klev,naero_tot,nbands_ rrtm) :: cg_ae73 REAL, DIMENSION(klon,klev,naero_tot) :: mass_temp 74 REAL, DIMENSION(klon,klev,naero_tot) :: mass_temp_pi 75 REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) :: tau_ae 76 REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) :: tau_ae_pi 77 REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) :: piz_ae 78 REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) :: cg_ae 81 79 82 80 … … 84 82 ! Proprietes optiques 85 83 ! 86 REAL:: alpha_aers_6bands(nbre_RH,nbands_ rrtm,naero_soluble) !--unit m2/g SO487 REAL:: alpha_aeri_6bands(nbands_ rrtm,naero_insoluble)88 REAL:: cg_aers_6bands(nbre_RH,nbands_ rrtm,naero_soluble) !--unit89 REAL:: cg_aeri_6bands(nbands_ rrtm,naero_insoluble)90 REAL:: piz_aers_6bands(nbre_RH,nbands_ rrtm,naero_soluble) !-- unit91 REAL:: piz_aeri_6bands(nbands_ rrtm,naero_insoluble) !-- unit84 REAL:: alpha_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble) !--unit m2/g SO4 85 REAL:: alpha_aeri_6bands(nbands_sw_rrtm,naero_insoluble) 86 REAL:: cg_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble) !--unit 87 REAL:: cg_aeri_6bands(nbands_sw_rrtm,naero_insoluble) 88 REAL:: piz_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble) !-- unit 89 REAL:: piz_aeri_6bands(nbands_sw_rrtm,naero_insoluble) !-- unit 92 90 93 91 INTEGER :: id … … 280 278 spsol = 0 281 279 spinsol = 0 282 if(NSW.NE.nbands_rrtm) then280 IF (NSW.NE.nbands_sw_rrtm) THEN 283 281 print *,'Erreur NSW doit etre egal a 6 pour cette routine' 284 282 stop 285 endif283 ENDIF 286 284 287 285 DO k=1, klev 288 286 DO i=1, klon 289 zrho=pplay(i,k)/t_seri(i,k)/RD ! kg/m3 290 !CDIR UNROLL=naero_spc 291 mass_temp(i,k,:) = m_allaer(i,k,:) / zrho / 1.e+9 292 !CDIR UNROLL=naero_spc 293 mass_temp_pi(i,k,:) = m_allaer_pi(i,k,:) / zrho / 1.e+9 287 !CDIR UNROLL=naero_tot 288 mass_temp(i,k,:) = m_allaer(i,k,:) / zrho(i,k) / 1.e+9 !--kg/kg 289 !CDIR UNROLL=naero_tot 290 mass_temp_pi(i,k,:) = m_allaer_pi(i,k,:) / zrho(i,k) / 1.e+9 294 291 zdp1(i,k)=pdel(i,k)/(RG*delt) ! air mass auxiliary variable --> zdp1 [kg/(m^2 *s)] 295 292 ENDDO … … 299 296 nb_aer = 2 300 297 ALLOCATE (aerosol_name(nb_aer)) 301 aerosol_name(1) = id_ASSO4M 302 aerosol_name(2) = id_CSSO4M 298 aerosol_name(1) = id_ASSO4M_phy 299 aerosol_name(2) = id_CSSO4M_phy 303 300 ELSEIF (flag_aerosol .EQ. 2) THEN 304 301 nb_aer = 2 305 302 ALLOCATE (aerosol_name(nb_aer)) 306 aerosol_name(1) = id_ASBCM 307 aerosol_name(2) = id_AIBCM 303 aerosol_name(1) = id_ASBCM_phy 304 aerosol_name(2) = id_AIBCM_phy 308 305 ELSEIF (flag_aerosol .EQ. 3) THEN 309 306 nb_aer = 2 310 307 ALLOCATE (aerosol_name(nb_aer)) 311 aerosol_name(1) = id_ASPOMM 312 aerosol_name(2) = id_AIPOMM 308 aerosol_name(1) = id_ASPOMM_phy 309 aerosol_name(2) = id_AIPOMM_phy 313 310 ELSEIF (flag_aerosol .EQ. 4) THEN 314 311 nb_aer = 3 315 312 ALLOCATE (aerosol_name(nb_aer)) 316 aerosol_name(1) = id_CSSSM 317 aerosol_name(2) = id_SSSSM 318 aerosol_name(3) = id_ASSSM 313 aerosol_name(1) = id_CSSSM_phy 314 aerosol_name(2) = id_SSSSM_phy 315 aerosol_name(3) = id_ASSSM_phy 319 316 ELSEIF (flag_aerosol .EQ. 5) THEN 320 317 nb_aer = 1 321 318 ALLOCATE (aerosol_name(nb_aer)) 322 aerosol_name(1) = id_CIDUSTM 319 aerosol_name(1) = id_CIDUSTM_phy 323 320 ELSEIF (flag_aerosol .EQ. 6) THEN 324 321 nb_aer = 10 325 322 ALLOCATE (aerosol_name(nb_aer)) 326 aerosol_name(1) = id_ASSO4M 327 aerosol_name(2) = id_ASBCM 328 aerosol_name(3) = id_AIBCM 329 aerosol_name(4) = id_ASPOMM 330 aerosol_name(5) = id_AIPOMM 331 aerosol_name(6) = id_CSSSM 332 aerosol_name(7) = id_SSSSM 333 aerosol_name(8) = id_ASSSM 334 aerosol_name(9) = id_CIDUSTM 335 aerosol_name(10)= id_CSSO4M 323 aerosol_name(1) = id_ASSO4M_phy 324 aerosol_name(2) = id_ASBCM_phy 325 aerosol_name(3) = id_AIBCM_phy 326 aerosol_name(4) = id_ASPOMM_phy 327 aerosol_name(5) = id_AIPOMM_phy 328 aerosol_name(6) = id_CSSSM_phy 329 aerosol_name(7) = id_SSSSM_phy 330 aerosol_name(8) = id_ASSSM_phy 331 aerosol_name(9) = id_CIDUSTM_phy 332 aerosol_name(10)= id_CSSO4M_phy 336 333 ENDIF 337 334 … … 354 351 DO i=1, KLON 355 352 rh(i,k)=MIN(RHcl(i,k)*100.,RH_MAX) 356 RH_num(i,k) = INT( 353 RH_num(i,k) = INT(rh(i,k)/10. + 1.) 357 354 IF (rh(i,k).GT.85.) RH_num(i,k)=10 358 355 IF (rh(i,k).GT.90.) RH_num(i,k)=11 … … 365 362 DO m=1,nb_aer ! tau is only computed for each mass 366 363 fac=1.0 367 IF (aerosol_name(m).EQ.id_ASBCM ) THEN364 IF (aerosol_name(m).EQ.id_ASBCM_phy) THEN 368 365 soluble=.TRUE. 369 366 spsol=1 370 ELSEIF (aerosol_name(m).EQ.id_ASPOMM ) THEN367 ELSEIF (aerosol_name(m).EQ.id_ASPOMM_phy) THEN 371 368 soluble=.TRUE. 372 369 spsol=2 373 ELSEIF (aerosol_name(m).EQ.id_ASSO4M ) THEN370 ELSEIF (aerosol_name(m).EQ.id_ASSO4M_phy) THEN 374 371 soluble=.TRUE. 375 372 spsol=3 376 373 fac=1.375 ! (NH4)2-SO4/SO4 132/96 mass conversion factor for OD 377 ELSEIF (aerosol_name(m).EQ.id_CSSO4M ) THEN374 ELSEIF (aerosol_name(m).EQ.id_CSSO4M_phy) THEN 378 375 soluble=.TRUE. 379 376 spsol=4 380 377 fac=1.375 ! (NH4)2-SO4/SO4 132/96 mass conversion factor for OD 381 ELSEIF (aerosol_name(m).EQ.id_SSSSM ) THEN378 ELSEIF (aerosol_name(m).EQ.id_SSSSM_phy) THEN 382 379 soluble=.TRUE. 383 380 spsol=5 384 ELSEIF (aerosol_name(m).EQ.id_CSSSM ) THEN381 ELSEIF (aerosol_name(m).EQ.id_CSSSM_phy) THEN 385 382 soluble=.TRUE. 386 383 spsol=6 387 ELSEIF (aerosol_name(m).EQ.id_ASSSM ) THEN384 ELSEIF (aerosol_name(m).EQ.id_ASSSM_phy) THEN 388 385 soluble=.TRUE. 389 386 spsol=7 390 ELSEIF (aerosol_name(m).EQ.id_CIDUSTM ) THEN387 ELSEIF (aerosol_name(m).EQ.id_CIDUSTM_phy) THEN 391 388 soluble=.FALSE. 392 389 spinsol=1 393 ELSEIF (aerosol_name(m).EQ.id_AIBCM ) THEN390 ELSEIF (aerosol_name(m).EQ.id_AIBCM_phy) THEN 394 391 soluble=.FALSE. 395 392 spinsol=2 396 ELSEIF (aerosol_name(m).EQ.id_AIPOMM ) THEN393 ELSEIF (aerosol_name(m).EQ.id_AIPOMM_phy) THEN 397 394 soluble=.FALSE. 398 395 spinsol=3 … … 471 468 DO i=1, KLON 472 469 !--anthropogenic aerosol 473 tau_allaer(i,k,2,inu)=tau_ae(i,k,id_ASSO4M ,inu)+tau_ae(i,k,id_CSSO4M,inu)+ &474 tau_ae(i,k,id_ASBCM ,inu)+tau_ae(i,k,id_AIBCM,inu)+ &475 tau_ae(i,k,id_ASPOMM ,inu)+tau_ae(i,k,id_AIPOMM,inu)+ &476 tau_ae(i,k,id_ASSSM ,inu)+tau_ae(i,k,id_CSSSM,inu)+ &477 tau_ae(i,k,id_SSSSM ,inu)+ tau_ae(i,k,id_CIDUSTM,inu)470 tau_allaer(i,k,2,inu)=tau_ae(i,k,id_ASSO4M_phy,inu)+tau_ae(i,k,id_CSSO4M_phy,inu)+ & 471 tau_ae(i,k,id_ASBCM_phy,inu)+tau_ae(i,k,id_AIBCM_phy,inu)+ & 472 tau_ae(i,k,id_ASPOMM_phy,inu)+tau_ae(i,k,id_AIPOMM_phy,inu)+ & 473 tau_ae(i,k,id_ASSSM_phy,inu)+tau_ae(i,k,id_CSSSM_phy,inu)+ & 474 tau_ae(i,k,id_SSSSM_phy,inu)+ tau_ae(i,k,id_CIDUSTM_phy,inu) 478 475 tau_allaer(i,k,2,inu)=MAX(tau_allaer(i,k,2,inu),1e-5) 479 476 480 piz_allaer(i,k,2,inu)=(tau_ae(i,k,id_ASSO4M ,inu)*piz_ae(i,k,id_ASSO4M,inu)+ &481 tau_ae(i,k,id_CSSO4M ,inu)*piz_ae(i,k,id_CSSO4M,inu)+ &482 tau_ae(i,k,id_ASBCM ,inu)*piz_ae(i,k,id_ASBCM,inu)+ &483 tau_ae(i,k,id_AIBCM ,inu)*piz_ae(i,k,id_AIBCM,inu)+ &484 tau_ae(i,k,id_ASPOMM ,inu)*piz_ae(i,k,id_ASPOMM,inu)+ &485 tau_ae(i,k,id_AIPOMM ,inu)*piz_ae(i,k,id_AIPOMM,inu)+ &486 tau_ae(i,k,id_ASSSM ,inu)*piz_ae(i,k,id_ASSSM,inu)+ &487 tau_ae(i,k,id_CSSSM ,inu)*piz_ae(i,k,id_CSSSM,inu)+ &488 tau_ae(i,k,id_SSSSM ,inu)*piz_ae(i,k,id_SSSSM,inu)+ &489 tau_ae(i,k,id_CIDUSTM ,inu)*piz_ae(i,k,id_CIDUSTM,inu)) &477 piz_allaer(i,k,2,inu)=(tau_ae(i,k,id_ASSO4M_phy,inu)*piz_ae(i,k,id_ASSO4M_phy,inu)+ & 478 tau_ae(i,k,id_CSSO4M_phy,inu)*piz_ae(i,k,id_CSSO4M_phy,inu)+ & 479 tau_ae(i,k,id_ASBCM_phy,inu)*piz_ae(i,k,id_ASBCM_phy,inu)+ & 480 tau_ae(i,k,id_AIBCM_phy,inu)*piz_ae(i,k,id_AIBCM_phy,inu)+ & 481 tau_ae(i,k,id_ASPOMM_phy,inu)*piz_ae(i,k,id_ASPOMM_phy,inu)+ & 482 tau_ae(i,k,id_AIPOMM_phy,inu)*piz_ae(i,k,id_AIPOMM_phy,inu)+ & 483 tau_ae(i,k,id_ASSSM_phy,inu)*piz_ae(i,k,id_ASSSM_phy,inu)+ & 484 tau_ae(i,k,id_CSSSM_phy,inu)*piz_ae(i,k,id_CSSSM_phy,inu)+ & 485 tau_ae(i,k,id_SSSSM_phy,inu)*piz_ae(i,k,id_SSSSM_phy,inu)+ & 486 tau_ae(i,k,id_CIDUSTM_phy,inu)*piz_ae(i,k,id_CIDUSTM_phy,inu)) & 490 487 /tau_allaer(i,k,2,inu) 491 488 piz_allaer(i,k,2,inu)=MAX(piz_allaer(i,k,2,inu),0.1) 492 489 493 cg_allaer(i,k,2,inu)=(tau_ae(i,k,id_ASSO4M ,inu)*piz_ae(i,k,id_ASSO4M,inu)*cg_ae(i,k,id_ASSO4M,inu)+ &494 tau_ae(i,k,id_CSSO4M ,inu)*piz_ae(i,k,id_CSSO4M,inu)*cg_ae(i,k,id_CSSO4M,inu)+ &495 tau_ae(i,k,id_ASBCM ,inu)*piz_ae(i,k,id_ASBCM,inu)*cg_ae(i,k,id_ASBCM,inu)+ &496 tau_ae(i,k,id_AIBCM ,inu)*piz_ae(i,k,id_AIBCM,inu)*cg_ae(i,k,id_AIBCM,inu)+ &497 tau_ae(i,k,id_ASPOMM ,inu)*piz_ae(i,k,id_ASPOMM,inu)*cg_ae(i,k,id_ASPOMM,inu)+ &498 tau_ae(i,k,id_AIPOMM ,inu)*piz_ae(i,k,id_AIPOMM,inu)*cg_ae(i,k,id_AIPOMM,inu)+ &499 tau_ae(i,k,id_ASSSM ,inu)*piz_ae(i,k,id_ASSSM,inu)*cg_ae(i,k,id_ASSSM,inu)+ &500 tau_ae(i,k,id_CSSSM ,inu)*piz_ae(i,k,id_CSSSM,inu)*cg_ae(i,k,id_CSSSM,inu)+ &501 tau_ae(i,k,id_SSSSM ,inu)*piz_ae(i,k,id_SSSSM,inu)*cg_ae(i,k,id_SSSSM,inu)+ &502 tau_ae(i,k,id_CIDUSTM ,inu)*piz_ae(i,k,id_CIDUSTM,inu)*cg_ae(i,k,id_CIDUSTM,inu))/ &490 cg_allaer(i,k,2,inu)=(tau_ae(i,k,id_ASSO4M_phy,inu)*piz_ae(i,k,id_ASSO4M_phy,inu)*cg_ae(i,k,id_ASSO4M_phy,inu)+ & 491 tau_ae(i,k,id_CSSO4M_phy,inu)*piz_ae(i,k,id_CSSO4M_phy,inu)*cg_ae(i,k,id_CSSO4M_phy,inu)+ & 492 tau_ae(i,k,id_ASBCM_phy,inu)*piz_ae(i,k,id_ASBCM_phy,inu)*cg_ae(i,k,id_ASBCM_phy,inu)+ & 493 tau_ae(i,k,id_AIBCM_phy,inu)*piz_ae(i,k,id_AIBCM_phy,inu)*cg_ae(i,k,id_AIBCM_phy,inu)+ & 494 tau_ae(i,k,id_ASPOMM_phy,inu)*piz_ae(i,k,id_ASPOMM_phy,inu)*cg_ae(i,k,id_ASPOMM_phy,inu)+ & 495 tau_ae(i,k,id_AIPOMM_phy,inu)*piz_ae(i,k,id_AIPOMM_phy,inu)*cg_ae(i,k,id_AIPOMM_phy,inu)+ & 496 tau_ae(i,k,id_ASSSM_phy,inu)*piz_ae(i,k,id_ASSSM_phy,inu)*cg_ae(i,k,id_ASSSM_phy,inu)+ & 497 tau_ae(i,k,id_CSSSM_phy,inu)*piz_ae(i,k,id_CSSSM_phy,inu)*cg_ae(i,k,id_CSSSM_phy,inu)+ & 498 tau_ae(i,k,id_SSSSM_phy,inu)*piz_ae(i,k,id_SSSSM_phy,inu)*cg_ae(i,k,id_SSSSM_phy,inu)+ & 499 tau_ae(i,k,id_CIDUSTM_phy,inu)*piz_ae(i,k,id_CIDUSTM_phy,inu)*cg_ae(i,k,id_CIDUSTM_phy,inu))/ & 503 500 (tau_allaer(i,k,2,inu)*piz_allaer(i,k,2,inu)) 504 501 505 502 !--natural aerosol 506 tau_allaer(i,k,1,inu)=tau_ae_pi(i,k,id_ASSO4M ,inu)+tau_ae_pi(i,k,id_CSSO4M,inu)+ &507 tau_ae_pi(i,k,id_ASBCM ,inu)+tau_ae_pi(i,k,id_AIBCM,inu)+ &508 tau_ae_pi(i,k,id_ASPOMM ,inu)+tau_ae_pi(i,k,id_AIPOMM,inu)+ &509 tau_ae_pi(i,k,id_ASSSM ,inu)+tau_ae_pi(i,k,id_CSSSM,inu)+ &510 tau_ae_pi(i,k,id_SSSSM ,inu)+ tau_ae_pi(i,k,id_CIDUSTM,inu)503 tau_allaer(i,k,1,inu)=tau_ae_pi(i,k,id_ASSO4M_phy,inu)+tau_ae_pi(i,k,id_CSSO4M_phy,inu)+ & 504 tau_ae_pi(i,k,id_ASBCM_phy,inu)+tau_ae_pi(i,k,id_AIBCM_phy,inu)+ & 505 tau_ae_pi(i,k,id_ASPOMM_phy,inu)+tau_ae_pi(i,k,id_AIPOMM_phy,inu)+ & 506 tau_ae_pi(i,k,id_ASSSM_phy,inu)+tau_ae_pi(i,k,id_CSSSM_phy,inu)+ & 507 tau_ae_pi(i,k,id_SSSSM_phy,inu)+ tau_ae_pi(i,k,id_CIDUSTM_phy,inu) 511 508 tau_allaer(i,k,1,inu)=MAX(tau_allaer(i,k,1,inu),1e-5) 512 509 513 piz_allaer(i,k,1,inu)=(tau_ae_pi(i,k,id_ASSO4M ,inu)*piz_ae(i,k,id_ASSO4M,inu)+ &514 tau_ae_pi(i,k,id_CSSO4M ,inu)*piz_ae(i,k,id_CSSO4M,inu)+ &515 tau_ae_pi(i,k,id_ASBCM ,inu)*piz_ae(i,k,id_ASBCM,inu)+ &516 tau_ae_pi(i,k,id_AIBCM ,inu)*piz_ae(i,k,id_AIBCM,inu)+ &517 tau_ae_pi(i,k,id_ASPOMM ,inu)*piz_ae(i,k,id_ASPOMM,inu)+ &518 tau_ae_pi(i,k,id_AIPOMM ,inu)*piz_ae(i,k,id_AIPOMM,inu)+ &519 tau_ae_pi(i,k,id_ASSSM ,inu)*piz_ae(i,k,id_ASSSM,inu)+ &520 tau_ae_pi(i,k,id_CSSSM ,inu)*piz_ae(i,k,id_CSSSM,inu)+ &521 tau_ae_pi(i,k,id_SSSSM ,inu)*piz_ae(i,k,id_SSSSM,inu)+ &522 tau_ae_pi(i,k,id_CIDUSTM ,inu)*piz_ae(i,k,id_CIDUSTM,inu)) &510 piz_allaer(i,k,1,inu)=(tau_ae_pi(i,k,id_ASSO4M_phy,inu)*piz_ae(i,k,id_ASSO4M_phy,inu)+ & 511 tau_ae_pi(i,k,id_CSSO4M_phy,inu)*piz_ae(i,k,id_CSSO4M_phy,inu)+ & 512 tau_ae_pi(i,k,id_ASBCM_phy,inu)*piz_ae(i,k,id_ASBCM_phy,inu)+ & 513 tau_ae_pi(i,k,id_AIBCM_phy,inu)*piz_ae(i,k,id_AIBCM_phy,inu)+ & 514 tau_ae_pi(i,k,id_ASPOMM_phy,inu)*piz_ae(i,k,id_ASPOMM_phy,inu)+ & 515 tau_ae_pi(i,k,id_AIPOMM_phy,inu)*piz_ae(i,k,id_AIPOMM_phy,inu)+ & 516 tau_ae_pi(i,k,id_ASSSM_phy,inu)*piz_ae(i,k,id_ASSSM_phy,inu)+ & 517 tau_ae_pi(i,k,id_CSSSM_phy,inu)*piz_ae(i,k,id_CSSSM_phy,inu)+ & 518 tau_ae_pi(i,k,id_SSSSM_phy,inu)*piz_ae(i,k,id_SSSSM_phy,inu)+ & 519 tau_ae_pi(i,k,id_CIDUSTM_phy,inu)*piz_ae(i,k,id_CIDUSTM_phy,inu)) & 523 520 /tau_allaer(i,k,1,inu) 524 521 piz_allaer(i,k,1,inu)=MAX(piz_allaer(i,k,1,inu),0.1) 525 522 526 cg_allaer(i,k,1,inu)=(tau_ae_pi(i,k,id_ASSO4M ,inu)*piz_ae(i,k,id_ASSO4M,inu)*cg_ae(i,k,id_ASSO4M,inu)+ &527 tau_ae_pi(i,k,id_CSSO4M ,inu)*piz_ae(i,k,id_CSSO4M,inu)*cg_ae(i,k,id_CSSO4M,inu)+ &528 tau_ae_pi(i,k,id_ASBCM ,inu)*piz_ae(i,k,id_ASBCM,inu)*cg_ae(i,k,id_ASBCM,inu)+ &529 tau_ae_pi(i,k,id_AIBCM ,inu)*piz_ae(i,k,id_AIBCM,inu)*cg_ae(i,k,id_AIBCM,inu)+ &530 tau_ae_pi(i,k,id_ASPOMM ,inu)*piz_ae(i,k,id_ASPOMM,inu)*cg_ae(i,k,id_ASPOMM,inu)+ &531 tau_ae_pi(i,k,id_AIPOMM ,inu)*piz_ae(i,k,id_AIPOMM,inu)*cg_ae(i,k,id_AIPOMM,inu)+ &532 tau_ae_pi(i,k,id_ASSSM ,inu)*piz_ae(i,k,id_ASSSM,inu)*cg_ae(i,k,id_ASSSM,inu)+ &533 tau_ae_pi(i,k,id_CSSSM ,inu)*piz_ae(i,k,id_CSSSM,inu)*cg_ae(i,k,id_CSSSM,inu)+ &534 tau_ae_pi(i,k,id_SSSSM ,inu)*piz_ae(i,k,id_SSSSM,inu)*cg_ae(i,k,id_SSSSM,inu)+ &535 tau_ae_pi(i,k,id_CIDUSTM ,inu)*piz_ae(i,k,id_CIDUSTM,inu)*cg_ae(i,k,id_CIDUSTM,inu))/ &523 cg_allaer(i,k,1,inu)=(tau_ae_pi(i,k,id_ASSO4M_phy,inu)*piz_ae(i,k,id_ASSO4M_phy,inu)*cg_ae(i,k,id_ASSO4M_phy,inu)+ & 524 tau_ae_pi(i,k,id_CSSO4M_phy,inu)*piz_ae(i,k,id_CSSO4M_phy,inu)*cg_ae(i,k,id_CSSO4M_phy,inu)+ & 525 tau_ae_pi(i,k,id_ASBCM_phy,inu)*piz_ae(i,k,id_ASBCM_phy,inu)*cg_ae(i,k,id_ASBCM_phy,inu)+ & 526 tau_ae_pi(i,k,id_AIBCM_phy,inu)*piz_ae(i,k,id_AIBCM_phy,inu)*cg_ae(i,k,id_AIBCM_phy,inu)+ & 527 tau_ae_pi(i,k,id_ASPOMM_phy,inu)*piz_ae(i,k,id_ASPOMM_phy,inu)*cg_ae(i,k,id_ASPOMM_phy,inu)+ & 528 tau_ae_pi(i,k,id_AIPOMM_phy,inu)*piz_ae(i,k,id_AIPOMM_phy,inu)*cg_ae(i,k,id_AIPOMM_phy,inu)+ & 529 tau_ae_pi(i,k,id_ASSSM_phy,inu)*piz_ae(i,k,id_ASSSM_phy,inu)*cg_ae(i,k,id_ASSSM_phy,inu)+ & 530 tau_ae_pi(i,k,id_CSSSM_phy,inu)*piz_ae(i,k,id_CSSSM_phy,inu)*cg_ae(i,k,id_CSSSM_phy,inu)+ & 531 tau_ae_pi(i,k,id_SSSSM_phy,inu)*piz_ae(i,k,id_SSSSM_phy,inu)*cg_ae(i,k,id_SSSSM_phy,inu)+ & 532 tau_ae_pi(i,k,id_CIDUSTM_phy,inu)*piz_ae(i,k,id_CIDUSTM_phy,inu)*cg_ae(i,k,id_CIDUSTM_phy,inu))/ & 536 533 (tau_allaer(i,k,1,inu)*piz_allaer(i,k,1,inu)) 537 534
Note: See TracChangeset
for help on using the changeset viewer.