Changeset 1347 for LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/aeropt_5wv.F90
- Timestamp:
- Apr 13, 2010, 5:12:56 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/aeropt_5wv.F90
r1299 r1347 11 11 USE DIMPHY 12 12 USE aero_mod 13 USE phys_local_var_mod, only: od550aer,od865aer,ec550aer,od550lt1aer 13 14 14 15 ! … … 141 142 REAL :: tau3d(KLON,KLEV), piz3d(KLON,KLEV), cg3d(KLON,KLEV) 142 143 REAL :: abs3d(KLON,KLEV) ! epaisseur optique d'absorption 143 144 REAL :: dh(KLON,KLEV) 144 145 145 146 REAL :: alpha_aers_5wv(nbre_RH,las,naero_soluble) ! ext. coeff. Soluble comp. units *** m2/g … … 621 622 DO k=1, klev 622 623 DO i=1, klon 624 ! IF (t_seri(i,k).EQ.0) stop 'stop aeropt_5wv T ' 625 ! IF (pplay(i,k).EQ.0) stop 'stop aeropt_5wv p ' 623 626 zrho=pplay(i,k)/t_seri(i,k)/RD ! kg/m3 627 dh(i,k)=pdel(i,k)/(gravit*zrho) 624 628 !CDIR UNROLL=naero_spc 625 629 mass_temp(i,k,:) = m_allaer(i,k,:) / zrho / 1.e+9 … … 846 850 ! ENDDO 847 851 852 DO i=1, klon 853 od550aer(i)=SUM(tausum(i,2,:)) 854 od865aer(i)=SUM(tausum(i,5,:)) 855 DO k=1, KLEV 856 ec550aer(i,k)=SUM(tau(i,k,2,:))/dh(i,k) 857 END DO 858 END DO 859 od550lt1aer(:)=tausum(:,2,id_ASSO4M)+tausum(:,2,id_ASBCM)+tausum(:,2,id_AIBCM)+ & 860 tausum(:,2,id_ASPOMM)+tausum(:,2,id_AIPOMM)+tausum(:,2,id_ASSSM)+ & 861 0.03*tausum(:,2,id_CSSSM)+0.4*tausum(:,2,id_CIDUSTM) 862 863 864 848 865 DEALLOCATE(aerosol_name) 849 866
Note: See TracChangeset
for help on using the changeset viewer.