Changeset 3052
- Timestamp:
- Oct 31, 2017, 12:24:32 PM (7 years ago)
- Location:
- LMDZ6/branches/IPSLCM6.0.13
- Files:
-
- 10 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/IPSLCM6.0.13/DefLists/CMIP6_ping_atmos.xml
r3021 r3052 222 222 <field id="CMIP6_rsdcs4co2" field_ref="rsdcs4co2" /> <!-- P1 (W m-2) downwelling_shortwave_flux_in_air_assuming_clear_sky : Downwelling clear-sky shortwave radiation calculated using carbon dioxide concentrations increased fourfold --> 223 223 <field id="CMIP6_rsdcsaf" field_ref="dummy_XYA" /> <!-- P1 (W m-2) downwelling_shortwave_flux_assuming_clean_clear_sky : Calculated in the absence of aerosols and clouds (following Ghan). This requires a double-call in the radiation code with precisely the same meteorology. --> 224 <field id="CMIP6_rsdcsafbnd" field_ref="dummy_ XYA" /> <!-- P1 (W m-2) band_downwelling_shortwave_flux_assuming_clean_clear_sky : Calculated in the absence of aerosols and clouds (following Ghan). This requires a double-call in the radiation code with precisely the same meteorology. -->225 <field id="CMIP6_rsdcsbnd" field_ref="dummy_ XYA" /> <!-- P1 (W m-2) band_downwelling_shortwave_flux_assuming_clear_sky : Calculated with aerosols but without clouds. This is a standard clear-sky calculation -->224 <field id="CMIP6_rsdcsafbnd" field_ref="dummy_not_provided" /> <!-- P1 (W m-2) band_downwelling_shortwave_flux_assuming_clean_clear_sky : Calculated in the absence of aerosols and clouds (following Ghan). This requires a double-call in the radiation code with precisely the same meteorology. --> 225 <field id="CMIP6_rsdcsbnd" field_ref="dummy_not_provided" /> <!-- P1 (W m-2) band_downwelling_shortwave_flux_assuming_clear_sky : Calculated with aerosols but without clouds. This is a standard clear-sky calculation --> 226 226 <field id="CMIP6_rsds" field_ref="SWdnSFC" /> <!-- P1 (W m-2) surface_downwelling_shortwave_flux_in_air : surface solar irradiance for UV calculations --> 227 227 <field id="CMIP6_rsdscs" field_ref="SWdnSFCclr" /> <!-- P1 (W m-2) surface_downwelling_shortwave_flux_in_air_assuming_clear_sky : surface solar irradiance clear sky for UV calculations --> 228 228 <field id="CMIP6_rsdscsaf" field_ref="dummy_XY" /> <!-- P1 (W m-2) surface_downwelling_shortwave_flux_in_air_assuming_clean_clear_sky : Calculated in the absence of aerosols and clouds. --> 229 <field id="CMIP6_rsdscsafbnd" field_ref="dummy_ XY" /> <!-- P1 (W m-2) surface_downwelling_shortwave_flux_in_air_assuming_clean_clear_sky : Calculated in the absence of aerosols and clouds, following Ghan (2013, ACP). This requires a double-call in the radiation code with precisely the same meteorology. -->230 <field id="CMIP6_rsdscsbnd" field_ref="dummy_ XY" /> <!-- P1 (W m-2) surface_downwelling_shortwave_flux_in_air_assuming_clear_sky : Calculated with aerosols but without clouds. This is a standard clear-sky calculation -->229 <field id="CMIP6_rsdscsafbnd" field_ref="dummy_not_provided" /> <!-- P1 (W m-2) surface_downwelling_shortwave_flux_in_air_assuming_clean_clear_sky : Calculated in the absence of aerosols and clouds, following Ghan (2013, ACP). This requires a double-call in the radiation code with precisely the same meteorology. --> 230 <field id="CMIP6_rsdscsbnd" field_ref="dummy_not_provided" /> <!-- P1 (W m-2) surface_downwelling_shortwave_flux_in_air_assuming_clear_sky : Calculated with aerosols but without clouds. This is a standard clear-sky calculation --> 231 231 <field id="CMIP6_rsdscsdiff" field_ref="dummy_XY" /> <!-- P1 (W m-2) surface_diffuse_downwelling_shortwave_flux_in_air_assuming_clear_sky : unset --> 232 232 <field id="CMIP6_rsdsdiff" field_ref="dummy_XY" /> <!-- P1 (W m-2) surface_diffuse_downwelling_shortwave_flux_in_air : unset --> … … 238 238 <field id="CMIP6_rsucs4co2" field_ref="rsucs4co2" /> <!-- P1 (W m-2) upwelling_shortwave_flux_in_air_assuming_clear_sky : Upwelling clear-sky shortwave radiation calculated using carbon dioxide concentrations increased fourfold --> 239 239 <field id="CMIP6_rsucsaf" field_ref="dummy_XYA" /> <!-- P1 (W m-2) upwelling_shortwave_flux_assuming_clean_clear_sky : Calculated in the absence of aerosols and clouds (following Ghan). This requires a double-call in the radiation code with precisely the same meteorology. --> 240 <field id="CMIP6_rsucsafbnd" field_ref="dummy_ XYA" /> <!-- P1 (W m-2) band_upwelling_shortwave_flux_assuming_clean_clear_sky : Calculated in the absence of aerosols and clouds (following Ghan). This requires a double-call in the radiation code with precisely the same meteorology. -->241 <field id="CMIP6_rsucsbnd" field_ref="dummy_ XYA" /> <!-- P1 (W m-2) band_upwelling_shortwave_flux_assuming_clear_sky : Calculated with aerosols but without clouds. This is a standard clear-sky calculation -->240 <field id="CMIP6_rsucsafbnd" field_ref="dummy_not_provided" /> <!-- P1 (W m-2) band_upwelling_shortwave_flux_assuming_clean_clear_sky : Calculated in the absence of aerosols and clouds (following Ghan). This requires a double-call in the radiation code with precisely the same meteorology. --> 241 <field id="CMIP6_rsucsbnd" field_ref="dummy_not_provided" /> <!-- P1 (W m-2) band_upwelling_shortwave_flux_assuming_clear_sky : Calculated with aerosols but without clouds. This is a standard clear-sky calculation --> 242 242 <field id="CMIP6_rsus" field_ref="SWupSFC" /> <!-- P1 (W m-2) surface_upwelling_shortwave_flux_in_air : Surface Upwelling Shortwave Radiation --> 243 243 <field id="CMIP6_rsuscs" field_ref="SWupSFCclr" /> <!-- P1 (W m-2) surface_upwelling_shortwave_flux_in_air_assuming_clear_sky : Surface Upwelling Clear-sky Shortwave Radiation --> 244 244 <field id="CMIP6_rsuscsaf" field_ref="dummy_XY" /> <!-- P1 (W m-2) surface_upwelling_shortwave_flux_in_air_assuming_clean_clear_sky : Surface Upwelling Clear-sky, Aerosol Free Shortwave Radiation --> 245 <field id="CMIP6_rsuscsafbnd" field_ref="dummy_ XY" /> <!-- P1 (W m-2) surface_upwelling_shortwave_flux_in_air_assuming_clean_clear_sky : Calculated in the absence of aerosols and clouds, following Ghan (ACP, 2013). This requires a double-call in the radiation code with precisely the same meteorology. -->246 <field id="CMIP6_rsuscsbnd" field_ref="dummy_ XY" /> <!-- P1 (W m-2) surface_upwelling_shortwave_flux_in_air_assuming_clear_sky : Calculated with aerosols but without clouds. This is a standard clear-sky calculation -->245 <field id="CMIP6_rsuscsafbnd" field_ref="dummy_not_provided" /> <!-- P1 (W m-2) surface_upwelling_shortwave_flux_in_air_assuming_clean_clear_sky : Calculated in the absence of aerosols and clouds, following Ghan (ACP, 2013). This requires a double-call in the radiation code with precisely the same meteorology. --> 246 <field id="CMIP6_rsuscsbnd" field_ref="dummy_not_provided" /> <!-- P1 (W m-2) surface_upwelling_shortwave_flux_in_air_assuming_clear_sky : Calculated with aerosols but without clouds. This is a standard clear-sky calculation --> 247 247 <field id="CMIP6_rsut" field_ref="SWupTOA" /> <!-- P1 (W m-2) toa_outgoing_shortwave_flux : at the top of the atmosphere --> 248 248 <field id="CMIP6_rsut4co2" field_ref="rsut4co2" /> <!-- P1 (W m-2) toa_outgoing_shortwave_flux : TOA Outgoing Shortwave Radiation calculated using carbon dioxide concentrations increased fourfold --> 249 249 <field id="CMIP6_rsutcs" field_ref="SWupTOAclr" /> <!-- P1 (W m-2) toa_outgoing_shortwave_flux_assuming_clear_sky : Calculated in the absence of clouds. --> 250 250 <field id="CMIP6_rsutcs4co2" field_ref="rsutcs4co2" /> <!-- P1 (W m-2) toa_outgoing_shortwave_flux_assuming_clear_sky : TOA Outgoing Clear-Sky Shortwave Radiation calculated using carbon dioxide concentrations increased fourfold --> 251 <field id="CMIP6_rsutcsafbnd" field_ref="dummy_ XY" /> <!-- P1 (W m-2) band_toa_outgoing_shortwave_flux_assuming_clean_clear_sky : Calculated in the absence of aerosols and clouds, following Ghan (2013, ACP). This requires a double-call in the radiation code with precisely the same meteorology. -->252 <field id="CMIP6_rsutcsbnd" field_ref="dummy_ XY" /> <!-- P1 (W m-2) toa_outgoing_shortwave_flux_assuming_clear_sky : Calculated with aerosols but without clouds. This is a standard clear-sky calculation -->251 <field id="CMIP6_rsutcsafbnd" field_ref="dummy_not_provided" /> <!-- P1 (W m-2) band_toa_outgoing_shortwave_flux_assuming_clean_clear_sky : Calculated in the absence of aerosols and clouds, following Ghan (2013, ACP). This requires a double-call in the radiation code with precisely the same meteorology. --> 252 <field id="CMIP6_rsutcsbnd" field_ref="dummy_not_provided" /> <!-- P1 (W m-2) toa_outgoing_shortwave_flux_assuming_clear_sky : Calculated with aerosols but without clouds. This is a standard clear-sky calculation --> 253 253 <field id="CMIP6_rtmt" field_ref="nettop" /> <!-- P1 (W m-2) net_downward_radiative_flux_at_top_of_atmosphere_model : Net Downward Radiative Flux at Top of Model : I.e., at the top of that portion of the atmosphere where dynamics are explicitly treated by the model. This is reported only if it differs from the net downward radiative flux at the top of the atmosphere. --> 254 254 <field id="CMIP6_rv" field_ref="dummy_XY" /> <!-- P1 (s-1) atmosphere_relative_vorticity : Relative Vorticity at 850 hPa --> … … 270 270 <field id="CMIP6_snrefr" field_ref="dummy_XY" /> <!-- P1 (kg m-2 s-1) surface_snow_and_ice_refreezing_flux : Re-freezing of water in the snow --> 271 271 <field id="CMIP6_snwc" field_ref="dummy_XY" /> <!-- P1 (kg m-2) canopy_snow_amount : Total water mass of the snowpack (liquid or frozen), averaged over a grid cell and interecepted by the canopy. --> 272 <field id="CMIP6_solbnd" field_ref=" dummy_XY" /> <!-- P1 (W m-2) band_solar_insolation : Top-of-Atmosphere Solar Insolation for each band -->272 <field id="CMIP6_solbnd" field_ref="solbnd" /> <!-- P1 (W m-2) band_solar_insolation : Top-of-Atmosphere Solar Insolation for each band --> 273 273 <field id="CMIP6_swsffluxaero" field_ref="dummy_XY" /> <!-- P2 (W m-2) shortwave__flux_due_to_volcanic_aerosols_at__the_surface : shortwave heating rate due to volcanic aerosols to be diagnosed through double radiation call, zonal average values required --> 274 274 <field id="CMIP6_swsrfasdust" field_ref="dummy_XY" /> <!-- P1 (W m-2) tendency_of_all_sky_surface_shortwave_flux_due_to_dust_ambient_aerosol_particles : Balkanski - LSCE --> -
LMDZ6/branches/IPSLCM6.0.13/DefLists/context_lmdz.xml
r3019 r3052 61 61 <axis id="bnds" standard_name="bounds" unit="1" > 62 62 </axis> 63 <axis id="spectband" standard_name="Sensor Band Central Radiation Wavenumber" unit="m-1"> 64 </axis> 63 65 64 66 <!-- Cosp axis definitions--> … … 109 111 <domain id="dom_glo" /> 110 112 <axis id="plev" /> 113 </grid> 114 115 <grid id="grid_glo_spectband"> 116 <domain id="dom_glo" /> 117 <axis id="spectband" /> 111 118 </grid> 112 119 -
LMDZ6/branches/IPSLCM6.0.13/DefLists/field_def_lmdz.xml
r3015 r3052 637 637 <field id="dtlif" long_name="Orography dT" unit="K/s" /> 638 638 <field id="dthin" long_name="Hines GWD dT" unit="K/s" /> 639 <field id="rsu" long_name="SW upward radiation" unit="W m-2" enabled="FALSE"/>640 <field id="rsd" long_name="SW downward radiation" unit="W m-2" enabled="FALSE"/>641 <field id="rlu" long_name="LW upward radiation" unit="W m-2" enabled="FALSE"/>642 <field id="rld" long_name="LW downward radiation" unit="W m-2" enabled="FALSE"/>643 <field id="rsucs" long_name="SW CS upward radiation" unit="W m-2" enabled="FALSE"/>644 <field id="rsdcs" long_name="SW CS downward radiation" unit="W m-2" enabled="FALSE"/>645 <field id="rlucs" long_name="LW CS upward radiation" unit="W m-2" enabled="FALSE"/>646 <field id="rldcs" long_name="LW CS downward radiation" unit="W m-2" enabled="FALSE"/>639 <field id="rsu" long_name="SW upward radiation" unit="W m-2" /> 640 <field id="rsd" long_name="SW downward radiation" unit="W m-2" /> 641 <field id="rlu" long_name="LW upward radiation" unit="W m-2" /> 642 <field id="rld" long_name="LW downward radiation" unit="W m-2" /> 643 <field id="rsucs" long_name="SW CS upward radiation" unit="W m-2" /> 644 <field id="rsdcs" long_name="SW CS downward radiation" unit="W m-2" /> 645 <field id="rlucs" long_name="LW CS upward radiation" unit="W m-2" /> 646 <field id="rldcs" long_name="LW CS downward radiation" unit="W m-2" /> 647 647 <field id="tnt" long_name="Tendency of air temperature" unit="K s-1" /> 648 648 <field id="tntc" long_name="Tendency of air temperature due to Moist Convection" unit="K s-1" /> … … 783 783 </field_group> 784 784 785 <field_group id="fields_solbnd_3D" grid_ref="grid_glo_spectband" > 786 <field id="solbnd" long_name="Top-of-Atmosphere Solar Insolation for each band" unit="W m-2" /> 787 </field_group> 788 785 789 <field_group id="fields_COSP_CALIPSO_2D" grid_ref="grid_glo" freq_op="3h"> 786 790 <field id="cllcalipso" long_name="Lidar Low-level Cloud Fraction" unit="1" detect_missing_value=".true." /> -
LMDZ6/branches/IPSLCM6.0.13/DefLists/file_def_histmth_lmdz.xml
r3003 r3052 660 660 </field_group> 661 661 662 <!-- VARS 3D --> 663 <field_group operation="average" freq_op="1ts" grid_ref="grid_glo_spectband"> 664 <field field_ref="solbnd" level="4" /> 665 </field_group> 666 662 667 </file> 663 668 -
LMDZ6/branches/IPSLCM6.0.13/libf/phylmd/clesphys.h
r2989 r3052 11 11 LOGICAL ok_limitvrai 12 12 LOGICAL ok_all_xml 13 LOGICAL ok_lwoff 13 14 INTEGER nbapp_rad, iflag_con, nbapp_cv, nbapp_wk, iflag_ener_conserv 14 15 REAL co2_ppm, co2_ppm0, solaire … … 137 138 & , iflag_ice_thermo, ok_gwd_rando, NSW, iflag_albedo & 138 139 & , ok_chlorophyll,ok_conserv_q, adjust_tropopause & 139 & , ok_daily_climoz, ok_all_xml 140 & , ok_daily_climoz, ok_all_xml, ok_lwoff & 140 141 & , iflag_phytrac 141 142 -
LMDZ6/branches/IPSLCM6.0.13/libf/phylmd/conf_phys_m.F90
r3023 r3052 91 91 LOGICAL, SAVE :: ok_newmicro_omp 92 92 LOGICAL, SAVE :: ok_all_xml_omp 93 LOGICAL, SAVE :: ok_lwoff_omp 93 94 LOGICAL, SAVE :: ok_journe_omp, ok_mensuel_omp, ok_instan_omp, ok_hf_omp 94 95 LOGICAL, SAVE :: ok_LES_omp … … 294 295 ok_all_xml_omp = .FALSE. 295 296 CALL getin('ok_all_xml', ok_all_xml_omp) 297 298 ! 299 !Config Key = ok_lwoff 300 !Config Desc = inhiber l effet radiatif LW des nuages 301 !Config Def = .FALSE. 302 !Config Help = 303 ! 304 ok_lwoff_omp = .FALSE. 305 CALL getin('ok_lwoff', ok_lwoff_omp) 296 306 ! 297 307 … … 2245 2255 2246 2256 ok_all_xml = ok_all_xml_omp 2257 ok_lwoff = ok_lwoff_omp 2247 2258 ok_newmicro = ok_newmicro_omp 2248 2259 ok_journe = ok_journe_omp -
LMDZ6/branches/IPSLCM6.0.13/libf/phylmd/phys_output_ctrlout_mod.F90
r3003 r3052 1417 1417 TYPE(ctrl_out), SAVE :: o_lcc3dstra = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11, 11/), & 1418 1418 'lcc3dstra', 'Stratiform cloud liquid fraction', '1', (/ ('', i=1, 10) /)) 1419 TYPE(ctrl_out), SAVE :: o_solbnd = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11, 11/), & 1420 'solbnd', 'Top-of-Atmosphere Solar Insolation for each band', 'W m-2', (/ ('', i=1, 10) /)) 1419 1421 TYPE(ctrl_out), SAVE :: o_stratomask = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11, 11/), & 1420 1422 'stratomask', 'Stratospheric fraction', '1', (/ ('', i=1, 10) /)) -
LMDZ6/branches/IPSLCM6.0.13/libf/phylmd/phys_output_mod.F90
r3003 r3052 137 137 REAL :: ini_heure 138 138 #endif 139 INTEGER :: ISW 140 REAL, DIMENSION(NSW) :: wl1_sun, wl2_sun !wavelength bounds (in um) for SW 141 REAL, DIMENSION(NSW) :: wn1_sun, wn2_sun !wavenumber bounds (in m-1) for SW 142 REAL, DIMENSION(NSW) :: spectband !mean wavenumb. of each sp.band 143 REAL, DIMENSION(NSW,2) :: spbnds_sun !bounds of spectband 139 144 140 145 WRITE(lunout,*) 'Debut phys_output_mod.F90' … … 354 359 CALL wxios_add_vaxis("Alt", & 355 360 levmax(iff) - levmin(iff) + 1, pseudoalt) 361 362 IF (NSW.EQ.6) THEN 363 ! 364 !wl1_sun: minimum bound of wavelength (in um) 365 ! 366 wl1_sun(1)=0.180 367 wl1_sun(2)=0.250 368 wl1_sun(3)=0.440 369 wl1_sun(4)=0.690 370 wl1_sun(5)=1.190 371 wl1_sun(6)=2.380 372 ! 373 !wl2_sun: maximum bound of wavelength (in um) 374 ! 375 wl2_sun(1)=0.250 376 wl2_sun(2)=0.440 377 wl2_sun(3)=0.690 378 wl2_sun(4)=1.190 379 wl2_sun(5)=2.380 380 wl2_sun(6)=4.000 381 ! 382 ELSE IF(NSW.EQ.2) THEN 383 ! 384 !wl1_sun: minimum bound of wavelength (in um) 385 ! 386 wl1_sun(1)=0.250 387 wl1_sun(2)=0.690 388 ! 389 !wl2_sun: maximum bound of wavelength (in um) 390 ! 391 wl2_sun(1)=0.690 392 wl2_sun(2)=4.000 393 ENDIF 394 395 DO ISW=1, NSW 396 wn1_sun(ISW)=1.e+6/wl1_sun(ISW) 397 wn2_sun(ISW)=1.e+6/wl2_sun(ISW) 398 spbnds_sun(ISW,1)=wn2_sun(ISW) 399 spbnds_sun(ISW,2)=wn1_sun(ISW) 400 spectband(ISW)=(wn1_sun(ISW)+wn2_sun(ISW))/2 401 ENDDO 402 ! 403 !!! ajout axe vertical spectband : solar band number 404 CALL wxios_add_vaxis("spectband", NSW, spectband, positif='down') 356 405 ELSE 357 406 ! NMC files -
LMDZ6/branches/IPSLCM6.0.13/libf/phylmd/phys_output_write_mod.F90
r3003 r3052 119 119 o_swtoacf_zero, o_swsrfcf_zero, & 120 120 o_topswai, o_solswai, o_toplwai, o_sollwai, o_scdnc, & 121 o_cldncl, o_reffclws, o_reffclwc, o_s tratomask,&121 o_cldncl, o_reffclws, o_reffclwc, o_solbnd, o_stratomask,& 122 122 o_cldnvi, o_lcc, o_lcc3d, o_lcc3dcon, & 123 123 o_lcc3dstra, o_reffclwtop, o_ec550aer, & … … 389 389 REAL, DIMENSION (klon,klev) :: zx_tmp_fi3d, zpt_conv 390 390 REAL, DIMENSION (klon,klev+1) :: zx_tmp_fi3d1 391 REAL, DIMENSION (klon,NSW) :: zx_tmp_fi3dsp 391 392 CHARACTER (LEN=4) :: bb2 392 393 INTEGER, DIMENSION(nbp_lon*nbp_lat) :: ndex2d … … 1383 1384 CALL histwrite_phy(o_zfull,zphi/RG) 1384 1385 1386 !solbnd begin 1387 #ifdef CPP_RRTM 1388 IF (iflag_rrtm.EQ.1) THEN 1389 IF (vars_defined) THEN 1390 DO ISW=1, NSW 1391 zx_tmp_fi3dsp(:,ISW) = swdn(:,klevp1)*swradcorr(:)*RSUN(ISW) 1392 ENDDO 1393 CALL histwrite_phy(o_solbnd, zx_tmp_fi3dsp) 1394 ENDIF 1395 ENDIF 1396 #endif 1397 1398 !solbnd end 1399 1385 1400 IF (flag_aerosol_strat.EQ.2) THEN 1386 1401 CALL histwrite_phy(o_stratomask, stratomask) -
LMDZ6/branches/IPSLCM6.0.13/libf/phylmd/physiq_mod.F90
r3028 r3052 1105 1105 REAL, dimension(klon, klev) :: cldfrarad ! fraction nuageuse 1106 1106 1107 !lwoff=y : offset LW CRE for radiation code and other schemes 1108 REAL, SAVE :: betalwoff 1109 !OMP THREADPRIVATE(betalwoff) 1110 ! 1107 1111 INTEGER :: nbtr_tmp ! Number of tracer inside concvl 1108 1112 REAL, dimension(klon,klev) :: sh_in ! Specific humidity entering in phytrac … … 1199 1203 pbase=0 1200 1204 !IM 180608 1201 1202 1205 1203 1206 itau_con=0 … … 1406 1409 !>jyg 1407 1410 ENDIF 1408 1409 1410 1411 !IM begin 1411 1412 print*,'physiq: clwcon rnebcon ratqs',clwcon(1,1),rnebcon(1,1) & … … 1791 1792 WRITE(*,*)'lat2_beta=',lat2_beta 1792 1793 WRITE(*,*)'mskocean_beta=',mskocean_beta 1794 1795 !lwoff=y : offset LW CRE for radiation code and other schemes 1796 !lwoff=y : betalwoff=1. 1797 betalwoff=0. 1798 IF (ok_lwoff) THEN 1799 betalwoff=1. 1800 ENDIF 1801 WRITE(*,*)'ok_lwoff=',ok_lwoff 1802 ! 1803 !lwoff=y to begin only sollw and sollwdown are set up to CS values 1804 sollw = sollw + betalwoff * (sollw0 - sollw) 1805 sollwdown(:)= sollwdown(:) + betalwoff *(-1.*ZFLDN0(:,1) - & 1806 sollwdown(:)) 1793 1807 ENDIF 1794 1808 ! … … 3606 3620 ! global 3607 3621 ! 3622 !IM 251017 begin 3623 print*,'physiq betaCRF global zdtime=',zdtime 3624 !IM 251017 end 3608 3625 DO k=1, klev 3609 3626 DO i=1, klon … … 3720 3737 print *,' ->radlwsw, number 1 ' 3721 3738 ENDIF 3739 3722 3740 ! 3723 3741 CALL radlwsw & … … 3736 3754 tau_aero_lw_rrtm, & 3737 3755 cldtaupirad,new_aod, & 3756 ! zqsat, flwcrad, fiwcrad, & 3738 3757 zqsat, flwc, fiwc, & 3739 3758 ref_liq, ref_ice, ref_liq_pi, ref_ice_pi, & … … 3758 3777 ZSWFT0_i, ZFSDN0, ZFSUP0) 3759 3778 3779 !lwoff=y, betalwoff=1. : offset LW CRE for radiation code and other 3780 !schemes 3781 toplw = toplw + betalwoff * (toplw0 - toplw) 3782 sollw = sollw + betalwoff * (sollw0 - sollw) 3783 lwdn = lwdn + betalwoff * (lwdn0 - lwdn) 3784 lwup = lwup + betalwoff * (lwup0 - lwup) 3785 sollwdown(:)= sollwdown(:) + betalwoff *(-1.*ZFLDN0(:,1) - & 3786 sollwdown(:)) 3787 cool = cool + betalwoff * (cool0 - cool) 3788 3760 3789 #ifndef CPP_XIOS 3761 3790 !--OB 30/05/2016 modified 21/10/2016 … … 3808 3837 tau_aero_lw_rrtm, & 3809 3838 cldtaupi,new_aod, & 3839 ! zqsat, flwcrad, fiwcrad, & 3810 3840 zqsat, flwc, fiwc, & 3811 3841 ref_liq, ref_ice, ref_liq_pi, ref_ice_pi, & … … 3829 3859 ZLWFT0_i, ZFLDN0, ZFLUP0, & 3830 3860 ZSWFT0_i, ZFSDN0, ZFSUP0) 3831 endif 3832 ! 3861 endif !ok_4xCO2atm 3833 3862 ENDIF ! aerosol_couple 3834 3863 itaprad = 0 … … 4413 4442 ! ------------------------------------------------- 4414 4443 ! 4415 4416 4444 include "calcul_STDlev.h" 4417 4445 !
Note: See TracChangeset
for help on using the changeset viewer.