--- trunk/libf/phylmd/Radlwsw/sw.f90 2011/10/07 13:11:58 53 +++ trunk/libf/phylmd/Radlwsw/sw.f90 2013/06/24 15:39:52 70 @@ -20,7 +20,7 @@ ! 3. COMPUTES FLUXES IN 2ND SPECTRAL INTERVAL (SW2S) ! REFERENCE. - ! SEE RADIATION'S PART OF THE ECMWF RESEARCH DEPARTMENT + ! SEE RADIATION PART OF THE ECMWF RESEARCH DEPARTMENT ! DOCUMENTATION, AND FOUQUART AND BONNEL (1980) ! AUTHOR. @@ -32,7 +32,7 @@ ! 03-11-27 J. QUAAS Introduce aerosol forcings (based on BOUCHER) USE clesphys, ONLY: bug_ozone - USE suphec_m, ONLY: rcpd, rday, rg, rmd, rmo3 + USE suphec_m, ONLY: rcpd, rday, rg, md, rmo3 USE raddim, ONLY: kdlon, kflev ! ARGUMENTS: @@ -101,7 +101,7 @@ DATA appel1er /.TRUE./ !jq-Introduced for aerosol forcings double precision flag_aer - logical ok_ade, ok_aie ! use aerosol forcings or not? + logical, intent(in):: ok_ade, ok_aie ! use aerosol forcings or not? double precision tauae(kdlon, kflev, 2) ! aerosol optical properties double precision pizae(kdlon, kflev, 2) ! aerosol optical properties(see aeropt.F) @@ -157,7 +157,7 @@ *PDP(JL, JK)*(101325.0/PPSOL(JL)) ELSE ! Correction MPL 100505 - ZOZ(JL, JK) = POZON(JL, JK)*RMD/RMO3*46.6968/RG*PDP(JL, JK) + ZOZ(JL, JK) = POZON(JL, JK)*MD/RMO3*46.6968/RG*PDP(JL, JK) ENDIF ENDDO ENDDO @@ -186,7 +186,7 @@ ENDDO ENDDO - flag_aer=0.0 + flag_aer=0. CALL SWU(PSCT, PCLDSW, PPMB, PPSOL, & PRMU0, PFRAC, PTAVE, PWV, & ZAKI, ZCLD, ZCLEAR, ZDSIG, ZFACT, ZRMU, ZSEC, ZUD) @@ -215,7 +215,7 @@ IF (ok_ade) THEN ! cloudy-sky + aerosol dir OB - flag_aer=1.0 + flag_aer=1. CALL SWU(PSCT, PCLDSW, PPMB, PPSOL, & PRMU0, PFRAC, PTAVE, PWV, & ZAKI, ZCLD, ZCLEAR, ZDSIG, ZFACT, ZRMU, ZSEC, ZUD)