Changeset 2303


Ignore:
Timestamp:
Jun 17, 2015, 8:06:14 PM (10 years ago)
Author:
jescribano
Message:

Bugs corrections, control vector is now fine mode+coarse mode and seasalt coarse+fine, change in emission scheme parameters, more outputs at 10h30 and 13h30 LT. (Pending correct optical and sedimentation parameters)

Location:
LMDZ5/branches/LMDZ5_SPLA/libf/phylmd
Files:
9 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/aeropt_spl.F

    r2217 r2303  
    200200     .                +ss_dust(2)*tr_seri(i,k,id_codu)*
    201201     .                zdz(i,k)*1.e6
    202       IF(id_scdu>0)  taue670_dustsco(i)=taue550_dustsco(i)+
     202      IF(id_scdu>0)  taue670_dustsco(i)=taue670_dustsco(i)+
    203203     .                ss_dustsco(2)*tr_seri(i,k,id_scdu)*
    204204     .                zdz(i,k)*1.e6
    … …  
    229229     .                +ss_dust(3)*tr_seri(i,k,id_codu)*
    230230     .                zdz(i,k)*1.e6
    231       IF(id_scdu>0)  taue865_dustsco(i)=taue550_dustsco(i)+
     231      IF(id_scdu>0)  taue865_dustsco(i)=taue865_dustsco(i)+
    232232     .                ss_dustsco(3)*tr_seri(i,k,id_scdu)*
    233233     .                zdz(i,k)*1.e6

  • LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/cv3a_uncompress.F90

    r1992 r2303  
    138138      evap1(idcum(i), k) = evap(i, k) !RomP
    139139      ep1(idcum(i), k) = ep(i, k) !RomP
    140       eplamm(idcum(i), k) = eplamm(i, k) !RomP+jyg
     140      !eplamm(idcum(i), k) = eplamm(i, k) !RomP+jyg
     141      eplamm1(idcum(i), k) = eplamm(i, k) !RomP+jyg _ 20150518
    141142      wdtraina1(idcum(i), k) = wdtraina(i, k) !RomP
    142143      wdtrainm1(idcum(i), k) = wdtrainm(i, k) !RomP
    … …  
    167168        sigij1(idcum(i), k, j) = sigij(i, k, j) !RomP
    168169        elij1(idcum(i), k, j) = elij(i, k, j) !RomP
    169         epmlmmm(idcum(i), k, j) = epmlmmm(i, k, j) !RomP+jyg
     170        !epmlmmm(idcum(i), k, j) = epmlmmm(i, k, j) !RomP+jyg
     171        epmlmmm1(idcum(i), k, j) = epmlmmm(i, k, j) !RomP+jyg 20150518
    170172      END DO
    171173    END DO

  • LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/dustemission_mod.F90

    r2217 r2303  
    3838!  real   , parameter :: cd=1.*roa/gravity
    3939! new values
    40 !  logical, parameter :: ok_splatunning=.true.
     40!  logical, parameter :: ok_splatuning=.true.
    4141! Div=3 from S. Alfaro (Sow et al ACPD 2011)
    4242!JE 20150206
    43   integer, parameter :: div1=3.
    44   integer, parameter :: div2=3.
    45   integer, parameter :: div3=3.
     43!  integer, parameter :: div1=3.
     44!  integer, parameter :: div2=3.
     45!  integer, parameter :: div3=3.
     46  integer, parameter :: div1=6.
     47  integer, parameter :: div2=6.
     48  integer, parameter :: div3=6.
    4649  real   , parameter :: e1=3.61/div1
    4750  real   , parameter :: e2=3.52/div2
    … …  
    161164  USE mod_grid_phy_lmdz
    162165  USE mod_phys_lmdz_para
     166  USE indice_sol_mod
    163167
    164168  IMPLICIT NONE
    … …  
    169173  REAL,DIMENSION(klon),     INTENT(IN)     :: xlat
    170174  REAL,DIMENSION(klon),     INTENT(IN)     :: xlon
    171   REAL,DIMENSION(klon),    INTENT(IN)     :: pctsrf
     175  REAL,DIMENSION(klon,nbsrf), INTENT(IN)     :: pctsrf
    172176  REAL,DIMENSION(klon),INTENT(IN)          :: zu10m   ! 10m zonal wind
    173177  REAL,DIMENSION(klon),INTENT(IN)          :: zv10m   ! meridional 10m wind
    … …  
    218222  ENDIF
    219223
    220   CALL adaptdustemission(debutphy,emisbinloc,emdustacc,emdustcoa,emdustsco)
     224  !CALL adaptdustemission(debutphy,emisbinloc,emdustacc,emdustcoa,emdustsco)
     225  CALL adaptdustemission(debutphy,emisbinloc,emdustacc,emdustcoa,emdustsco,maskdust,pctsrf)
    221226  ! output in kg/m2/s
    222227
    … …  
    258263
    259264 SUBROUTINE adaptdustemission(debutphy,emisbinlocal, &
    260                                emdustacc,emdustcoa,emdustsco)
     265                emdustacc,emdustcoa,emdustsco,maskdust,pctsrf)
     266!                               emdustacc,emdustcoa,emdustsco)
    261267
    262268  USE dimphy
    … …  
    265271  USE mod_grid_phy_lmdz
    266272  USE mod_phys_lmdz_para
     273  USE indice_sol_mod
    267274
    268275  IMPLICIT NONE
    … …  
    276283  INTEGER,SAVE ::iminacclow,iminacchigh,imincoalow
    277284  INTEGER,SAVE ::imincoahigh,iminscohigh,iminscolow
     285  INTEGER,DIMENSION(klon) :: maskdust ! where the emissions were calculated
     286  REAL,DIMENSION(klon,nbsrf),     INTENT(IN)     :: pctsrf
    278287!  real,parameter :: sizeacclow=0.03
    279288!  real,parameter :: sizeacchigh=0.5
    … …  
    286295  real,parameter :: sizescolow=6.
    287296  real,parameter :: sizescohigh=30.  ! in micrometers
    288 
    289   real,parameter :: tunningfactor=4.5  ! factor for fine bins!!! important!!
     297!--------------------------------
     298  real,parameter :: tuningfactorfine=1.0  ! factor for fine bins!!! important!!
     299!  real,parameter :: tuningfactorfine=4.5  ! factor for fine bins!!! important!!
     300  real,parameter :: tuningfactorcoa=4.0  ! factor for coarse bins!!! important!!
     301!  real,parameter :: tuningfactorcoa=4.5  ! factor for coarse bins!!! important!!
     302  real,parameter :: tuningfactorsco=4.0  ! factor for supercoarse bins!!! important!!
     303!  real,parameter :: tuningfactorsco=4.5  ! factor for supercoarse bins!!! important!!
     304  real,parameter :: basesumemission= 0.0  !1.e-6  ! emissions to SUM to each land pixel FOR ASSIMILATION ONLY important!!  in mg/m2/s, per bin
     305 !basesumemission = 1.e-6 increase the AOD in about 12%  (0.03 of AOD) ,
     306 !while 1e-8 increase in about 0.12%  (0.003 of AOD)
     307
     308  real,dimension(klon) :: basesumacc,basesumcoa,basesumsco
     309!--------------------------------
    290310!JE20140915  real,parameter :: sizeacclow=0.06
    291311!JE20140915  real,parameter :: sizeacchigh=1.0
    … …  
    317337   call abort_gcm('adaptdustemission', 'Dust range problem',1)
    318338  endif
    319   print *,'ALL DUST BIN: tunning EMISSION factor= ',tunningfactor
     339  print *,'FINE DUST BIN: tuning EMISSION factor= ',tuningfactorfine
     340  print *,'COA DUST BIN: tuning EMISSION factor= ',tuningfactorcoa
     341  print *,'SCO DUST BIN: tuning EMISSION factor= ',tuningfactorsco
     342  print *,'ALL DUST BIN: SUM to the emissions (mg/m2/s) = ',basesumemission
    320343  auxr1=9999.
    321344  auxr2=9999.
    … …  
    386409
    387410! estimate and integrate bins into only accumulation and coarse
     411do k=1,klon
     412  basesumacc(k)=basesumemission*(pctsrf(k,is_ter))*1.e-6 ! from mg/m2/s
     413  basesumcoa(k)=basesumemission*(pctsrf(k,is_ter))*1.e-6
     414  basesumsco(k)=basesumemission*(pctsrf(k,is_ter))*1.e-6
     415enddo
    388416
    389417
    … …  
    395423   auxr1=auxr1+emisbinlocal(k,i)
    396424  enddo
    397   emdustacc(k)= auxr1*tunningfactor
     425  emdustacc(k)=(auxr1 + basesumacc(k))*tuningfactorfine
    398426  do i=imincoalow,imincoahigh-1
    399427    auxr2=auxr2+emisbinlocal(k,i)
    400428  enddo
    401   emdustcoa(k)=auxr2*tunningfactor
     429  emdustcoa(k)=(auxr2 + basesumcoa(k))*tuningfactorcoa
    402430  do i=iminscolow,iminscohigh-1
    403431    auxr3=auxr3+emisbinlocal(k,i)
    404432  enddo
    405   emdustsco(k)=auxr3*tunningfactor
     433  emdustsco(k)=(auxr3 + basesumsco(k))*tuningfactorsco
    406434enddo
     435
     436
     437!do k=1,klon
     438!auxr1=0.0
     439!auxr2=0.0
     440!auxr3=0.0
     441!  do i=iminacclow,iminacchigh-1
     442!   auxr1=auxr1+emisbinlocal(k,i)
     443!  enddo
     444!  emdustacc(k)= auxr1*tuningfactor
     445!  do i=imincoalow,imincoahigh-1
     446!    auxr2=auxr2+emisbinlocal(k,i)
     447!  enddo
     448!  emdustcoa(k)=auxr2*tuningfactor
     449!  do i=iminscolow,iminscohigh-1
     450!    auxr3=auxr3+emisbinlocal(k,i)
     451!  enddo
     452!  emdustsco(k)=auxr3*tuningfactor
     453!enddo
     454!
     455
     456
     457
    407458
    408459!JEdbg<<
    … …  
    462513  USE mod_grid_phy_lmdz
    463514  USE mod_phys_lmdz_para
     515  USE indice_sol_mod
    464516
    465517  IMPLICIT NONE
    … …  
    467519  REAL,DIMENSION(klon),     INTENT(IN)     :: xlat
    468520  REAL,DIMENSION(klon),     INTENT(IN)     :: xlon
    469   REAL,DIMENSION(klon),     INTENT(IN)     :: pctsrf
     521  ! JE20150605<< better to read
     522  ! REAL,DIMENSION(klon),     INTENT(IN)     :: pctsrf
     523  REAL,DIMENSION(klon,nbsrf),     INTENT(IN)     :: pctsrf
     524  ! JE20150605>>
    470525
    471526  !Local
    … …  
    613668      !&      xlat(i).ge.latmin.and.xlat(i).le.latmax    &
    614669      !&      .and.pctsrf(i)>0.5.and.Pini(i,nts)>0.)THEN
    615       IF(pctsrf(i)>0.5.and.Pini(i,nts)>0.)THEN
    616          
     670  ! JE20150605<< easier to read
     671  !    IF(pctsrf(i)>0.5.and.Pini(i,nts)>0.)THEN
     672      IF(pctsrf(i,is_ter)>0.5.and.Pini(i,nts)>0.)THEN
     673  ! JE20150605>>
    617674           sol(i,nts) = solini(i,nts)
    618675             P(i,nts) = Pini(i,nts)

  • LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/phys_output_write_mod.F90

    r2196 r2303  
    2424       diag_aod550_dust, diag_aod670_dust, diag_aod865_dust , &
    2525       diag_aod550_dustsco, diag_aod670_dustsco, diag_aod865_dustsco, &
    26        aod550_aqua, aod670_aqua, aod865_aqua, &
    27        aod550_terra, aod670_terra, aod865_terra, &
     26!       aod550_aqua, aod670_aqua, aod865_aqua, &
     27!       aod550_terra, aod670_terra, aod865_terra, &
     28       aod550_aqua,aod550_tr2_aqua,aod550_ss_aqua,aod550_dust_aqua,aod550_dustsco_aqua,&
     29       aod670_aqua,aod670_tr2_aqua,aod670_ss_aqua,aod670_dust_aqua,aod670_dustsco_aqua,&
     30       aod865_aqua,aod865_tr2_aqua,aod865_ss_aqua,aod865_dust_aqua,aod865_dustsco_aqua,&
     31       aod550_terra,aod550_tr2_terra,aod550_ss_terra,aod550_dust_terra,aod550_dustsco_terra,&
     32       aod670_terra,aod670_tr2_terra,aod670_ss_terra,aod670_dust_terra,aod670_dustsco_terra,&
     33       aod865_terra,aod865_tr2_terra,aod865_ss_terra,aod865_dust_terra,aod865_dustsco_terra,&
    2834       trm01,trm02,trm03,trm04,trm05, &
    2935       sconc01,sconc02,sconc03,sconc04,sconc05, &
    … …  
    204210         o_m1dflux,o_m2dflux,o_m3dflux, &
    205211         o_taue550, &
    206         o_taue670,o_taue865, &
     212         o_taue670,o_taue865, &
    207213         o_taue550_tr2, o_taue670_tr2, o_taue865_tr2, &
    208214         o_taue550_ss,o_taue670_ss, o_taue865_ss, &
    … …  
    211217         o_taue550_aqua, o_taue670_aqua, o_taue865_aqua, &
    212218         o_taue550_terra, o_taue670_terra, o_taue865_terra, &
     219         o_taue550_fine_aqua     ,         o_taue670_fine_aqua     ,  &
     220         o_taue865_fine_aqua     ,         o_taue550_coss_aqua      ,  &
     221         o_taue670_coss_aqua      ,         o_taue865_coss_aqua      ,  &
     222         o_taue550_codu_aqua    ,         o_taue670_codu_aqua    ,  &
     223         o_taue865_codu_aqua    ,         o_taue670_scdu_aqua ,  &
     224         o_taue550_scdu_aqua ,         o_taue865_scdu_aqua ,  &
     225         o_taue550_fine_terra     ,         o_taue670_fine_terra     ,&
     226         o_taue865_fine_terra     ,         o_taue550_coss_terra      ,&
     227         o_taue670_coss_terra      ,         o_taue865_coss_terra      ,&
     228         o_taue550_codu_terra    ,         o_taue670_codu_terra    ,&
     229         o_taue865_codu_terra    ,         o_taue670_scdu_terra ,&
     230         o_taue550_scdu_terra ,         o_taue865_scdu_terra ,&
    213231         o_trm01,o_trm02,o_trm03,o_trm04,o_trm05,&
    214232         o_sconc01,o_sconc02,o_sconc03,o_sconc04,o_sconc05, &

  • LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/physiq.F90

    r2217 r2303  
    36423642  nbreg_wstardust=1
    36433643  fileregionsdimsind='regions_ind_meta'
    3644 !  fileregionsdimsdust='regions_dustacc_meta'
    3645   fileregionsdimsdust='regions_dust_meta'
     3644  fileregionsdimsdust='regions_dustacc_meta'
     3645!  fileregionsdimsdust='regions_dust_meta'
    36463646  fileregionsdimsbb='regions_bb_meta'
    36473647  fileregionsdimswstar='regions_pwstarwake_meta'
    … …  
    37053705  paramname_dustcoa='dustcoasco'
    37063706  paramname_dustsco='dustcoasco'
    3707   paramname_dustacc='dust'
    3708   paramname_dustcoa='dust'
    3709   paramname_dustsco='dust'
     3707!  paramname_dustacc='dust'
     3708!  paramname_dustcoa='dust'
     3709!  paramname_dustsco='dust'
    37103710  paramname_wstarBL='pwstarbl'
    37113711  paramname_wstarWAKE='pwstarwake'
    3712   paramname_ssacc='ss'
    3713   paramname_sscoa='ss'
     3712  paramname_ssacc='ssacc'
     3713  paramname_sscoa='sscoa'
    37143714
    37153715  filescaleparams='modvalues.nc'
    … …  
    37393739  print *, 'wstardustBL', param_wstarBLperregion
    37403740  print *, 'wstardustWAKE', param_wstarWAKEperregion
    3741   print *, 'dustsco', scale_param_dustsco
    37423741  print *, 'ff', scale_param_ff 
    37433742  print *, 'bb', scale_param_bb 

  • LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/phytracr_spl_mod.F90

    r2217 r2303  
    7373
    7474
     75      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod550_terra  ! AOD at terra overpass time ( 10.30 local hour)
     76      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod550_tr2_terra  ! AOD at terra overpass time ( 10.30 local hour)
     77      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod550_ss_terra  ! AOD at terra overpass time ( 10.30 local hour)
     78      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod550_dust_terra  ! AOD at terra overpass time ( 10.30 local hour)
     79      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod550_dustsco_terra  ! AOD at terra overpass time ( 10.30 local hour)
     80      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod670_terra  ! AOD at terra overpass time ( 10.30 local hour)
     81      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod670_tr2_terra  ! AOD at terra overpass time ( 10.30 local hour)
     82      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod670_ss_terra  ! AOD at terra overpass time ( 10.30 local hour)
     83      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod670_dust_terra  ! AOD at terra overpass time ( 10.30 local hour)
     84      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod670_dustsco_terra  ! AOD at terra overpass time ( 10.30 local hour)
     85      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod865_terra  ! AOD at terra overpass time ( 10.30 local hour)
     86      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod865_tr2_terra  ! AOD at terra overpass time ( 10.30 local hour)
     87      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod865_ss_terra  ! AOD at terra overpass time ( 10.30 local hour)
     88      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod865_dust_terra  ! AOD at terra overpass time ( 10.30 local hour)
     89      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod865_dustsco_terra  ! AOD at terra overpass time ( 10.30 local hour)
     90
     91
    7592      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod550_aqua  ! AOD at aqua overpass time ( 13.30 local hour)
    76       REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod550_terra ! AOD at terra overpass time ( 10.30 local hour)
     93      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod550_tr2_aqua  ! AOD at aqua overpass time ( 13.30 local hour)
     94      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod550_ss_aqua  ! AOD at aqua overpass time ( 13.30 local hour)
     95      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod550_dust_aqua  ! AOD at aqua overpass time ( 13.30 local hour)
     96      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod550_dustsco_aqua  ! AOD at aqua overpass time ( 13.30 local hour)
    7797      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod670_aqua  ! AOD at aqua overpass time ( 13.30 local hour)
    78       REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod670_terra ! AOD at terra overpass time ( 10.30 local hour)
     98      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod670_tr2_aqua  ! AOD at aqua overpass time ( 13.30 local hour)
     99      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod670_ss_aqua  ! AOD at aqua overpass time ( 13.30 local hour)
     100      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod670_dust_aqua  ! AOD at aqua overpass time ( 13.30 local hour)
     101      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod670_dustsco_aqua  ! AOD at aqua overpass time ( 13.30 local hour)
    79102      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod865_aqua  ! AOD at aqua overpass time ( 13.30 local hour)
    80       REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod865_terra ! AOD at terra overpass time ( 10.30 local hour)
    81 
    82 !$OMP THREADPRIVATE(aod550_aqua,aod550_terra,aod670_aqua,aod670_terra)
    83 !$OMP THREADPRIVATE(aod865_aqua,aod865_terra)
     103      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod865_tr2_aqua  ! AOD at aqua overpass time ( 13.30 local hour)
     104      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod865_ss_aqua  ! AOD at aqua overpass time ( 13.30 local hour)
     105      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod865_dust_aqua  ! AOD at aqua overpass time ( 13.30 local hour)
     106      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: aod865_dustsco_aqua  ! AOD at aqua overpass time ( 13.30 local hour)
     107
     108!$OMP THREADPRIVATE(aod550_aqua,aod550_tr2_aqua,aod550_ss_aqua,aod550_dust_aqua,aod550_dustsco_aqua)
     109!$OMP THREADPRIVATE(aod670_aqua,aod670_tr2_aqua,aod670_ss_aqua,aod670_dust_aqua,aod670_dustsco_aqua)
     110!$OMP THREADPRIVATE(aod865_aqua,aod865_tr2_aqua,aod865_ss_aqua,aod865_dust_aqua,aod865_dustsco_aqua)
     111!$OMP THREADPRIVATE(aod550_terra,aod550_tr2_terra,aod550_ss_terra,aod550_dust_terra,aod550_dustsco_terra)
     112!$OMP THREADPRIVATE(aod670_terra,aod670_tr2_terra,aod670_ss_terra,aod670_dust_terra,aod670_dustsco_terra)
     113!$OMP THREADPRIVATE(aod865_terra,aod865_tr2_terra,aod865_ss_terra,aod865_dust_terra,aod865_dustsco_terra)
     114
     115
    84116      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: sconc01 ! surface concentration
    85117      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: trm01   ! burden
    … …  
    378410
    379411!JE20141224 >>
    380       REAL, DIMENSION(klon_glo) :: aod550_aqua_glo  ! AOD at aqua overpass time ( 13.30 local hour)
    381       REAL, DIMENSION(klon_glo) :: aod550_terra_glo ! AOD at terra overpass time ( 10.30 local hour)
    382       REAL, DIMENSION(klon_glo) :: aod670_aqua_glo  ! AOD at aqua overpass time ( 13.30 local hour)
    383       REAL, DIMENSION(klon_glo) :: aod670_terra_glo ! AOD at terra overpass time ( 10.30 local hour)
    384       REAL, DIMENSION(klon_glo) :: aod865_aqua_glo  ! AOD at aqua overpass time ( 13.30 local hour)
    385       REAL, DIMENSION(klon_glo) :: aod865_terra_glo ! AOD at terra overpass time ( 10.30 local hour)
     412!JE20150518<<
     413      REAL, DIMENSION(klon_glo)  :: aod550_terra_glo  ! AOD at terra overpass time ( 10.30 local hour)
     414      REAL, DIMENSION(klon_glo)  :: aod550_tr2_terra_glo  ! AOD at terra overpass time ( 10.30 local hour)
     415      REAL, DIMENSION(klon_glo)  :: aod550_ss_terra_glo  ! AOD at terra overpass time ( 10.30 local hour)
     416      REAL, DIMENSION(klon_glo)  :: aod550_dust_terra_glo  ! AOD at terra overpass time ( 10.30 local hour)
     417      REAL, DIMENSION(klon_glo)  :: aod550_dustsco_terra_glo  ! AOD at terra overpass time ( 10.30 local hour)
     418      REAL, DIMENSION(klon_glo)  :: aod670_terra_glo  ! AOD at terra overpass time ( 10.30 local hour)
     419      REAL, DIMENSION(klon_glo)  :: aod670_tr2_terra_glo  ! AOD at terra overpass time ( 10.30 local hour)
     420      REAL, DIMENSION(klon_glo)  :: aod670_ss_terra_glo  ! AOD at terra overpass time ( 10.30 local hour)
     421      REAL, DIMENSION(klon_glo)  :: aod670_dust_terra_glo  ! AOD at terra overpass time ( 10.30 local hour)
     422      REAL, DIMENSION(klon_glo)  :: aod670_dustsco_terra_glo  ! AOD at terra overpass time ( 10.30 local hour)
     423      REAL, DIMENSION(klon_glo)  :: aod865_terra_glo  ! AOD at terra overpass time ( 10.30 local hour)
     424      REAL, DIMENSION(klon_glo)  :: aod865_tr2_terra_glo  ! AOD at terra overpass time ( 10.30 local hour)
     425      REAL, DIMENSION(klon_glo)  :: aod865_ss_terra_glo  ! AOD at terra overpass time ( 10.30 local hour)
     426      REAL, DIMENSION(klon_glo)  :: aod865_dust_terra_glo  ! AOD at terra overpass time ( 10.30 local hour)
     427      REAL, DIMENSION(klon_glo)  :: aod865_dustsco_terra_glo  ! AOD at terra overpass time ( 10.30 local hour)
     428
     429      REAL, DIMENSION(klon_glo)  :: aod550_aqua_glo  ! AOD at aqua overpass time ( 13.30 local hour)
     430      REAL, DIMENSION(klon_glo)  :: aod550_tr2_aqua_glo  ! AOD at aqua overpass time ( 13.30 local hour)
     431      REAL, DIMENSION(klon_glo)  :: aod550_ss_aqua_glo  ! AOD at aqua overpass time ( 13.30 local hour)
     432      REAL, DIMENSION(klon_glo)  :: aod550_dust_aqua_glo  ! AOD at aqua overpass time ( 13.30 local hour)
     433      REAL, DIMENSION(klon_glo)  :: aod550_dustsco_aqua_glo  ! AOD at aqua overpass time ( 13.30 local hour)
     434      REAL, DIMENSION(klon_glo)  :: aod670_aqua_glo  ! AOD at aqua overpass time ( 13.30 local hour)
     435      REAL, DIMENSION(klon_glo)  :: aod670_tr2_aqua_glo  ! AOD at aqua overpass time ( 13.30 local hour)
     436      REAL, DIMENSION(klon_glo)  :: aod670_ss_aqua_glo  ! AOD at aqua overpass time ( 13.30 local hour)
     437      REAL, DIMENSION(klon_glo)  :: aod670_dust_aqua_glo  ! AOD at aqua overpass time ( 13.30 local hour)
     438      REAL, DIMENSION(klon_glo)  :: aod670_dustsco_aqua_glo  ! AOD at aqua overpass time ( 13.30 local hour)
     439      REAL, DIMENSION(klon_glo)  :: aod865_aqua_glo  ! AOD at aqua overpass time ( 13.30 local hour)
     440      REAL, DIMENSION(klon_glo)  :: aod865_tr2_aqua_glo  ! AOD at aqua overpass time ( 13.30 local hour)
     441      REAL, DIMENSION(klon_glo)  :: aod865_ss_aqua_glo  ! AOD at aqua overpass time ( 13.30 local hour)
     442      REAL, DIMENSION(klon_glo)  :: aod865_dust_aqua_glo  ! AOD at aqua overpass time ( 13.30 local hour)
     443      REAL, DIMENSION(klon_glo)  :: aod865_dustsco_aqua_glo  ! AOD at aqua overpass time ( 13.30 local hour)
     444!!!!!!!!!!!!!
     445!JE20150518>>
     446
     447
    386448
    387449
    … …  
    9931055
    9941056
    995       lmt_dms(:)=0.0
    996       aux_var2(:)=0.0
    997       aux_var3(:,:)=0.0
    998       source_tr(:,:)=0.0
    999       flux_tr(:,:)=0.0
    1000       flux_sparam_bb(:)=0.0
    1001       flux_sparam_ff(:)=0.0
    1002       flux_sparam_ind(:)=0.0
    1003       flux_sparam_ddfine(:)=0.0
    1004       flux_sparam_ddcoa(:)=0.0
    1005       flux_sparam_ddsco(:)=0.0
    1006       flux_sparam_ssfine(:)=0.0
    1007       flux_sparam_sscoa(:)=0.0
    10081057     
    10091058!JE20140507      d_tr_dry(:,:)=0.0
    … …  
    17401789        ALLOCATE(iregion_wstardust(klon))
    17411790
     1791!JE20150518<<
    17421792       ALLOCATE(masque_aqua(klon)) 
    17431793       ALLOCATE(masque_terra(klon)) 
    1744        ALLOCATE(aod550_aqua(klon)) 
    1745        ALLOCATE(aod550_terra(klon)) 
    1746        ALLOCATE(aod670_aqua(klon)) 
    1747        ALLOCATE(aod670_terra(klon)) 
    1748        ALLOCATE(aod865_aqua(klon)) 
    1749        ALLOCATE(aod865_terra(klon)) 
     1794!       ALLOCATE(aod550_aqua(klon)) 
     1795!       ALLOCATE(aod550_terra(klon)) 
     1796!       ALLOCATE(aod670_aqua(klon)) 
     1797!       ALLOCATE(aod670_terra(klon)) 
     1798!       ALLOCATE(aod865_aqua(klon)) 
     1799!       ALLOCATE(aod865_terra(klon))
     1800
     1801      ALLOCATE( aod550_terra(klon)) 
     1802      ALLOCATE( aod550_tr2_terra(klon)) 
     1803      ALLOCATE( aod550_ss_terra(klon))   
     1804      ALLOCATE( aod550_dust_terra(klon))   
     1805      ALLOCATE( aod550_dustsco_terra(klon))   
     1806      ALLOCATE( aod670_terra(klon))   
     1807      ALLOCATE( aod670_tr2_terra(klon)) 
     1808      ALLOCATE( aod670_ss_terra(klon)) 
     1809      ALLOCATE( aod670_dust_terra(klon)) 
     1810      ALLOCATE( aod670_dustsco_terra(klon)) 
     1811      ALLOCATE( aod865_terra(klon))   
     1812      ALLOCATE( aod865_tr2_terra(klon)) 
     1813      ALLOCATE( aod865_ss_terra(klon)) 
     1814      ALLOCATE( aod865_dust_terra(klon)) 
     1815      ALLOCATE( aod865_dustsco_terra(klon)) 
     1816
     1817      ALLOCATE( aod550_aqua(klon)) 
     1818      ALLOCATE( aod550_tr2_aqua(klon)) 
     1819      ALLOCATE( aod550_ss_aqua(klon))   
     1820      ALLOCATE( aod550_dust_aqua(klon))   
     1821      ALLOCATE( aod550_dustsco_aqua(klon))   
     1822      ALLOCATE( aod670_aqua(klon))   
     1823      ALLOCATE( aod670_tr2_aqua(klon)) 
     1824      ALLOCATE( aod670_ss_aqua(klon)) 
     1825      ALLOCATE( aod670_dust_aqua(klon)) 
     1826      ALLOCATE( aod670_dustsco_aqua(klon)) 
     1827      ALLOCATE( aod865_aqua(klon))   
     1828      ALLOCATE( aod865_tr2_aqua(klon)) 
     1829      ALLOCATE( aod865_ss_aqua(klon)) 
     1830      ALLOCATE( aod865_dust_aqua(klon)) 
     1831      ALLOCATE( aod865_dustsco_aqua(klon)) 
     1832 
    17501833
    17511834       masque_aqua(:)=0
    17521835       masque_terra(:)=0
    1753        aod550_aqua(:)=0.
    1754        aod550_terra(:)=0.
    1755        aod670_aqua(:)=0.
    1756        aod670_terra(:)=0.
    1757        aod865_aqua(:)=0.
    1758        aod865_terra(:)=0.
     1836!       aod550_aqua(:)=0.
     1837!       aod550_terra(:)=0.
     1838!       aod670_aqua(:)=0.
     1839!       aod670_terra(:)=0.
     1840!       aod865_aqua(:)=0.
     1841!       aod865_terra(:)=0.
     1842
     1843      aod550_terra(:)=0. 
     1844      aod550_tr2_terra(:)=0. 
     1845      aod550_ss_terra(:)=0.   
     1846      aod550_dust_terra(:)=0.   
     1847      aod550_dustsco_terra(:)=0.   
     1848      aod670_terra(:)=0.   
     1849      aod670_tr2_terra(:)=0. 
     1850      aod670_ss_terra(:)=0. 
     1851      aod670_dust_terra(:)=0. 
     1852      aod670_dustsco_terra(:)=0. 
     1853      aod865_terra(:)=0.   
     1854      aod865_tr2_terra(:)=0. 
     1855      aod865_ss_terra(:)=0. 
     1856      aod865_dust_terra(:)=0. 
     1857      aod865_dustsco_terra(:)=0. 
     1858      aod550_aqua(:)=0. 
     1859      aod550_tr2_aqua(:)=0. 
     1860      aod550_ss_aqua(:)=0.   
     1861      aod550_dust_aqua(:)=0.   
     1862      aod550_dustsco_aqua(:)=0.   
     1863      aod670_aqua(:)=0.   
     1864      aod670_tr2_aqua(:)=0. 
     1865      aod670_ss_aqua(:)=0. 
     1866      aod670_dust_aqua(:)=0. 
     1867      aod670_dustsco_aqua(:)=0. 
     1868      aod865_aqua(:)=0.   
     1869      aod865_tr2_aqua(:)=0. 
     1870      aod865_ss_aqua(:)=0. 
     1871      aod865_dust_aqua(:)=0. 
     1872      aod865_dustsco_aqua(:)=0. 
     1873!JE20150518>>
     1874
     1875
     1876
     1877
     1878
    17591879!
    17601880!Config Key  = iflag_lscav
    … …  
    17851905!$OMP END MASTER
    17861906!$OMP BARRIER
    1787 
     1907       ENDIF ! debutphy
     1908     
    17881909      lmt_dms(:)=0.0
    17891910      aux_var2(:)=0.0
    … …  
    17991920      flux_sparam_ssfine(:)=0.0
    18001921      flux_sparam_sscoa(:)=0.0
    1801        ENDIF ! debutphy
    1802      
    18031922
    18041923! initialiation for time computation
    … …  
    20542173!      ENDIF
    20552174!
    2056       !c_FullName1='regions_dustacc'
    2057       c_FullName1='regions_dust'
     2175      c_FullName1='regions_dustacc'
     2176      !c_FullName1='regions_dust'
    20582177      call readregions_spl(iregion_dust,c_FullName1)
    20592178      c_FullName1='regions_ind'
    … …  
    20812200      IF (is_mpi_root .AND. is_omp_root) THEN
    20822201      DO k=1,klon_glo
    2083         WRITE(25,'(i)') iauxklon_glo(k)
     2202        WRITE(25,'(i10)') iauxklon_glo(k)
    20842203      ENDDO
    20852204      ENDIF ! mpi root
    … …  
    20902209      IF (is_mpi_root .AND. is_omp_root) THEN
    20912210      DO k=1,klon_glo
    2092         WRITE(55,'(i)') iauxklon_glo(k)
     2211        WRITE(55,'(i10)') iauxklon_glo(k)
    20932212      ENDDO
    20942213      ENDIF ! mpi root
    … …  
    20992218      IF (is_mpi_root .AND. is_omp_root) THEN
    21002219      DO k=1,klon_glo
    2101         WRITE(75,'(i)') iauxklon_glo(k)
     2220        WRITE(75,'(i10)') iauxklon_glo(k)
    21022221      ENDDO
    21032222      ENDIF ! mpi root
    … …  
    21082227      IF (is_mpi_root .AND. is_omp_root) THEN
    21092228      DO k=1,klon_glo
    2110         WRITE(95,'(i)') iauxklon_glo(k)
     2229        WRITE(95,'(i10)') iauxklon_glo(k)
    21112230      ENDDO
    21122231      ENDIF ! mpi root
    … …  
    32913410                        id_prec, id_fine, id_coss, id_codu, id_scdu,     &
    32923411                        ok_chimeredust,                                 &
    3293                     diff_aod550_tot,diag_aod670_tot,diag_aod865_tot,     &
    3294                     diff_aod550_tr2,diag_aod670_tr2,diag_aod865_tr2,     &
    3295                     diag_aod550_ss,diag_aod670_ss,diag_aod865_ss,        &
     3412                    diff_aod550_tot, diag_aod670_tot, diag_aod865_tot,     &
     3413                    diff_aod550_tr2, diag_aod670_tr2, diag_aod865_tr2,     &
     3414                    diag_aod550_ss,  diag_aod670_ss,  diag_aod865_ss,        &
    32963415                    diag_aod550_dust,diag_aod670_dust,diag_aod865_dust,  &
    32973416           diag_aod550_dustsco,diag_aod670_dustsco,diag_aod865_dustsco) 
    … …  
    33163435      CALL satellite_out_spla(jD_cur,jH_cur,pdtphys,xlat,xlon,   &
    33173436                              masque_aqua_cur, masque_terra_cur )
     3437      IF (jH_cur-pdtphys/86400. .LT. 0.) THEN
     3438       !new utc day: put in 0 everything
     3439!JE20150518<<
     3440!       aod550_aqua(:) =0.
     3441!       aod550_terra(:) =0.
     3442!       aod670_aqua(:) =0.
     3443!       aod670_terra(:) =0.
     3444!       aod865_aqua(:) =0.
     3445!       aod865_terra(:) =0.
     3446       masque_aqua(:) =0
     3447       masque_terra(:) =0
     3448       aod550_terra(:)=0. 
     3449       aod550_tr2_terra(:)=0. 
     3450       aod550_ss_terra(:)=0.   
     3451       aod550_dust_terra(:)=0.   
     3452       aod550_dustsco_terra(:)=0.   
     3453       aod670_terra(:)=0.   
     3454       aod670_tr2_terra(:)=0. 
     3455       aod670_ss_terra(:)=0. 
     3456       aod670_dust_terra(:)=0. 
     3457       aod670_dustsco_terra(:)=0. 
     3458       aod865_terra(:)=0.   
     3459       aod865_tr2_terra(:)=0. 
     3460       aod865_ss_terra(:)=0. 
     3461       aod865_dust_terra(:)=0. 
     3462       aod865_dustsco_terra(:)=0. 
     3463       aod550_aqua(:)=0. 
     3464       aod550_tr2_aqua(:)=0. 
     3465       aod550_ss_aqua(:)=0.   
     3466       aod550_dust_aqua(:)=0.   
     3467       aod550_dustsco_aqua(:)=0.   
     3468       aod670_aqua(:)=0.   
     3469       aod670_tr2_aqua(:)=0. 
     3470       aod670_ss_aqua(:)=0. 
     3471       aod670_dust_aqua(:)=0. 
     3472       aod670_dustsco_aqua(:)=0. 
     3473       aod865_aqua(:)=0.   
     3474       aod865_tr2_aqua(:)=0. 
     3475       aod865_ss_aqua(:)=0. 
     3476       aod865_dust_aqua(:)=0. 
     3477       aod865_dustsco_aqua(:)=0. 
     3478!JE20150518>>
     3479      ENDIF
    33183480
    33193481      DO i=1,klon
    3320          aod550_aqua(i)=aod550_aqua(i)+   &
     3482!         aod550_aqua(i)=aod550_aqua(i)+   &
     3483!                       masque_aqua_cur(i)*diff_aod550_tot(i)
     3484!         aod670_aqua(i)=aod670_aqua(i)+   &
     3485!                        masque_aqua_cur(i)*diag_aod670_tot(i)
     3486!         aod865_aqua(i)=aod865_aqua(i)+   &
     3487!                       masque_aqua_cur(i)*diag_aod865_tot(i)
     3488
     3489       aod550_terra(i)=aod550_terra(i)+   &
     3490                       masque_terra_cur(i)*diff_aod550_tot(i)
     3491       aod550_tr2_terra(i)= aod550_tr2_terra(i)+ &
     3492                       masque_terra_cur(i)*diff_aod550_tr2(i)
     3493       aod550_ss_terra(i)=aod550_ss_terra(i) + &
     3494                       masque_terra_cur(i)*diag_aod550_ss(i)
     3495       aod550_dust_terra(i)=  aod550_dust_terra(i) + &
     3496                       masque_terra_cur(i)*diag_aod550_dust(i)
     3497       aod550_dustsco_terra(i)= aod550_dustsco_terra(i) + &
     3498                       masque_terra_cur(i)*diag_aod550_dustsco(i)
     3499       aod670_terra(i)=aod670_terra(i)+   &
     3500                       masque_terra_cur(i)*diag_aod670_tot(i)
     3501       aod670_tr2_terra(i)= aod670_tr2_terra(i)+ &
     3502                       masque_terra_cur(i)*diag_aod670_tr2(i)
     3503       aod670_ss_terra(i)=aod670_ss_terra(i) + &
     3504                       masque_terra_cur(i)*diag_aod670_ss(i)
     3505       aod670_dust_terra(i)=  aod670_dust_terra(i) + &
     3506                       masque_terra_cur(i)*diag_aod670_dust(i)
     3507       aod670_dustsco_terra(i)= aod670_dustsco_terra(i) + &
     3508                       masque_terra_cur(i)*diag_aod670_dustsco(i)
     3509       aod865_terra(i)=aod865_terra(i)+   &
     3510                       masque_terra_cur(i)*diag_aod865_tot(i)
     3511       aod865_tr2_terra(i)= aod865_tr2_terra(i)+ &
     3512                       masque_terra_cur(i)*diag_aod865_tr2(i)
     3513       aod865_ss_terra(i)=aod865_ss_terra(i) + &
     3514                       masque_terra_cur(i)*diag_aod865_ss(i)
     3515       aod865_dust_terra(i)=  aod865_dust_terra(i) + &
     3516                       masque_terra_cur(i)*diag_aod865_dust(i)
     3517       aod865_dustsco_terra(i)= aod865_dustsco_terra(i) + &
     3518                       masque_terra_cur(i)*diag_aod865_dustsco(i)
     3519
     3520
     3521
     3522       aod550_aqua(i)=aod550_aqua(i)+   &
    33213523                       masque_aqua_cur(i)*diff_aod550_tot(i)
    3322          aod670_aqua(i)=aod670_aqua(i)+   &
    3323                         masque_aqua_cur(i)*diag_aod670_tot(i)
    3324          aod865_aqua(i)=aod865_aqua(i)+   &
     3524       aod550_tr2_aqua(i)= aod550_tr2_aqua(i)+ &
     3525                       masque_aqua_cur(i)*diff_aod550_tr2(i)
     3526       aod550_ss_aqua(i)=aod550_ss_aqua(i) + &
     3527                       masque_aqua_cur(i)*diag_aod550_ss(i)
     3528       aod550_dust_aqua(i)=  aod550_dust_aqua(i) + &
     3529                       masque_aqua_cur(i)*diag_aod550_dust(i)
     3530       aod550_dustsco_aqua(i)= aod550_dustsco_aqua(i) + &
     3531                       masque_aqua_cur(i)*diag_aod550_dustsco(i)
     3532       aod670_aqua(i)=aod670_aqua(i)+   &
     3533                       masque_aqua_cur(i)*diag_aod670_tot(i)
     3534       aod670_tr2_aqua(i)= aod670_tr2_aqua(i)+ &
     3535                       masque_aqua_cur(i)*diag_aod670_tr2(i)
     3536       aod670_ss_aqua(i)=aod670_ss_aqua(i) + &
     3537                       masque_aqua_cur(i)*diag_aod670_ss(i)
     3538       aod670_dust_aqua(i)=  aod670_dust_aqua(i) + &
     3539                       masque_aqua_cur(i)*diag_aod670_dust(i)
     3540       aod670_dustsco_aqua(i)= aod670_dustsco_aqua(i) + &
     3541                       masque_aqua_cur(i)*diag_aod670_dustsco(i)
     3542       aod865_aqua(i)=aod865_aqua(i)+   &
    33253543                       masque_aqua_cur(i)*diag_aod865_tot(i)
     3544       aod865_tr2_aqua(i)= aod865_tr2_aqua(i)+ &
     3545                       masque_aqua_cur(i)*diag_aod865_tr2(i)
     3546       aod865_ss_aqua(i)=aod865_ss_aqua(i) + &
     3547                       masque_aqua_cur(i)*diag_aod865_ss(i)
     3548       aod865_dust_aqua(i)=  aod865_dust_aqua(i) + &
     3549                       masque_aqua_cur(i)*diag_aod865_dust(i)
     3550       aod865_dustsco_aqua(i)= aod865_dustsco_aqua(i) + &
     3551                       masque_aqua_cur(i)*diag_aod865_dustsco(i)
     3552!         aod550_terra(i)=aod550_terra(i)+  &
     3553!                       masque_terra_cur(i)*diff_aod550_tot(i)
     3554!         aod670_terra(i)=aod670_terra(i)+  &
     3555!                       masque_terra_cur(i)*diag_aod670_tot(i)
     3556!         aod865_terra(i)=aod865_terra(i)+   &
     3557!                       masque_terra_cur(i)*diag_aod865_tot(i)
    33263558         masque_aqua(i)=masque_aqua(i)+masque_aqua_cur(i)
    3327          aod550_terra(i)=aod550_terra(i)+  &
    3328                        masque_terra_cur(i)*diff_aod550_tot(i)
    3329          aod670_terra(i)=aod670_terra(i)+  &
    3330                        masque_terra_cur(i)*diag_aod670_tot(i)
    3331          aod865_terra(i)=aod865_terra(i)+   &
    3332                        masque_terra_cur(i)*diag_aod865_tot(i)
    33333559         masque_terra(i)=masque_terra(i)+masque_terra_cur(i)
    33343560      ENDDO
    … …  
    33413567                   aod670_aqua(i)=aod670_aqua(i)/masque_aqua(i)
    33423568                   aod865_aqua(i)=aod865_aqua(i)/masque_aqua(i)
     3569                   aod550_tr2_aqua(i)=aod550_tr2_aqua(i)/masque_aqua(i)
     3570                   aod670_tr2_aqua(i)=aod670_tr2_aqua(i)/masque_aqua(i)
     3571                   aod865_tr2_aqua(i)=aod865_tr2_aqua(i)/masque_aqua(i)
     3572                   aod550_ss_aqua(i)=aod550_ss_aqua(i)/masque_aqua(i)
     3573                   aod670_ss_aqua(i)=aod670_ss_aqua(i)/masque_aqua(i)
     3574                   aod865_ss_aqua(i)=aod865_ss_aqua(i)/masque_aqua(i)
     3575                   aod550_dust_aqua(i)=aod550_dust_aqua(i)/masque_aqua(i)
     3576                   aod670_dust_aqua(i)=aod670_dust_aqua(i)/masque_aqua(i)
     3577                   aod865_dust_aqua(i)=aod865_dust_aqua(i)/masque_aqua(i)
     3578                   aod550_dustsco_aqua(i)=aod550_dustsco_aqua(i)/masque_aqua(i)
     3579                   aod670_dustsco_aqua(i)=aod670_dustsco_aqua(i)/masque_aqua(i)
     3580                   aod865_dustsco_aqua(i)=aod865_dustsco_aqua(i)/masque_aqua(i)
    33433581               ELSE
    33443582                   aod550_aqua(i) = -999.
    33453583                   aod670_aqua(i) = -999.
    33463584                   aod865_aqua(i) = -999.
     3585                   aod550_tr2_aqua(i)= -999.
     3586                   aod670_tr2_aqua(i)= -999.
     3587                   aod865_tr2_aqua(i)= -999.
     3588                   aod550_ss_aqua(i)= -999.
     3589                   aod670_ss_aqua(i)= -999.
     3590                   aod865_ss_aqua(i)= -999.
     3591                   aod550_dust_aqua(i)= -999.
     3592                   aod670_dust_aqua(i)= -999.
     3593                   aod865_dust_aqua(i)= -999.
     3594                   aod550_dustsco_aqua(i)= -999.
     3595                   aod670_dustsco_aqua(i)= -999.
     3596                   aod865_dustsco_aqua(i)= -999.
    33473597               ENDIF
    33483598               IF (masque_terra(i).GT. 0) THEN
    3349                    aod550_terra(i) = aod550_terra(i)/masque_terra(i)
     3599                   aod550_terra(i)=aod550_terra(i)/masque_terra(i)
    33503600                   aod670_terra(i)=aod670_terra(i)/masque_terra(i)
    33513601                   aod865_terra(i)=aod865_terra(i)/masque_terra(i)
    3352 
    3353                ELSE
     3602                   aod550_tr2_terra(i)=aod550_tr2_terra(i)/masque_terra(i)
     3603                   aod670_tr2_terra(i)=aod670_tr2_terra(i)/masque_terra(i)
     3604                   aod865_tr2_terra(i)=aod865_tr2_terra(i)/masque_terra(i)
     3605                   aod550_ss_terra(i)=aod550_ss_terra(i)/masque_terra(i)
     3606                   aod670_ss_terra(i)=aod670_ss_terra(i)/masque_terra(i)
     3607                   aod865_ss_terra(i)=aod865_ss_terra(i)/masque_terra(i)
     3608                   aod550_dust_terra(i)=aod550_dust_terra(i)/masque_terra(i)
     3609                   aod670_dust_terra(i)=aod670_dust_terra(i)/masque_terra(i)
     3610                   aod865_dust_terra(i)=aod865_dust_terra(i)/masque_terra(i)
     3611                   aod550_dustsco_terra(i)=aod550_dustsco_terra(i)/masque_terra(i)
     3612                   aod670_dustsco_terra(i)=aod670_dustsco_terra(i)/masque_terra(i)
     3613                   aod865_dustsco_terra(i)=aod865_dustsco_terra(i)/masque_terra(i)
     3614               ELSE
    33543615                   aod550_terra(i) = -999.
    33553616                   aod670_terra(i) = -999.
    33563617                   aod865_terra(i) = -999.
     3618                   aod550_tr2_terra(i)= -999.
     3619                   aod670_tr2_terra(i)= -999.
     3620                   aod865_tr2_terra(i)= -999.
     3621                   aod550_ss_terra(i)= -999.
     3622                   aod670_ss_terra(i)= -999.
     3623                   aod865_ss_terra(i)= -999.
     3624                   aod550_dust_terra(i)= -999.
     3625                   aod670_dust_terra(i)= -999.
     3626                   aod865_dust_terra(i)= -999.
     3627                   aod550_dustsco_terra(i)= -999.
     3628                   aod670_dustsco_terra(i)= -999.
     3629                   aod865_dustsco_terra(i)= -999.
    33573630               ENDIF
     3631!              IF (masque_terra(i).GT. 0) THEN
     3632!                   aod550_terra(i) = aod550_terra(i)/masque_terra(i)
     3633!                   aod670_terra(i)=aod670_terra(i)/masque_terra(i)
     3634!                   aod865_terra(i)=aod865_terra(i)/masque_terra(i)
     3635!
     3636!               ELSE
     3637!                   aod550_terra(i) = -999.
     3638!                   aod670_terra(i) = -999.
     3639!                   aod865_terra(i) = -999.
     3640!               ENDIF
    33583641          ENDDO         
    33593642!      !write  dbg
    … …  
    34003683!$OMP BARRIER
    34013684      ENDIF
    3402        !put in 0 everything
    3403        aod550_aqua(:) =0.
    3404        aod550_terra(:) =0.
    3405        aod670_aqua(:) =0.
    3406        aod670_terra(:) =0.
    3407        aod865_aqua(:) =0.
    3408        aod865_terra(:) =0.
    3409        masque_aqua(:) =0
    3410        masque_terra(:) =0
     3685!       !put in 0 everything
     3686!       aod550_aqua(:) =0.
     3687!       aod550_terra(:) =0.
     3688!       aod670_aqua(:) =0.
     3689!       aod670_terra(:) =0.
     3690!       aod865_aqua(:) =0.
     3691!       aod865_terra(:) =0.
     3692!       masque_aqua(:) =0
     3693!       masque_terra(:) =0
    34113694      ENDIF
    34123695
    … …  
    46924975
    46934976      tia_spla=tia_spla+REAL(ti_spla)/REAL(clock_rate)
    4694   print *,'---times for this timestep: time proc, time proc/time pytracr_spl-'
     4977  print *,'times for this timestep:timeproc,timeproc/time_pytracr_spl-'
    46954978     print *,'time spla',REAL(ti_spla)/REAL(clock_rate)                &
    46964979      ,REAL(ti_spla)/REAL(ti_spla)                                     
    … …  
    47735056      OPEN (UNIT=1,FILE=trim(adjustl(fileregions)))
    47745057      READ(1,'(a)') auxstr
    4775       READ(1,'(i)') nbreg
     5058      READ(1,'(i10)') nbreg
    47765059      CLOSE(UNIT=1)
    47775060      ENDIF
    … …  
    47975080      OPEN (UNIT=1,FILE=trim(adjustl(fileregionsdimsind)))
    47985081      READ(1,'(a)') auxstr
    4799       READ(1,'(i)') nbreg_ind
     5082      READ(1,'(i10)') nbreg_ind
    48005083      CLOSE(UNIT=1)
    48015084
    48025085      OPEN (UNIT=1,FILE=trim(adjustl(fileregionsdimsdust)))
    48035086      READ(1,'(a)') auxstr
    4804       READ(1,'(i)') nbreg_dust
     5087      READ(1,'(i10)') nbreg_dust
    48055088      CLOSE(UNIT=1)
    48065089
    48075090      OPEN (UNIT=1,FILE=trim(adjustl(fileregionsdimsbb)))
    48085091      READ(1,'(a)') auxstr
    4809       READ(1,'(i)') nbreg_bb
     5092      READ(1,'(i10)') nbreg_bb
    48105093      CLOSE(UNIT=1)
    48115094     
    … …  
    48345117      OPEN(1,file=trim(adjustl(filenameregion)))
    48355118      DO k=1,klon_glo
    4836       READ(1,'(i)') iregion_glo(k)
     5119      READ(1,'(i10)') iregion_glo(k)
    48375120      ENDDO
    48385121      CLOSE(UNIT=1)

  • LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/read_surface.F90

    r2175 r2303  
    2121!
    2222       real tmp_dyn(iip1,jjp1)
     23! JE no se       real tmp_dyn_glo(nbp_lon+1,nbp_lat)
    2324       real tmp_dyn_glo(nbp_lon,nbp_lat)
    2425       REAL tmp_dyn_invers(iip1,jjp1)
     26! JE no se       real tmp_dyn_invers_glo(nbp_lon+1,nbp_lat)
    2527       real tmp_dyn_invers_glo(nbp_lon,nbp_lat)
    2628       real tmp_fi(klon)
    … …  
    6365!       print *,'stat,i',status,i,outcycle,aux4s
    6466!       print *,'ifclause',status.NE. NF_NOERR ,outcycle == .false.
    65        IF ((.not.(status.NE. NF_NOERR) ).and.( outcycle == .false.)) THEN
     67       IF ((.not.(status.NE. NF_NOERR) ).and.( .not. outcycle )) THEN
    6668         outcycle=.true.
    6769         latstr=aux4s

  • LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/spla_output_dat.h

    r2176 r2303  
    4242  'taue865_dustsco','Tau ext 865dustsco','', (/ ('', i=1, 9) /))
    4343
     44
     45  type(ctrl_out),save :: o_taue550_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     46  'taue550_terra','Tau ext 550 terra','', &
     47      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     48         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     49  type(ctrl_out),save :: o_taue550_fine_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     50  'taue550_fine_terra','Tau ext fine 550 terra','', &
     51      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     52         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     53  type(ctrl_out),save :: o_taue550_coss_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     54  'taue550_coss_terra','Tau ext coss 550 terra','', &
     55      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     56         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     57  type(ctrl_out),save :: o_taue550_codu_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     58  'taue550_codu_terra','Tau ext codu 550 terra','', &
     59      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     60         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     61  type(ctrl_out),save :: o_taue550_scdu_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     62  'taue550_scdu_terra','Tau ext scdu 550 terra','', &
     63      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     64         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     65
     66  type(ctrl_out),save :: o_taue670_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     67  'taue670_terra','Tau ext 670 terra','', &
     68      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     69         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     70  type(ctrl_out),save :: o_taue670_fine_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     71  'taue670_fine_terra','Tau ext fine 670 terra','', &
     72      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     73         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     74  type(ctrl_out),save :: o_taue670_coss_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     75  'taue670_coss_terra','Tau ext coss 670 terra','', &
     76      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     77         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     78  type(ctrl_out),save :: o_taue670_codu_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     79  'taue670_codu_terra','Tau ext codu 670 terra','', &
     80      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     81         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     82  type(ctrl_out),save :: o_taue670_scdu_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     83  'taue670_scdu_terra','Tau ext scdu 670 terra','', &
     84      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     85         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     86
     87  type(ctrl_out),save :: o_taue865_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     88  'taue865_terra','Tau ext 865 terra','', &
     89      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     90         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     91  type(ctrl_out),save :: o_taue865_fine_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     92  'taue865_fine_terra','Tau ext fine 865 terra','', &
     93      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     94         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     95  type(ctrl_out),save :: o_taue865_coss_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     96  'taue865_coss_terra','Tau ext coss 865 terra','', &
     97      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     98         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     99  type(ctrl_out),save :: o_taue865_codu_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     100  'taue865_codu_terra','Tau ext codu 865 terra','', &
     101      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     102         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     103  type(ctrl_out),save :: o_taue865_scdu_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     104  'taue865_scdu_terra','Tau ext scdu 865 terra','', &
     105      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     106         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     107
    44108  type(ctrl_out),save :: o_taue550_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    45109  'taue550_aqua','Tau ext 550 aqua','', &
    46110      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    47111         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    48   type(ctrl_out),save :: o_taue550_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    49   'taue550_terra','Tau ext 550 terra','', &
     112  type(ctrl_out),save :: o_taue550_fine_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     113  'taue550_fine_aqua','Tau ext fine 550 aqua','', &
     114      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     115         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     116  type(ctrl_out),save :: o_taue550_coss_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     117  'taue550_coss_aqua','Tau ext coss 550 aqua','', &
     118      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     119         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     120  type(ctrl_out),save :: o_taue550_codu_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     121  'taue550_codu_aqua','Tau ext codu 550 aqua','', &
     122      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     123         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     124  type(ctrl_out),save :: o_taue550_scdu_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     125  'taue550_scdu_aqua','Tau ext scdu 550 aqua','', &
    50126      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    51127         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    … …  
    55131      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    56132         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    57 
    58   type(ctrl_out),save :: o_taue670_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    59   'taue670_terra','Tau ext 670 terra','', &
     133  type(ctrl_out),save :: o_taue670_fine_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     134  'taue670_fine_aqua','Tau ext fine 670 aqua','', &
     135      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     136         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     137  type(ctrl_out),save :: o_taue670_coss_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     138  'taue670_coss_aqua','Tau ext coss 670 aqua','', &
     139      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     140         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     141  type(ctrl_out),save :: o_taue670_codu_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     142  'taue670_codu_aqua','Tau ext codu 670 aqua','', &
     143      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     144         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     145  type(ctrl_out),save :: o_taue670_scdu_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     146  'taue670_scdu_aqua','Tau ext scdu 670 aqua','', &
    60147      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    61148         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    … …  
    65152      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    66153         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    67 
    68   type(ctrl_out),save :: o_taue865_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
    69   'taue865_terra','Tau ext 865 terra','', &
    70       (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    71          't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     154  type(ctrl_out),save :: o_taue865_fine_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     155  'taue865_fine_aqua','Tau ext fine 865 aqua','', &
     156      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     157         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     158  type(ctrl_out),save :: o_taue865_coss_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     159  'taue865_coss_aqua','Tau ext coss 865 aqua','', &
     160      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     161         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     162  type(ctrl_out),save :: o_taue865_codu_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     163  'taue865_codu_aqua','Tau ext codu 865 aqua','', &
     164      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     165         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     166  type(ctrl_out),save :: o_taue865_scdu_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     167  'taue865_scdu_aqua','Tau ext scdu 865 aqua','', &
     168      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     169         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     170
     171
     172
     173!  type(ctrl_out),save :: o_taue550_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     174!  'taue550_terra','Tau ext 550 terra','', &
     175!      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     176!         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     177!
     178!  type(ctrl_out),save :: o_taue670_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     179!  'taue670_aqua','Tau ext 670 aqua','', &
     180!      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     181!         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     182!
     183!  type(ctrl_out),save :: o_taue670_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     184!  'taue670_terra','Tau ext 670 terra','', &
     185!      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     186!         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     187!
     188!  type(ctrl_out),save :: o_taue865_aqua     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     189!  'taue865_aqua','Tau ext 865 aqua','', &
     190!      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     191!         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
     192!
     193!  type(ctrl_out),save :: o_taue865_terra     = ctrl_out((/ 4, 4, 4, 10, 10, 10, 10, 10, 10 /), &
     194!  'taue865_terra','Tau ext 865 terra','', &
     195!      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
     196!         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    72197
    73198       

  • LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/spla_output_write.h

    r2196 r2303  
    2727      CALL histwrite_phy( o_taue865_aqua     , aod865_aqua    )
    2828      CALL histwrite_phy( o_taue865_terra    , aod865_terra    )
     29
     30      IF(id_fine>0) CALL histwrite_phy( o_taue550_fine_aqua     ,aod550_tr2_aqua   )
     31      IF(id_fine>0) CALL histwrite_phy( o_taue670_fine_aqua     ,aod670_tr2_aqua   )
     32      IF(id_fine>0) CALL histwrite_phy( o_taue865_fine_aqua     ,aod865_tr2_aqua   )
     33      IF(id_coss>0) CALL histwrite_phy( o_taue550_coss_aqua      ,aod550_ss_aqua    )
     34      IF(id_coss>0) CALL histwrite_phy( o_taue670_coss_aqua      ,aod670_ss_aqua    )
     35      IF(id_coss>0) CALL histwrite_phy( o_taue865_coss_aqua      ,aod865_ss_aqua    )
     36      IF(id_codu>0) CALL histwrite_phy( o_taue550_codu_aqua    ,aod550_dust_aqua    )
     37      IF(id_codu>0) CALL histwrite_phy( o_taue670_codu_aqua    ,aod670_dust_aqua    )
     38      IF(id_codu>0) CALL histwrite_phy( o_taue865_codu_aqua    ,aod865_dust_aqua    )
     39      IF(id_scdu>0) CALL histwrite_phy( o_taue670_scdu_aqua ,aod670_dustsco_aqua    )
     40      IF(id_scdu>0) CALL histwrite_phy( o_taue550_scdu_aqua ,aod550_dustsco_aqua    )
     41      IF(id_scdu>0) CALL histwrite_phy( o_taue865_scdu_aqua ,aod865_dustsco_aqua    )
     42     
     43      IF(id_fine>0) CALL histwrite_phy( o_taue550_fine_terra     ,aod550_tr2_terra   )
     44      IF(id_fine>0) CALL histwrite_phy( o_taue670_fine_terra     ,aod670_tr2_terra   )
     45      IF(id_fine>0) CALL histwrite_phy( o_taue865_fine_terra     ,aod865_tr2_terra   )
     46      IF(id_coss>0) CALL histwrite_phy( o_taue550_coss_terra      ,aod550_ss_terra    )
     47      IF(id_coss>0) CALL histwrite_phy( o_taue670_coss_terra      ,aod670_ss_terra    )
     48      IF(id_coss>0) CALL histwrite_phy( o_taue865_coss_terra      ,aod865_ss_terra    )
     49      IF(id_codu>0) CALL histwrite_phy( o_taue550_codu_terra    ,aod550_dust_terra    )
     50      IF(id_codu>0) CALL histwrite_phy( o_taue670_codu_terra    ,aod670_dust_terra    )
     51      IF(id_codu>0) CALL histwrite_phy( o_taue865_codu_terra    ,aod865_dust_terra    )
     52      IF(id_scdu>0) CALL histwrite_phy( o_taue670_scdu_terra ,aod670_dustsco_terra    )
     53      IF(id_scdu>0) CALL histwrite_phy( o_taue550_scdu_terra ,aod550_dustsco_terra    )
     54      IF(id_scdu>0) CALL histwrite_phy( o_taue865_scdu_terra ,aod865_dustsco_terra    )
     55
     56
     57
    2958      IF(id_prec>0) CALL histwrite_phy( o_trm01 , trm01 )
    3059      IF(id_fine>0) CALL histwrite_phy( o_trm02 , trm02 )
Note: See TracChangeset for help on using the changeset viewer.