- Timestamp:
- Dec 26, 2013, 6:34:40 PM (11 years ago)
- Location:
- trunk/LMDZ.GENERIC/libf/phystd
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.GENERIC/libf/phystd/aeropacity.F90
r1026 r1132 381 381 ELSEIF (pplev(ig,l) .lt. pres_top_tropo .and. pplev(ig,l) .gt. pres_bottom_strato) THEN 382 382 expfactor=log(obs_tau_col_strato/obs_tau_col_tropo)/log(pres_bottom_strato/pres_top_tropo) 383 aerosol(ig,l,iaer)= obs_tau_col_tropo*((pplev(ig,l)/pres_top_tropo)**expfactor)*aerosol(ig,l,iaer) 383 aerosol(ig,l,iaer)= obs_tau_col_tropo*((pplev(ig,l)/pres_top_tropo)**expfactor)*aerosol(ig,l,iaer)/1.5 384 384 !! 4. strato layer 385 385 ELSEIF (pplev(ig,l) .le. pres_bottom_strato .and. pplev(ig,l) .gt. pres_top_strato) THEN -
trunk/LMDZ.GENERIC/libf/phystd/iniaerosol.F
r1026 r1132 55 55 pres_bottom_tropo= 66000.0 !! In Pa ! 56 56 pres_top_tropo= 18000.0 !! In Pa 57 pres_bottom_strato= 10000.0 !! In Pa57 pres_bottom_strato= 2000.0 !! In Pa 58 58 pres_top_strato= 100.0 !! In Pa 59 59 size_tropo=2.e-6 !! particle size in the tropospheric layer, in meters … … 62 62 write(*,*) '--- Two-layer aerosol = ', iaero_back2lay 63 63 64 65 64 write(*,*) '=== Number of aerosols= ', ia 66 65
Note: See TracChangeset
for help on using the changeset viewer.