Changeset 2217

Timestamp:

Feb 25, 2015, 2:43:34 PM (9 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:

• aeropt_spl.F (modified) (2 diffs)
• aeropt_spl_F_20140902 (deleted)
• dustemission_mod.F90 (modified) (11 diffs)
• dustemission_mod_F90_bck (deleted)
• gastoparticle_dot_F_orig20140120 (deleted)
• physiq.F90 (modified) (2 diffs)
• phytracr_sp_dot_F_20140121_seq (deleted)
• phytracr_spl.F_20140506_seq (deleted)
• phytracr_spl.F_20140507_seq (deleted)
• phytracr_spl_dot_F_20140210_seq (deleted)
• phytracr_spl_dot_F_20140210b_seq (deleted)
• phytracr_spl_mod.F90 (modified) (1 diff)
• read_newemissions_dot_F_orig_je20140113_seq (deleted)
• sediment_mod.F (modified) (1 diff)

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

@@ -50 +50 @@
       PARAMETER (rh_int=12)
       REAL auxreal
-c      REAL ss_a(nb_lambda,rh_int,nbtr-1)
+c      REAL ss_a(nb_lambda,int,nbtr-1)
 c      DATA ss_a/72*1./
       REAL ss_dust(nb_lambda), ss_acc550(rh_int), alpha_acc
@@ -90 +90 @@
 c
       IF (ok_chimeredust) THEN 
-      ss_dust=(/0.5167768,0.5684330,0.6531643/)
-      ss_dustsco=(/0.1003391,0.1012288,0.1024651/)
+!JE20150212<< : changes in ustar in dustmod changes emission distribution
+!      ss_dust=(/0.5167768,0.5684330,0.6531643/)
+!      ss_dustsco=(/0.1003391,0.1012288,0.1024651/)
+      ss_dust=(/0.4670522, 0.5077308, 0.5745184/)
+      ss_dustsco=(/0.099858, 0.1007395, 0.1019673/)
+!JE20150212>>
+
       ELSE
       ss_dust=(/0.564, 0.614, 0.700/)

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

@@ -6 +6 @@
 !  INTEGER, PARAMETER     :: nbins=800  ! number of aerosol bins: for spla
 !  INTEGER, PARAMETER     :: nbins=8000  ! number of aerosol bins: for spla
+  
+  INTEGER, PARAMETER     :: flag_feff=1 ! 0: deactivate feff (drag partition scheme)
   INTEGER, PARAMETER     :: nbins=800  ! number of aerosol bins: for spla
   INTEGER, PARAMETER     :: nmode=3   ! number of soil-dust modes 
@@ -30 +32 @@
 !  real   , parameter :: e2=3.52/div2
 !  real   , parameter :: e3=3.46/div3
-!  real   , parameter :: factorflux=10.
 !  real   , parameter :: rop=2.65 ! particle density g/m3
 !  real   , parameter :: roa=0.001227  ! air density g/m3
@@ -46 +47 @@
   real   , parameter :: e2=3.52/div2
   real   , parameter :: e3=3.46/div3
-  real   , parameter :: factorflux=1.
   real   , parameter :: rop=2.65 ! particle density g/m3
   real   , parameter :: roa=0.001227  ! air density g/m3
@@ -53 +53 @@
 ! C=2.61 from Marticorena and Bergametti 1995 instead of Gillete and Chen 2001
 ! (recommended C=1.1  in supply-limited dust source area.. )
-  real   , parameter :: cd=2.61*roa/gravity
+   real   , parameter :: cd=2.61*roa/gravity
+!  real   , parameter :: cd=1.0*roa/gravity
 !JE20150202>>>>
   real,parameter     :: beta=16300.
@@ -286 +287 @@
   real,parameter :: sizescohigh=30.  ! in micrometers
 
-  real,parameter :: tunningfine=1.0  ! factor for fine bins!!! important!!
+  real,parameter :: tunningfactor=4.5  ! factor for fine bins!!! important!!
 !JE20140915  real,parameter :: sizeacclow=0.06
 !JE20140915  real,parameter :: sizeacchigh=1.0
@@ -316 +317 @@
    call abort_gcm('adaptdustemission', 'Dust range problem',1)
   endif
-  print *,'FINE DUST BIN: tunning EMISSION factor= ',tunningfine
+  print *,'ALL DUST BIN: tunning EMISSION factor= ',tunningfactor
   auxr1=9999.
   auxr2=9999.
@@ -394 +395 @@
    auxr1=auxr1+emisbinlocal(k,i)
   enddo
-  emdustacc(k)= auxr1*tunningfine
+  emdustacc(k)= auxr1*tunningfactor
   do i=imincoalow,imincoahigh-1
     auxr2=auxr2+emisbinlocal(k,i)
   enddo
-  emdustcoa(k)=auxr2
+  emdustcoa(k)=auxr2*tunningfactor
   do i=iminscolow,iminscohigh-1
     auxr3=auxr3+emisbinlocal(k,i)
   enddo
-  emdustsco(k)=auxr3
+  emdustsco(k)=auxr3*tunningfactor
 enddo
 
@@ -822 +823 @@
           enddo
         enddo
+! JE20150120<<
+  if (flag_feff .eq. 0) then
+    print *,'JE_dbg FORCED deactivated feff'
+    do i=1,klon
+      do k=1,ntyp
+        feff(i,k)=1.
+      enddo
+    enddo
+  endif
+! JE20150120>>
+
+
 if (1==1) then
 !  !  CALL writefield_phy("AA",tmp1(1:klon,1:5),5)
@@ -1062 +1075 @@
   REAL :: pdfcum,U10mMOD,pdfu,weilambda
   REAL :: z0salt,ceff,cerod,cpcent
-!JE20150202  !REAL :: cdnms,ustarns,modwm,utmin
-  REAL :: cdnms,ustarns,modwm
+  REAL :: cdnms,ustarns,modwm,utmin
+!JE20150202   REAL :: cdnms,ustarns,modwm
   REAL :: fdp1,fdp2,ad1,ad2,ad3,flux_diam
   REAL :: dfec1,dfec2,dfec3,t1,t2,t3,p1,p2,p3,dec,ec
@@ -1259 +1272 @@
              ENDDO !n=1,ntyp
 70 CONTINUE
-!factorflux
         fluxdust(i,1)=fluxdust(i,1)+flux1*probu(kwb)
         fluxdust(i,2)=fluxdust(i,2)+flux2*probu(kwb)
         fluxdust(i,3)=fluxdust(i,3)+flux3*probu(kwb)
    ENDDO !kwb=1,nwb
-!JE20150202 <<
-!      m1dflux(i)=10.*fluxdust(i,1)
-!      m2dflux(i)=10.*fluxdust(i,2)          ! tous en Kg/m2/s
-!      m3dflux(i)=10.*fluxdust(i,3)
-      m1dflux(i)=factorflux*10.*fluxdust(i,1)
-      m2dflux(i)=factorflux*10.*fluxdust(i,2)          ! tous en Kg/m2/s
-      m3dflux(i)=factorflux*10.*fluxdust(i,3)
-!JE20150202 >>
+      m1dflux(i)=10.*fluxdust(i,1)
+      m2dflux(i)=10.*fluxdust(i,2)          ! tous en Kg/m2/s
+      m3dflux(i)=10.*fluxdust(i,3)
 
 
@@ -1285 +1292 @@
     DO j=1,nbins
 !JE20150202 <<
-!    emisbin(k,j) = emisbin(k,j)+10*fluxdust(k,i)*massfrac(i,j)  
-    emisbin(k,j) = emisbin(k,j)+factorflux*fluxdust(k,i)*massfrac(i,j)  
+    emisbin(k,j) = emisbin(k,j)+10*fluxdust(k,i)*massfrac(i,j)  
 !JE20150202 >>
     ENDDO !j, nbind

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

@@ -3642 +3642 @@
   nbreg_wstardust=1
   fileregionsdimsind='regions_ind_meta'
-  fileregionsdimsdust='regions_dustacc_meta'
+!  fileregionsdimsdust='regions_dustacc_meta'
+  fileregionsdimsdust='regions_dust_meta'
   fileregionsdimsbb='regions_bb_meta'
   fileregionsdimswstar='regions_pwstarwake_meta'
@@ -3704 +3705 @@
   paramname_dustcoa='dustcoasco'
   paramname_dustsco='dustcoasco'
+  paramname_dustacc='dust'
+  paramname_dustcoa='dust'
+  paramname_dustsco='dust'
   paramname_wstarBL='pwstarbl'
   paramname_wstarWAKE='pwstarwake'

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

@@ -2054 +2054 @@
 !      ENDIF
 !
-      c_FullName1='regions_dustacc'
+      !c_FullName1='regions_dustacc'
+      c_FullName1='regions_dust'
       call readregions_spl(iregion_dust,c_FullName1)
       c_FullName1='regions_ind'

LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/sediment_mod.F

@@ -76 +76 @@
 !JE20140911        mmd_dustsco=15.56554  !micrometer for bin 0.5-20 and 0.5-10 um
         IF (ok_chimeredust) THEN 
-        mmd_dust=3.761212  !micrometer for bin 0.5-3 and 0.5-10 um
-        mmd_dustsco=15.06167  !micrometer for bin 3-20 and 0.5-10 um
+!JE20150212<< : changes in ustar in dustmod changes emission distribution
+!        mmd_dust=3.761212  !micrometer for bin 0.5-3 and 0.5-10 um
+!        mmd_dustsco=15.06167  !micrometer for bin 3-20 and 0.5-10 um
+!JE20150212>>
+        mmd_dust=3.983763
+        mmd_dustsco=15.10854 
         ELSE 
         mmd_dust=11.0  !micrometer for bin 0.5-20 and 0.5-10 um