- Timestamp:
- Sep 17, 2012, 1:45:13 PM (12 years ago)
- Location:
- LMDZ5/trunk/libf/phylmd
- Files:
-
- 3 edited
-
phys_output_mod.F90 (modified) (2 diffs)
-
phys_output_write.h (modified) (2 diffs)
-
sw_aeroAR4.F90 (modified) (21 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/phys_output_mod.F90
r1648 r1653 388 388 389 389 type(ctrl_out),save :: o_topswad = ctrl_out((/ 2, 10, 10, 10, 10, 10 /),'topswad') 390 type(ctrl_out),save :: o_topswad0 = ctrl_out((/ 2, 10, 10, 10, 10, 10 /),'topswad0') 390 391 type(ctrl_out),save :: o_topswai = ctrl_out((/ 2, 10, 10, 10, 10, 10 /),'topswai') 391 392 type(ctrl_out),save :: o_solswad = ctrl_out((/ 2, 10, 10, 10, 10, 10 /),'solswad') 393 type(ctrl_out),save :: o_solswad0 = ctrl_out((/ 2, 10, 10, 10, 10, 10 /),'solswad0') 392 394 type(ctrl_out),save :: o_solswai = ctrl_out((/ 2, 10, 10, 10, 10, 10 /),'solswai') 393 395 … … 1128 1130 o_topswad%flag,o_topswad%name, "ADE at TOA", "W/m2") 1129 1131 CALL histdef2d(iff,clef_stations(iff), & 1132 o_topswad0%flag,o_topswad0%name, "ADE clear-sky at TOA", "W/m2") 1133 CALL histdef2d(iff,clef_stations(iff), & 1130 1134 o_solswad%flag,o_solswad%name, "ADE at SRF", "W/m2") 1135 CALL histdef2d(iff,clef_stations(iff), & 1136 o_solswad0%flag,o_solswad0%name, "ADE clear-sky at SRF", "W/m2") 1131 1137 1132 1138 CALL histdef2d(iff,clef_stations(iff), & -
LMDZ5/trunk/libf/phylmd/phys_output_write.h
r1648 r1653 1258 1258 $ topswad_aero) 1259 1259 ENDIF 1260 IF (o_topswad0%flag(iff)<=lev_files(iff)) THEN 1261 CALL histwrite_phy(nid_files(iff), 1262 $clef_stations(iff), 1263 $o_topswad0%name,itau_w, 1264 $ topswad0_aero) 1265 ENDIF 1260 1266 IF (o_solswad%flag(iff)<=lev_files(iff)) THEN 1261 1267 CALL histwrite_phy(nid_files(iff), … … 1263 1269 $o_solswad%name,itau_w, 1264 1270 $ solswad_aero) 1271 ENDIF 1272 IF (o_solswad0%flag(iff)<=lev_files(iff)) THEN 1273 CALL histwrite_phy(nid_files(iff), 1274 $clef_stations(iff), 1275 $o_solswad0%name,itau_w, 1276 $ solswad0_aero) 1265 1277 ENDIF 1266 1278 -
LMDZ5/trunk/libf/phylmd/sw_aeroAR4.F90
r1641 r1653 56 56 ! -------------- 57 57 ! ORIGINAL : 89-07-14 58 ! 95-01-01 J.-J. MORCRETTE Direct/Diffuse Albedo 59 ! 03-11-27 J. QUAAS Introduce aerosol forcings (based on BOUCHER) 60 ! 09-04 A. COZIC - C.DEANDREIS Indroduce NAT/BC/POM/DUST/SS aerosol forcing 58 ! 1995-01-01 J.-J. MORCRETTE Direct/Diffuse Albedo 59 ! 2003-11-27 J. QUAAS Introduce aerosol forcings (based on BOUCHER) 60 ! 2009-04 A. COZIC - C.DEANDREIS Indroduce NAT/BC/POM/DUST/SS aerosol forcing 61 ! 2012-09 O. BOUCHER - reorganise aerosol cases with ok_ade and ok_aie 61 62 ! ------------------------------------------------------------------ 62 63 ! … … 82 83 83 84 REAL(KIND=8) PCLDSW(KDLON,KFLEV) ! CLOUD FRACTION 84 REAL(KIND=8) PTAU(KDLON,2,KFLEV) ! CLOUD OPTICAL THICKNESS 85 REAL(KIND=8) PTAU(KDLON,2,KFLEV) ! CLOUD OPTICAL THICKNESS (pre-industrial value) 85 86 REAL(KIND=8) PCG(KDLON,2,KFLEV) ! ASYMETRY FACTOR 86 87 REAL(KIND=8) POMEGA(KDLON,2,KFLEV) ! SINGLE SCATTERING ALBEDO … … 140 141 REAL(KIND=8) pizaero(kdlon,kflev,9,2) ! (see aeropt.F) 141 142 REAL(KIND=8) cgaero(kdlon,kflev,9,2) ! -"- 142 REAL(KIND=8) PTAUA(KDLON,2,KFLEV) ! CLOUD OPTICAL THICKNESS (pre -industrialvalue)143 REAL(KIND=8) PTAUA(KDLON,2,KFLEV) ! CLOUD OPTICAL THICKNESS (present-day value) 143 144 REAL(KIND=8) POMEGAA(KDLON,2,KFLEV) ! SINGLE SCATTERING ALBEDO 144 145 REAL(KIND=8) PTOPSWADAERO(KDLON) ! SHORTWAVE FLUX AT T.O.A.(+AEROSOL DIR) 145 146 REAL(KIND=8) PSOLSWADAERO(KDLON) ! SHORTWAVE FLUX AT SURFACE(+AEROSOL DIR) 146 REAL(KIND=8) PTOPSWAD0AERO(KDLON) 147 REAL(KIND=8) PSOLSWAD0AERO(KDLON) 147 REAL(KIND=8) PTOPSWAD0AERO(KDLON) ! SHORTWAVE FLUX AT T.O.A.(+AEROSOL DIR) 148 REAL(KIND=8) PSOLSWAD0AERO(KDLON) ! SHORTWAVE FLUX AT SURFACE(+AEROSOL DIR) 148 149 REAL(KIND=8) PTOPSWAIAERO(KDLON) ! SHORTWAVE FLUX AT T.O.A.(+AEROSOL IND) 149 150 REAL(KIND=8) PSOLSWAIAERO(KDLON) ! SHORTWAVE FLUX AT SURFACE(+AEROSOL IND) 150 REAL(KIND=8) PTOPSWAERO(KDLON,9) 151 REAL(KIND=8) PTOPSW0AERO(KDLON,9) 152 REAL(KIND=8) PSOLSWAERO(KDLON,9) 153 REAL(KIND=8) PSOLSW0AERO(KDLON,9) 151 REAL(KIND=8) PTOPSWAERO(KDLON,9) ! SW TOA AS DRF nat & ant 152 REAL(KIND=8) PTOPSW0AERO(KDLON,9) ! SW SRF AS DRF nat & ant 153 REAL(KIND=8) PSOLSWAERO(KDLON,9) ! SW TOA CS DRF nat & ant 154 REAL(KIND=8) PSOLSW0AERO(KDLON,9) ! SW SRF CS DRF nat & ant 154 155 REAL(KIND=8) PTOPSWCFAERO(KDLON,3) ! SW TOA AS cloudRF nat & ant 155 156 REAL(KIND=8) PSOLSWCFAERO(KDLON,3) ! SW SRF AS cloudRF nat & ant … … 179 180 180 181 ! Key to define the aerosol effect acting on climate 181 ! 0: aerosol feedback active according to ok_ade, ok_aie DEFAULT 182 ! 1: no feedback , zero aerosol fluxes are used for climate, diagnostics according to ok_ade_ok_aie 183 ! 2: feedback according to total aerosol direct effect used for climate, diagnostics according to ok_ade, ok_aie 184 ! 3: feedback according to natural aerosol direct effect used for climate, diagnostics according to ok_ade_ok_aie 185 186 INTEGER,SAVE :: AEROSOLFEEDBACK_ACTIVE = 0 182 ! OB: AEROSOLFEEDBACK_ACTIVE is now a LOGICAL 183 ! TRUE: fluxes use natural and/or anthropogenic aerosols according to ok_ade and ok_aie, DEFAULT 184 ! FALSE: fluxes use no aerosols (case 1) 185 186 LOGICAL,SAVE :: AEROSOLFEEDBACK_ACTIVE = .TRUE. 187 187 !$OMP THREADPRIVATE(AEROSOLFEEDBACK_ACTIVE) 188 188 189 189 CHARACTER (LEN=20) :: modname='sw_aeroAR4' 190 190 CHARACTER (LEN=80) :: abort_message 191 192 IF ((.not. ok_ade) .and. (AEROSOLFEEDBACK_ACTIVE .ge. 2)) THEN193 abort_message ='Error: direct effect is not activated but assumed to be active - see sw_aeroAR4.F90'194 CALL abort_gcm (modname,abort_message,1)195 ENDIF196 AEROSOLFEEDBACK_ACTIVE=MIN(MAX(AEROSOLFEEDBACK_ACTIVE,0),3)197 IF (AEROSOLFEEDBACK_ACTIVE .gt. 3) THEN198 abort_message ='Error: AEROSOLFEEDBACK_ACTIVE options go only until 3'199 CALL abort_gcm (modname,abort_message,1)200 ENDIF201 191 202 192 IF(.NOT.initialized) THEN … … 209 199 ALLOCATE(ZFSUPAI_AERO(KDLON,KFLEV+1)) 210 200 ALLOCATE(ZFSDNAI_AERO(KDLON,KFLEV+1)) 211 ALLOCATE(ZFSUP_AERO (KDLON,KFLEV+1,9)) 212 ALLOCATE(ZFSDN_AERO (KDLON,KFLEV+1,9)) 213 ALLOCATE(ZFSUP0_AERO(KDLON,KFLEV+1,9)) 214 ALLOCATE(ZFSDN0_AERO(KDLON,KFLEV+1,9)) 201 !-OB decrease size of these arrays to what is needed 202 ! | direct effect 203 !ind effect | no aerosol natural total 204 !natural (PTAU) | 1 3 2 --ZFSUP/ZFSDN 205 !total (PTAUA) | 5 4 --ZFSUP/ZFSDN 206 !no cloud | 1 3 2 --ZFSUP0/ZFSDN0 207 ! 208 ALLOCATE(ZFSUP_AERO (KDLON,KFLEV+1,5)) 209 ALLOCATE(ZFSDN_AERO (KDLON,KFLEV+1,5)) 210 ALLOCATE(ZFSUP0_AERO(KDLON,KFLEV+1,3)) 211 ALLOCATE(ZFSDN0_AERO(KDLON,KFLEV+1,3)) 212 ! end OB modif 215 213 ZFSUPAD_AERO(:,:)=0. 216 214 ZFSDNAD_AERO(:,:)=0. … … 226 224 227 225 IF (appel1er) THEN 228 WRITE(lunout,*) 226 WRITE(lunout,*)'SW calling frequency : ', swpas 229 227 WRITE(lunout,*) " In general, it should be 1" 230 228 appel1er = .FALSE. … … 241 239 ENDDO 242 240 243 ! clear sky is either computed IF no direct effect is asked for, or for extended diag)244 IF ( swaero_diag .or. .not. ok_ade) THEN241 ! clear sky with no aerosols at all is computed IF ACTIVEFEEDBACK_ACTIVE is false or for extended diag 242 IF ( swaero_diag .or. .not. AEROSOLFEEDBACK_ACTIVE ) THEN 245 243 246 244 ! clear-sky: zero aerosol effect … … 268 266 ENDDO 269 267 ENDDO 270 ENDIF 271 272 ! cloudy sky is either computed IF no indirect effect is asked for, or for extended diag)273 IF ( swaero_diag .or. .not. ok_aie) THEN268 ENDIF ! swaero_diag .or. .not. AEROSOLFEEDBACK_ACTIVE 269 270 ! cloudy sky with no aerosols at all is either computed IF no indirect effect is asked for, or for extended diag 271 IF ( swaero_diag .or. .not. AEROSOLFEEDBACK_ACTIVE ) THEN 274 272 ! cloudy-sky: zero aerosol effect 275 273 flag_aer=0.0 … … 297 295 ENDDO 298 296 ENDDO 299 ENDIF 300 297 ENDIF ! swaero_diag .or. .not. AEROSOLFEEDBACK_ACTIVE 301 298 302 299 IF (ok_ade) THEN 303 300 304 ! clear sky (Anne Cozic 03/07/2007)direct effect of total aerosol301 ! clear sky direct effect of total aerosol 305 302 ! CAS AER (2) 306 303 flag_aer=1.0 … … 329 326 ENDDO 330 327 331 ! cloudy sky is either computed IF no indirect effect is asked for, or for extended diag) 332 IF (swaero_diag .or. .not. ok_aie) THEN 333 ! cloudy-sky aerosol direct effect of total aerosol 328 IF ( swaero_diag ) THEN 329 ! clear sky direct effect natural aerosol 330 ! CAS AER (3) 331 flag_aer=1.0 332 CALL SWU_LMDAR4(PSCT,ZCLDSW0,PPMB,PPSOL,& 333 PRMU0,PFRAC,PTAVE,PWV,& 334 ZAKI,ZCLD,ZCLEAR,ZDSIG,ZFACT,ZRMU,ZSEC,ZUD) 335 INU = 1 336 CALL SW1S_LMDAR4(INU, PAER, flag_aer,& 337 tauaero(:,:,3,:), pizaero(:,:,3,:), cgaero(:,:,3,:),& 338 PALBD, PALBP, PCG, ZCLD, ZCLEAR, PCLDSW,& 339 ZDSIG, POMEGA, ZOZ, ZRMU, ZSEC, PTAU, ZUD,& 340 ZFD, ZFU) 341 INU = 2 342 CALL SW2S_LMDAR4(INU, PAER, flag_aer,& 343 tauaero(:,:,3,:), pizaero(:,:,3,:), cgaero(:,:,3,:),& 344 ZAKI, PALBD, PALBP, PCG, ZCLD, ZCLEAR, PCLDSW,& 345 ZDSIG, POMEGA, ZOZ, ZRMU, ZSEC, PTAU, ZUD,& 346 PWV, PQS,& 347 ZFDOWN, ZFUP) 348 349 DO JK = 1 , KFLEV+1 350 DO JL = 1, KDLON 351 ZFSUP0_AERO(JL,JK,3) = (ZFUP(JL,JK) + ZFU(JL,JK)) * ZFACT(JL) 352 ZFSDN0_AERO(JL,JK,3) = (ZFDOWN(JL,JK) + ZFD(JL,JK)) * ZFACT(JL) 353 ENDDO 354 ENDDO 355 ENDIF !swaero_diag 356 357 ! cloudy-sky with natural aerosols for indirect effect 358 ! but total aerosols for direct effect 359 ! PTAU 360 ! CAS AER (2) 334 361 flag_aer=1.0 335 362 CALL SWU_LMDAR4(PSCT,PCLDSW,PPMB,PPSOL,& … … 356 383 ENDDO 357 384 ENDDO 358 ENDIF 359 360 ! natural aeroosl clear sky is computed for extended diag) 361 IF ( swaero_diag ) THEN 362 ! clear sky direct effect natural aerosol 363 flag_aer=1.0 364 CALL SWU_LMDAR4(PSCT,ZCLDSW0,PPMB,PPSOL,& 365 PRMU0,PFRAC,PTAVE,PWV,& 366 ZAKI,ZCLD,ZCLEAR,ZDSIG,ZFACT,ZRMU,ZSEC,ZUD) 367 INU = 1 368 CALL SW1S_LMDAR4(INU, PAER, flag_aer,& 369 tauaero(:,:,3,:), pizaero(:,:,3,:), cgaero(:,:,3,:),& 370 PALBD, PALBP, PCG, ZCLD, ZCLEAR, PCLDSW,& 371 ZDSIG, POMEGA, ZOZ, ZRMU, ZSEC, PTAU, ZUD,& 372 ZFD, ZFU) 373 INU = 2 374 CALL SW2S_LMDAR4(INU, PAER, flag_aer,& 375 tauaero(:,:,3,:), pizaero(:,:,3,:), cgaero(:,:,3,:),& 376 ZAKI, PALBD, PALBP, PCG, ZCLD, ZCLEAR, PCLDSW,& 377 ZDSIG, POMEGA, ZOZ, ZRMU, ZSEC, PTAU, ZUD,& 378 PWV, PQS,& 379 ZFDOWN, ZFUP) 380 381 DO JK = 1 , KFLEV+1 382 DO JL = 1, KDLON 383 ZFSUP0_AERO(JL,JK,3) = (ZFUP(JL,JK) + ZFU(JL,JK)) * ZFACT(JL) 384 ZFSDN0_AERO(JL,JK,3) = (ZFDOWN(JL,JK) + ZFD(JL,JK)) * ZFACT(JL) 385 ENDDO 386 ENDDO 387 ENDIF 388 389 ! cloud sky natural is for extended diagnostics 385 ! ENDIF 386 387 ENDIF !-end ok_ade 388 389 IF ( ok_ade .or. (.not.ok_ade .and. ok_aie) ) THEN 390 391 ! cloudy-sky with natural aerosols for indirect effect 392 ! and natural aerosols for direct effect 393 ! PTAU 394 ! CAS AER (3) 390 395 IF ( swaero_diag ) THEN 391 396 ! cloudy-sky direct effect natural aerosol … … 414 419 ENDDO 415 420 ENDDO 416 ENDIF 417 418 ENDIF ! ok_ade 419 420 ! cloudy sky needs to be computed in all cases IF ok_aie is activated 421 IF (ok_aie) THEN 422 !jq cloudy-sky + aerosol direct + aerosol indirect of total aerosol 421 ENDIF !swaero_diag 422 423 ENDIF !--true/false or false/true 424 425 IF (ok_ade .and. ok_aie) THEN 426 427 ! cloudy-sky with total aerosols for indirect effect 428 ! and total aerosols for direct effect 429 ! PTAUA 430 ! CAS AER (2) 423 431 flag_aer=1.0 424 432 CALL SWU_LMDAR4(PSCT,PCLDSW,PPMB,PPSOL,& … … 438 446 PWV, PQS,& 439 447 ZFDOWN, ZFUP) 448 440 449 DO JK = 1 , KFLEV+1 441 450 DO JL = 1, KDLON … … 444 453 ENDDO 445 454 ENDDO 455 456 ENDIF ! ok_ade .and. ok_aie 457 458 IF (ok_aie) THEN 459 ! cloudy-sky with total aerosols for indirect effect 460 ! and natural aerosols for direct effect 461 ! PTAUA 462 ! CAS AER (3) 463 flag_aer=1.0 464 CALL SWU_LMDAR4(PSCT,PCLDSW,PPMB,PPSOL,& 465 PRMU0,PFRAC,PTAVE,PWV,& 466 ZAKI,ZCLD,ZCLEAR,ZDSIG,ZFACT,ZRMU,ZSEC,ZUD) 467 INU = 1 468 CALL SW1S_LMDAR4(INU, PAER, flag_aer,& 469 tauaero(:,:,3,:), pizaero(:,:,3,:), cgaero(:,:,3,:),& 470 PALBD, PALBP, PCG, ZCLD, ZCLEAR, PCLDSW,& 471 ZDSIG, POMEGAA, ZOZ, ZRMU, ZSEC, PTAUA, ZUD,& 472 ZFD, ZFU) 473 INU = 2 474 CALL SW2S_LMDAR4(INU, PAER, flag_aer,& 475 tauaero(:,:,3,:), pizaero(:,:,3,:), cgaero(:,:,3,:),& 476 ZAKI, PALBD, PALBP, PCG, ZCLD, ZCLEAR, PCLDSW,& 477 ZDSIG, POMEGAA, ZOZ, ZRMU, ZSEC, PTAUA, ZUD,& 478 PWV, PQS,& 479 ZFDOWN, ZFUP) 480 481 DO JK = 1 , KFLEV+1 482 DO JL = 1, KDLON 483 ZFSUP_AERO(JL,JK,5) = (ZFUP(JL,JK) + ZFU(JL,JK)) * ZFACT(JL) 484 ZFSDN_AERO(JL,JK,5) = (ZFDOWN(JL,JK) + ZFD(JL,JK)) * ZFACT(JL) 485 ENDDO 486 ENDDO 487 446 488 ENDIF ! ok_aie 447 489 … … 450 492 itapsw = itapsw + 1 451 493 452 IF ( AEROSOLFEEDBACK_ACTIVE .eq. 0) THEN494 IF ( AEROSOLFEEDBACK_ACTIVE ) THEN 453 495 IF ( ok_ade .and. ok_aie ) THEN 454 496 ZFSUP(:,:) = ZFSUP_AERO(:,:,4) … … 457 499 ZFSDN0(:,:) = ZFSDN0_AERO(:,:,2) 458 500 ENDIF 501 459 502 IF ( ok_ade .and. (.not. ok_aie) ) THEN 460 503 ZFSUP(:,:) = ZFSUP_AERO(:,:,2) … … 465 508 466 509 IF ( (.not. ok_ade) .and. ok_aie ) THEN 467 print*,'Warning: indirect effect in cloudy regions includes direct aerosol effect'468 ZFS UP(:,:) = ZFSUP_AERO(:,:,4)469 ZFS DN(:,:) = ZFSDN_AERO(:,:,4)470 ZFS UP0(:,:) = ZFSUP0_AERO(:,:,1)471 ZFSDN0(:,:) = ZFSDN0_AERO(:,:,1)472 ENDIF 510 ZFSUP(:,:) = ZFSUP_AERO(:,:,5) 511 ZFSDN(:,:) = ZFSDN_AERO(:,:,5) 512 ZFSUP0(:,:) = ZFSUP0_AERO(:,:,3) 513 ZFSDN0(:,:) = ZFSDN0_AERO(:,:,3) 514 ENDIF 515 473 516 IF ((.not. ok_ade) .and. (.not. ok_aie)) THEN 517 ZFSUP(:,:) = ZFSUP_AERO(:,:,3) 518 ZFSDN(:,:) = ZFSDN_AERO(:,:,3) 519 ZFSUP0(:,:) = ZFSUP0_AERO(:,:,3) 520 ZFSDN0(:,:) = ZFSDN0_AERO(:,:,3) 521 ENDIF 522 523 ! MS the following allows to compute the forcing diagostics without 524 ! letting the aerosol forcing act on the meteorology 525 ! SEE logic above 526 ELSE 474 527 ZFSUP(:,:) = ZFSUP_AERO(:,:,1) 475 528 ZFSDN(:,:) = ZFSDN_AERO(:,:,1) … … 478 531 ENDIF 479 532 480 ! MS the following allows to compute the forcing diagostics without 481 ! letting the aerosol forcing act on the meteorology 482 ! SEE logic above 483 ELSEIF ( AEROSOLFEEDBACK_ACTIVE .gt. 0) THEN 484 ZFSUP(:,:) = ZFSUP_AERO(:,:,AEROSOLFEEDBACK_ACTIVE) 485 ZFSDN(:,:) = ZFSDN_AERO(:,:,AEROSOLFEEDBACK_ACTIVE) 486 ZFSUP0(:,:) = ZFSUP0_AERO(:,:,AEROSOLFEEDBACK_ACTIVE) 487 ZFSDN0(:,:) = ZFSDN0_AERO(:,:,AEROSOLFEEDBACK_ACTIVE) 488 ENDIF 489 490 533 ! Now computes heating rates 491 534 DO k = 1, KFLEV 492 535 kpl1 = k+1 … … 511 554 PTOPSW(i) = ZFSDN(i,KFLEV+1) - ZFSUP(i,KFLEV+1) 512 555 513 514 556 ! net anthropogenic forcing direct and 1st indirect effect diagnostics 515 557 ! requires a natural aerosol field read and used 516 558 ! Difference of net fluxes from double call to radiation 517 559 518 519 560 IF (ok_ade) THEN 520 561 521 562 ! indices 1: natural; 2 anthropogenic 563 522 564 ! TOA/SRF all sky natural forcing 523 565 PSOLSWAERO(i,1) = (ZFSDN_AERO(i,1,3) - ZFSUP_AERO(i,1,3))-(ZFSDN_AERO(i,1,1) - ZFSUP_AERO(i,1,1)) 524 566 PTOPSWAERO(i,1) = (ZFSDN_AERO(i,KFLEV+1,3) - ZFSUP_AERO(i,KFLEV+1,3))- (ZFSDN_AERO(i,KFLEV+1,1) - ZFSUP_AERO(i,KFLEV+1,1)) 525 567 568 ! TOA/SRF clear sky natural forcing 569 PSOLSW0AERO(i,1) = (ZFSDN0_AERO(i,1,3) - ZFSUP0_AERO(i,1,3))-(ZFSDN0_AERO(i,1,1) - ZFSUP0_AERO(i,1,1)) 570 PTOPSW0AERO(i,1) = (ZFSDN0_AERO(i,KFLEV+1,3) - ZFSUP0_AERO(i,KFLEV+1,3))-(ZFSDN0_AERO(i,KFLEV+1,1) - ZFSUP0_AERO(i,KFLEV+1,1)) 571 572 IF (ok_aie) THEN 573 574 ! TOA/SRF all sky anthropogenic forcing 575 PSOLSWAERO(i,2) = (ZFSDN_AERO(i,1,4) - ZFSUP_AERO(i,1,4))-(ZFSDN_AERO(i,1,5) - ZFSUP_AERO(i,1,5)) 576 PTOPSWAERO(i,2) = (ZFSDN_AERO(i,KFLEV+1,4) - ZFSUP_AERO(i,KFLEV+1,4))- (ZFSDN_AERO(i,KFLEV+1,5) - ZFSUP_AERO(i,KFLEV+1,5)) 577 578 ELSE 579 526 580 ! TOA/SRF all sky anthropogenic forcing 527 581 PSOLSWAERO(i,2) = (ZFSDN_AERO(i,1,2) - ZFSUP_AERO(i,1,2))-(ZFSDN_AERO(i,1,3) - ZFSUP_AERO(i,1,3)) 528 582 PTOPSWAERO(i,2) = (ZFSDN_AERO(i,KFLEV+1,2) - ZFSUP_AERO(i,KFLEV+1,2))- (ZFSDN_AERO(i,KFLEV+1,3) - ZFSUP_AERO(i,KFLEV+1,3)) 529 583 530 ! TOA/SRF clear sky natural forcing 531 PSOLSW0AERO(i,1) = (ZFSDN0_AERO(i,1,3) - ZFSUP0_AERO(i,1,3))-(ZFSDN0_AERO(i,1,1) - ZFSUP0_AERO(i,1,1)) 532 PTOPSW0AERO(i,1) = (ZFSDN0_AERO(i,KFLEV+1,3) - ZFSUP0_AERO(i,KFLEV+1,3))-(ZFSDN0_AERO(i,KFLEV+1,1) - ZFSUP0_AERO(i,KFLEV+1,1)) 584 ENDIF 533 585 534 586 ! TOA/SRF clear sky anthropogenic forcing … … 536 588 PTOPSW0AERO(i,2) = (ZFSDN0_AERO(i,KFLEV+1,2) - ZFSUP0_AERO(i,KFLEV+1,2))-(ZFSDN0_AERO(i,KFLEV+1,3) - ZFSUP0_AERO(i,KFLEV+1,3)) 537 589 590 ! direct anthropogenic forcing , as in old LMDzT, however differences of net fluxes 591 PSOLSWADAERO(i) = PSOLSWAERO(i,2) 592 PTOPSWADAERO(i) = PTOPSWAERO(i,2) 593 PSOLSWAD0AERO(i) = PSOLSW0AERO(i,2) 594 PTOPSWAD0AERO(i) = PTOPSW0AERO(i,2) 595 596 ! OB: these diagnostics may not always work but who need them 538 597 ! Cloud forcing indices 1: natural; 2 anthropogenic; 3: zero aerosol direct effect 539 598 ! Instantaneously computed cloudy sky direct aerosol effect, cloud forcing due to aerosols above clouds … … 552 611 PTOPSWCFAERO(i,3) = (ZFSDN_AERO(i,KFLEV+1,1) - ZFSUP_AERO(i,KFLEV+1,1))- (ZFSDN0_AERO(i,KFLEV+1,1) - ZFSUP0_AERO(i,KFLEV+1,1)) 553 612 554 ! direct anthropogenic forcing , as in old LMDzT, however differences of net fluxes555 PSOLSWADAERO(i) = PSOLSWAERO(i,2)556 PTOPSWADAERO(i) = PTOPSWAERO(i,2)557 PSOLSWAD0AERO(i) = PSOLSW0AERO(i,2)558 PTOPSWAD0AERO(i) = PTOPSW0AERO(i,2)559 560 613 ENDIF 561 614 562 563 615 IF (ok_aie) THEN 616 IF (ok_ade) THEN 564 617 PSOLSWAIAERO(i) = (ZFSDN_AERO(i,1,4) - ZFSUP_AERO(i,1,4))-(ZFSDN_AERO(i,1,2) - ZFSUP_AERO(i,1,2)) 565 618 PTOPSWAIAERO(i) = (ZFSDN_AERO(i,KFLEV+1,4) - ZFSUP_AERO(i,KFLEV+1,4))-(ZFSDN_AERO(i,KFLEV+1,2) - ZFSUP_AERO(i,KFLEV+1,2)) 619 ELSE 620 PSOLSWAIAERO(i) = (ZFSDN_AERO(i,1,5) - ZFSUP_AERO(i,1,5))-(ZFSDN_AERO(i,1,3) - ZFSUP_AERO(i,1,3)) 621 PTOPSWAIAERO(i) = (ZFSDN_AERO(i,KFLEV+1,5) - ZFSUP_AERO(i,KFLEV+1,5))-(ZFSDN_AERO(i,KFLEV+1,3) - ZFSUP_AERO(i,KFLEV+1,3)) 622 ENDIF 566 623 ENDIF 567 624 568 ENDDO 625 ENDDO 626 569 627 END SUBROUTINE SW_AEROAR4
Note: See TracChangeset
for help on using the changeset viewer.
