Ignore:
Timestamp:
Feb 25, 2015, 2:43:34 PM (10 years ago)
Author:
jescribano
Message:

Bugs corrections. Included a correction/tunning factor for the Chimere-dust emissions, Constant of MB95 equal to 2.61 as in MB95. No spurious increase of u* before horizontal flux calculations in the dust emission scheme. Values of AG00 binding energies fixed as the original AG00 divided by 3 as is Sow et al 2011 ACPD.

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

Legend:

Unmodified
Added
Removed
  • LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/aeropt_spl.F

    r2196 r2217  
    5050      PARAMETER (rh_int=12)
    5151      REAL auxreal
    52 c      REAL ss_a(nb_lambda,rh_int,nbtr-1)
     52c      REAL ss_a(nb_lambda,int,nbtr-1)
    5353c      DATA ss_a/72*1./
    5454      REAL ss_dust(nb_lambda), ss_acc550(rh_int), alpha_acc
     
    9090c
    9191      IF (ok_chimeredust) THEN
    92       ss_dust=(/0.5167768,0.5684330,0.6531643/)
    93       ss_dustsco=(/0.1003391,0.1012288,0.1024651/)
     92!JE20150212<< : changes in ustar in dustmod changes emission distribution
     93!      ss_dust=(/0.5167768,0.5684330,0.6531643/)
     94!      ss_dustsco=(/0.1003391,0.1012288,0.1024651/)
     95      ss_dust=(/0.4670522, 0.5077308, 0.5745184/)
     96      ss_dustsco=(/0.099858, 0.1007395, 0.1019673/)
     97!JE20150212>>
     98
    9499      ELSE
    95100      ss_dust=(/0.564, 0.614, 0.700/)
  • LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/dustemission_mod.F90

    r2196 r2217  
    66!  INTEGER, PARAMETER     :: nbins=800  ! number of aerosol bins: for spla
    77!  INTEGER, PARAMETER     :: nbins=8000  ! number of aerosol bins: for spla
     8 
     9  INTEGER, PARAMETER     :: flag_feff=1 ! 0: deactivate feff (drag partition scheme)
    810  INTEGER, PARAMETER     :: nbins=800  ! number of aerosol bins: for spla
    911  INTEGER, PARAMETER     :: nmode=3   ! number of soil-dust modes
     
    3032!  real   , parameter :: e2=3.52/div2
    3133!  real   , parameter :: e3=3.46/div3
    32 !  real   , parameter :: factorflux=10.
    3334!  real   , parameter :: rop=2.65 ! particle density g/m3
    3435!  real   , parameter :: roa=0.001227  ! air density g/m3
     
    4647  real   , parameter :: e2=3.52/div2
    4748  real   , parameter :: e3=3.46/div3
    48   real   , parameter :: factorflux=1.
    4949  real   , parameter :: rop=2.65 ! particle density g/m3
    5050  real   , parameter :: roa=0.001227  ! air density g/m3
     
    5353! C=2.61 from Marticorena and Bergametti 1995 instead of Gillete and Chen 2001
    5454! (recommended C=1.1  in supply-limited dust source area.. )
    55   real   , parameter :: cd=2.61*roa/gravity
     55   real   , parameter :: cd=2.61*roa/gravity
     56!  real   , parameter :: cd=1.0*roa/gravity
    5657!JE20150202>>>>
    5758  real,parameter     :: beta=16300.
     
    286287  real,parameter :: sizescohigh=30.  ! in micrometers
    287288
    288   real,parameter :: tunningfine=1.0  ! factor for fine bins!!! important!!
     289  real,parameter :: tunningfactor=4.5  ! factor for fine bins!!! important!!
    289290!JE20140915  real,parameter :: sizeacclow=0.06
    290291!JE20140915  real,parameter :: sizeacchigh=1.0
     
    316317   call abort_gcm('adaptdustemission', 'Dust range problem',1)
    317318  endif
    318   print *,'FINE DUST BIN: tunning EMISSION factor= ',tunningfine
     319  print *,'ALL DUST BIN: tunning EMISSION factor= ',tunningfactor
    319320  auxr1=9999.
    320321  auxr2=9999.
     
    394395   auxr1=auxr1+emisbinlocal(k,i)
    395396  enddo
    396   emdustacc(k)= auxr1*tunningfine
     397  emdustacc(k)= auxr1*tunningfactor
    397398  do i=imincoalow,imincoahigh-1
    398399    auxr2=auxr2+emisbinlocal(k,i)
    399400  enddo
    400   emdustcoa(k)=auxr2
     401  emdustcoa(k)=auxr2*tunningfactor
    401402  do i=iminscolow,iminscohigh-1
    402403    auxr3=auxr3+emisbinlocal(k,i)
    403404  enddo
    404   emdustsco(k)=auxr3
     405  emdustsco(k)=auxr3*tunningfactor
    405406enddo
    406407
     
    822823          enddo
    823824        enddo
     825! JE20150120<<
     826  if (flag_feff .eq. 0) then
     827    print *,'JE_dbg FORCED deactivated feff'
     828    do i=1,klon
     829      do k=1,ntyp
     830        feff(i,k)=1.
     831      enddo
     832    enddo
     833  endif
     834! JE20150120>>
     835
     836
    824837if (1==1) then
    825838!  !  CALL writefield_phy("AA",tmp1(1:klon,1:5),5)
     
    10621075  REAL :: pdfcum,U10mMOD,pdfu,weilambda
    10631076  REAL :: z0salt,ceff,cerod,cpcent
    1064 !JE20150202  !REAL :: cdnms,ustarns,modwm,utmin
    1065   REAL :: cdnms,ustarns,modwm
     1077  REAL :: cdnms,ustarns,modwm,utmin
     1078!JE20150202   REAL :: cdnms,ustarns,modwm
    10661079  REAL :: fdp1,fdp2,ad1,ad2,ad3,flux_diam
    10671080  REAL :: dfec1,dfec2,dfec3,t1,t2,t3,p1,p2,p3,dec,ec
     
    12591272             ENDDO !n=1,ntyp
    1260127370 CONTINUE
    1261 !factorflux
    12621274        fluxdust(i,1)=fluxdust(i,1)+flux1*probu(kwb)
    12631275        fluxdust(i,2)=fluxdust(i,2)+flux2*probu(kwb)
    12641276        fluxdust(i,3)=fluxdust(i,3)+flux3*probu(kwb)
    12651277   ENDDO !kwb=1,nwb
    1266 !JE20150202 <<
    1267 !      m1dflux(i)=10.*fluxdust(i,1)
    1268 !      m2dflux(i)=10.*fluxdust(i,2)          ! tous en Kg/m2/s
    1269 !      m3dflux(i)=10.*fluxdust(i,3)
    1270       m1dflux(i)=factorflux*10.*fluxdust(i,1)
    1271       m2dflux(i)=factorflux*10.*fluxdust(i,2)          ! tous en Kg/m2/s
    1272       m3dflux(i)=factorflux*10.*fluxdust(i,3)
    1273 !JE20150202 >>
     1278      m1dflux(i)=10.*fluxdust(i,1)
     1279      m2dflux(i)=10.*fluxdust(i,2)          ! tous en Kg/m2/s
     1280      m3dflux(i)=10.*fluxdust(i,3)
    12741281
    12751282
     
    12851292    DO j=1,nbins
    12861293!JE20150202 <<
    1287 !    emisbin(k,j) = emisbin(k,j)+10*fluxdust(k,i)*massfrac(i,j) 
    1288     emisbin(k,j) = emisbin(k,j)+factorflux*fluxdust(k,i)*massfrac(i,j) 
     1294    emisbin(k,j) = emisbin(k,j)+10*fluxdust(k,i)*massfrac(i,j) 
    12891295!JE20150202 >>
    12901296    ENDDO !j, nbind
  • LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/physiq.F90

    r2196 r2217  
    36423642  nbreg_wstardust=1
    36433643  fileregionsdimsind='regions_ind_meta'
    3644   fileregionsdimsdust='regions_dustacc_meta'
     3644!  fileregionsdimsdust='regions_dustacc_meta'
     3645  fileregionsdimsdust='regions_dust_meta'
    36453646  fileregionsdimsbb='regions_bb_meta'
    36463647  fileregionsdimswstar='regions_pwstarwake_meta'
     
    37043705  paramname_dustcoa='dustcoasco'
    37053706  paramname_dustsco='dustcoasco'
     3707  paramname_dustacc='dust'
     3708  paramname_dustcoa='dust'
     3709  paramname_dustsco='dust'
    37063710  paramname_wstarBL='pwstarbl'
    37073711  paramname_wstarWAKE='pwstarwake'
  • LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/phytracr_spl_mod.F90

    r2196 r2217  
    20542054!      ENDIF
    20552055!
    2056       c_FullName1='regions_dustacc'
     2056      !c_FullName1='regions_dustacc'
     2057      c_FullName1='regions_dust'
    20572058      call readregions_spl(iregion_dust,c_FullName1)
    20582059      c_FullName1='regions_ind'
  • LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/sediment_mod.F

    r2196 r2217  
    7676!JE20140911        mmd_dustsco=15.56554  !micrometer for bin 0.5-20 and 0.5-10 um
    7777        IF (ok_chimeredust) THEN
    78         mmd_dust=3.761212  !micrometer for bin 0.5-3 and 0.5-10 um
    79         mmd_dustsco=15.06167  !micrometer for bin 3-20 and 0.5-10 um
     78!JE20150212<< : changes in ustar in dustmod changes emission distribution
     79!        mmd_dust=3.761212  !micrometer for bin 0.5-3 and 0.5-10 um
     80!        mmd_dustsco=15.06167  !micrometer for bin 3-20 and 0.5-10 um
     81!JE20150212>>
     82        mmd_dust=3.983763
     83        mmd_dustsco=15.10854
    8084        ELSE
    8185        mmd_dust=11.0  !micrometer for bin 0.5-20 and 0.5-10 um
Note: See TracChangeset for help on using the changeset viewer.