Changeset 2324
- Timestamp:
- Jul 8, 2015, 5:20:22 PM (10 years ago)
- Location:
- LMDZ5/trunk/libf/phylmd
- Files:
-
- 3 edited
-
aeropt_2bands.F90 (modified) (12 diffs)
-
aeropt_5wv.F90 (modified) (5 diffs)
-
readaerosol_optic.F90 (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/aeropt_2bands.F90
r2311 r2324 29 29 REAL, INTENT(in) :: delt 30 30 REAL, DIMENSION(klon,klev,naero_tot), INTENT(in) :: m_allaer 31 !RAF32 31 REAL, DIMENSION(klon,klev,naero_tot), INTENT(in) :: m_allaer_pi 33 32 REAL, DIMENSION(klon,klev), INTENT(in) :: RHcl ! humidite relative ciel clair 34 !RAF REAL, DIMENSION(klon,naero_tot),INTENT(in) :: fractnat_allaer35 33 INTEGER, INTENT(in) :: flag_aerosol 36 34 REAL, DIMENSION(klon,klev), INTENT(in) :: pplay … … 48 46 ! 49 47 REAL, DIMENSION(klon,klev,naero_tot,nbands) :: tau_ae 50 !RAF51 48 REAL, DIMENSION(klon,klev,naero_tot,nbands) :: tau_ae_pi 52 49 REAL, DIMENSION(klon,klev,naero_tot,nbands) :: piz_ae … … 594 591 DO i=1, klon 595 592 zrho=pplay(i,k)/t_seri(i,k)/RD ! kg/m3 596 !CDIR UNROLL=naero_spc597 593 mass_temp(i,k,:) = m_allaer(i,k,:) / zrho / 1.e+9 598 !RAF zrho599 !CDIR UNROLL=naero_spc600 594 mass_temp_pi(i,k,:) = m_allaer_pi(i,k,:) / zrho / 1.e+9 601 595 zdp1(i,k)=pdel(i,k)/(gravit*delt) ! air mass auxiliary variable --> zdp1 [kg/(m^2 *s)] … … 653 647 654 648 655 !!CDIR ON_ADB(RH_tab)656 649 !CDIR ON_ADB(fact_RH) 657 650 !CDIR SHORTLOOP … … 661 654 662 655 DO k=1, KLEV 663 !!CDIR ON_ADB(RH_tab)664 656 !CDIR ON_ADB(fact_RH) 665 657 DO i=1, KLON … … 944 936 DO i=1, KLON 945 937 tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSO4M_phy,inu)+tau_ae(i,k,id_CSSO4M_phy,inu)+ & 946 tau_ae(i,k,id_ASBCM_phy,inu)+tau_ae(i,k,id_AIBCM_phy,inu)+ & 947 tau_ae(i,k,id_ASPOMM_phy,inu)+tau_ae(i,k,id_AIPOMM_phy,inu)+ & 938 tau_ae(i,k,id_ASBCM_phy,inu)+tau_ae(i,k,id_AIBCM_phy,inu)+ & 939 tau_ae(i,k,id_ASPOMM_phy,inu)+tau_ae(i,k,id_AIPOMM_phy,inu)+ & 948 940 tau_ae(i,k,id_ASSSM_phy,inu)+tau_ae(i,k,id_CSSSM_phy,inu)+ & 949 941 tau_ae(i,k,id_SSSSM_phy,inu)+ tau_ae(i,k,id_CIDUSTM_phy,inu) 950 942 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5) 951 943 952 953 954 tau_ae(i,k,id_ASBCM_phy,inu)*piz_ae(i,k,id_ASBCM_phy,inu)+&955 tau_ae(i,k,id_AIBCM_phy,inu)*piz_ae(i,k,id_AIBCM_phy,inu)+&956 957 tau_ae(i,k,id_AIPOMM_phy,inu)*piz_ae(i,k,id_AIPOMM_phy,inu)+ &958 tau_ae(i,k,id_ASSSM_phy,inu)*piz_ae(i,k,id_ASSSM_phy,inu)+&959 tau_ae(i,k,id_CSSSM_phy,inu)*piz_ae(i,k,id_CSSSM_phy,inu)+&960 tau_ae(i,k,id_SSSSM_phy,inu)*piz_ae(i,k,id_SSSSM_phy,inu)+&961 944 piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSO4M_phy,inu)*piz_ae(i,k,id_ASSO4M_phy,inu)+ & 945 tau_ae(i,k,id_CSSO4M_phy,inu)*piz_ae(i,k,id_CSSO4M_phy,inu)+ & 946 tau_ae(i,k,id_ASBCM_phy,inu)*piz_ae(i,k,id_ASBCM_phy,inu)+ & 947 tau_ae(i,k,id_AIBCM_phy,inu)*piz_ae(i,k,id_AIBCM_phy,inu)+ & 948 tau_ae(i,k,id_ASPOMM_phy,inu)*piz_ae(i,k,id_ASPOMM_phy,inu)+ & 949 tau_ae(i,k,id_AIPOMM_phy,inu)*piz_ae(i,k,id_AIPOMM_phy,inu)+ & 950 tau_ae(i,k,id_ASSSM_phy,inu)*piz_ae(i,k,id_ASSSM_phy,inu)+ & 951 tau_ae(i,k,id_CSSSM_phy,inu)*piz_ae(i,k,id_CSSSM_phy,inu)+ & 952 tau_ae(i,k,id_SSSSM_phy,inu)*piz_ae(i,k,id_SSSSM_phy,inu)+ & 953 tau_ae(i,k,id_CIDUSTM_phy,inu)*piz_ae(i,k,id_CIDUSTM_phy,inu)) & 962 954 /tau_allaer(i,k,mrfspecies,inu) 963 964 965 955 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),0.1) 956 957 cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSO4M_phy,inu)* & 966 958 piz_ae(i,k,id_ASSO4M_phy,inu)*cg_ae(i,k,id_ASSO4M_phy,inu)+ & 967 959 tau_ae(i,k,id_CSSO4M_phy,inu)*piz_ae(i,k,id_CSSO4M_phy,inu)* & … … 991 983 DO k=1, KLEV 992 984 DO i=1, KLON 993 !RAF 994 985 986 tau_allaer(i,k,mrfspecies,inu)=tau_ae_pi(i,k,id_ASSO4M_phy,inu)+ & 995 987 tau_ae_pi(i,k,id_CSSO4M_phy,inu)+ & 996 tau_ae_pi(i,k,id_ASBCM_phy,inu)+ &997 tau_ae_pi(i,k,id_AIBCM_phy,inu)+ &988 tau_ae_pi(i,k,id_ASBCM_phy,inu)+ & 989 tau_ae_pi(i,k,id_AIBCM_phy,inu)+ & 998 990 tau_ae_pi(i,k,id_ASPOMM_phy,inu)+ & 999 tau_ae_pi(i,k,id_AIPOMM_phy,inu)+ & 1000 tau_ae_pi(i,k,id_ASSSM_phy,inu)+ &1001 tau_ae_pi(i,k,id_CSSSM_phy,inu)+ &1002 tau_ae_pi(i,k,id_SSSSM_phy,inu)+ &991 tau_ae_pi(i,k,id_AIPOMM_phy,inu)+ & 992 tau_ae_pi(i,k,id_ASSSM_phy,inu)+ & 993 tau_ae_pi(i,k,id_CSSSM_phy,inu)+ & 994 tau_ae_pi(i,k,id_SSSSM_phy,inu)+ & 1003 995 tau_ae_pi(i,k,id_CIDUSTM_phy,inu) 1004 1005 1006 996 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5) 997 998 piz_allaer(i,k,mrfspecies,inu)=(tau_ae_pi(i,k,id_ASSO4M_phy,inu)*piz_ae(i,k,id_ASSO4M_phy,inu)+ & 1007 999 tau_ae_pi(i,k,id_CSSO4M_phy,inu)*piz_ae(i,k,id_CSSO4M_phy,inu)+ & 1008 tau_ae_pi(i,k,id_ASBCM_phy,inu)*piz_ae(i,k,id_ASBCM_phy,inu)+ &1009 tau_ae_pi(i,k,id_AIBCM_phy,inu)*piz_ae(i,k,id_AIBCM_phy,inu)+ &1000 tau_ae_pi(i,k,id_ASBCM_phy,inu)*piz_ae(i,k,id_ASBCM_phy,inu)+ & 1001 tau_ae_pi(i,k,id_AIBCM_phy,inu)*piz_ae(i,k,id_AIBCM_phy,inu)+ & 1010 1002 tau_ae_pi(i,k,id_ASPOMM_phy,inu)*piz_ae(i,k,id_ASPOMM_phy,inu)+ & 1011 tau_ae_pi(i,k,id_AIPOMM_phy,inu)*piz_ae(i,k,id_AIPOMM_phy,inu)+ & 1012 tau_ae_pi(i,k,id_ASSSM_phy,inu)*piz_ae(i,k,id_ASSSM_phy,inu)+ &1013 tau_ae_pi(i,k,id_CSSSM_phy,inu)*piz_ae(i,k,id_CSSSM_phy,inu)+ &1014 tau_ae_pi(i,k,id_SSSSM_phy,inu)*piz_ae(i,k,id_SSSSM_phy,inu)+ &1003 tau_ae_pi(i,k,id_AIPOMM_phy,inu)*piz_ae(i,k,id_AIPOMM_phy,inu)+ & 1004 tau_ae_pi(i,k,id_ASSSM_phy,inu)*piz_ae(i,k,id_ASSSM_phy,inu)+ & 1005 tau_ae_pi(i,k,id_CSSSM_phy,inu)*piz_ae(i,k,id_CSSSM_phy,inu)+ & 1006 tau_ae_pi(i,k,id_SSSSM_phy,inu)*piz_ae(i,k,id_SSSSM_phy,inu)+ & 1015 1007 tau_ae_pi(i,k,id_CIDUSTM_phy,inu)*piz_ae(i,k,id_CIDUSTM_phy,inu)) & 1016 1008 /tau_allaer(i,k,mrfspecies,inu) 1017 1018 1019 1009 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),0.1) 1010 1011 cg_allaer(i,k,mrfspecies,inu)=(& 1020 1012 tau_ae_pi(i,k,id_ASSO4M_phy,inu)*piz_ae(i,k,id_ASSO4M_phy,inu)*cg_ae(i,k,id_ASSO4M_phy,inu)+ & 1021 1013 tau_ae_pi(i,k,id_CSSO4M_phy,inu)*piz_ae(i,k,id_CSSO4M_phy,inu)*cg_ae(i,k,id_CSSO4M_phy,inu)+ & … … 1036 1028 DO k=1, KLEV 1037 1029 DO i=1, KLON 1038 1039 1040 1030 tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASBCM_phy,inu)+tau_ae(i,k,id_AIBCM_phy,inu) 1031 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5) 1032 piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASBCM_phy,inu)*piz_ae(i,k,id_ASBCM_phy,inu) & 1041 1033 +tau_ae(i,k,id_AIBCM_phy,inu)*piz_ae(i,k,id_AIBCM_phy,inu))/ & 1042 1034 tau_allaer(i,k,mrfspecies,inu) 1043 1035 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),0.1) 1044 1036 cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASBCM_phy,inu)*piz_ae(i,k,id_ASBCM_phy,inu) *cg_ae(i,k,id_ASBCM_phy,inu)& 1045 1037 +tau_ae(i,k,id_AIBCM_phy,inu)*piz_ae(i,k,id_AIBCM_phy,inu)*cg_ae(i,k,id_AIBCM_phy,inu))/ & … … 1053 1045 DO i=1, KLON 1054 1046 tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSO4M_phy,inu)+tau_ae(i,k,id_CSSO4M_phy,inu) 1055 1047 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5) 1056 1048 piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_CSSO4M_phy,inu)*piz_ae(i,k,id_CSSO4M_phy,inu) & 1057 1049 +tau_ae(i,k,id_ASSO4M_phy,inu)*piz_ae(i,k,id_ASSO4M_phy,inu))/ & 1058 1050 tau_allaer(i,k,mrfspecies,inu) 1059 1051 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),0.1) 1060 1052 cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_CSSO4M_phy,inu)* & 1061 1053 piz_ae(i,k,id_CSSO4M_phy,inu) *cg_ae(i,k,id_CSSO4M_phy,inu)& … … 1072 1064 tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASPOMM_phy,inu)+tau_ae(i,k,id_AIPOMM_phy,inu) 1073 1065 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5) 1074 1066 piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASPOMM_phy,inu)*piz_ae(i,k,id_ASPOMM_phy,inu) & 1075 1067 +tau_ae(i,k,id_AIPOMM_phy,inu)*piz_ae(i,k,id_AIPOMM_phy,inu))/ & 1076 1068 tau_allaer(i,k,mrfspecies,inu) 1077 1078 cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASPOMM_phy,inu)*piz_ae(i,k,id_ASPOMM_phy,inu)*cg_ae(i,k,id_ASPOMM_phy,inu)&1069 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),0.1) 1070 cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASPOMM_phy,inu)*piz_ae(i,k,id_ASPOMM_phy,inu)*cg_ae(i,k,id_ASPOMM_phy,inu)& 1079 1071 +tau_ae(i,k,id_AIPOMM_phy,inu)*piz_ae(i,k,id_AIPOMM_phy,inu)*cg_ae(i,k,id_AIPOMM_phy,inu))/ & 1080 1072 (tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu)) … … 1087 1079 DO i=1, KLON 1088 1080 tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_CIDUSTM_phy,inu) 1089 1081 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5) 1090 1082 piz_allaer(i,k,mrfspecies,inu)=piz_ae(i,k,id_CIDUSTM_phy,inu) 1091 1083 cg_allaer(i,k,mrfspecies,inu)=cg_ae(i,k,id_CIDUSTM_phy,inu) 1092 1084 ENDDO 1093 1085 ENDDO … … 1135 1127 ENDDO 1136 1128 ENDDO 1137 1138 1129 1139 1130 inu=1 1140 1131 DO i=1, KLON 1141 1132 absvisaer(i)=SUM((1-piz_allaer(i,:,:,inu))*tau_allaer(i,:,:,inu)) 1142 END DO1133 ENDDO 1143 1134 1144 1135 DEALLOCATE(aerosol_name) -
LMDZ5/trunk/libf/phylmd/aeropt_5wv.F90
r2311 r2324 67 67 ! 68 68 REAL, DIMENSION(klon), INTENT(out) :: ai ! POLDER aerosol index 69 ! REAL, DIMENSION(klon,nwave,naero_spc), INTENT(out) :: tausum70 ! REAL, DIMENSION(klon,klev,nwave,naero_spc), INTENT(out) :: tau71 69 REAL, DIMENSION(klon,nwave,naero_tot), INTENT(out) :: tausum 72 70 REAL, DIMENSION(klon,klev,nwave,naero_tot), INTENT(out) :: tau … … 626 624 DO k=1, klev 627 625 DO i=1, klon 628 ! IF (t_seri(i,k).EQ.0) stop 'stop aeropt_5wv T '629 ! IF (pplay(i,k).EQ.0) stop 'stop aeropt_5wv p '630 626 zrho=pplay(i,k)/t_seri(i,k)/RD ! kg/m3 631 627 dh(i,k)=pdel(i,k)/(gravit*zrho) 632 !CDIR UNROLL=naero_spc633 628 mass_temp(i,k,:) = m_allaer(i,k,:) / zrho / 1.e+9 634 629 zdp1(i,k)=pdel(i,k)/(gravit*delt) ! air mass auxiliary variable --> zdp1 [kg/(m^2 *s)] … … 760 755 ENDIF 761 756 762 !Bug 21 12 10 AI763 ! used_tau(spsol)=.TRUE.764 757 IF (soluble) then 765 758 used_tau(spsol)=.TRUE. 766 759 ELSE 767 760 used_tau(naero_soluble+spinsol)=.TRUE. 768 761 ENDIF … … 866 859 DO m=1,naero_spc 867 860 od550aer(i)=od550aer(i)+tausum(i,2,m) 868 END DO 869 END DO 861 ENDDO 862 ENDDO 863 870 864 DO i=1, klon 871 865 od865aer(i)=0. 872 866 DO m=1,naero_spc 873 867 od865aer(i)=od865aer(i)+tausum(i,5,m) 874 END DO 875 END DO 868 ENDDO 869 ENDDO 870 876 871 DO i=1, klon 877 872 DO k=1, KLEV … … 879 874 DO m=1,naero_spc 880 875 ec550aer(i,k)=ec550aer(i,k)+tau(i,k,2,m)/dh(i,k) 881 END 882 END 883 END 876 ENDDO 877 ENDDO 878 ENDDO 884 879 885 od550lt1aer(:)=tausum(:,2,id_ASSO4M_phy)+tausum(:,2,id_ASBCM_phy)+tausum(:,2,id_AIBCM_phy)+ & 886 tausum(:,2,id_ASPOMM_phy)+tausum(:,2,id_AIPOMM_phy)+tausum(:,2,id_ASSSM_phy)+ & 887 0.03*tausum(:,2,id_CSSSM_phy)+0.4*tausum(:,2,id_CIDUSTM_phy) 888 889 880 od550lt1aer(:)=tausum(:,2,id_ASSO4M_phy)+tausum(:,2,id_ASBCM_phy)+tausum(:,2,id_AIBCM_phy)+ & 881 tausum(:,2,id_ASPOMM_phy)+tausum(:,2,id_AIPOMM_phy)+tausum(:,2,id_ASSSM_phy)+ & 882 0.03*tausum(:,2,id_CSSSM_phy)+0.4*tausum(:,2,id_CIDUSTM_phy) 890 883 891 884 DEALLOCATE(aerosol_name) -
LMDZ5/trunk/libf/phylmd/readaerosol_optic.F90
r2323 r2324 68 68 REAL, DIMENSION(klon,klev) :: ssacu_pi 69 69 REAL, DIMENSION(klon,klev) :: pdel 70 REAL, DIMENSION(klon,klev,naero_ spc) :: m_allaer71 REAL, DIMENSION(klon,klev,naero_ spc) :: m_allaer_pi !RAF70 REAL, DIMENSION(klon,klev,naero_tot) :: m_allaer 71 REAL, DIMENSION(klon,klev,naero_tot) :: m_allaer_pi !RAF 72 72 ! REAL, DIMENSION(klon,naero_tot) :: fractnat_allaer !RAF delete?? 73 73 … … 156 156 m_allaer(:,:,id_AIBCM_phy) = bcins(:,:) ! AIBCM 157 157 m_allaer(:,:,id_AIPOMM_phy) = pomins(:,:) ! AIPOMM 158 m_allaer(:,:,id_ASNO3M_phy) = 0.0 159 m_allaer(:,:,id_CSNO3M_phy) = 0.0 160 m_allaer(:,:,id_CINO3M_phy) = 0.0 158 161 159 162 !RAF 160 m_allaer_pi(:,:,1) = bcsol_pi(:,:) ! ASBCM pre-ind 161 m_allaer_pi(:,:,2) = pomsol_pi(:,:) ! ASPOMM pre-ind 162 m_allaer_pi(:,:,3) = sulfate_pi(:,:) ! ASSO4M (= SO4) pre-ind 163 m_allaer_pi(:,:,4) = 0. ! CSSO4M pre-ind 164 m_allaer_pi(:,:,5) = sssupco_pi(:,:) ! SSSSM pre-ind 165 m_allaer_pi(:,:,6) = sscoarse_pi(:,:) ! CSSSM pre-ind 166 m_allaer_pi(:,:,7) = ssacu_pi(:,:) ! ASSSM pre-ind 167 m_allaer_pi(:,:,8) = cidust_pi(:,:) ! CIDUSTM pre-ind 168 m_allaer_pi(:,:,9) = bcins_pi(:,:) ! AIBCM pre-ind 169 m_allaer_pi(:,:,10) = pomins_pi(:,:) ! AIPOMM pre-ind 163 m_allaer_pi(:,:,id_ASBCM_phy) = bcsol_pi(:,:) ! ASBCM pre-ind 164 m_allaer_pi(:,:,id_ASPOMM_phy) = pomsol_pi(:,:) ! ASPOMM pre-ind 165 m_allaer_pi(:,:,id_ASSO4M_phy) = sulfate_pi(:,:) ! ASSO4M (= SO4) pre-ind 166 m_allaer_pi(:,:,id_CSSO4M_phy) = 0. ! CSSO4M pre-ind 167 m_allaer_pi(:,:,id_SSSSM_phy) = sssupco_pi(:,:) ! SSSSM pre-ind 168 m_allaer_pi(:,:,id_ASSSM_phy) = sscoarse_pi(:,:) ! CSSSM pre-ind 169 m_allaer_pi(:,:,id_CIDUSTM_phy) = ssacu_pi(:,:) ! ASSSM pre-ind 170 m_allaer_pi(:,:,id_AIBCM_phy) = cidust_pi(:,:) ! CIDUSTM pre-ind 171 m_allaer_pi(:,:,id_AIBCM_phy) = bcins_pi(:,:) ! AIBCM pre-ind 172 m_allaer_pi(:,:,id_AIPOMM_phy) = pomins_pi(:,:) ! AIPOMM pre-ind 173 m_allaer_pi(:,:,id_ASNO3M_phy) = 0.0 174 m_allaer_pi(:,:,id_CSNO3M_phy) = 0.0 175 m_allaer_pi(:,:,id_CINO3M_phy) = 0.0 170 176 171 177 ! … … 212 218 ! Diagnostics calculation for CMIP5 protocol 213 219 sconcso4(:)=m_allaer(:,1,id_ASSO4M_phy)*1.e-9 214 !sconcno3(:)=m_allaer(:,1,id_ASNO3M_phy)*1.e-9220 sconcno3(:)=m_allaer(:,1,id_ASNO3M_phy)*1.e-9 215 221 sconcoa(:)=(m_allaer(:,1,id_ASPOMM_phy)+m_allaer(:,1,id_AIPOMM_phy))*1.e-9 216 222 sconcbc(:)=(m_allaer(:,1,id_ASBCM_phy)+m_allaer(:,1,id_AIBCM_phy))*1.e-9 … … 218 224 sconcdust(:)=m_allaer(:,1,id_CIDUSTM_phy)*1.e-9 219 225 concso4(:,:)=m_allaer(:,:,id_ASSO4M_phy)*1.e-9 220 !concno3(:,:)=m_allaer(:,:,id_ASNO3M_phy)*1.e-9226 concno3(:,:)=m_allaer(:,:,id_ASNO3M_phy)*1.e-9 221 227 concoa(:,:)=(m_allaer(:,:,id_ASPOMM_phy)+m_allaer(:,:,id_AIPOMM_phy))*1.e-9 222 228 concbc(:,:)=(m_allaer(:,:,id_ASBCM_phy)+m_allaer(:,:,id_AIBCM_phy))*1.e-9
Note: See TracChangeset
for help on using the changeset viewer.
