Changeset 2196 for LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/sediment_mod.F
- Timestamp:
- Feb 6, 2015, 6:53:57 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/sediment_mod.F
r2175 r2196 3 3 SUBROUTINE sediment_mod(t_seri,pplay,zrho,paprs,time_step,RHcl, 4 4 . id_coss,id_codu,id_scdu, 5 . ok_chimeredust, 5 6 . sed_ss,sed_dust,sed_dustsco,tr_seri) 6 7 cnhl . xlon,xlat, … … 29 30 REAL paprs(klon, klev+1) !pressure at interface of layers Z (Pa) 30 31 REAL time_step !time step (sec) 31 32 LOGICAL ok_chimeredust 32 33 REAL xlat(klon) ! latitudes pour chaque point 33 34 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 75 76 !JE20140911 mmd_dustsco=15.56554 !micrometer for bin 0.5-20 and 0.5-10 um 77 IF (ok_chimeredust) THEN 76 78 mmd_dust=3.761212 !micrometer for bin 0.5-3 and 0.5-10 um 77 79 mmd_dustsco=15.06167 !micrometer for bin 3-20 and 0.5-10 um 78 80 ELSE 81 mmd_dust=11.0 !micrometer for bin 0.5-20 and 0.5-10 um 82 mmd_dustsco=100. ! absurd value, bin not used in this scheme 83 ENDIF 79 84 80 85 … … 103 108 c 104 109 c 110 IF(id_coss>0) THEN 105 111 DO k=1, klev 106 112 DO i=1,klon … … 159 165 sed_ss(i)=sed_flux(i,1)*1.e6*1.e3 !--unit mg/m2/s 160 166 ENDDO !klon 167 ELSE 168 DO i=1, klon 169 sed_ss(i)=0. 170 ENDDO 171 ENDIF 161 172 c 162 173 c … … 165 176 c 166 177 c 178 IF(id_codu>0) THEN 167 179 DO k=1, klev 168 180 DO i=1,klon … … 206 218 sed_dust(i)=sed_flux(i,1)*1.e6*1.e3 !--unit mg/m2/s 207 219 ENDDO !klon 220 ELSE 221 DO i=1, klon 222 sed_dust(i)=0. 223 ENDDO 224 ENDIF 225 c 226 208 227 209 228 c--------- For scoarse dust ------------------ 210 229 c 211 230 c 231 IF(id_scdu>0) THEN 212 232 DO k=1, klev 213 233 DO i=1,klon … … 252 272 sed_dustsco(i)=sed_flux(i,1)*1.e6*1.e3 !--unit mg/m2/s 253 273 ENDDO !klon 274 ELSE 275 DO i=1, klon 276 sed_dustsco(i)=0. 277 ENDDO 278 ENDIF 279 c 254 280 255 281
Note: See TracChangeset
for help on using the changeset viewer.