Changeset 2685 for trunk/LMDZ.MARS/libf
- Timestamp:
- May 27, 2022, 4:57:03 PM (2 years ago)
- Location:
- trunk/LMDZ.MARS/libf/phymars
- Files:
-
- 4 edited
-
callradite_mod.F (modified) (5 diffs)
-
physiq_mod.F (modified) (20 diffs)
-
rocketduststorm_mod.F90 (modified) (2 diffs)
-
topmons_mod.F90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/libf/phymars/callradite_mod.F
r2678 r2685 7 7 SUBROUTINE callradite(icount,ngrid,nlayer,nq,zday,ls,pq,albedo, 8 8 $ emis,mu0,pplev,pplay,pt,tsurf,fract,dist_sol,igout, 9 $ dtlw,dtsw,fluxsurf_lw,fluxsurf_sw,fluxtop_lw, 10 $ fluxtop_sw,tau_pref_scenario,tau_pref_gcm, 9 $ dtlw,dtsw,fluxsurf_lw,fluxsurf_dn_sw,fluxsurf_up_sw, 10 $ fluxtop_lw,fluxtop_dn_sw,fluxtop_up_sw, 11 & tau_pref_scenario,tau_pref_gcm, 11 12 & tau,aerosol,dsodust,tauscaling,dust_rad_adjust,IRtoVIScoef, 12 13 $ taucloudtes,rdust,rice,nuice,riceco2,nuiceco2,co2ice, … … 141 142 c nueffrad(ngrid,nlayer,naerkind) Aerosol effective variance 142 143 143 c144 c output:145 c -------146 c dtlw (ngrid,nlayer) longwave (IR) heating rate (K/s)147 c dtsw(ngrid,nlayer) shortwave (Solar) heating rate (K/s)148 c fluxsurf_lw(ngrid) surface downward flux tota LW (thermal IR) (W.m-2)149 c fluxsurf_sw(ngrid,1) surface downward flux SW for solar band#1 (W.m-2)150 c fluxsurf_sw(ngrid,2) surface downward flux SW for solar band#2 (W.m-2)151 c152 c fluxtop_lw(ngrid) outgoing upward flux tota LW (thermal IR) (W.m-2)153 c fluxtop_sw(ngrid,1) outgoing upward flux SW for solar band#1 (W.m-2)154 c fluxtop_sw(ngrid,2) outgoing upward flux SW for solar band#2 (W.m-2)155 156 c tau Column total visible dust optical depth at each point157 c aerosol(ngrid,nlayer,naerkind) aerosol extinction optical depth158 c at reference wavelength "longrefvis" set159 c in dimradmars_h , in each layer, for one of160 c the "naerkind" kind of aerosol optical161 c properties.162 144 c======================================================================= 163 145 c … … 189 171 REAL,INTENT(IN) :: tsurf(ngrid) 190 172 REAL,INTENT(IN) :: dist_sol,mu0(ngrid),fract(ngrid) 191 REAL,INTENT(OUT) :: dtlw(ngrid,nlayer),dtsw(ngrid,nlayer) 192 REAL,INTENT(OUT) :: fluxsurf_lw(ngrid), fluxtop_lw(ngrid) 193 REAL,INTENT(OUT) :: fluxsurf_sw(ngrid,2), fluxtop_sw(ngrid,2) 173 REAL,INTENT(OUT) :: dtlw(ngrid,nlayer) ! longwave (IR) heating rate (K/s) 174 REAL,INTENT(OUT) :: dtsw(ngrid,nlayer) ! shortwave (Solar) heating rate (K/s) 175 REAL,INTENT(OUT) :: fluxsurf_lw(ngrid) ! total LW (thermal IR) downward flux 176 ! (W.m-2) at the surface 177 REAL,INTENT(OUT) :: fluxtop_lw(ngrid) ! outgoing total LW (thermal IR) 178 ! upward flux (W.m-2) at the top of the atm. 179 REAL,INTENT(OUT) :: fluxsurf_dn_sw(ngrid,2) ! surface downward SW flux for 180 ! solar bands #1 and #2 (W.m-2) 181 REAL,INTENT(OUT) :: fluxsurf_up_sw(ngrid,2) ! surface upward SW flux for 182 ! solar bands #1 and #2 (W.m-2) 183 REAL,INTENT(OUT) :: fluxtop_dn_sw(ngrid,2) ! incoming downward SW flux for 184 ! solar bands #1 and #2 (W.m-2) at top of atm. 185 REAL,INTENT(OUT) :: fluxtop_up_sw(ngrid,2) ! outgoing upward SW flux for 186 ! solar bands #1 and #2 (W.m-2) at top of atm. 194 187 REAL,INTENT(OUT) :: tau_pref_scenario(ngrid) ! prescribed dust column 195 188 ! visible opacity at odpref from scenario 196 189 REAL,INTENT(OUT) :: tau_pref_gcm(ngrid) ! computed dust column 197 190 ! visible opacity at odpref in the GCM 198 REAL,INTENT(OUT) :: tau(ngrid,naerkind) 191 REAL,INTENT(OUT) :: tau(ngrid,naerkind) ! Column visible optical depth 192 ! for each aerosol 199 193 REAL,INTENT(OUT) :: taucloudtes(ngrid)! Cloud opacity at infrared 200 194 ! reference wavelength using 201 195 ! Qabs instead of Qext 202 196 ! (direct comparison with TES) 203 REAL,INTENT(OUT) :: aerosol(ngrid,nlayer,naerkind) 197 REAL,INTENT(OUT) :: aerosol(ngrid,nlayer,naerkind) ! aerosol extinction 198 ! optical depth at reference wavelength "longrefvis", 199 ! set in dimradmars_h, for each kind of aerosol 204 200 REAL,INTENT(INOUT) :: dsodust(ngrid,nlayer) 205 201 REAL,INTENT(OUT) :: rdust(ngrid,nlayer) ! Dust geometric mean radius (m) … … 585 581 enddo 586 582 583 ! copy SW fluxes at surface and TOA 587 584 do ig = 1,nd 588 fluxsurf_sw(ig0+ig,1) = zfluxd_sw(ig,1,1) 589 fluxsurf_sw(ig0+ig,2) = zfluxd_sw(ig,1,2) 590 fluxtop_sw(ig0+ig,1) = zfluxu_sw(ig,nlaylte+1,1) 591 fluxtop_sw(ig0+ig,2) = zfluxu_sw(ig,nlaylte+1,2) 585 ! surface downward SW flux 586 fluxsurf_dn_sw(ig0+ig,1) = zfluxd_sw(ig,1,1) 587 fluxsurf_dn_sw(ig0+ig,2) = zfluxd_sw(ig,1,2) 588 ! surface upward SW flux 589 fluxsurf_up_sw(ig0+ig,1) = zfluxu_sw(ig,1,1) 590 fluxsurf_up_sw(ig0+ig,2) = zfluxu_sw(ig,1,2) 591 ! downward SW flux at top of atmosphere 592 fluxtop_dn_sw(ig0+ig,1) = zfluxd_sw(ig,nlaylte+1,1) 593 fluxtop_dn_sw(ig0+ig,2) = zfluxd_sw(ig,nlaylte+1,2) 594 ! upward SW flux at top of atmosphere 595 fluxtop_up_sw(ig0+ig,1) = zfluxu_sw(ig,nlaylte+1,1) 596 fluxtop_up_sw(ig0+ig,2) = zfluxu_sw(ig,nlaylte+1,2) 592 597 enddo 593 598 … … 613 618 PRINT*,albedo(igout,1),emis(igout),mu0(igout), 614 619 s fract(igout), fluxsurf_lw(igout), 615 $ fluxsurf_ sw(igout,1)+fluxsurf_sw(igout,2)620 $ fluxsurf_dn_sw(igout,1)+fluxsurf_dn_sw(igout,2) 616 621 PRINT*,'Tlay Tlev Play Plev dT/dt SW dT/dt LW (K/s)' 617 622 PRINT*,'daysec',daysec -
trunk/LMDZ.MARS/libf/phymars/physiq_mod.F
r2678 r2685 284 284 REAL sl_ra, sl_di0 285 285 REAL sky 286 REAL direct_flux(ngrid)! Incident direct solar flux on Mars at surface (W.m-2)286 REAL fluxsurf_dir_dn_sw(ngrid) ! Incident direct solar flux on Mars at surface (W.m-2) 287 287 288 288 REAL,PARAMETER :: stephan = 5.67e-08 ! Stephan Boltzman constant … … 298 298 299 299 REAL fluxsurf_lw(ngrid) !incident LW (IR) surface flux (W.m-2) 300 REAL fluxsurf_sw(ngrid,2) !incident SW (solar) surface flux (W.m-2) 300 REAL fluxsurf_dn_sw(ngrid,2) ! Incident SW (solar) surface flux (W.m-2) 301 REAL fluxsurf_up_sw(ngrid,2) ! Reflected SW (solar) surface flux (W.m-2) 301 302 REAL fluxtop_lw(ngrid) !Outgoing LW (IR) flux to space (W.m-2) 302 REAL fluxtop_sw(ngrid,2) !Outgoing SW (solar) flux to space (W.m-2) 303 REAL fluxtop_dn_sw(ngrid,2) ! Incoming SW (solar) flux from space (W.m-2) 304 REAL fluxtop_up_sw(ngrid,2) ! Outgoing SW (solar) flux to space (W.m-2) 303 305 REAL tau_pref_scenario(ngrid) ! prescribed dust column visible opacity 304 306 ! at odpref … … 386 388 REAL zq(ngrid,nlayer,nq) 387 389 388 REAL fluxtop_sw_tot(ngrid), fluxsurf_sw_tot(ngrid) 390 REAL fluxtop_dn_sw_tot(ngrid), fluxtop_up_sw_tot(ngrid) 391 REAL fluxsurf_dn_sw_tot(ngrid), fluxsurf_up_sw_tot(ngrid) 389 392 character*2 str2 390 393 ! character*5 str5 … … 502 505 real zdtlwclf(ngrid,nlayer) 503 506 real fluxsurf_lwclf(ngrid) 504 real fluxsurf_ swclf(ngrid,2)507 real fluxsurf_dn_swclf(ngrid,2),fluxsurf_up_swclf(ngrid,2) 505 508 real fluxtop_lwclf(ngrid) 506 real fluxtop_ swclf(ngrid,2)509 real fluxtop_dn_swclf(ngrid,2),fluxtop_up_swclf(ngrid,2) 507 510 real taucloudtesclf(ngrid) 508 511 real tf_clf, ntf_clf ! tf: fraction of clouds, ntf: fraction without clouds … … 976 979 CALL callradite(icount,ngrid,nlayer,nq,zday,zls,pq,albedo, 977 980 & emis,mu0,zplev,zplay,pt,tsurf,fract,dist_sol,igout, 978 & zdtlw,zdtsw,fluxsurf_lw,fluxsurf_sw,fluxtop_lw, 979 & fluxtop_sw,tau_pref_scenario,tau_pref_gcm, 981 & zdtlw,zdtsw,fluxsurf_lw,fluxsurf_dn_sw,fluxsurf_up_sw, 982 & fluxtop_lw,fluxtop_dn_sw,fluxtop_up_sw, 983 & tau_pref_scenario,tau_pref_gcm, 980 984 & tau,aerosol,dsodust,tauscaling,dust_rad_adjust,IRtoVIScoef, 981 985 & taucloudtes,rdust,rice,nuice,riceco2,nuiceco2,co2ice, … … 992 996 CALL callradite(icount,ngrid,nlayer,nq,zday,zls,pq, 993 997 & albedo,emis,mu0,zplev,zplay,pt,tsurf,fract, 994 & dist_sol,igout,zdtlwclf,zdtswclf,fluxsurf_lwclf, 995 & fluxsurf_swclf,fluxtop_lwclf,fluxtop_swclf, 998 & dist_sol,igout,zdtlwclf,zdtswclf, 999 & fluxsurf_lwclf,fluxsurf_dn_swclf,fluxsurf_up_swclf, 1000 & fluxtop_lwclf,fluxtop_dn_swclf,fluxtop_up_swclf, 996 1001 & tau_pref_scenario,tau_pref_gcm,tau,aerosol, 997 1002 & dsodust,tauscaling,dust_rad_adjust,IRtoVIScoef, … … 1009 1014 fluxsurf_lw(ig) = ntf_clf*fluxsurf_lwclf(ig) 1010 1015 & + tf_clf*fluxsurf_lw(ig) 1011 fluxsurf_sw(ig,1) = ntf_clf*fluxsurf_swclf(ig,1) 1012 & + tf_clf*fluxsurf_sw(ig,1) 1013 fluxsurf_sw(ig,2) = ntf_clf*fluxsurf_swclf(ig,2) 1014 & + tf_clf*fluxsurf_sw(ig,2) 1016 fluxsurf_dn_sw(ig,1:2) = 1017 & ntf_clf*fluxsurf_dn_swclf(ig,1:2) 1018 & + tf_clf*fluxsurf_dn_sw(ig,1:2) 1019 fluxsurf_up_sw(ig,1:2) = 1020 & ntf_clf*fluxsurf_up_swclf(ig,1:2) 1021 & + tf_clf*fluxsurf_up_sw(ig,1:2) 1015 1022 fluxtop_lw(ig) = ntf_clf*fluxtop_lwclf(ig) 1016 1023 & + tf_clf*fluxtop_lw(ig) 1017 fluxtop_ sw(ig,1) = ntf_clf*fluxtop_swclf(ig,1)1018 & + tf_clf*fluxtop_sw(ig,1)1019 fluxtop_ sw(ig,2) = ntf_clf*fluxtop_swclf(ig,2)1020 & + tf_clf*fluxtop_sw(ig,2)1024 fluxtop_dn_sw(ig,1:2)=ntf_clf*fluxtop_dn_swclf(ig,1:2) 1025 & + tf_clf*fluxtop_dn_sw(ig,1:2) 1026 fluxtop_up_sw(ig,1:2)=ntf_clf*fluxtop_up_swclf(ig,1:2) 1027 & + tf_clf*fluxtop_up_sw(ig,1:2) 1021 1028 taucloudtes(ig) = ntf_clf*taucloudtesclf(ig) 1022 1029 & + tf_clf*taucloudtes(ig) … … 1067 1074 sl_the = theta_sl(ig) 1068 1075 IF (sl_the .ne. 0.) THEN 1069 ztim1=fluxsurf_ sw(ig,1)+fluxsurf_sw(ig,2)1076 ztim1=fluxsurf_dn_sw(ig,1)+fluxsurf_dn_sw(ig,2) 1070 1077 DO l=1,2 1071 1078 sl_lct = ptime*24. + 180.*longitude(ig)/pi/15. … … 1075 1082 sl_alb = albedo(ig,l) 1076 1083 sl_psi = psi_sl(ig) 1077 sl_fl0 = fluxsurf_ sw(ig,l)1084 sl_fl0 = fluxsurf_dn_sw(ig,l) 1078 1085 sl_di0 = 0. 1079 1086 if (mu0(ig) .gt. 0.) then … … 1081 1088 sl_di0 = sl_di0*flux_1AU/dist_sol/dist_sol 1082 1089 sl_di0 = sl_di0/ztim1 1083 sl_di0 = fluxsurf_ sw(ig,l)*sl_di01090 sl_di0 = fluxsurf_dn_sw(ig,l)*sl_di0 1084 1091 endif 1085 1092 ! you never know (roundup concern...) … … 1090 1097 & sl_di0, sl_fl0, sl_flu ) 1091 1098 !!!!!!!!!!!!!!!!!!!!!!!!!! 1092 fluxsurf_ sw(ig,l) = sl_flu1099 fluxsurf_dn_sw(ig,l) = sl_flu 1093 1100 ENDDO 1094 1101 !!! compute correction on IR flux as well … … 1111 1118 DO ig=1,ngrid 1112 1119 fluxrad_sky(ig)=emis(ig)*fluxsurf_lw(ig) 1113 $ +fluxsurf_ sw(ig,1)*(1.-albedo(ig,1))1114 $ +fluxsurf_ sw(ig,2)*(1.-albedo(ig,2))1120 $ +fluxsurf_dn_sw(ig,1)*(1.-albedo(ig,1)) 1121 $ +fluxsurf_dn_sw(ig,2)*(1.-albedo(ig,2)) 1115 1122 ENDDO 1116 1123 … … 2040 2047 $ co2ice,albedo,emis,rdust, 2041 2048 $ zdtc,zdtsurfc,pdpsrf,zduc,zdvc,zdqc, 2042 $ fluxsurf_ sw,zls,2049 $ fluxsurf_dn_sw,zls, 2043 2050 $ zdqssed_co2,zcondicea_co2microp, 2044 2051 & zdqsc) … … 2278 2285 2279 2286 c Sum of fluxes in solar spectral bands (for output only) 2280 DO ig=1,ngrid 2281 fluxtop_sw_tot(ig)=fluxtop_sw(ig,1) + fluxtop_sw(ig,2) 2282 fluxsurf_sw_tot(ig)=fluxsurf_sw(ig,1) + fluxsurf_sw(ig,2) 2283 ENDDO 2287 fluxtop_dn_sw_tot(1:ngrid)=fluxtop_dn_sw(1:ngrid,1) + 2288 & fluxtop_dn_sw(1:ngrid,2) 2289 fluxtop_up_sw_tot(1:ngrid)=fluxtop_up_sw(1:ngrid,1) + 2290 & fluxtop_up_sw(1:ngrid,2) 2291 fluxsurf_dn_sw_tot(1:ngrid)=fluxsurf_dn_sw(1:ngrid,1) + 2292 & fluxsurf_dn_sw(1:ngrid,2) 2293 fluxsurf_up_sw_tot(1:ngrid)=fluxsurf_up_sw(1:ngrid,1) + 2294 & fluxsurf_up_sw(1:ngrid,2) 2295 2284 2296 c ******* TEST ****************************************************** 2285 2297 ztim1 = 999 … … 2452 2464 do ig=1,ngrid 2453 2465 if(mu0(ig).le.0.01) then 2454 direct_flux(ig) = 0.2466 fluxsurf_dir_dn_sw(ig) = 0. 2455 2467 else 2456 direct_flux(ig) = flux_1AU/dist_sol/dist_sol*mu0(ig)*2468 fluxsurf_dir_dn_sw(ig) = flux_1AU/dist_sol/dist_sol*mu0(ig)* 2457 2469 & exp(-(tau(ig,iaer_dust_doubleq)+ 2458 2470 & tau(ig,iaer_h2o_ice))/mu0(ig)) … … 2657 2669 & "Thermal IR radiative flux to surface","W.m-2",2, 2658 2670 & fluxsurf_lw) 2659 call wstats(ngrid,"fluxsurf_sw", 2660 & "Solar radiative flux to surface","W.m-2",2, 2661 & fluxsurf_sw_tot) 2671 call wstats(ngrid,"fluxsurf_dn_sw", 2672 & "Incoming Solar radiative flux to surface","W.m-2",2, 2673 & fluxsurf_dn_sw_tot) 2674 call wstats(ngrid,"fluxsurf_up_sw", 2675 & "Reflected Solar radiative flux from surface","W.m-2",2, 2676 & fluxsurf_up_sw_tot) 2662 2677 call wstats(ngrid,"fluxtop_lw", 2663 2678 & "Thermal IR radiative flux to space","W.m-2",2, 2664 2679 & fluxtop_lw) 2665 call wstats(ngrid,"fluxtop_sw", 2666 & "Solar radiative flux to space","W.m-2",2, 2667 & fluxtop_sw_tot) 2680 call wstats(ngrid,"fluxtop_dn_sw", 2681 & "Incoming Solar radiative flux from space","W.m-2",2, 2682 & fluxtop_dn_sw_tot) 2683 call wstats(ngrid,"fluxtop_up_sw", 2684 & "Outgoing Solar radiative flux to space","W.m-2",2, 2685 & fluxtop_up_sw_tot) 2668 2686 call wstats(ngrid,"temp","Atmospheric temperature","K",3,zt) 2669 2687 call wstats(ngrid,"u","Zonal (East-West) wind","m.s-1",3,zu) … … 2685 2703 c call wstats(ngrid,"lw_htrt","lw heat.rate", 2686 2704 c & "W.m-2",3,zdtlw) 2687 call wstats(ngrid," direct_flux",2688 & "Direct incoming flux at surface",2689 & "W.m-2",2, direct_flux)2705 call wstats(ngrid,"fluxsurf_dir_dn_sw", 2706 & "Direct incoming SW flux at surface", 2707 & "W.m-2",2,fluxsurf_dir_dn_sw) 2690 2708 2691 2709 if (calltherm) then … … 2921 2939 comm_HR_LW(1:ngrid,1:nlayer) = zdtlw(1:ngrid,1:nlayer) 2922 2940 !state real SWDOWNZ ij misc 1 - h "SWDOWNZ" "DOWNWARD SW FLUX AT SURFACE" "W m-2" 2923 comm_SWDOWNZ(1:ngrid) = fluxsurf_ sw_tot(1:ngrid)2941 comm_SWDOWNZ(1:ngrid) = fluxsurf_dn_sw_tot(1:ngrid) 2924 2942 !state real TAU_DUST ij misc 1 - h "TAU_DUST" "REFERENCE VISIBLE DUST OPACITY" "" 2925 2943 comm_TAU_DUST(1:ngrid) = tau_pref_gcm(1:ngrid) … … 2995 3013 call WRITEDIAGFI(ngrid,"fluxsurf_lw","fluxsurf_lw","W.m-2",2, 2996 3014 & fluxsurf_lw) 2997 call WRITEDIAGFI(ngrid,"fluxsurf_ sw","fluxsurf_sw","W.m-2",2,2998 & fluxsurf_sw_tot)3015 call WRITEDIAGFI(ngrid,"fluxsurf_dn_sw","fluxsurf_dn_sw", 3016 & "W.m-2",2,fluxsurf_dn_sw_tot) 2999 3017 call WRITEDIAGFI(ngrid,"fluxtop_lw","fluxtop_lw","W.m-2",2, 3000 3018 & fluxtop_lw) 3001 call WRITEDIAGFI(ngrid,"fluxtop_ sw","fluxtop_sw","W.m-2",2,3002 & fluxtop_sw_tot)3019 call WRITEDIAGFI(ngrid,"fluxtop_up_sw","fluxtop_up_sw","W.m-2", 3020 & 2,fluxtop_up_sw_tot) 3003 3021 call WRITEDIAGFI(ngrid,"temp","temperature","K",3,zt) 3004 3022 call WRITEDIAGFI(ngrid,"Sols","Time","sols",0,[zday]) … … 4016 4034 CALL send_xios_field("area",cell_area) 4017 4035 4018 ! CALL send_xios_field("ISR",fluxtop_ sw_tot)4036 ! CALL send_xios_field("ISR",fluxtop_dn_sw_tot) 4019 4037 CALL send_xios_field("OLR",fluxtop_lw) 4020 4038 -
trunk/LMDZ.MARS/libf/phymars/rocketduststorm_mod.F90
r2643 r2685 167 167 ! variables for the radiative transfer 168 168 REAL fluxsurf_lw1(ngrid) 169 REAL fluxsurf_ sw1(ngrid,2)169 REAL fluxsurf_dn_sw1(ngrid,2),fluxsurf_up_sw1(ngrid,2) 170 170 REAL fluxtop_lw1(ngrid) 171 REAL fluxtop_ sw1(ngrid,2)171 REAL fluxtop_dn_sw1(ngrid,2),fluxtop_up_sw1(ngrid,2) 172 172 REAL tau(ngrid,naerkind) 173 173 REAL aerosol(ngrid,nlayer,naerkind) … … 255 255 CALL callradite(icount,ngrid,nlayer,nq,zday,zls,pq,albedo, & 256 256 emis,mu0,pplev,pplay,pt,tsurf,fract,dist_sol,igout, & 257 zdtlw1,zdtsw1,fluxsurf_lw1,fluxsurf_sw1,fluxtop_lw1, & 258 fluxtop_sw1,tau_pref_scenario,tau_pref_gcm, & 257 zdtlw1,zdtsw1,fluxsurf_lw1,fluxsurf_dn_sw1,fluxsurf_up_sw1, & 258 fluxtop_lw1,fluxtop_dn_sw1,fluxtop_up_sw1, & 259 tau_pref_scenario,tau_pref_gcm, & 259 260 tau,aerosol,dsodust,tauscaling,dust_rad_adjust,IRtoVIScoef, & 260 261 taucloudtes,rdust,rice,nuice,riceco2,nuiceco2,co2ice,rstormdust,rtopdust, & -
trunk/LMDZ.MARS/libf/phymars/topmons_mod.F90
r2643 r2685 120 120 ! Variables for radiative transfer 121 121 REAL fluxsurf_lw1(ngrid) 122 REAL fluxsurf_ sw1(ngrid,2)122 REAL fluxsurf_dn_sw1(ngrid,2),fluxsurf_up_sw1(ngrid,2) 123 123 REAL fluxtop_lw1(ngrid) 124 REAL fluxtop_ sw1(ngrid,2)124 REAL fluxtop_dn_sw1(ngrid,2),fluxtop_up_sw1(ngrid,2) 125 125 REAL tau(ngrid,naerkind) 126 126 REAL taucloudtes(ngrid) … … 286 286 CALL callradite(icount,ngrid,nlayer,nq,zday,zls,zq,albedo, & 287 287 emis,mu0,pplev,pplay,pt,tsurf,fract,dist_sol,igout, & 288 zdtlw1,zdtsw1,fluxsurf_lw1,fluxsurf_sw1,fluxtop_lw1, & 289 fluxtop_sw1,tau_pref_scenario,tau_pref_gcm, & 288 zdtlw1,zdtsw1,fluxsurf_lw1,fluxsurf_dn_sw1,fluxsurf_up_sw1, & 289 fluxtop_lw1,fluxtop_dn_sw1,fluxtop_up_sw1, & 290 tau_pref_scenario,tau_pref_gcm, & 290 291 tau,aerosol,dsodust,tauscaling,dust_rad_adjust,IRtoVIScoef, & 291 292 taucloudtes,rdust,rice,nuice,riceco2,nuiceco2,co2ice,rstormdust,rtopdust, &
Note: See TracChangeset
for help on using the changeset viewer.
