Index: LMDZ6/trunk/libf/phylmd/radlwsw_m.F90 =================================================================== --- LMDZ6/trunk/libf/phylmd/radlwsw_m.F90 (revision 5664) +++ LMDZ6/trunk/libf/phylmd/radlwsw_m.F90 (revision 5665) @@ -10,8 +10,8 @@ SUBROUTINE radlwsw( & debut, dist, rmu0, fract, & - !albedo SB >>> - ! paprs, pplay,tsol,alb1, alb2, & + !albedo SB >>> + ! paprs, pplay,tsol,alb1, alb2, & paprs, pplay,tsol,SFRWL,alb_dir, alb_dif, & - !albedo SB <<< + !albedo SB <<< t,q,wo,& cldfra, cldemi, cldtaupd,& @@ -38,9 +38,9 @@ solsw_aero, solsw0_aero, & topswcf_aero, solswcf_aero,& - !-C. Kleinschmitt for LW diagnostics + !-C. Kleinschmitt for LW diagnostics toplwad_aero, sollwad_aero,& toplwai_aero, sollwai_aero, & toplwad0_aero, sollwad0_aero, & - !-end + !-end ZLWFT0_i, ZFLDN0, ZFLUP0, & ZSWFT0_i, ZFSDN0, ZFSUP0, & @@ -310,10 +310,10 @@ REAL(KIND=8) cloud_cover_sw(klon) REAL(KIND=8), dimension(klon,klev+1) :: ZFLUX_DIR_i, & ! Direct compt of surf flux into horizontal plane - ZFLUX_DIR_CLEAR_i ! CS Direct + ZFLUX_DIR_CLEAR_i ! CS Direct REAL(KIND=8), dimension(klon,klev+1) :: ZFLUX_DIR, & ! Direct compt of surf flux into horizontal plane - ZFLUX_DIR_CLEAR ! CS Direct + ZFLUX_DIR_CLEAR ! CS Direct REAL(KIND=8), dimension(klon) :: ZFLUX_DIR_INTO_SUN -!!!!!!! Declarations specifiques pour ECRAD !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! + ! Declarations specifiques pour ECRAD ! ! AI 02.2021 #ifdef CPP_ECRAD @@ -341,8 +341,8 @@ ! Interm REAL(KIND=8), dimension(klon) :: ZFLUX_UV, & ! UV flux - ZFLUX_PAR, & ! photosynthetically active radiation similarly - ZFLUX_PAR_CLEAR, & ! CS photosynthetically - ZFLUX_SW_DN_TOA(klon), & ! DN SW flux at TOA - ZEMIS_OUT ! effective broadband emissivity + ZFLUX_PAR, & ! photosynthetically active radiation similarly + ZFLUX_PAR_CLEAR, & ! CS photosynthetically + ZFLUX_SW_DN_TOA(klon), & ! DN SW flux at TOA + ZEMIS_OUT ! effective broadband emissivity REAL(KIND=8), dimension(klon,klev+1) :: ZLWDERIVATIVE ! LW derivatives @@ -364,10 +364,10 @@ ZO2_off(klon,klev) ! O2#endif #endif -!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! + ! REAL(kind=8) POZON(kdlon, kflev, size(wo, 3)) ! mass fraction of ozone ! "POZON(:, :, 1)" is for the average day-night field, ! "POZON(:, :, 2)" is for daylight time. -!!!!! Modif MPL 6.01.09 avec RRTM, on passe de 5 a 6 + ! Modif MPL 6.01.09 avec RRTM, on passe de 5 a 6 REAL(KIND=8) PAER(kdlon,kflev,6) REAL(KIND=8) PCLDLD(kdlon,kflev) @@ -424,5 +424,5 @@ ! "POZON(:, :, 1)" is for the average day-night field, ! "POZON(:, :, 2)" is for daylight time. -!!!!! Modif MPL 6.01.09 avec RRTM, on passe de 5 a 6 + ! Modif MPL 6.01.09 avec RRTM, on passe de 5 a 6 REAL(KIND=8) PAER_i(kdlon,kflev,6) REAL(KIND=8) PDP_i(klon,klev) @@ -562,10 +562,10 @@ IF (type_trac == 'repr') THEN -IF (CPPKEY_REPROBUS) THEN - IF (iflag_rrtm==0) THEN - IF (ok_SUNTIME) PSCT = solaireTIME/zdist/zdist - print*,'Constante solaire: ',PSCT*zdist*zdist - ENDIF -END IF + IF (CPPKEY_REPROBUS) THEN + IF (iflag_rrtm==0) THEN + IF (ok_SUNTIME) PSCT = solaireTIME/zdist/zdist + print*,'Constante solaire: ',PSCT*zdist*zdist + ENDIF + END IF ENDIF @@ -597,5 +597,5 @@ !albedo SB <<< - PEMIS(i) = 1.0 !!!!! A REVOIR (MPL) + PEMIS(i) = 1.0 ! A REVOIR (MPL) PVIEW(i) = 1.66 PPSOL(i) = paprs(iof+i,1) @@ -646,8 +646,8 @@ IF (type_trac == 'repr') THEN -IF (CPPKEY_REPROBUS) THEN - ndimozon = size(wo, 3) - CALL RAD_INTERACTIF(POZON,iof) -END IF + IF (CPPKEY_REPROBUS) THEN + ndimozon = size(wo, 3) + CALL RAD_INTERACTIF(POZON,iof) + END IF ENDIF ! @@ -658,9 +658,9 @@ ENDDO ! -!!!!! Modif MPL 6.01.09 avec RRTM, on passe de 5 a 6 + ! Modif MPL 6.01.09 avec RRTM, on passe de 5 a 6 DO kk = 1, 6 DO k = 1, kflev DO i = 1, kdlon - PAER(i,k,kk) = 1.0E-15 !!!!! A REVOIR (MPL) + PAER(i,k,kk) = 1.0E-15 ! A REVOIR (MPL) ENDDO ENDDO @@ -679,5 +679,5 @@ !===== iflag_rrtm ================================================ ! - IF (iflag_rrtm == 0) THEN !!!! remettre 0 juste pour tester l'ancien rayt via rrtm + IF (iflag_rrtm == 0) THEN ! remettre 0 juste pour tester l'ancien rayt via rrtm ! !--- Mise a zero des tableaux output du rayonnement LW-AR4 ---------- @@ -906,7 +906,7 @@ DO k=1,kflev POZON_i(1:klon,k,:)=POZON(1:klon,kflev+1-k,:) -!!! POZON_i(1:klon,k)=POZON(1:klon,k) !!! on laisse 1=sol et klev=top + ! POZON_i(1:klon,k)=POZON(1:klon,k) ! on laisse 1=sol et klev=top ! print *,'Juste avant RECMWFL: k tsol temp',k,tsol,t(1,k) -!!!!!!! Modif MPL 6.01.09 avec RRTM, on passe de 5 a 6 + ! Modif MPL 6.01.09 avec RRTM, on passe de 5 a 6 DO i=1,6 PAER_i(1:klon,k,i)=PAER(1:klon,kflev+1-k,i) @@ -989,10 +989,10 @@ ! & PTOPSWCFAERO,PSOLSWCFAERO,& ! & PSWADAERO,& !--NL - !!--LW diagnostics CK + !--LW diagnostics CK ! & PTOPLWADAERO,PSOLLWADAERO,& ! & PTOPLWAD0AERO,PSOLLWAD0AERO,& ! & PTOPLWAIAERO,PSOLLWAIAERO,& ! & PLWADAERO,& !--NL - !!..end + !..end ! & ok_ade, ok_aie, ok_volcan, flag_aerosol,flag_aerosol_strat,& ! & flag_aer_feedback) @@ -1118,5 +1118,5 @@ ! IF (fract(i) == 0.) THEN -!!!!! A REVOIR MPL (20090630) ca n a pas de sens quand fract=0 + ! A REVOIR MPL (20090630) ca n a pas de sens quand fract=0 ! pas plus que dans le sw_AR4 zalbpla(i) = 1.0e+39 @@ -1124,6 +1124,6 @@ zalbpla(i) = ZFSUP(i,klev+1)/ZFSDN(i,klev+1) ENDIF -!!! 5 juin 2015 -!!! Correction MP bug RRTM + ! 5 juin 2015 + ! Correction MP bug RRTM zsollwdown(i)= -1.*ZFLDN(i,1) ENDDO @@ -1317,7 +1317,7 @@ ! AI mars 2022 SOLARIRAD = solaire/zdist/zdist - !! diagnos pour la comparaison a la version offline -!!! - Gas en VMR pour offline et MMR pour online -!!! - on utilise pour solarirrad une valeur constante + ! diagnos pour la comparaison a la version offline + ! - Gas en VMR pour offline et MMR pour online + ! - on utilise pour solarirrad une valeur constante if (lldebug_for_offline) then SOLARIRAD = 1366.0896 @@ -1366,37 +1366,37 @@ & debut, ok_volcan, flag_aerosol_strat, & & day_cur, current_time, & - ! Cste solaire/(d_Terre-Soleil)**2 + ! Cste solaire/(d_Terre-Soleil)**2 & SOLARIRAD, & - ! Cos(angle zin), temp sol + ! Cos(angle zin), temp sol & rmu0, tsol, & - ! Albedo diffuse et directe + ! Albedo diffuse et directe & PALBD_NEW,PALBP_NEW, & - ! Emessivite : PEMIS_WINDOW (???), & + ! Emessivite : PEMIS_WINDOW (???), & & ZEMIS, ZEMISW, & - ! longitude(rad), sin(latitude), PMASQ_ ??? + ! longitude(rad), sin(latitude), PMASQ_ ??? & ZGELAM, ZGEMU, & - ! Temp et pres aux interf, vapeur eau, Satur spec humid + ! Temp et pres aux interf, vapeur eau, Satur spec humid & paprs_i, ZTH_i, q_i, qsat_i, & - ! Gas + ! Gas & ZCO2, ZCH4, ZN2O, ZNO2, ZCFC11, ZCFC12, ZHCFC22, & & ZCCL4, POZON_i(:,:,1), ZO2, & - ! nuages : + ! nuages : & cldfra_i, flwc_i, fiwc_i, ZQ_SNOW, & - ! rayons effectifs des gouttelettes + ! rayons effectifs des gouttelettes & ref_liq_i, ref_ice_i, & - ! aerosols + ! aerosols & ZAEROSOL_OLD, ZAEROSOL, & - ! Outputs - ! Net flux : + ! Outputs + ! Net flux : & ZSWFT_i, ZLWFT_i, ZSWFT0_ii, ZLWFT0_ii, & - ! DWN flux : + ! DWN flux : & ZFSDWN_i, ZFLUX_i(:,2,:), ZFCDWN_i, ZFLUC_i(:,2,:), & - ! UP flux : + ! UP flux : & ZFSUP_i, ZFLUX_i(:,1,:), ZFCUP_i, ZFLUC_i(:,1,:), & - ! Surf Direct flux : ATTENTION + ! Surf Direct flux : ATTENTION & ZFLUX_DIR_i, ZFLUX_DIR_CLEAR_i, ZFLUX_DIR_INTO_SUN, & - ! UV and para flux + ! UV and para flux & ZFLUX_UV, ZFLUX_PAR, ZFLUX_PAR_CLEAR, & - ! & ZFLUX_SW_DN_TOA, + ! & ZFLUX_SW_DN_TOA, & ZEMIS_OUT, ZLWDERIVATIVE, & & PSFSWDIF, PSFSWDIR, & @@ -1410,37 +1410,37 @@ & debut, ok_volcan, flag_aerosol_strat, & & day_cur, current_time, & - ! Cste solaire/(d_Terre-Soleil)**2 + ! Cste solaire/(d_Terre-Soleil)**2 & SOLARIRAD, & - ! Cos(angle zin), temp sol + ! Cos(angle zin), temp sol & rmu0, tsol, & - ! Albedo diffuse et directe + ! Albedo diffuse et directe & PALBD_NEW,PALBP_NEW, & - ! Emessivite : PEMIS_WINDOW (???), & + ! Emessivite : PEMIS_WINDOW (???), & & ZEMIS, ZEMISW, & - ! longitude(rad), sin(latitude), PMASQ_ ??? + ! longitude(rad), sin(latitude), PMASQ_ ??? & ZGELAM, ZGEMU, & - ! Temp et pres aux interf, vapeur eau, Satur spec humid + ! Temp et pres aux interf, vapeur eau, Satur spec humid & paprs_i, ZTH_i, q_i, qsat_i, & - ! Gas + ! Gas & ZCO2, ZCH4, ZN2O, ZNO2, ZCFC11, ZCFC12, ZHCFC22, & & ZCCL4, POZON_i(:,:,1), ZO2, & - ! nuages : + ! nuages : & cldfra_i, flwc_i, fiwc_i, ZQ_SNOW, & - ! rayons effectifs des gouttelettes + ! rayons effectifs des gouttelettes & ref_liq_i, ref_ice_i, & - ! aerosols + ! aerosols & ZAEROSOL_OLD, ZAEROSOL, & - ! Outputs - ! Net flux : + ! Outputs + ! Net flux : & ZSWFT_i, ZLWFT_i, ZSWFT0_ii, ZLWFT0_ii, & - ! DWN flux : + ! DWN flux : & ZFSDWN_i, ZFLUX_i(:,2,:), ZFCDWN_i, ZFLUC_i(:,2,:), & - ! UP flux : + ! UP flux : & ZFSUP_i, ZFLUX_i(:,1,:), ZFCUP_i, ZFLUC_i(:,1,:), & - ! Surf Direct flux : ATTENTION + ! Surf Direct flux : ATTENTION & ZFLUX_DIR_i, ZFLUX_DIR_CLEAR_i, ZFLUX_DIR_INTO_SUN, & - ! UV and para flux + ! UV and para flux & ZFLUX_UV, ZFLUX_PAR, ZFLUX_PAR_CLEAR, & - ! & ZFLUX_SW_DN_TOA, + ! & ZFLUX_SW_DN_TOA, & ZEMIS_OUT, ZLWDERIVATIVE, & & PSFSWDIF, PSFSWDIR, &