Changeset 2196 for LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/finemission.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/finemission.F (modified) (4 diffs)

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

@@ -68 +68 @@
       DO i=1,klon
          IF (iregion_ind(i).GT.0) THEN
-           source_tr(i,id_fine)=source_tr(i,id_fine)+
+       IF(id_fine>0)    source_tr(i,id_fine)=source_tr(i,id_fine)+
      .                (scale_param_ff(iregion_ind(i))*lmt_bcff(i)+     !g/m2/s
      .                 scale_param_ff(iregion_ind(i))*lmt_omff(i)
      .                 )     * 1.e4                        !g/m2/s
 c
-           flux_tr(i,id_fine)=flux_tr(i,id_fine)+
+      IF(id_fine>0)     flux_tr(i,id_fine)=flux_tr(i,id_fine)+
      .                (scale_param_ff(iregion_ind(i))*lmt_bcff(i)+     !mg/m2/s
      .                 scale_param_ff(iregion_ind(i))*lmt_omff(i)
@@ -84 +84 @@
          ENDIF
          IF (iregion_bb(i).GT.0) THEN
-           source_tr(i,id_fine)=source_tr(i,id_fine)+
+       IF(id_fine>0)    source_tr(i,id_fine)=source_tr(i,id_fine)+
      .                (scale_param_bb(iregion_bb(i))*lmt_bcbb_l(i)+   !g/m2/s
      .                 scale_param_bb(iregion_bb(i))*lmt_ombb_l(i)    !g/m2/s
      .                 )     * 1.e4                        !g/m2/s
 c
-           flux_tr(i,id_fine)=flux_tr(i,id_fine)+
+       IF(id_fine>0)    flux_tr(i,id_fine)=flux_tr(i,id_fine)+
      .                (scale_param_bb(iregion_bb(i))*lmt_bcbb_l(i)+   !mg/m2/s
      .                 scale_param_bb(iregion_bb(i))*lmt_ombb_l(i)+   !mg/m2/s
@@ -101 +101 @@
      .                   *1.e4*1.e3
          ENDIF
-         source_tr(i,id_fine)=source_tr(i,id_fine)+
+       IF(id_fine>0)  source_tr(i,id_fine)=source_tr(i,id_fine)+
      .                (lmt_bcnff(i)+lmt_bcba(i)+lmt_omnff(i)+
      .                 lmt_omnat(i)+lmt_omba(i))     * 1.e4           !g/m2/s
 c
-         flux_tr(i,id_fine)=flux_tr(i,id_fine)+
+       IF(id_fine>0)  flux_tr(i,id_fine)=flux_tr(i,id_fine)+
      .                (lmt_bcnff(i)+lmt_omnff(i)+lmt_omnat(i)+
      .                 lmt_omba(i)+lmt_bcba(i))     * 1.e4 *1.e3      !mg/m2/s
@@ -135 +135 @@
 c
          IF (iregion_bb(i) .GT.0) THEN
-           tr_seri(i,k,id_fine)=tr_seri(i,k,id_fine)+
+        IF(id_fine>0)   tr_seri(i,k,id_fine)=tr_seri(i,k,id_fine)+
      .              (scale_param_bb(iregion_bb(i))*lmt_bcbb_h(i)+
      .               scale_param_bb(iregion_bb(i))*lmt_ombb_h(i))