Changeset 5605


Ignore:
Timestamp:
Apr 7, 2025, 4:05:50 PM (2 months ago)
Author:
lebasn
Message:

StratAer?: add new diags (4 new longwaves for ext+AOD and SO2 chem loss) + update field/file_def.

Location:
LMDZ6/trunk
Files:
8 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/DefLists/field_def_lmdz.xml

    r5570 r5605  
    10311031
    10321032    <field_group id="fields_strataer_3D" grid_ref="grid_glo_presnivs" operation="average" >
     1033        <field id="ext_strat_443"        long_name="Strat. aerosol extinction coefficient at 443 nm"       unit="1/m" />
    10331034        <field id="ext_strat_550"        long_name="Strat. aerosol extinction coefficient at 550 nm"       unit="1/m" />
    1034         <field id="ext_strat_1020"       long_name="Strat. aerosol extinction coefficient at 1020 nm"      unit="1/m" />
     1035        <field id="ext_strat_670"        long_name="Strat. aerosol extinction coefficient at 670 nm"       unit="1/m" />
     1036        <field id="ext_strat_765"        long_name="Strat. aerosol extinction coefficient at 765 nm"       unit="1/m" />
     1037        <field id="ext_strat_1020"      long_name="Strat. aerosol extinction coefficient at 1020 nm"     unit="1/m" />
     1038        <field id="ext_strat_10um"     long_name="Strat. aerosol extinction coefficient at 10 um"         unit="1/m" />
    10351039        <field id="budg_3D_nucl"         long_name="H2SO4 nucleation mass flux"                            unit="kg(S)/m2/layer/s" />
    10361040        <field id="budg_3D_cond_evap"    long_name="H2SO4 net condensation/evaporation mass flux"          unit="kg(S)/m2/layer/s" />
    … …  
    10421046        <field id="SAD_sulfate"                long_name="SAD WET sulfate aerosols"                        unit="cm2/cm3" />
    10431047        <field id="reff_sulfate"                long_name="Effective    radius of WET sulfate aerosols"              unit="cm" />
    1044         <field id="sulfMMR"                long_name="Sulfate aerosol concentration (dry mass mixing ratio)"                        unit="kg(H2SO4)/kg(air)" />
     1048        <field id="sulfMMR"                long_name="Sulfate aerosol concentration (dry mass mixing ratio)"  unit="kg(H2SO4)/kg(air)" />
     1049        <field id="SO2_CHLM"            long_name="SO2 chemical loss rate"                        unit="mole/cm3/s" />
    10451050        <field id="OCS_lifetime"         long_name="OCS lifetime"                                          unit="s" />
    10461051        <field id="SO2_lifetime"         long_name="SO2 lifetime"                                          unit="s" />
    … …  
    11651170
    11661171    <field_group id="fields_strataer_2D" grid_ref="grid_glo" operation="average">
     1172        <field id="OD443_strat_only"   long_name="Stratospheric Aerosol Optical depth at 443 nm "        unit="1" />
    11671173        <field id="OD550_strat_only"   long_name="Stratospheric Aerosol Optical depth at 550 nm "        unit="1" />
     1174        <field id="OD670_strat_only"   long_name="Stratospheric Aerosol Optical depth at 670 nm "        unit="1" />
     1175        <field id="OD765_strat_only"   long_name="Stratospheric Aerosol Optical depth at 765 nm "        unit="1" />
    11681176        <field id="OD1020_strat_only"  long_name="Stratospheric Aerosol Optical depth at 1020 nm "       unit="1" />
     1177        <field id="OD10um_strat_only"  long_name="Stratospheric Aerosol Optical depth at 10 um "       unit="1" />
    11691178        <field id="surf_PM25_sulf"     long_name="Sulfate PM2.5 concentration at the surface"            unit="ug/m3" />
    11701179        <field id="budg_dep_dry_ocs"   long_name="OCS dry deposition flux"                               unit="kg(S)/m2/s" />

  • LMDZ6/trunk/DefLists/file_def_histdaystrataer_lmdz.xml

    r5150 r5605  
    44     
    55      <field_group grid_ref="grid_out"  level="3">
    6         <field field_ref="OD550_strat_only" level="1" />
    7         <field field_ref="OD1020_strat_only" level="1" />
     6        <field field_ref="OD443_strat_only" level="2" />
     7        <field field_ref="OD550_strat_only" level="1" />
     8        <field field_ref="OD670_strat_only" level="2" />
     9        <field field_ref="OD765_strat_only" level="2" />
     10        <field field_ref="OD1020_strat_only" level="1" />
     11        <field field_ref="OD10um_strat_only" level="2" />
    812        <field field_ref="surf_PM25_sulf" level="2" />
    913        <field field_ref="budg_dep_dry_ocs" level="3" />
    … …  
    2933     
    3034      <field_group grid_ref="grid_out_presnivs" level="10">
    31         <field field_ref="ext_strat_550" level="1" />
    32         <field field_ref="ext_strat_1020" level="5" />
     35        <field field_ref="ext_strat_443"  level="5" />
     36        <field field_ref="ext_strat_550"  level="1" />
     37        <field field_ref="ext_strat_670"  level="5" />
     38        <field field_ref="ext_strat_765"  level="5" />
     39        <field field_ref="ext_strat_1020" level="5" />
     40        <field field_ref="ext_strat_10um" level="5" />
    3341        <field field_ref="budg_3D_nucl" level="10" />
    3442        <field field_ref="budg_3D_cond_evap" level="10" />
    … …  
    4149        <field field_ref="reff_sulfate" level="5" />
    4250        <field field_ref="sulfMMR" level="1" />
     51        <field field_ref="SO2_CHLM" level="3" />
    4352        <field field_ref="OCS_lifetime" level="10" />
    4453        <field field_ref="SO2_lifetime" level="10" />

  • LMDZ6/trunk/DefLists/file_def_histstrataer_lmdz.xml

    r5150 r5605  
    77       
    88        <field_group grid_ref="grid_out"  level="3">
    9           <field field_ref="OD550_strat_only" level="1" />
    10           <field field_ref="OD1020_strat_only" level="1" />
     9          <field field_ref="OD443_strat_only" level="2" />
     10          <field field_ref="OD550_strat_only" level="1" />
     11          <field field_ref="OD670_strat_only" level="2" />
     12          <field field_ref="OD765_strat_only" level="2" />
     13          <field field_ref="OD1020_strat_only" level="1" />
     14          <field field_ref="OD10um_strat_only" level="2" />
    1115          <field field_ref="surf_PM25_sulf" level="1" />
    1216          <field field_ref="budg_dep_dry_ocs" level="3" />
    … …  
    3236       
    3337        <field_group grid_ref="grid_out_presnivs" level="5">
    34           <field field_ref="ext_strat_550"  level="1" />
    35           <field field_ref="ext_strat_1020" level="1" />
     38          <field field_ref="ext_strat_443"  level="2" />
     39          <field field_ref="ext_strat_550"  level="1" />
     40          <field field_ref="ext_strat_670"  level="2" />
     41          <field field_ref="ext_strat_765"  level="2" />
     42          <field field_ref="ext_strat_1020" level="1" />
     43          <field field_ref="ext_strat_10um" level="2" />
    3644          <field field_ref="budg_3D_nucl" level="1" />
    3745          <field field_ref="budg_3D_cond_evap" level="1" />
    … …  
    4452          <field field_ref="reff_sulfate" level="1" />
    4553          <field field_ref="sulfMMR" level="1" />
     54          <field field_ref="SO2_CHLM" level="1" />
    4655          <field field_ref="OCS_lifetime" level="1" />
    4756          <field field_ref="SO2_lifetime" level="1" />
    4857          <field field_ref="vsed_aer" level="2" />
    4958          <field field_ref="f_r_wet" level="1" />
    50           <field field_ref="mass" level="2" />
    51           <field field_ref="temp" level="2" />
    52           <field field_ref="pres" level="2" />
     59          <field field_ref="mass" level="1" />
     60          <field field_ref="temp" level="1" />
     61          <field field_ref="pres" level="1" />
    5362          <field field_ref="h2o"  level="1" />
    5463          <field field_ref="dqch4" level="1" />

  • LMDZ6/trunk/libf/phylmd/StratAer/calcaerosolstrato_rrtm.f90

    r5559 r5605  
    77  USE iniprint_mod_h
    88  USE phys_state_var_mod, ONLY: tau_aero_sw_rrtm, piz_aero_sw_rrtm, cg_aero_sw_rrtm, tau_aero_lw_rrtm
    9   USE phys_local_var_mod, ONLY: mdw, tausum_aero, tausum_strat, tau_strat_550, tau_strat_1020, stratomask
     9  USE phys_local_var_mod, ONLY: mdw, tausum_aero, tausum_strat, tau_strat_443, tau_strat_550, &
     10       tau_strat_670, tau_strat_765, tau_strat_1020, tau_strat_10um, stratomask
    1011  USE aero_mod
    1112  USE dimphy
    … …  
    8687  DO k=1,klev
    8788    IF (stratomask(i,k).GT.0.5) THEN
    88       tausum_strat(i,1)=tausum_strat(i,1)+tau_strat_wave(i,k,2)  !--550 nm
    89       tausum_strat(i,2)=tausum_strat(i,2)+tau_strat_wave(i,k,5)  !--1020 nm
    90       tausum_strat(i,3)=tausum_strat(i,3)+tau_strat_wave(i,k,6)  !--10 um
     89      tausum_strat(i,1)=tausum_strat(i,1)+tau_strat_wave(i,k,1)  !--443 nm
     90      tausum_strat(i,2)=tausum_strat(i,2)+tau_strat_wave(i,k,2)  !--550 nm
     91      tausum_strat(i,3)=tausum_strat(i,3)+tau_strat_wave(i,k,3)  !--670 nm
     92      tausum_strat(i,4)=tausum_strat(i,4)+tau_strat_wave(i,k,4)  !--765 nm
     93      tausum_strat(i,5)=tausum_strat(i,5)+tau_strat_wave(i,k,5)  !--1020 nm
     94      tausum_strat(i,6)=tausum_strat(i,6)+tau_strat_wave(i,k,6)  !--10 um
    9195    ENDIF
    9296  ENDDO
    … …  
    97101    zrho=pplay(i,k)/t_seri(i,k)/RD            !air density in kg/m3
    98102    zdz=(paprs(i,k)-paprs(i,k+1))/zrho/RG     !thickness of layer in m
    99     tau_strat_550(i,k)=tau_strat_wave(i,k,2)/zdz
    100     tau_strat_1020(i,k)=tau_strat_wave(i,k,5)/zdz
     103        tau_strat_443(i,k)=tau_strat_wave(i,k,1)/zdz
     104        tau_strat_550(i,k)=tau_strat_wave(i,k,2)/zdz
     105        tau_strat_670(i,k)=tau_strat_wave(i,k,3)/zdz
     106        tau_strat_765(i,k)=tau_strat_wave(i,k,4)/zdz
     107        tau_strat_1020(i,k)=tau_strat_wave(i,k,5)/zdz
     108        tau_strat_10um(i,k)=tau_strat_wave(i,k,6)/zdz
    101109  ENDDO
    102110  ENDDO

  • LMDZ6/trunk/libf/phylmd/StratAer/so2_to_h2so4.f90

    r5268 r5605  
    99  USE yomcst_mod_h, ONLY : RG, RD
    1010  ! lifetime (sec) et O3_clim (VMR)
    11   USE phys_local_var_mod, ONLY : SO2_lifetime, H2SO4_lifetime, O3_clim, budg_3D_so2_to_h2so4, budg_so2_to_h2so4
     11  USE phys_local_var_mod, ONLY : SO2_lifetime,H2SO4_lifetime,O3_clim,budg_3D_so2_to_h2so4,budg_so2_to_h2so4,SO2_chlm
    1212  USE strataer_local_var_mod, ONLY : flag_OH_reduced, flag_H2SO4_photolysis, flag_min_rreduce
    1313 
    … …  
    3232  budg_3D_so2_to_h2so4(:,:)=0.0
    3333  budg_so2_to_h2so4(:)=0.0
    34 
     34  SO2_chlm(:,:)=0.0
     35 
    3536  DO ilon=1, klon
    3637     DO ilev=1, klev
    … …  
    108109           ! IF (SO2_lifetime(ilon,ilev).GT.0.0 .AND. SO2_lifetime(ilon,ilev).LT.1.E10) THEN
    109110           
     111           SO2_chlm(ilon,ilev) = tr_seri(ilon,ilev,id_SO2_strat)*(1.0-exp(-pdtphys/rreduce)) * &
     112                pplay(ilon,ilev)/(t_seri(ilon,ilev)*1.38e-19)  / pdtphys  !SO2 loss rate [mole.cm-3.s-1]
    110113           
    111114           IF (flag_H2SO4_photolysis) THEN
    … …  
    145148                (paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG/pdtphys
    146149           budg_so2_to_h2so4(ilon)=budg_so2_to_h2so4(ilon)+budg_3D_so2_to_h2so4(ilon,ilev)
     150        ELSE
     151           ! troposphere
     152           ! SO2 tropospheric lifetime (in sec) set to 5 days
     153           rreduce =  5.0*24.0*60.0*60.0
     154           rrate =tr_seri(ilon,ilev,id_SO2_strat)*(1.0-exp(-pdtphys/rreduce))
     155           tr_seri(ilon,ilev,id_SO2_strat)=tr_seri(ilon,ilev,id_SO2_strat) - rrate
     156           SO2_chlm(ilon,ilev) = rrate * &
     157                pplay(ilon,ilev)/(t_seri(ilon,ilev)*1.38e-19)  / pdtphys  !SO2 loss rate [moleccm-3s-1]
    147158        ENDIF
    148159        ! IF (is_strato(ilon,ilev)) THEN

  • LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90

    r5569 r5605  
    745745      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: cg_bin
    746746!$OMP THREADPRIVATE(cg_bin)
     747      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: SO2_chlm
     748!$OMP THREADPRIVATE(SO2_chlm)
     749      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tau_strat_443
     750!$OMP THREADPRIVATE(tau_strat_443)
    747751      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tau_strat_550
    748752!$OMP THREADPRIVATE(tau_strat_550)
     753      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tau_strat_670
     754!$OMP THREADPRIVATE(tau_strat_670)
     755      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tau_strat_765
     756!$OMP THREADPRIVATE(tau_strat_765)
    749757      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tau_strat_1020
    750758!$OMP THREADPRIVATE(tau_strat_1020)
     759      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tau_strat_10um
     760!$OMP THREADPRIVATE(tau_strat_10um)
    751761      REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tausum_strat
    752762!$OMP THREADPRIVATE(tausum_strat)
    … …  
    12541264      ALLOCATE (piz_bin(nbands_sw_rrtm+nbands_lw_rrtm+nwave,nbtr))
    12551265      ALLOCATE (cg_bin(nbands_sw_rrtm+nbands_lw_rrtm+nwave,nbtr))
     1266      ALLOCATE (SO2_chlm(klon,klev))
     1267      ALLOCATE (tau_strat_443(klon,klev))
    12561268      ALLOCATE (tau_strat_550(klon,klev))
     1269      ALLOCATE (tau_strat_670(klon,klev))
     1270      ALLOCATE (tau_strat_765(klon,klev))
    12571271      ALLOCATE (tau_strat_1020(klon,klev))
     1272      ALLOCATE (tau_strat_10um(klon,klev))
    12581273      ALLOCATE (tausum_strat(klon,3))
    12591274      ALLOCATE (budg_dep_dry_ocs(klon))
    … …  
    16521667      DEALLOCATE (piz_bin)
    16531668      DEALLOCATE (cg_bin)
     1669      DEALLOCATE (SO2_chlm)
     1670      DEALLOCATE (tau_strat_443)
    16541671      DEALLOCATE (tau_strat_550)
     1672      DEALLOCATE (tau_strat_670)
     1673      DEALLOCATE (tau_strat_765)
    16551674      DEALLOCATE (tau_strat_1020)
     1675      DEALLOCATE (tau_strat_10um)
    16561676      DEALLOCATE (tausum_strat)
    16571677      DEALLOCATE (surf_PM25_sulf)

  • LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

    r5491 r5605  
    14471447
    14481448!--extinction coefficient
     1449  TYPE(ctrl_out), SAVE :: o_ext_strat_443 = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
     1450    'ext_strat_443', 'Strat. aerosol extinction coefficient at 443 nm', '1/m', (/ ('', i=1, 10) /))
    14491451  TYPE(ctrl_out), SAVE :: o_ext_strat_550 = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
    14501452    'ext_strat_550', 'Strat. aerosol extinction coefficient at 550 nm', '1/m', (/ ('', i=1, 10) /))
     1453  TYPE(ctrl_out), SAVE :: o_ext_strat_670 = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
     1454    'ext_strat_670', 'Strat. aerosol extinction coefficient at 670 nm', '1/m', (/ ('', i=1, 10) /))
     1455  TYPE(ctrl_out), SAVE :: o_ext_strat_765 = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
     1456    'ext_strat_765', 'Strat. aerosol extinction coefficient at 765 nm', '1/m', (/ ('', i=1, 10) /))
    14511457  TYPE(ctrl_out), SAVE :: o_ext_strat_1020 = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
    14521458    'ext_strat_1020', 'Strat. aerosol extinction coefficient at 1020 nm', '1/m', (/ ('', i=1, 10) /))
     1459  TYPE(ctrl_out), SAVE :: o_ext_strat_10um = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
     1460    'ext_strat_10um', 'Strat. aerosol extinction coefficient at 10 um', '1/m', (/ ('', i=1, 10) /))
    14531461!--strat aerosol optical depth
     1462  TYPE(ctrl_out), SAVE :: o_tau_strat_443 = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
     1463    'OD443_strat_only', 'Stratospheric Aerosol Optical depth at 443 nm ', '1', (/ ('', i=1, 10) /))
    14541464  TYPE(ctrl_out), SAVE :: o_tau_strat_550 = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
    14551465    'OD550_strat_only', 'Stratospheric Aerosol Optical depth at 550 nm ', '1', (/ ('', i=1, 10) /))
     1466  TYPE(ctrl_out), SAVE :: o_tau_strat_670 = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
     1467    'OD670_strat_only', 'Stratospheric Aerosol Optical depth at 670 nm ', '1', (/ ('', i=1, 10) /))
     1468  TYPE(ctrl_out), SAVE :: o_tau_strat_765 = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
     1469    'OD765_strat_only', 'Stratospheric Aerosol Optical depth at 765 nm ', '1', (/ ('', i=1, 10) /))
    14561470  TYPE(ctrl_out), SAVE :: o_tau_strat_1020 = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
    14571471    'OD1020_strat_only', 'Stratospheric Aerosol Optical depth at 1020 nm ', '1', (/ ('', i=1, 10) /))
     1472  TYPE(ctrl_out), SAVE :: o_tau_strat_10um = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
     1473    'OD10um_strat_only', 'Stratospheric Aerosol Optical depth at 10 um ', '1', (/ ('', i=1, 10) /))
    14581474  TYPE(ctrl_out), SAVE :: o_SAD_sulfate = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
    14591475    'SAD_sulfate', 'SAD WET sulfate aerosols', 'cm2/cm3', (/ ('', i=1, 10) /))
    … …  
    14671483  TYPE(ctrl_out), SAVE :: o_R2SO4 = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
    14681484    'R2SO4', 'H2SO4 mass fraction in aerosol', '%', (/ ('', i=1, 10) /))
     1485  TYPE(ctrl_out), SAVE :: o_SO2_chlm = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
     1486    'SO2_CHLM', 'SO2 chemical loss rate', 'part/cm3/s', (/ ('', i=1, 10) /))
    14691487  TYPE(ctrl_out), SAVE :: o_OCS_lifetime = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
    14701488    'OCS_lifetime', 'OCS lifetime', 's', (/ ('', i=1, 10) /))

  • LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

    r5585 r5605  
    264264         o_budg_emi_ocs, o_budg_emi_so2, o_budg_emi_h2so4, o_budg_emi_part, &
    265265         o_budg_ocs_to_so2, o_budg_so2_to_h2so4, o_budg_h2so4_to_part, &
    266          o_surf_PM25_sulf, o_ext_strat_550, o_tau_strat_550, &
    267          o_vsed_aer, o_tau_strat_1020, o_ext_strat_1020, o_f_r_wet, &
    268          o_SAD_sulfate, o_reff_sulfate, o_sulfmmr, o_nd_mode, o_sulfmmr_mode
     266         o_surf_PM25_sulf, o_ext_strat_443, o_tau_strat_443, o_ext_strat_550, o_tau_strat_550, &
     267         o_ext_strat_670, o_tau_strat_670, o_ext_strat_765, o_tau_strat_765, o_vsed_aer, &
     268         o_tau_strat_1020, o_ext_strat_1020, o_tau_strat_10um, o_ext_strat_10um, o_f_r_wet, &
     269         o_SAD_sulfate, o_reff_sulfate, o_sulfmmr, o_nd_mode, o_sulfmmr_mode,o_SO2_chlm
    269270
    270271    USE lmdz_lscp_ini, ONLY: ok_poprecip
    … …  
    424425         budg_emi_ocs, budg_emi_so2, budg_emi_h2so4, budg_emi_part, &
    425426         budg_ocs_to_so2, budg_so2_to_h2so4, budg_h2so4_to_part, &
    426          surf_PM25_sulf, tau_strat_550, tausum_strat, &
    427          vsed_aer, tau_strat_1020, f_r_wet, &
    428          SAD_sulfate, reff_sulfate, sulfmmr, nd_mode, sulfmmr_mode
     427         surf_PM25_sulf,  tau_strat_443, tau_strat_550, tau_strat_670, tau_strat_765, &
     428         tausum_strat, vsed_aer, tau_strat_1020, tau_strat_10um, f_r_wet, &
     429         SAD_sulfate, reff_sulfate, sulfmmr, nd_mode, sulfmmr_mode, SO2_chlm
    429430
    430431    USE carbon_cycle_mod, ONLY: fco2_ff, fco2_bb, fco2_land, fco2_ocean
    … …  
    18711872          CALL histwrite_phy(o_vsed_aer, vsed_aer)
    18721873          CALL histwrite_phy(o_f_r_wet, f_r_wet)
     1874          CALL histwrite_phy(o_SO2_chlm, SO2_chlm)
     1875          CALL histwrite_phy(o_ext_strat_443, tau_strat_443)
    18731876          CALL histwrite_phy(o_ext_strat_550, tau_strat_550)
     1877          CALL histwrite_phy(o_ext_strat_670, tau_strat_670)
     1878          CALL histwrite_phy(o_ext_strat_765, tau_strat_765)
    18741879          CALL histwrite_phy(o_ext_strat_1020, tau_strat_1020)
    1875           CALL histwrite_phy(o_tau_strat_550, tausum_strat(:,1))
    1876           CALL histwrite_phy(o_tau_strat_1020, tausum_strat(:,2))
     1880          CALL histwrite_phy(o_ext_strat_10um, tau_strat_10um)
     1881          CALL histwrite_phy(o_tau_strat_443, tausum_strat(:,1))
     1882          CALL histwrite_phy(o_tau_strat_550, tausum_strat(:,2))
     1883          CALL histwrite_phy(o_tau_strat_670, tausum_strat(:,3))
     1884          CALL histwrite_phy(o_tau_strat_765, tausum_strat(:,4))
     1885          CALL histwrite_phy(o_tau_strat_1020, tausum_strat(:,5))
     1886          CALL histwrite_phy(o_tau_strat_10um, tausum_strat(:,6))
    18771887          CALL histwrite_phy(o_SAD_sulfate, SAD_sulfate)
    18781888          CALL histwrite_phy(o_reff_sulfate, reff_sulfate)
Note: See TracChangeset for help on using the changeset viewer.