Changeset 2196 for LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/sediment_mod.F

Timestamp:

Feb 6, 2015, 6:53:57 PM (10 years ago)

Author:

jescribano

Message:

Dust emission scheme changes: (1) Included possibility of use previous dust emission scheme (with 2 dust bins). (2) Parameter of Marticorena and Bergametti 1995 set to it's original value (2.61) for testing purposes with pdtphys=15min. (3) Cleaning ustar calculation.

File:

• LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/sediment_mod.F (modified) (8 diffs)

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

@@ -3 +3 @@
       SUBROUTINE sediment_mod(t_seri,pplay,zrho,paprs,time_step,RHcl,
      .                                       id_coss,id_codu,id_scdu,
+     .                                        ok_chimeredust,
      .                           sed_ss,sed_dust,sed_dustsco,tr_seri)
 cnhl     .                                       xlon,xlat,
@@ -29 +30 @@
        REAL paprs(klon, klev+1)    !pressure at interface of layers Z (Pa)
        REAL time_step            !time step (sec)
-
+       LOGICAL ok_chimeredust
        REAL xlat(klon)       ! latitudes pour chaque point 
        REAL xlon(klon)       ! longitudes pour chaque point 
@@ -74 +75 @@
 !JE20140911        mmd_dust=5.156346  !micrometer for bin 0.5-20 and 0.5-10 um
 !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
-
+        ELSE 
+        mmd_dust=11.0  !micrometer for bin 0.5-20 and 0.5-10 um
+        mmd_dustsco=100. ! absurd value, bin not used in this scheme
+        ENDIF
 
 
@@ -103 +108 @@
 c
 c
+       IF(id_coss>0) THEN
        DO k=1, klev
        DO i=1,klon
@@ -159 +165 @@
          sed_ss(i)=sed_flux(i,1)*1.e6*1.e3    !--unit mg/m2/s
        ENDDO          !klon
+       ELSE
+        DO i=1, klon
+          sed_ss(i)=0.
+        ENDDO
+       ENDIF
 c
 c
@@ -165 +176 @@
 c
 c
+       IF(id_codu>0) THEN
        DO k=1, klev
        DO i=1,klon
@@ -206 +218 @@
          sed_dust(i)=sed_flux(i,1)*1.e6*1.e3    !--unit mg/m2/s
        ENDDO          !klon
+       ELSE
+        DO i=1, klon
+          sed_dust(i)=0.
+        ENDDO
+       ENDIF
+c
+
 
 c--------- For scoarse  dust ------------------
 c
 c
+       IF(id_scdu>0) THEN
        DO k=1, klev
        DO i=1,klon
@@ -252 +272 @@
          sed_dustsco(i)=sed_flux(i,1)*1.e6*1.e3    !--unit mg/m2/s
        ENDDO          !klon
+       ELSE
+        DO i=1, klon
+          sed_dustsco(i)=0.
+        ENDDO
+       ENDIF
+c